We investigate the possibility of a supernova in supermassive ($5 \times 10^4 \;M_\odot$) population III stars induced by a general relativistic instability occurring in the helium burning phase. This explosion could occur via rapid helium burning during an early contraction of the isentropic core. Such an explosion would be visible to future telescopes and could disrupt the proposed direct collapse formation channel for early universe supermassive black holes. We simulate first the stellar evolution from hydrogen burning using a 1D stellar evolution code with a post Newtonian approximation; at the point of dynamical collapse, we switch to a 1D (general relativistic) hydrodynamics code with the Misner-Sharpe metric. In opposition to a previous study, we do not find an explosion in the non rotating case, although our model is close to exploding for a similar mass to the explosion in the previous study. When we include slow rotation, we find one exploding model, and we conclude that there likely exist additional exploding models, though they may be rare.

We have devised fitting formulae for evolution tracks of massive stars with $8 \lesssim M/M_\odot \lesssim 160$ under extreme metal poor (EMP) environments for $\log (Z/Z_\odot) = -2, -4, -5, -6$, and $-8$, where $M_\odot$ and $Z_\odot$ are the solar mass and metallicity, respectively. Our fitting formulae are based on reference stellar models which we have newly obtained by simulating the time evolutions of EMP stars. Our fitting formulae take into account stars ending with blue supergiant (BSG) stars, and stars skipping Hertzsprung gap (HG) phases and blue loops, which are characteristics of massive EMP stars. In our fitting formulae, stars may remain BSG stars when they finish their core Helium burning (CHeB) phases. Our fitting formulae are in good agreement with our stellar evolution models. We can use these fitting formulae on the SSE, BSE, NBODY4, and NBODY6 codes, which are widely used for population synthesis calculations and star cluster simulations. These fitting formulae should be useful to make theoretical templates of binary black holes formed under EMP environments.

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the $9.6\, \rm M_\odot$ progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the $9.6\, \rm M_\odot$ progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

We report, for the first time, the detection of the Mn-K$\alpha$ line in the Type IIb supernova (SN IIb) remnant, Cassiopeia A. Manganese ($^{55}$Mn after decay of $^{55}$Co), a neutron-rich element, together with chromium ($^{52}$Cr after decay of $^{52}$Fe), is mainly synthesized at the explosive incomplete Si burning regime. Therefore, the Mn/Cr mass ratio with its neutron excess reflects the neutronization at the relevant burning layer during the explosion. Chandra's archival X-ray data of Cassiopeia A indicate a low Mn/Cr mass ratio with values in the range 0.10--0.66, which, when compared to one-dimensional SN explosion models, requires that the electron fraction be 0.4990 $\lesssim Y_{\rm e} \lesssim$ 0.5 at the incomplete Si burning layer. An explosion model assuming a solar-metallicity progenitor with a typical explosion energy ($1 \times 10^{51}$ erg) fails to reproduce such a high electron fraction. In such models, the explosive Si-burning regime extends only to the Si/O layer established during the progenitor's hydrostatic evolution; the $Y_e$ in the Si/O layer is lower than the value required by our observational constraints. We can satisfy the observed Mn/Cr mass ratio if the explosive Si-burning regime were to extend into the O/Ne hydrostatic layer, which has a higher $Y_{\rm e}$. This would require an energetic ($> 2 \times 10^{51}$ erg) and/or asymmetric explosion of a sub-solar metallicity progenitor ($Z \lesssim 0.5Z_{\odot}$) for Cassiopeia A. The low initial metallicity can be used to rule out a single-star progenitor, leaving the possibility of a binary progenitor with a compact companion (white dwarf, neutron star or black hole). We discuss the detectability of X-rays from Bondi accretion onto such a compact companion around the explosion site. We also discuss other possible mass-loss scenarios for the progenitor system of Cassiopeia A.

We report on robust measurements of elemental abundances of the Type IIn supernova SN 1978K, based on the high-resolution X-ray spectrum obtained with the Reflection Grating Spectrometer (RGS) onboard XMM-Newton. The RGS clearly resolves a number of emission lines, including N Ly$\alpha$, O Ly$\alpha$, O Ly$\beta$, Fe XVII, Fe XVIII, Ne He$\alpha$ and Ne Ly$\alpha$ for the first time from SN 1978K. The X-ray spectrum can be represented by an absorbed, two-temperature thermal emission model, with temperatures of $kT \sim 0.6$ keV and $2.7$ keV. The elemental abundances are obtained to be N $=$ $2.36_{-0.80}^{+0.88}$, O $=$ $0.20 \pm{0.05}$, Ne $=$ $0.47 \pm{0.12}$, Fe $=$ $0.15_{-0.02}^{+0.01}$ times the solar values. The low metal abundances except for N show that the X-ray emitting plasma originates from the circumstellar medium blown by the progenitor star. The abundances of N and O are far from CNO-equilibrium abundances expected for the surface composition of a luminous blue variable, and resemble the H-rich envelope of less-massive stars with masses of 10-25 M$_\odot$. Together with other peculiar properties of SN 1978K, i.e., a low expansion velocity of 500-1000 km s$^{-1}$ and SN IIn-like optical spectra, we propose that SN 1978K is a result of either an electron-capture SN from a super asymptotic giant branch star, or a weak Fe core-collapse explosion of a relatively low-mass ($\sim$10 M$_\odot$) or high-mass ($\sim$20-25 M$_\odot$) red supergiant star. However, these scenarios can not naturally explain the high mass-loss rate of the order of $\dot{M} \sim 10^{-3} \rm{M_{\odot}\ yr^{-1}}$ over $\gtrsim$1000 yr before the explosion, which is inferred by this work as well as many other earlier studies. Further theoretical studies are required to explain the high mass-loss rates at the final evolutionary stages of massive stars.

(Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 years, accounting for the pre-SN structure of the progenitor star. The simulations include all physical processes relevant for the complex phases of SN evolution and for the interaction of the SNR with the highly inhomogeneous ambient environment around SN 1987A. From the simulations, we synthesize observables to be compared with observations. By comparing the model results with observations, we constrained the initial SN anisotropy causing Doppler shifts observed in emission lines of heavy elements from ejecta, and leading to the remnant evolution observed in the X-ray band in the last 30 years. In particular, we found that the high mixing of ejecta unveiled by high redshifts and broadenings of [FeII] and $^{44}$Ti lines require a highly asymmetric SN explosion channeling a significant fraction of energy along an axis almost lying in the plane of the central equatorial ring around SN 1987A, roughly along the line-of-sight but with an offset of 40 deg, with the lobe propagating away from the observer slightly more energetic than the other. We found unambiguously that the observed distribution of ejecta and the dynamical and radiative properties of the SNR can be best reproduced if the structure of the progenitor star was that of a blue supergiant resulted from the merging of two massive stars.

We perform three-dimensional hydrodynamic simulations of aspherical core-collapse supernovae focusing on the matter mixing in SN 1987A. The impacts of four progenitor (pre-supernova) models and parameterized aspherical explosions are investigated. The four pre-supernova models include a blue supergiant (BSG) model based on a slow merger scenario developed recently for the progenitor of SN 1987A (Urushibata et al. 2018). The others are a BSG model based on a single star evolution and two red supergiant (RSG) models. Among the investigated explosion (simulation) models, a model with the binary merger progenitor model and with an asymmetric bipolar-like explosion, which invokes a jetlike explosion, best reproduces constraints on the mass of high velocity $^{56}$Ni, as inferred from the observed [Fe II] line profiles. The advantage of the binary merger progenitor model for the matter mixing is the flat and less extended $\rho \,r^3$ profile of the C+O core and the helium layer, which may be characterized by the small helium core mass. From the best explosion model, the direction of the bipolar explosion axis (the strongest explosion direction), the neutron star (NS) kick velocity, and its direction are predicted. Other related implications and future prospects are also given.

Double neutron star (DNS) systems are produced from massive binaries. A supernova (SN) explosion of an extremely stripped star is expected to occur at the final stage of DNS formation. This type of SN is called ultra-stripped SN (USSN). A recent research reveals that a type Ic SN, iPTF 14gqr (SN 2014ft) has low ejecta mass ($\approx0.2~\mathrm{M}_\odot$) and its progenitor has a helium envelope with its mass $\sim0.01~\mathrm{M}_\odot$. Then this SN is interpreted as an USSN, and this is the first discovery of the USSN. Furthermore, the observation of iPTF 14gqr provides us with some information about its formation history. Here, we perform population synthesis calculations so as to estimate the detection rate of iPTF 14gqr like USSNe by optical transient surveys, intermediate Palomar Transient Factory (iPTF), Zwicky Transient Facility (ZTF), and Large Synoptic Survey Telescope (LSST). We find that iPTF, ZTF, and LSST can observe iPTF 14gqr like USSNe at the rates of 1, 40, and 5 $\mathrm{yr}^{-1}$, respectively. The iPTF may detect 4 iPTF 14gqr like USSNe during four years observation. Additionally, we investigate effects of some binary physics on the orbital parameters of DNS systems and formation channels.

We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution code, we first calculate the evolution of massive stars with an initial mass of 9-40 $M_\odot$. Four different overshoot parameters are applied, and CO core mass trend similar to previous works is obtained in the 1D models. Selecting eleven 1D models that have a silicon and oxygen coexisting layer, we perform 2D hydrodynamics simulations of the evolution $\sim$100 s until the onset of core-collapse. We find that convection with large-scale eddies and the turbulent Mach number $\sim$0.1 is obtained in the models having a Si/O layer with a scale of 10$^8$ cm, whereas most models that have an extended O/Si layer up to a few $\times 10^9$ cm exhibit lower turbulent velocity. Our results indicate that the supernova progenitors that possess a thick Si/O layer could provide a preferable condition for perturbation-aided explosions. We perform 3D simulation of a 25 $M_\odot$ model, which exhibits large-scale convection in the 2D models. The 3D model develops large ($\ell = 2$) convection similar to the 2D model, however, the turbulent velocity is lower. By estimating the neutrino emission properties of the 3D model, we point out that a time modulation of the event rates, if observed in KamLAND and Hyper-Kamiokande, would provide an important information about structural changes in the presupernova convective layer.

We investigate remnant neutron star masses (in particular, the minimum allowed mass) by performing advanced stellar evolution calculations and neutrino-radiation hydrodynamics simulations for core-collapse supernova explosions. We find that, based on standard astrophysical scenarios, low-mass carbon-oxygen cores can have sufficiently massive iron cores that eventually collapse, explode as supernovae, and give rise to remnant neutron stars that have a minimum mass of 1.17 M$_\odot$ --- compatible with the lowest mass of the neutron star precisely measured in a binary system of PSR J0453+1559.

A critical mass ONe core with a high ignition density is considered to end in gravitational collapse leading to neutron star formation. Being distinct from a Fe core collapse, the final evolution involves combustion flame propagation, in which complex phase transition from ONe elements into the nuclear-statistical-equilibrium (NSE) state takes place. We simulate the core evolution from the O+Ne ignition until the bounce shock penetrates the whole core, using a state-of-the-art 1D Lagrangian neutrino-radiation-hydrodynamic code, in which important nuclear burning, electron capture, and neutrino reactions are taken into account. Special care is also taken in making a stable initial condition by importing the stellar EOS, which is used for the progenitor evolution calculation, and by improving the remapping process. We find that the central ignition leads to intense $\nu_e$ radiation with $L_{\nu_e} \gtrsim 10^{51}$ erg s$^{-1}$ powered by fast electron captures onto NSE isotopes. This pre-bounce $\nu_e$ radiation heats the surroundings by the neutrino-electron scattering, which acts as a new driving mechanism of the flame propagation together with the adiabatic contraction. The resulting flame velocity of $\sim10^8$ cm s$^{-1}$ will be more than one-order-of-magnitude faster than that of laminar flame driven by heat conduction. We also find that the duration of the pre-bounce $\nu_e$ radiation phase depends on the degree of the core hydrostatic/dynamical stability. Therefore, the future detection of the pre-bounce neutrino is important not only to discriminate the ONe core collapse from the Fe core collapse but also to constrain the progenitor hydrodynamical stability.

We investigate the collective neutrino oscillations under the three flavor multi-angle approximation in a spherically symmetric simulation of failed supernovae. A failed supernova emits high neutrino fluxes in a short time, while intense accretion proceeds with high electron density enough to experience re-collapse into a black hole. Our results show that matter-induced effects completely dominate over neutrino self-interaction effects and multi-angle matter suppression occurs at all time snapshots we studied. These facts suggest us that only MSW resonances affect the neutrino flavor conversions in failed supernovae and simple spectra will be observed at neutrino detectors. We also estimate the neutrino event rate in current and future neutrino detectors from a source at 10 kpc as a Galactic event. The time evolution of neutrino detection could provide information about the dense and hot matter and constrain the neutrino mass ordering problem.

Recent theory predicts that a first star is born with a massive initial mass of $\gtrsim$ 100 $M_\odot$. Pair instability supernova (PISN) is a common fate for such a massive star. Our final goal is to prove the existence of PISN and thus the high mass nature of the initial mass function in the early universe by conducting {\it abundance profiling}, in which properties of a hypothetical first star is constrained by metal-poor star abundances. In order to determine reliable and useful abundances, we investigate the PISN nucleosynthesis taking both rotating and non-rotating progenitors for the first time. We show that the initial and CO core mass ranges for PISNe depend on the envelope structures: non-magnetic rotating models developing inflated envelopes have a lower-shifted CO mass range of $\sim$ 70--125 $M_\odot$, while non-rotating and magnetic rotating models with deflated envelopes have a range of $\sim$ 80--135 $M_\odot$. However, we find no significant difference in explosive yields from rotating and non-rotating progenitors, except for large nitrogen production in non-magnetic rotating models. Furthermore, we conduct the first systematic comparison between theoretical yields and a large sample of metal-poor star abundances. We find that the predicted low [Na/Mg] $\sim$ $-1.5$ and high [Ca/Mg] $\sim$ $0.5$--$1.3$ abundance ratios are the most important to discriminate PISN signatures from normal metal-poor star abundances, and confirm that no currently observed metal-poor star matches with the PISN abundance. Extensive discussion on the non-detection is finally made.

Even after elaborate investigations for 30 years, we still do not know well how the progenitor of SN 1987A has evolved. To explain unusual red-to-blue evolution, previous studies suggest that in a red giant stage either the increase of surface He abundance or the envelope mass was necessary. It is usually supposed that the He enhancement is caused by the rotational mixing, and the mass increase is by a binary merger. Thus, we have investigated these scenarios thoroughly. The obtained findings are that rotating single star models do not satisfy all the observational constraints and that the enhancement of envelope mass alone does not explain observations. Here, we consider a slow merger scenario in which both the He abundance and the envelope mass enhancements are expected to occur. We indeed show that most observational constraints such as the red-to-blue evolution, lifetime, total mass, position in the HR diagram at collapse, and the chemical anomalies are well reproduced by the merger model of 14 and 9 M$_{\odot}$ stars. We also discuss the effects of the added envelope spin in the merger scenarios.

This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\bar{\nu}_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($\nu_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $\nu_e$ and $\bar{\nu}_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\bar{\nu}_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $\nu_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\bar{\nu}_e$ events will be largest for JUNO, 200-900 $\bar{\nu}_e$'s, depending on the mass hierarchy whereas the number of $\nu_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\bar{\nu}_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $\nu_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.

We investigate the explosive nucleosynthesis during two dimensional neutrino-driven explosion of ultra-stripped Type Ic supernovae evolved from 1.45 and 1.5 M$_\odot$ CO stars. These supernovae explode with the explosion energy of $\sim 10^{50}$ erg and release $\sim 0.1$ M$_\odot$ ejecta. The light trans-iron elements Ga-Zr are produced in the neutrino-irradiated ejecta. The abundance distribution of these elements has a large uncertainty because of the uncertainty of the electron fraction of the neutrino-irradiated ejecta. The yield of the elements will be less than 0.01 M$_\odot$. Ultra-stripped supernova and core-collapse supernova evolved from a light CO core can be main sources of the light trans-iron elements. They could also produce neutron-rich nuclei $^{48}$Ca. The ultra-stripped supernovae eject $^{56}$Ni of $\sim$ 0.006 - 0.01 M$_\odot$. If most of neutrino-irradiated ejecta is proton-rich, $^{56}$Ni will be produced more abundantly. The light curves of these supernovae indicate sub-luminous fast decaying explosion with the peak magnitude of about $-15$ - $-16$. Future observations of ultra-stripped supernovae could give a constraint to the event rate of a class of neutron star mergers.

Aiming to study GRB engine duration, we present numerical simulations to investigate collapsar jets. We consider typical explosion energy ($10^{52}$ erg) but different engine durations, in the widest domain to date from 0.1 to 100 s. We employ an AMR 2D hydrodynamical code. Our results show that engine duration strongly influences jet nature. We show that the efficiency of launching and collimating relativistic outflow increases with engine duration, until the intermediate engine range where it is the highest, past this point to long engine range, the trend is slightly reversed; we call this point where acceleration and collimation are the highest "sweet spot" ($\sim$10 -- 30 s). Moreover, jet energy flux shows that variability is also high in this duration domain. We argue that not all engine durations can produce the collimated, relativistic and variable long GRB-jets. Considering a typical progenitor and engine energy, we conclude that the ideal engine duration to reproduce a long GRB is $\sim$10 -- 30 s, where the launch of relativistic, collimated and variable jets is favored. We note that this duration domain makes a good link with Lazzati et al. (2013), which suggested that the bulk of BATSE's long GRBs is powered by $\sim$10 -- 20 s collapsar engines.

Supermassive stars (SMSs) of mass $\gtrsim 10^5 M_\odot$ are candidates for seeds of supermassive black holes found in the center of many massive galaxies. We simulate the gravitational collapse of a rigidly rotating SMS core including nuclear burning effects in axisymmetric numerical-relativity simulation. We find that for realistic initial conditions, the nuclear burning does not play an important role. After the collapse, a torus surrounding a rotating black hole is formed and a fraction of the torus material is ejected. We quantitatively study the relation between the properties of these objects and rotation. We find that if a SMS core is sufficiently rapidly rotating, the torus and outflow mass have approximately $6\%$ and $1\%$ of the initial mass, respectively. The typical average velocity and the total kinetic energy of the outflow are $0.2~c$ and $10^{54-56}$ erg where $c$ is the speed of light. Finally, we briefly discuss the possibility for observing the torus and outflow.

The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $\gtrsim 10^5~M_\odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $\dot M \lesssim 0.1~M_\odot~{\rm yr}^{-1}$. With $\dot{M} \simeq 0.3 - 1~M_\odot~{\rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M \simeq 2 - 3.5~\times 10^5~M_\odot$. In an extreme case with $10~M_\odot~{\rm yr}^{-1}$, the star does not collapse until the mass reaches $\simeq 8.0\times 10^5~M_\odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.

When a supernova explosion occurs in neighbors around hundreds pc, current and future neutrino detectors are expected to observe neutrinos from the presupernova star before the explosion. We show a possibility for obtaining the evidence for burning processes in the central region of presupernova stars though the observations of neutrino signals by current and future neutrino detectors such as KamLAND, JUNO, and Hyper-Kamiokande. We also investigate supernova alarms using neutrinos from presupernova stars in neighbors. If a supernova explodes at ~ 200 pc, future 20 kton size liquid scintillation detectors are expected to observe hundreds neutrino events. We also propose a possibility of the detection of neutrino events by Gd-loaded Hyper-Kamiokande using delayed $\gamma$-ray signals. These detectors could observe detailed time variation of neutrino events. The neutrino emission rate increases by the core contraction in the final evolution stage. However, the O and Si shell burnings suppress the neutrino emission for a moment. The observed decrease in the neutrino event rate before hours to the explosion is possibly evidence for the shell burnings. The observations of detailed time evolution of presupernova neutrino events could reveal properties of burning processes in the central region of presupernova stars.

We derive the gravitational waveform from the collapse of a rapidly rotating supermassive star (SMS) core leading directly to a seed of a supermassive black hole (SMBH) in axisymmetric numerical-relativity simulations. We find that the peak strain amplitude of gravitational waves emitted during the black-hole formation is $\approx 5 \times 10^{-21}$ at the frequency $f \approx 5$\,mHz for an event at the cosmological redshift $z=3$, if the collapsing SMS core is in the hydrogen-burning phase. Such gravitational waves will be detectable by space laser interferometric detectors like eLISA with signal-to-noise ratio $\approx 10$, if the sensitivity is as high as LISA for $f=1$--10\,mHz. The detection of the gravitational-wave signal will provide a potential opportunity for testing the direct-collapse scenario for the formation of a seed of SMBHs.

Massive stars having a CO core of $\sim$40-60 M$_\odot$ experience pulsational pair-instability (PPI) after carbon-burning. This instability induces strong pulsations of the whole star and a part of outer envelope is ejected. We investigate the evolution and mass ejection of metal-poor very massive stars which experience PPI. We use stellar models with initial masses of 140, 200, and 250 M$_\odot$ and the metallicity Z=0.004. Their masses decrease to 54.09, 58.65, and 61.03 M$_\odot$ before the neon-burning owing to mass-loss and He mass fraction at the surface becomes about 20%. During the PPI period of $\sim$1-2000 yr, they experience six, four, and three pulsations, respectively. The larger CO-core model has the longer PPI period and ejects the larger amount of mass. Since almost all surface He has been lost by the pulsations, these stars become Type Ic supernovae if they explode. Light curves during the PPI stage and supernovae are investigated and are implicated in luminous supernovae. The luminosity created by the interaction of different PPI ejecta becomes $M_{\rm bol} \sim -16$ to $-20$. The interaction between the circumstellar shell ejected by PPI and the supernova ejecta can be more luminous. These luminous transients could be an origin of Type I superluminous supernovae and supernovae with precursor.

Energetics of nuclear reaction is fundamentally important to understand the mechanism of pair instability supernovae (PISNe). Based on the hydrodynamic equations and thermodynamic relations, we derive exact expressions for energy conservation suitable to be solved in simulation. We also show that some formulae commonly used in the literature are obtained as approximations of the exact expressions. We simulate the evolution of very massive stars of ~100-320 Msun with zero- and 1/10 Zsun, and calculate further explosions as PISNe, applying each of the exact and approximate formulae. The calculations demonstrate that the explosion properties of PISN, such as the mass range, the 56Ni yield, and the explosion energy, are significantly affected by applying the different energy generation rates. We discuss how these results affect the estimate of the PISN detection rate, which depends on the theoretical predictions of such explosion properties.

We study explosion characteristics of ultra-stripped supernovae (SNe), which are candidates of SNe generating binary neutron stars (NSs). As a first step, we perform stellar evolutionary simulations of bare carbon-oxygen cores of mass from 1.45 to 2.0 $M_\odot$ until the iron cores become unstable and start collapsing. We then perform axisymmetric hydrodynamics simulations with spectral neutrino transport using these stellar evolution outcomes as initial conditions. All models exhibit successful explosions driven by neutrino heating. The diagnostic explosion energy, ejecta mass, Ni mass, and NS mass are typically $\sim 10^{50}$ erg, $\sim 0.1 M_\odot$, $\sim 0.01M_\odot$, and $\approx 1.3 M_\odot$, which are compatible with observations of rapidly-evolving and luminous transient such as SN 2005ek. We also find that the ultra-stripped SN is a candidate for producing the secondary low-mass NS in the observed compact binary NSs like PSR J0737-3039.

Aiming to distinguish two types of progenitors of core collapse supernovae, i.e., one with a core composed mainly of oxygen and neon (abbreviated as ONe core) and the other with an iron core (or Fe core), we calculated the luminosities and spectra of neutrinos emitted from these cores prior to gravitational collapse, taking neutrino oscillation into account. We found that the total energies emitted as $\bar{\nu}_e$ from the ONe core are $\lesssim 10^{46}\ {\rm erg}$, which is much smaller than $\sim 10^{47}\ {\rm erg}$ for Fe cores. The average energy, on the other hand, is twice as large for the ONe core as those for the Fe cores. The neutrinos produced by the plasmon decays in the ONe core are more numerous than those from the electron-positron annihilation in both cores but they have much lower average energies $\lesssim 1\ {\rm MeV}$. Although it is difficult to detect the pre-supernova neutrinos from the ONe core even if it is located within 200$\ $pc from the earth, we expect $ \sim 9 - 43$ and $\sim 7 - 61$ events for Fe cores at KamLAND and Super-Kamiokande, respectively, depending on the progenitor mass and neutrino-mass hierarchy. These numbers might be increased by an order of magnitude if we envisage next-generation detectors such as JUNO. We will hence be able to distinguish the two types of progenitors by the detection or non-detection of the pre-supernova neutrinos if they are close enough ($\lesssim1\ {\rm kpc}$).

We perform stellar evolution simulation of first stars and calculate stellar yields from the first supernovae. The initial masses are taken from 12 to 140 Msun to cover the whole range of core-collapse supernova progenitors, and stellar rotation is included, which results in efficient internal mixing. A weak explosion is assumed in supernova yield calculations, thus only outer distributed matter, which is not affected by the explosive nucleosynthesis, is ejected in the models. We show that the initial mass and the rotation affect the explosion yield. All the weak explosion models have abundances of [C/O] larger than unity. Stellar yields from massive progenitors of > 40-60 Msun show enhancement of Mg and Si. Rotating models yield abundant Na and Al. And Ca is synthesized in non-rotating heavy massive models of > 80 Msun. We fit the stellar yields to the three most iron-deficient stars, and constrain the initial parameters of the mother progenitor stars. The abundance pattern in SMSS 0313-6708 is well explained by 50-80 Msun non-rotating models, rotating 30-40 Msun models well fit the abundance of HE 0107-5240, and both non-rotating and rotating 15-40 Msun models explain HE 1327-2326. The presented analysis will be applicable to other carbon-enhanced hyper metal poor stars observed in the future. The abundance analyses will give valuable information about the characteristics about the first stars.

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

We investigate the possibility of a super-luminous Type Ic core-collapse supernovae producing a large amount of 56Ni. Very massive stars with a main-sequence mass larger than 100 Msun and a metallicity 0.001 < Z < 0.004 are expected to explode as super-luminous Type Ic supernovae. Stars with ~ 110 - 150 Msun and Z < 0.001 would explode as Type Ic pulsational pair-instability supernovae if the whole H and He layers has been lost by the mass loss during pulsational pair-instability. We evaluate the total ejecta mass and the yields of 56Ni, O, and Si in core-collapse supernovae evolved from very massive stars. We adopt 43.1 and 61.1 Msun WO stars with Z=0.004 as supernova progenitors expected to explode as Type Ic core-collapse supernovae. These progenitors have masses of 110 and 250 Msun at the zero-age main sequence. Spherical explosions with an explosion energy larger than 2e52 erg produce more than 3.5 Msun 56Ni, enough to reproduce the light curve of SN 2007bi. Asphericity of the explosion affects the total ejecta mass as well as the yields of 56Ni, O, and Si. Aspherical explosions of the 110 and 250 Msun models reproduce the 56Ni yield of SN 2007bi. These explosions will also show large velocity dispersion. An aspherical core-collapse supernova evolved from a very massive star is a possibility of the explosion of SN 2007bi.

We provide progenitor models for electron capture supernovae (ECSNe) with detailed evolutionary calculation. We include minor electron capture nuclei using a large nuclear reaction network with updated reaction rates. For electron captures, the Coulomb correction on the rates is treated and contribution of neutron-rich nuclei is taken into account in the nuclear statistical equilibrium (NSE) composition. We calculate the evolution of the most massive super asymptotic giant branch stars and show that these stars undergo off-center carbon burnings and form ONe cores at the center. These cores get heavier up to the critical mass of 1.367 Msun and keep contracting even after the initiation of O+Ne deflagration. Though inclusion of minor electron capture nuclei causes convective URCA process at the contraction phase, such process will have minor effect on the evolution. On the other hand, electron captures by neutron-rich nuclei in the NSE region have more significant effect. Also, we discuss the uniqueness of the critical core mass for ECSNe and the effect of mass loss on the plausibility of our models for ECSN progenitors.

We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.

In order to explore various aspects of stellar evolution, supernovae, gamma ray bursts and nucleosynthesis, we have developed a new efficient stellar evolution code. In this paper we describe this new code and compare the results with the ones calculated by the previous code. Specifically we focus on the progenitor evolution of lower end of the Fe-core collapse supernovae, and mass distribution of remnant neutron stars. We describe how different assumptions will lead different neutron star mass distribution. We also review recent works of our research group.

Very massive stars are thought to be formed in the early Universe because of a lack of cooling process by heavy elements, and might have been responsible for the later evolution of the Universe. We had an interest in vibrational stability of their evolution and carried out the linear nonadiabatic analysis of radial and nonradial oscillations for population III very massive main-sequence stars with $500-3000M_{\sun}$. We found that only the radial fundamental mode becomes unstable due to the $\varepsilon$-mechanism for these stars. The instability appears just after the CNO cycle is activated and the nuclear energy generation rate becomes large enough to stop the pre--main-sequence contraction, and continues during the early stage of the core hydrogen burning. Besides, we roughly estimated amount of mass loss due to the instability to know its significance.

We calculate the evolution of fluorine in the solar neighborhood with the \nu-process of core-collapse supernovae, the results of which are in good agreement with the observations of field stars. The \nu-process operating in supernovae causes the [F/O] ratio to plateau at [O/H] < -1.2, followed by a rapid increase toward [O/H] ~ -0.5 from the contribution of Asymptotic Giant Branch stars. The plateau value of [F/O] depends on the neutrino luminosity released by core-collapse supernovae and may be constrained by using future observations of field stars at low metallicities. For globular clusters, the handful of [F/O] measurements suggest that the relative contribution from low-mass supernovae is smaller in these systems than in the field.

We study the role of the unknown microphysical properties of carbonaceous dust particles in determining the amount and size distribution of carbonaceous dust condensed in type II supernova explosions. We parametrize the microphysical properties in terms of the shape factor of the grain and the sticking coefficient of gas-phase carbon atoms onto the grain surfaces. We find that the amount of dust formed is fairly independent of these properties, within the parameter range considered, though limited by the available amount of carbon atoms not locked in CO molecules. However, we find that the condensation times and size distributions of dust grains depend sensitively on the microphysical parameters, with the mass distributions being weighted toward larger effective radii for conditions considering grains with higher sticking coefficients and/or more aspherical shapes. We discuss that this leads to important consequences on the predicted extinction law of SN dust and on the survival rate of the formed grains as they pass through the reverse shock of the SN. We conclude that a more detailed understanding of the dust formation process and of the microphysical properties of each dust species needs to be achieved before robust prediction on the SN dust yields can be performed.

Primordial stars formed in the early universe are thought to be hosted by compact dark matter (DM) halos. If DM consists of Weakly Interacting Massive Particles (WIMPs), such stars may be powered by DM annihilation during the early phases of their evolutions. We study the pre-main sequence evolutions of the primordial star using a detailed stellar evolution code under the assumption that the annihilation of adiabatically contracted WIMPs DM within the star provides a sufficient energy to sustain the stellar equilibrium. We follow the evolution of accreting stars using several gas mass accretion rates derived from cosmological simulations. We show that the stellar mass becomes very large, up to 900 - 1000 M_sun when the star reaches the main-sequence phase for a reasonable set of model parameters such as DM particle mass and the annihilation cross section. During the dark star phase, the star expands over a thousand solar-radii, while the surface temperature remains below 10^4 K. The energy generated by nuclear reactions is not dominant during this phase. We also study models with different gas mass accretion rates and the DM particle masses. All our models for different DM particle masses pass the dark star phase. The final mass of the dark stars is essentially unchanged for DM mass of m_DM <= 10 GeV. Gravitational collapse of the massive dark stars will leave massive black holes with mass as large as 1000 M_sun in the early universe.

We investigate the formation of dust grains in the ejecta of Type Ia supernovae (SNe Ia), adopting the carbon-deflagration W7 model. In the calculations of dust formation, we apply the nucleation and grain growth theory and consider the two cases with and without formation of CO and SiO molecules. The results of the calculations show that for the sticking probability of alpha_j=1, C, silicate, Si, and FeS grains can condense at early times of ~100--300 days after the explosion, whereas Fe and SiC grains cannot form substantially. Due to the low gas density in SNe Ia with no H-envelope, the average radii of the newly formed grains are generally below 0.01 micron, being much smaller than those in Type II-P SNe. This supports our previous conclusion that the radius of dust formed in the ejecta is smaller in SNe with less massive envelopes. The total dust mass ranges from 3x10^{-4} M_sun to 0.2 M_sun for alpha_j=0.1--1, depending on whether or not CO and SiO molecules are formed. We also estimate the optical depths and thermal emission by the newly formed dust and compare to the relevant observations of SNe Ia. We find that the formation of C grains in SNe Ia is suppressed to be consistent with the observational constraints. This implies that energetic photons and electrons heavily depress the formation efficiency of C grains or that the outermost C-O layer of SNe Ia is almost fully burned. Finally, we perform the calculations of dust destruction in the SN remnants and find that dust grains formed in the ejecta of SNe Ia are almost completely destroyed in the shocked gas before being injected into the interstellar medium. This indicates that SNe Ia are unlikely to be the major sources of interstellar dust.

We develop one-zone galaxy formation models in the early Universe, taking into account dust formation and evolution by supernova (SN) explosions. Especially we focus on the time evolution of dust size distribution, because ${\rm H}_{2}$ formation on the dust surface plays a critical role in the star formation process in the early Universe. In the model we assume that star formation rate (SFR) is proportional to the total amount of ${\rm H}_{2}$. We consistently treat (i) the formation and size evolution of dust, (ii) the chemical reaction networks including ${\rm H}_{2}$ formation both on the surface of dust and in gas phase, and (iii) the SFR in the model. First, we find that, because of dust destruction due to both reverse and forward shocks driven by SNe, H$_{2}$ formation is more suppressed than that without dust destruction. At the galaxy age of $sim0.8 {\rm Gyr}$, for galaxy models with virial mass $M_{\rm vir}=10^{9} M_{\odot}$ and formation redshift $z_{\rm vir}=10$, the molecular fraction is 2.5 orders of magnitude less in the model with dust destruction by both shocks than that in the model without dust destruction. Second, we show that the H$_{2}$ formation rate strongly depends on the ISM density around SN progenitors. The SFR in higher ISM density is lower, since dust destruction by reverse shocks is more effective in higher ISM density. We conclude that not only the amount but also the size distribution of dust being related with the star formation activity strongly affects the evolution of galaxies in the early Universe.

Isotope ratios have opened a new window into the study of the details of stellar evolution, supernovae, and galactic chemical evolution. We present the evolution of the isotope ratios of elemental abundances (from C to Zn) in the solar neighbourhood, bulge, halo, and thick disk, using chemical evolution models with updated yields of Asymptotic Giant Branch (AGB) stars and core-collapse supernovae. The evolutionary history of each element is different owing to the effects of the initial progenitor mass and metallicity on element production. In the bulge and thick disk the star formation timescale is shorter than in the solar neighbourhood, leading to higher [alpha/Fe] ratios. Likewise, the smaller contribution from Type Ia supernovae in these regions leads to lower [Mn/Fe] ratios. Also in the bulge, the abundances of [(Na, Al, P, Cl, K, Sc, Cu, Zn)/Fe] are higher because of the effect of metallicity on element production from core-collapse supernovae. According to our predictions, it is possible to find metal-rich stars ([Fe/H]>-1) that formed in the early Universe as a result of rapid star formation. The chemical enrichment timescale of the halo is longer than in the solar neighbourhood, and consequently the ratios of [(C, F)/Fe] and 12C/13C are higher owing to a significant contribution from low-mass AGB stars. While the [alpha/Fe] and [Mn/Fe] ratios are the same as in the solar neighbourhood, the [(Na, Al, P, Cl, K, Sc, Cu, Zn)/Fe] ratios are predicted to be lower. Furthermore, we predict that isotope ratios such as 24Mg/25,26Mg are larger because of the contribution from low-metallicity supernovae. Using isotopic ratios it is possible to select stars that formed in a system with a low chemical enrichment efficiency such as the satellite galaxies that were accreted onto our own Milky Way Galaxy.

SN 2007bi is an extremely luminous Type Ic supernova. This supernova is thought to be evolved from a very massive star, and two possibilities have been proposed for the explosion mechanism. One possibility is a pair-instability supernova with an M_{CO} ~ 100 M_sun CO core progenitor. Another possibility is a core-collapse supernova with M_{CO} ~ 40 M_sun. We investigate the evolution of very massive stars with main-sequence mass M_{MS} = 100 - 500 M_sun and Z_0 = 0.004, which is in the metallicity range of the host galaxy of SN 2007bi, to constrain the progenitor of SN 2007bi. The supernova type relating to the surface He abundance is also discussed. The main-sequence mass of the progenitor exploding as a pair-instability supernova could be M_{MS} ~ 515 - 575 M_sun. The minimum main-sequence mass could be 310 M_sun when uncertainties in the mass-loss rate are considered. A star with M_{MS} ~ 110 - 280 M_sun evolves to a CO star, appropriate for the core-collapse supernova of SN 2007bi. Arguments based on the probability of pair-instability and core-collapse supernovae favour the hypothesis that SN 2007bi originated from a core-collapse supernova event.

We introduce our newly developed two different, three dimensional magneto hydrodynamical codes in detail. One of our codes is written in the Newtonian limit (NMHD) and the other is in the fully general relativistic code (GRMHD). Both codes employ adaptive mesh refinement and, in GRMHD, the metric is evolved with the "Baumgarte-Shapiro-Shibata-Nakamura" formalism known as the most stable method at present. We did several test problems and as for the first practical test, we calculated gravitational collapse of a $15M_\odot$ star. Main features found from our calculations are; (1) High velocity bipolar outflow is driven from the proto-neutronstar and moves through along the rotational axis in strongly magnetized models; (2) A one-armed spiral structure appears which is originated from the low-$|T/W|$ instability; (3) By comparing GRMHD and NMHD models, the maximum density increases about $\sim30%$ in GRMHD models due to the stronger gravitational effect. These features agree very well with previous studies and our codes are thus reliable to numerical simulation of gravitational collapse of massive stars.

There have been suggestions that the abundance of Extremely Metal-Poor (EMP) stars can be reproduced by Hypernovae (HNe), not by normal supernovae (SNe). However, recently it was also suggested that if the innermost neutron-rich or proton-rich matter is ejected, the abundance patterns of ejected matter are changed, and normal SNe may also reproduce the observations of EMP stars. In this letter, we calculate explosive nucleosynthesis with various Ye and entropy, and investigate whether normal SNe with this innermost matter, which we call "hot-bubble" component, can reproduce the abundance of EMP stars. We find that neutron-rich (Ye = 0.45-0.49) and proton-rich (Ye = 0.51-0.55) matter can increase Zn/Fe and Co/Fe ratios as observed, but tend to overproduce other Fe-peak elements. In addition to it, we find that if slightly proton-rich matter with 0.50 <= Ye < 0.501 with s/kb ~ 15-40 is ejected as much as ~ 0.06 Msolar, even normal SNe can reproduce the abundance of EMP stars, though it requires fine-tuning of Ye. On the other hand, HNe can more easily reproduce the observations of EMP stars without fine-tuning. Our results imply that HNe are the most possible origin of the abundance pattern of EMP stars.

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.

We propose a model in which intermediate-mass black holes (IMBHs) with mass of ~10000 Msun are formed in early dark matter halos. We carry out detailed stellar evolution calculations for accreting primordial stars including annihilation energy of dark matter particles. We follow the stellar core evolution consistently up to gravitational collapse. We show that very massive stars, as massive as 10000 Msun, can be formed in an early dark matter halo. Such stars are extremely bright with Log L/Lsun > 8.2. They gravitationally collapse to form IMBHs. These black holes could have seeded the formation of early super-massive blackholes.

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 calculate the evolution of zero-metallicity Population III (Pop III) stars whose mass grows from the initial mass of $\sim 1M_{\odot}$ by accreting the surrounding gases. Our calculations cover a whole evolutionary stages from the pre-main sequence, via various nuclear burning stages, through the final core collapse or pair-creation instability phases. We adopt the following stellar mass-dependent accretion rates which are derived from cosmological simulations of early structure formation based on the low mass dark matter halos at redshifts $z \sim 20$: (1) the accretion rates for the first generation (Pop III.1) stars and (2) the rates for zero-metallicity but the second generation (Pop III.2) stars which are affected by radiation from the Pop III.1 stars. For comparison, we also study the evolution with the mass-dependent accretion rates which are affected by radiatibe feedback. We show that the final mass of Pop III.1 stars can be as large as $\sim 1000M_{\odot}$, beyond the mass range ($140 - 300M_{\odot}$) for the pair-instability supernovae. Such massive stars undergo core-collapse to form intermediate-mass black holes, which may be the seeds for merger trees to supermassive black holes. On the other hand, Pop III.2 stars become less massive ($\lsim 40 - 60M_{\odot}$), being in the mass range of ordinary iron core-collapse stars. Such stars explode and eject heavy elements to contribute to chemical enrichment of the early universe as observed in the abundance patterns of extremely metal-poor stars in the Galactic halo.

Analyzing the solar system abundances, we have found two empirical abundance scaling laws concerning the p- and s-nuclei with the same atomic number. The first scaling is s/p ratios are almost constant over a wide range of the atomic number, where the p-nculei are lighter than the s-nuclei by two or four neutrons. The second scaling is p/p ratios are almost constant, where the second $p$-nuclei are lighter than the first p-nucleus by two neutrons. These scalings are a piece of evidence that most p-nuclei are dominantly synthesized by the gamma-process in supernova explosions. The scalings lead to a novel concept of "universality of gamma-process" that the s/p and p/p ratios of nuclei produced by individual gamma-processes are almost constant, respectively. We have calculated the ratios by gamma-process based on core-collapse supernova explosion models under various astrophysical conditions and found that the scalings hold for materials produced by individual gamma-processes independent of the astrophysical conditions assumed. The universality originates from three mechanisms: the shifts of the gamma-process layers to keep their peak temperature, the weak s-process in pre-supernovae, and the independence of the s/p ratios of the nuclear reactions. The results further suggest an extended universality that the s/p ratios in the gamma-process layers are not only constant but also centered on a specific value of 3. With this specific value and the first scaling, we estimate that the ratios of $s$-process abundance contributions from the AGB stars to the massive stars are almost 6.7 for the $s$-nuclei of A > 90. We find that large enhancements of s/p ratios for Ce, Er, and W are a piece of evidence that the weak s-process actually occurred before SNe.

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

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.

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}$.

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.

Motivated by the discovery of extremely bright supernovae SNe1999as and 2006gy, we have investigated how much 56Ni mass can be synthesized in core-collapse massive supernovae (SNe). We calculate the evolution of several very massive stars with initial masses M <~ 100 Msun from the main-sequence to the beginning of the Fe-core collapse and simulate their explosions and nucleosynthesis. In order to avoid complications associated with strong mass-loss, we only consider metal-poor stars with initial metallicity Z = Zsun/200. However, our results are applicable to higher metallicity models with similar C+O core masses. The C+O core mass for the 100Msun model is M_CO = 42.6Msun and this is the heaviest model in the literature for which Fe-core collapse SN is explored. The synthesized 56Ni mass increases with the increasing explosion energy and progenitor mass. For the explosion energy of E_51 = E/10^{51} erg =30, for example, the 56Ni masses of M(56Ni) = 2.2, 2.3, 5.0, and 6.6 Msun can be produced for the progenitors with initial masses of 30, 50, 80 and 100 Msun (or C+O core masses M_CO = 11.4, 19.3, 34.0 and 42.6 Msun), respectively. We find that producing M(56Ni) ~ 4Msun as seen in SN1999as is possible for M_CO >~ 34 Msun and E_{51} >~ 20. Producing M(56Ni) ~ 13Msun as suggested for SN2006gy requires a too large explosion energy for M_CO <~ 43Msun, but it may be possible with a reasonable explosion energy if M_CO >~ 60Msun.

We investigate the effects of neutrino-nucleus interactions (the nu-process) on the production of iron-peak elements in Population III core-collapse supernovae. The nu-process and the following proton and neutron capture reactions produce odd-Z iron-peak elements in complete and incomplete Si burning region. This reaction sequence enhances the abundances of Sc, Mn, and Co in the supernova ejecta. The supernova explosion models of 15 M_sol and 25 M_sol stars with the nu-process well reproduce the averaged Mn/Fe ratio observed in extremely metal-poor halo stars. In order to reproduce the observed Mn/Fe ratio, the total neutrino energy in the supernovae should be 3 - 9 x 10^{53} ergs. Stronger neutrino irradiation and other production sites are necessary to reproduce the observed Sc/Fe and Co/Fe ratios, although these ratios increase by the nu-process.

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.

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.

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

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.

We calculate evolution, collapse, explosion, and nucleosynthesis of Population III very-massive stars with 500$M_{\odot}$ and 1000$M_{\odot}$. Presupernova evolution is calculated in spherical symmetry. Collapse and explosion are calculated by a two-dimensional code, based on the bipolar jet models. We compare the results of nucleosynthesis with the abundance patterns of intracluster matter, hot gases in M82, and extremely metal-poor stars in the Galactic halo. It was found that both 500$M_{\odot}$ and 1000$M_{\odot}$ models enter the region of pair-instability but continue to undergo core collapse. In the presupernova stage, silicon burning regions occupy a large fraction, more than 20% of the total mass. For moderately aspherical explosions, the patterns of nucleosynthesis match the observational data of both intracluster medium and M82. Our results suggest that explosions of Population III core-collapse very-massive stars contribute significantly to the chemical evolution of gases in clusters of galaxies. For Galactic halo stars, our [O/Fe] ratios are smaller than the observational abundances. However, our proposed scenario is naturally consistent with this outcome. The final black hole masses are $\sim 230M_{\odot}$ and $\sim 500M_{\odot}$ for the $500M_{\odot}$ and 1000$M_{\odot}$ models, respectively. This result may support the view that Population III very massive stars are responsible for the origin of intermediate mass black holes which were recently reported to be discovered.

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.

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.

We calculate nucleosynthesis in Population (Pop) III supernovae (SNe) and compare the yields with various abundance patterns of extremely metal-poor (EMP) stars. We assume that the observed EMP stars are the second generation stars, which have the metal-abundance patterns of Pop III SNe. Previous theoretical yields of Pop III SNe cannot explain the trends in the abundance ratios among iron-peak elements (Mn, Co, Ni, Zn)/Fe as well as the large C/Fe ratio observed in certain EMP stars with [Fe/H] <~ -2.5. In the present paper, we show that if we introduce higher explosion energies and mixing-fallback in the core-collapse SN models of M ~ 20 - 130 Msun, the above abundance features of both typical and C-rich EMP stars can be much better explained. We suggest that the abundance patterns of the [Fe/H] ~ -2.5 stars correspond to supernova yields with normal explosion energies, while those of the carbon un-enhanced ([C/Fe] < 1) stars with [Fe/H] =~ -4 ~ - 3 correspond to high-energy supernova yields. The abundance patterns of the C-rich ([C/Fe]>~ 2) and low [Fe/H] (=~ -5 \~ -3.5) stars can be explained with the yields of faint SNe that eject little 56Ni as observed in SN1997D. In the supernova-induced star formation model, we can qualitatively explain why the EMP stars formed by the faint or energetic supernovae have lower [Fe/H] than the EMP stars formed by normal supernovae. We also examine how the abundance ratios among iron-peak elements depend on the electron mole fraction Ye, and conclude that a large explosion energy is still needed to realize the large Co/Fe and Zn/Fe ratios observed in typical EMP stars with [Fe/H] <~ -3.5.

With analyzing the solar system abundance, we have found two universal scaling laws concerning the p- and s-nuclei. They indicate that the gamma-process in supernova (SN) explosions is the most promising origin of the p-nuclei that has been discussed with many possible nuclear reactions and sites for about fifty years. In addition the scalings provide new concepts: an universality of the gamma-process and a new nuclear cosmochronometer. We have carried out gamma-process nucleosynthesis calculations for typical core-collapse SN explosion models and the results satisfy the observed scalings.

We investigate the formation of dust grains in the ejecta of population III supernovae including pair--instability supernovae, applying a theory of non-- steady state nucleation and grain growth. In the calculations, the time evolution of gas temperature in theejecta, which strongly affects the number density and size of newly formed grains, is calculated by solving the radiative transfer equation taking account of the energy deposition of radio active elements. Two extreme cases are considered for the mixing of elements in the ejecta; unmixed and uniformly mixed cases within the He--core. The results of calculations are summarized as the followings; in the unmixed ejecta, a variety of grain species condense, reflecting the difference of the elemental composition at the formation site in the ejecta, otherwise only oxide grains condense in the uniformly mixed ejecta. The average size of newly formed grains spans the range of three orders of magnitude, depending on the grain species and the formation condition, and the maximum radius is limited to less than 1 $\mu$m, which does not depend on the progenitor mass. The size distribution function summed up over all grain species is approximated by a power--law formula whose index is -3.5 for the larger radius and -2.5 for the smaller one; the radius at the crossover point ranges from 0.004 to 0.1 $\mu$m, depending on the model of supernovae. The fraction of mass locked into dust grains increases with increasing the progenitor mass; 2--5 % of the progenitor mass for core collapse supernovae and 15--30 % for pair--instability supernovae whose progenitor mass ranges from 140 to 260 $M_{\odot}$. Thus, if the very massive stars populate the first generation stars, a large amount of dust grains would be produced in the early universe.

A key question for supernova cosmology is whether the peak luminosities of Type Ia supernovae (SNe Ia) are sufficiently free from the effects of cosmic and galactic evolution. To answer this question, we review the currently popular scenario of SN Ia progenitors, i.e., the single degenerate scenario for the Chandrasekhar mass white dwarf (WD) models. We identify the progenitor's evolution with two channels: (1) the WD+RG (red-giant) and (2) the WD+MS (near main-sequence He-rich star) channels. The strong wind from accreting WDs plays a key role, which yields important age and metallicity effects on the evolution. We suggest that the variation of the carbon mass fraction $X$(C) in the C+O WD (or the variation of the initial WD mass) causes the diversity of SN Ia brightness. This model can explain the observed dependence of SNe Ia brightness on the galaxy types. We then predict how SN Ia brightness evolves along the redshift (with changing metallicity and age) for elliptical and spiral galaxies. Such evolutionary effects along the redshift can be corrected as has been made for local SNe Ia. We also touch on several related issues: (1) the abundance pattern of stars in dwarf spheroidal galaxies in relation to the metallicity effect on SNe Ia, (2) effects of angular momentum brought into the WD in relation to the diversities and the fate of double degenerates, and (3) possible presence of helium in the peculiar SN Ia 2000cx in relation to the sub-Chandrasekhar mass model.

During the last few years, a number of exceptional core-collapse supernovae (SNe) have been discovered. Their kinetic energy of the explosions are larger by more than an order of magnitude than the typical values for this type of SNe, so that these SNe have been called `Hypernovae'. We first describe how the basic properties of hypernovae can be derived from observations and modeling. These hypernovae seem to come from rather massive stars, thus forming black holes. On the other hand, there are some examples of massive SNe with only a small kinetic energy. We suggest that 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. Hypernova nucleosynthesis is characterized by larger abundance ratios (Zn,Co,V,Ti)/Fe and smaller (Mn,Cr)/Fe. Nucleosynthesis in Faint supernovae is characterized by a large amount of fall-back. We show that the abundance pattern of the most Fe deficient 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.

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.

It has been proposed theoretically that the first generation of stars in the Universe (population III) would be as massive as 100 solar masses (100Mo), because of inefficient cooling of the precursor gas clouds. Recently, the most iron-deficient (but still carbon-rich) low-mass star -- HE0107-5240 -- was discovered. If this is a population III that gained its metals (elements heavier than helium) after its formation, it would challenge the theoretical picture of the formation of the first stars. Here we report that the patterns of elemental abundance in HE0107-5240 (and other extremely metal-poor stars) are in good accord with the nucleosynthesis that occurs in stars with masses of 20-130Mo when they become supernovae if, during the explosions, the ejecta undergo substantial mixing and fall-back to form massive black holes. Such supernovae have been observed. The abundance patterns are not, however, consistent with enrichment by supernovae from stars in the range 130-300 Mo. We accordingly infer that the first-generation supernovae came mostly from explosions of ~ 20-130Mo stars; some of these produced iron-poor but carbon- and oxygen-rich ejecta. Low-mass second-generation stars, like HE0107-5240, could form because the carbon and oxygen provided pathways for gas to cool.

Supernova (SN) 2002ap in M74 was observed in the $UBVRIJHK$ bands for the first 40 days following its discovery (2002 January 29) until it disappeared because of solar conjunction, and then in June after it reappeared. The magnitudes and dates of peak brightness in each band were determined. While the rate of increase of the brightness before the peak is almost independent of wavelength, the subsequent rate of decrease becomes smaller with wavelength from the $U$ to the $R$ band, and is constant at wavelengths beyond $I$. The photometric evolution is faster than in the well-known ``hypernovae'' SNe~1998bw and 1997ef, indicating that SN 2002ap ejected less mass. The bolometric light curve of SN 2002ap for the full period of observations was constructed. The absolute magnitude is found to be much fainter than that of SN 1998bw, but is similar to that of SN 1997ef, which lies at the faint end of the hypernova population. The bolometric light curve at the early epochs was best reproduced with the explosion of a C+O star that ejects $2.5~M_\sun$ with kinetic energy $E_{\rm K}=4\times 10^{51}~{\rm ergs}$. A comparison of the predicted brightness of SN 2002ap with that observed after solar conjunction may imply that $\gamma$-ray deposition at the later epochs was more efficient than in the model. This may be due to an asymmetric explosion.

We review the characteristics of nucleosynthesis in 'Hypernovae', i.e., core-collapse supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs). The hypernova yields show the following characteristics: 1) The mass ratio between the complete and incomplete Si burning regions is larger in hypernovae than normal supernovae. As a result, higher energy explosions tend to produce larger [(Zn, Co, V)/Fe] and smaller [(Mn, Cr)/Fe], which could explain the trend observed in very metal-poor stars. 2) Because of enhanced $\alpha$-rich freezeout, $^{44}$Ca, $^{48}$Ti, and $^{64}$Zn are produced more abundantly than in normal supernovae. The large [(Ti, Zn)/Fe] ratios observed in very metal poor stars strongly suggest a significant contribution of hypernovae. 3) Oxygen burning takes place in more extended regions in hypernovae to synthesize a larger amount of Si, S, Ar, and Ca ("Si"), which makes the "Si"/O ratio larger. The abundance pattern of the starburst galaxy M82 may be attributed to hypernova explosions. We thus suggest that hypernovae make important contribution to the early Galactic (and cosmic) chemical evolution.

X-ray observations have shown that the chemical abundance in the starburst galaxy M82 is quite rich in Si and S compared with oxygen. Such an abundance pattern cannot be explained with any combination of conventional Type I and II supernova yields. Also the energy to heavy element mass ratio of the observed hot plasma is much higher than the value resulted from normal supernovae. We calculate explosive nucleosynthesis in core-collapse hypernovae and show that the abundance pattern and the large ratio between the energy and the heavy element mass can be explained with the hypernova nucleosynthesis. Such hypernova explosions are expected to occur for stars more massive than >~ 20-25 Msun, and likely dominating the starburst, because the age after the starburst in M82 is estimated to be as short as ~ 10^6 - 10^7 yr. We also investigate pair-instability supernovae (~ 150-300 Msun) and conclude that the energy to heavy element mass ratio in these supernovae is too small to explain the observation.

We calculate evolution and nucleosynthesis in massive Pop III stars with M = 13 ~ 270Msun, and compare the results with abundances of very metal-poor halo stars. The observed abundances can be explained by the energetic core-collapse supernovae with M < ~ 130Msun (``hypernovae'') but not by pair-instability supernovae (PISNe) with M ~ 140-270 Msun. This result constrain the IMF for the Pop III and very metal-poor Pop II stars.

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are presented, and the properties of the SN are investigated through models of its spectral evolution and its light curve. The SN is spectroscopically similar to the "hypernova" SN 1997ef. However, its kinetic energy [$\sim (4-10) \times 10^{51}$ erg] and the mass ejected (2.5-5 $M_{\odot}$) are smaller, resulting in a faster-evolving light curve. The SN synthesized $\sim 0.07 M_{\odot}$ of $^{56}$Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone should not be used to define a hypernova, whose defining character, namely very broad spectral features, is the result of a high kinetic energy. The likely main-sequence mass of the progenitor star was 20-25 $M_{\odot}$, which is also lower than that of both hypernovae SNe 1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the hypernova sequence as it is known so far. Observations of the nebular spectrum, which is expected to dominate by summer 2002, are necessary to confirm these values.

We calculate nucleosynthesis in core-collapse explosions of massive Pop III stars, and compare the results with abundances of metal-poor halo stars to constrain the parameters of Pop III supernovae. We focus on iron-peak elements and, in particular, we try to reproduce the large [Zn/Fe] observed in extremely metal-poor stars. The interesting trends of the observed ratios [Zn, Co, Mn, Cr, V/Fe] can be related to the variation of the relative mass of the complete and incomplete Si-burning regions in supernova ejecta. We find that [Zn/Fe] is larger for deeper mass-cuts, smaller neutron excess, and larger explosion energies. The large [Zn/Fe] and [O/Fe] observed in the very metal-poor halo stars suggest deep mixing of complete Si-burning material and a significant amount of fall-back in Type II supernovae. Furthermore, large explosion energies (E_51 >~ 2 for M ~ 13 Msun and E_51 >~ 20 for M >~ 20 Msun) are required to reproduce [Zn/Fe] ~ 0.5. The observed trends of the abundance ratios among the iron-peak elements are better explained with this high energy (``Hypernova'') models rather than the simple ``deep'' mass-cut effect, because the overabundance of Ni can be avoided in the hypernova models. We also present the yields of pair-instability supernova explosions of M = 130 - 300 Msun stars, and discuss that the abundance features of very metal-poor stars cannot be explained by pair-instability supernovae.

We review the characteristics of nucleosynthesis in 'Hypernovae', i.e., supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs). The hypernova yields compared to those of ordinary core-collapse supernovae show the following characteristics: 1) Complete Si-burning takes place in more extended region, so that the mass ratio between the complete and incomplete Si burning regions is generally larger in hypernovae than normal supernovae. As a result, higher energy explosions tend to produce larger [(Zn, Co)/Fe], smaller [(Mn, Cr)/Fe], and larger [Fe/O], which could explain the trend observed in very metal-poor stars. 2) Si-burning takes place in lower density regions, so that the effects of $\alpha$-rich freezeout is enhanced. Thus $^{44}$Ca, $^{48}$Ti, and $^{64}$Zn are produced more abundantly than in normal supernovae. The large [(Ti, Zn)/Fe] ratios observed in very metal-poor stars strongly suggest a significant contribution of hypernovae. 3) Oxygen burning also takes place in more extended regions for the larger explosion energy. Then a larger amount of Si, S, Ar, and Ca ("Si") are synthesized, which makes the "Si"/O ratio larger. The abundance pattern of the starburst galaxy M82 may be attributed to hypernova explosions. Asphericity in the explosions strengthens the nucleosynthesis properties of hypernovae except for "Si"/O. We thus suggest that hypernovae make important contribution to the early Galactic (and cosmic) chemical evolution.

We review the characteristics of nucleosynthesis and radioactivities in 'Hypernovae', i.e., supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs) and their $\gamma$-ray line signatures. We also discuss the $^{44}$Ti line $\gamma$-rays from SN1987A and the detectability with INTEGRAL. Signatures of hypernova nucleosynthesis are seen in the large [(Ti, Zn)/Fe] ratios in very metal poor stars. Radioactivities in hypernovae compared to those of ordinary core-collapse supernovae show the following characteristics: 1) The complete Si burning region is more extended, so that the ejected mass of $^{56}$Ni can be much larger. 2) Si-burning takes place in higher entropy and more $\alpha$-rich environment. Thus the $^{44}$Ti abundance relative to $^{56}$Ni is much larger. In aspherical explosions, $^{44}$Ti is even more abundant and ejected with velocities as high as $\sim$ 15,000 km s$^{-1}$, which could be observed in $\gamma$-ray line profiles. 3) The abundance of $^{26}$Al is not so sensitive to the explosion energy, while the $^{60}$Fe abundance is enhanced by a factor of $\sim$ 3.

We study nucleosynthesis in 'hypernovae', i.e., supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs) for both spherical and aspherical explosions. The hypernova yields compared to those of ordinary core-collapse supernovae show the following characteristics: 1) Complete Si-burning takes place in more extended region, so that the mass ratio between the complete and incomplete Si burning regions is generally larger in hypernovae than normal supernovae. As a result, higher energy explosions tend to produce larger [(Zn, Co)/Fe], small [(Mn, Cr)/Fe], and larger [Fe/O], which could explain the trend observed in very metal-poor stars. 2) Si-burning takes place in lower density regions, so that the effects of $\alpha$-rich freezeout is enhanced. Thus $^{44}$Ca, $^{48}$Ti, and $^{64}$Zn are produced more abundantly than in normal supernovae. The large [(Ti, Zn)/Fe] ratios observed in very metal poor stars strongly suggest a significant contribution of hypernovae. 3) Oxygen burning also takes place in more extended regions for the larger explosion energy. Then a larger amount of Si, S, Ar, and Ca ("Si") are synthesized, which makes the "Si"/O ratio larger. The abundance pattern of the starburst galaxy M82 may be attributed to hypernova explosions. Asphericity in the explosions strengthens the nucleosynthesis properties of hypernovae except for "Si"/O. We thus suggest that hypernovae make important contribution to the early Galactic (and cosmic) chemical evolution.

We examine the characteristics of nucleosynthesis in 'hypernovae', i.e., supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs). We carry out detailed nucleosynthesis calculations for these energetic explosions and compare the yields with those of ordinary core-collapse supernovae. We find that both complete and incomplete Si-burning takes place over more extended, lower density regions, so that the alpha-rich freezeout is enhanced and produces more Ti in comparison with ordinary supernova nucleosynthesis. In addition, oxygen and carbon burning takes place in more extended, lower density regions than in ordinary supernovae. Therefore, the fuel elements O, C, Al are less abundant while a larger amount of Si, S, Ar, and Ca ("Si") are synthesized by oxygen burning; this leads to larger ratios of "Si"/O in the ejecta. Enhancement of the mass ratio between complete and incomplete Si-burning regions in the ejecta may explain the abundance ratios among iron-peak elements in metal-poor stars. Also the enhanced "Si"/O ratio may explain the abundance ratios observed in star burst galaxies. We also discuss other implications of enhanced [Ti/Fe] and [Fe/O] for Galactic chemical evolution and the abundances of low mass black hole binaries.

We examine the characteristics of nucleosynthesis in 'hypernovae', i.e., supernovae with very large explosion energies ($\gsim 10^{52}$ ergs). Implications for the cosmic chemical evolution and the abundances in M82 are discussed.

The use of Type Ia supernovae as distance indicators relies on the determination of their brightness. This is not constant, but it can be calibrated using an observed relation between the brightness and the properties of the optical light curve (decline rate, width, shape), which indicates that brighter SNe have broader, slower light curves. However, the physical basis for this relation is not yet fully understood. Among possible causes are different masses of the progenitor white dwarfs or different opacities in Chandrasekhar-mass explosions. We parametrise the Chandrasekhar-mass models presented by Iwamoto et al (1999), which synthesize different amounts of Ni, and compute bolometric light curves and spectra at various epochs. Since opacity in SNe Ia is due mostly to spectral lines, it should depend on the mass of Fe-peak elements synthesized in the explosion, and on the temperature in the ejecta. Bolometric light curves computed using these prescriptions for the optical opacity reproduce the relation between brightness and decline rate. Furthermore, when spectra are calculated, the change in colour between maximum and two weeks later allows the observed relation between M_B(Max) and Dm_{15}(B) to be reproduced quite nicely. Spectra computed at various epochs compare well with corresponding spectra of spectroscopically normal SNeIa selected to cover a similar range of Dm_{15}(B) values.

Very recently the Chandra First Light Observation discovered a point-like source in Cassiopeia A (Cas A) supernova remnant. This detection was subsequently confirmed by the analyses of the archival data from both ROSAT and Einstein observations. Here we compare the results from these observations with the scenarios involving both black holes (BH) and neutron stars (NS). If this point source is a BH we offer as a promising model a disk-corona type model with low accretion rate where a soft photon source at ~ 0.1 keV is Comptonized by higher energy electrons in the corona. If it is a NS the dominant radiation observed by Chandra most likely originates from smaller, hotter regions of the stellar surface, but we argue that it is still worthwhile to compare the cooler component from the rest of the surface with cooling theories. We emphasize that the detection of this point source itself should potentially provide enormous impacts on the theories of supernova explosion, progenitor scenario, compact remnant formation, accretion to compact objects, and neutron star thermal evolution.

Relatively uniform light curves and spectral evolution of Type Ia supernovae (SNe Ia) have led to the use of SNe Ia as a ``standard candle'' to determine cosmological parameters. Whether a statistically significant value of the cosmological constant can be obtained depends on whether the peak luminosities of SNe Ia are sufficiently free from the effects of cosmic and galactic evolutions. Here we first review the single degenerate scenario for the Chandrasekhar mass white dwarf (WD) models of SNe Ia. We identify the progenitor's evolution and population with two channels: (1) the WD+RG (red-giant) and (2) the WD+MS (near main-sequence He-rich star) channels. In these channels, the strong wind from accreting WDs plays a key role, which yields important age and metallicity effects on the evolution. We then address the questions whether the nature of SNe Ia depends systematically on environmental properties such as metallicity and age of the progenitor system and whether significant evolutionary effects exist. We suggest that the variation of the carbon mass fraction $X$(C) in the C+O WD (or the variation of the initial WD mass) causes the diversity of the brightness of SNe Ia. This model can explain the observed dependences of SNe Ia brightness on the galaxy types and the distance from the galactic center. Finally, applying the metallicity effect on the evolution of SN Ia progenitors, we make a prediction of the cosmic supernova rate history as a composite of the supernova rates in different types of galaxies.

We calculate explosive nucleosynthesis in Chandrasekhar mass models for Type Ia Supernovae(SNe Ia) to obtain new constraints on the rate of matter accretion onto the progenitor white dwarf and on the ignition density of central carbon deflagration. The calculated abundance of the Fe-group neutron-rich nuclei is highly sensitive to the electron captures taking place in the central layers. The yields obtained from a slow central deflagration, and from a fast deflagration or delayed detonation in the outer layers, are combined and put to comparison with solar isotopic abundances. We found that (1) to avoid too large ratios of $^{54}$Cr/$^{56}$Fe and $^{50}$Ti/$^{56}$Fe, the ignition density should be as low as \ltsim 2 \e9 \gmc, and that (2) to avoid the overproduction of $^{58}$Ni and $^{54}$Fe, either the flame speed should not exceed a few % of the sound speed in the central low $Y_e$ layers, or the progenitor star has to be metal-poor compared with solar. Such low central densities can be realized by a rapid accretion as fast as $\dot M$ \gtsim 1 $\times$ 10$^{-7}$M$_\odot$ yr$^{-1}$. In order to reproduce the solar abundance of $^{48}$Ca, one also needs progenitor systems that undergo ignition at higher densities. We also found that the total amount of $^{56}$Ni, the Si-Ca/Fe ratio, and the abundance of elements like Mn and Cr (incomplete Si-burning ashes), depend on the density of the deflagration-detonation transition in delayed detonations. Our nucleosynthesis results favor transition densities slightly below 2.2$\times 10^7$ g cm$^{-3}$.

The Chandrasekhar mass model for Type Ia Supernovae (SNe Ia) has received increasing support from recent comparisons of observations with light curve predictions and modeling of synthetic spectra. It explains SN Ia events via thermonuclear explosions of accreting white dwarfs in binary stellar systems, being caused by central carbon ignition when the white dwarf approaches the Chandrasekhar mass. As the electron gas in white dwarfs is degenerate, characterized by high Fermi energies for the high density regions in the center, electron capture on intermediate mass and Fe-group nuclei plays an important role in explosive burning. Electron capture affects the central electron fraction Y_e, which determines the composition of the ejecta from such explosions. Up to the present, astrophysical tabulations based on shell model matrix elements were only available for light nuclei in the sd-shell. Recently new Shell Model Monte Carlo (SMMC) and large-scale shell model diagonalization calculations have also been performed for pf-shell nuclei. These lead in general to a reduction of electron capture rates in comparison with previous, more phenomenological, approaches. Making use of these new shell model based rates, we present the first results for the composition of Fe-group nuclei produced in the central regions of SNe Ia and possible changes in the constraints on model parameters like ignition densities and burning front speeds.

We calculate presupernova evolutions and supernova explosions of massive stars (M=13-25 Mo) for various metallicities. We find the following characteristic abundance patterns of nucleosynthesis in the metal-free (Pop III) stars. (1) The alpha-nuclei (from C to Zn) are more efficiently produced than other isotopes, and the abundance pattern of alpha-nuclei can be similar to the solar abundance. In particular, near solar ratios of alpha elements/Fe might be a signature of Pop III which could produce a large amount of Fe. (2) The abundance ratios of odd Z to even Z elements such as Na/Mg and Al/Mg become smaller for lower metallicity. However, these ratios almost saturate below Z <~ 10^{-5}, and [Na, Al/Mg] ~ - 1 for Pop III and low metal Pop II nucleosynthesis. This result is consistent with abundance pattern of metal poor stars, in which these ratios also saturate around -1. We suggest that these stars with the lowest [Na/Mg] or [Al/Mg] may contain the abundance pattern of Pop III nucleosynthesis. (3) Metal poor stars show interesting trends in the ratios of [Cr, Mn, Co/Fe]. We discuss that these trends are not explained by the differences in metallicity, but by the relative thickness between the complete and the incomplete Si burning layers. Large [Co/Fe] and small [Cr, Mn/Fe] values found in the observations are explained if mass cut is deep or if matter is ejected from complete Si burning layer in a form of a jet or bullets. (4) We also find that primary ^{14}N production occurs in the massive Pop III stars, because these stars have radiative H-rich envelopes so that the convective layer in the He-shell burning region can reach the H-rich region.

It is generally believed that cosmological Gamma Ray Bursts (GRBs) are produced by the deceleration of relativistic objects with Lorentz factor (Gamma) >~ 100. We study the possibility that some GRBs are produced along with relativistic matter ejection from supernovae. In this model, it is quite likely that the matter has to travel through the progenitor's thick envelope before generating GRBs. Under the assumption that the ejected matter is described as a single collective matter, we obtain constraints on the matter to have Gamma >~ 100 at the breakout of the progenitor. One advantage of considering this type of model is that the expected GRB energy is sufficiently large, in contrast to the GRB generation model by the shock breakout in the energetic supernova explosion. We find that in general the cross section of the matter has to be very small compared with the progenitor's radius and thus the matter has to be bullet (or jet)-like rather than shell-like.

The influence of the metallicity at the main sequence on the chemical structure of the exploding white dwarf, the nucleosynthesis during the explosion and the light curves of an individual Type Ia supernovae have been studied. Detailed calculations of the stellar evolution, the explosion, and light curves of delayed detonation models are presented. Detailed stellar evolution calculations with a main sequence mass of 7 M_o have been performed to test the influence of the metallicity Z on the structure of the progenitor. A change of Z influences the central helium burning and, consequently, the size of the C/O core which becomes a C/O white dwarf and its C/O ratio. Subsequently, the white dwarf may grow to the Chandrasekhar mass and explode as a Type Ia supernovae. Consequently, the C/O structure of the exploding white dwarf depends on Z. Since C and O are the fuel for the thermonuclear explosion, Z indirectly changes the energetics of the explosion. In our example, changing Z from Pop I to Pop II causes a change in the zero point of the maximum brightness/decline relation by about $0.1^m$ and a change in the rise time by about 1 day. Combined with previous studies, We provide a relation between a change of the rise time in SNeIa and the off-set in the brightness/decline relation. We find that dM = dt where dM is the offset in the relation and dt is the change of the rise time in days. This relation provides a way to detect and test for evolutionary effects, and to estimate the possible consequences for the determination of the cosmological constant. However, due to uncertainties in the time evolution of the properties of the progenitors, we cannot predict whether evolution becomes important at z=0.5, 1 or 2. This question has to be answered by observations.

Relatively uniform light curves and spectral evolution of Type Ia supernovae (SNe Ia) have led to the use of SNe Ia as a ``standard candle'' to determine cosmological parameters, such as the Hubble constant, the density parameter, and the cosmological constant. Whether a statistically significant value of the cosmological constant can be obtained depends on whether the peak luminosities of SNe Ia are sufficiently free from the effects of cosmic and galactic evolutions. Here we first review the single degenerate scenario for the Chandrasekhar mass white dwarf (WD) models of SNe Ia. We identify the progenitor's evolution and population with two channels: (1) the WD+RG (red-giant) and (2) the WD+MS (near main-sequence He-rich star) channels. In these channels, the strong wind from accreting white dwarfs plays a key role, which yields important age and metallicity effects on the evolution. We then address the questions whether the nature of SNe Ia depends systematically on environmental properties such as metallicity and age of the progenitor system and whether significant evolutionary effects exist. We suggest that the variation of the carbon mass fraction $X$(C) in the C+O WD (or the variation of the initial WD mass) causes the diversity of the brightness of SNe Ia. This model can explain the observed dependence of SNe Ia brighness on the galaxy types. Finally, applying the metallicity effect on the evolution of SN Ia progenitors, we make a prediction of the cosmic supernova rate history as a composite of the supernova rates in different types of galaxies.

Observations suggest that the properties of Type Ia supernovae (SNe Ia) may depend on environmental characteristics, such as morphology, metallicity, and age of host galaxies. The influence of these environmental properties on the resulting SNe Ia is studied in this paper. First it is shown that the carbon mass fraction X(C) in the C+O white dwarf SN Ia progenitors tends to be smaller for lower metallicity and older the binary system age. It is then suggested that the variation of X(C) causes the diversity in the brightness of SNe Ia: a smaller X(C) leads to a dimmer SN Ia. Further studies of the propagation of the turbulent flame are necessary to confirm this relation. Our model for the SN Ia progenitors then predicts that when the progenitors belong to an older population or to a low metallicity environment, the number of bright SNe Ia is reduced, so that the variation in brightness among the SNe Ia is also smaller. Thus our model can explain why the mean SN Ia brightness and its dispersion depend on the morphology of the host galaxies and on the distance of the SN from the center of the galaxy. It is further predicted that at higher redshift (z >~ 1) both the the mean brightness of SNe Ia and its variation should be smaller in spiral galaxies than in elliptical galaxies. These variations are within the range observed in nearby SNe Ia. In so far as the variation in X(C) is the most important cause for the diversity among SNe Ia, the light curve shape method currently used to determine the absolute magnitude of SNe Ia can be applied also to high redshift SNe Ia.

We have found a new evolutionary path to Type Ia supernovae (SNe Ia) which has been overlooked in previous work. In this scenario, a carbon-oxygen white dwarf (C+O WD) is originated, not from an asymptotic giant branch star with a C+O core, but from a red-giant star with a helium core of $\sim 0.8-2.0 M_\odot$. The helium star, which is formed after the first common envelope evolution, evolves to form a C+O WD of $\sim 0.8-1.1 M_\odot$ with transferring a part of the helium envelope onto the secondary main-sequence star. This new evolutionary path, together with the optically thick wind from mass-accreting white dwarf, provides a much wider channel to SNe Ia than previous scenarios. A part of the progenitor systems are identified as the luminous supersoft X-ray sources or the recurrent novae like U Sco, which are characterized by the accretion of helium-rich matter. The white dwarf accretes hydrogen-rich, helium-enhanced matter from a lobe-filling, slightly evolved companion at a critical rate and blows excess matter in the wind. The white dwarf grows in mass to the Chandrasekhar mass limit and explodes as an SN Ia. A theoretical estimate indicates that this channel contributes a considerable part of the inferred rate of SNe Ia in our Galaxy, i.e., the rate is about ten times larger than the previous theoretical estimates for white dwarfs with slightly evolved companions.

We explore the effects on nucleosynthesis in Type II supernovae of various parameters (mass cut, neutron excess, explosion energy, progenitor mass) in order to explain the observed trends of the iron-peak element abundance ratios ([Cr/Fe], [Mn/Fe], [Co/Fe] and [Ni/Fe]) in halo stars as a function of metallicity for the range $ -4 \le$ [Fe/H] $\le -2.5$. [Cr/Fe] and [Mn/Fe] decrease with decreasing [Fe/H], while [Co/Fe] behaves the opposite way and increases. We show that such a behavior can be explained by a variation of mass cuts in Type II supernovae as a function of progenitor mass, which provides a changing mix of nucleosynthesis from an alpha-rich freeze-out of Si-burning and incomplete Si-burning. This explanation is consistent with the amount of ejected $^{56}$Ni determined from modeling the early light curves of individual supernovae. We also suggest that the ratio [H/Fe] of halo stars is mainly determined by the mass of interstellar hydrogen mixed with the ejecta of a single supernova which is larger for larger explosion energy and the larger Str\"omgren radius of the progenitor.

Recent observations have revealed that Type Ia supernovae (SNe Ia) are not perfect standard candles but show some variations in their absolute magnitudes, light curve shapes, and spectra. The C/O ratio in the SNe Ia progenitors (C-O white dwarfs) may be related to this variation. In this work, we systematically investigate the effects of stellar mass (M) and metallicity (Z) on the C/O ratio and its distribution in the C-O white dwarfs by calculating stellar evolution from the main-sequence through the end of the second dredge-up for M=3-9 Mo and Z=0.001-0.03. We find that the total carbon mass fraction just before SN Ia explosion varies in the range 0.36 -- 0.5. We also calculate the metallicity dependence of the main-sequence-mass range of the SN Ia progenitor white dwarfs. Our results show that the maximum main-sequence mass to form C-O white dwarfs decreases significantly toward lower metallicity, and the number of SN Ia progenitors may be underestimated if metallicity effectis neglected. We discuss the implications of these results on the variation of SNe Ia, determination of cosmological parameters, luminosity function of white dwarfs, and the galactic chemical evolution.

The discovery of the peculiar supernova (SN) 1998bw and its possible association with the gamma-ray burst (GRB) 980425$^{1,2,3}$ provide new clues to the understanding of the explosion mechanism of very massive stars and to the origin of some classes of gamma-ray bursts. Its spectra indicate that SN~1998bw is a type Ic supernova$^{3,4}$, but its peak luminosity is unusually high compared with typical type Ic supernovae$^3$. Here we report our findings that the optical spectra and the light curve of SN 1998bw can be well reproduced by an extremely energetic explosion of a massive carbon+oxygen (C+O) star. The kinetic energy is as large as $\sim 2-5 \times 10^{52}$ ergs, more than ten times the previously known energy of supernovae. For this reason, the explosion may be called a `hypernova'. Such a C+O star is the stripped core of a very massive star that has lost its H and He envelopes. The extremely large energy, suggesting the existence of a new mechanism of massive star explosion, can cause a relativistic shock that may be linked to the gamma-ray burst.

The thermal evolution of a neutron star is studied by including the energy loss due to axion emission. Two axion models and three types of neutron-star matter equation of state are used with the effects of nucleon superfluidity properly taken into account. In comparison with the observational data of PSR0656+14 from ROSAT, the upper limits on the axion mass are found to be $m_a < 0.06 - 0.3$ eV and 0.08$ -$ 0.8 eV for the KSVZ and DFSZ axion models, respectively, with the soft equation of state giving the most stringent limits.

The discovery of X-ray, optical and radio afterglows of gamma-ray bursts (GRBs) and the measurements of the distances to some of them have established that these events come from Gpc distances and are the most powerful photon emitters known in the Universe, with peak luminosities up to 10^52 erg/s. We here report the discovery of an optical transient, in the BeppoSAX Wide Field Camera error box of GRB980425, which occurred within about a day of the gamma-ray burst. Its optical light curve, spectrum and location in a spiral arm of the galaxy ESO 184-G82, at a redshift z = 0.0085, show that the transient is a very luminous type Ic supernova, SN1998bw. The peculiar nature of SN1998bw is emphasized by its extraordinary radio properties which require that the radio emitter expand at relativistical speed. Since SN1998bw is very different from all previously observed afterglows of GRBs, our discovery raises the possibility that very different mechanisms may give rise to GRBs, which differ little in their gamma-ray properties.

It is tempting to inflate along one of the many flat directions that arise in supersymmetric theories. The required flatness of the potential to obtain sufficient inflation and to not overproduce density fluctuations occurs naturally. However, the density perturbations (in the case of a single moduli field) that arise from inflaton quantum fluctuations are too small for structure formation. Here we propose that topological defects (such as cosmic strings), which arise during a phase transition near the end of moduli inflation can provide an alternative source of structure. The strings produced will be `fat', yet have the usual evolution by the time of nucleosynthesis. Possible models are discussed.