We consider black hole formation in failed supernovae when a dense circumstellar medium (CSM) is present around the massive star progenitor. By utilizing radiation hydrodynamical simulations, we calculate the mass ejection of blue supergiants and Wolf-Rayet stars in the collapsing phase and the radiative shock occurring between the ejecta and the ambient CSM. We find that the resultant emission is redder and dimmer than normal supernovae (bolometric luminosity of $\sim 10^{40}-10^{41}\ {\rm erg\ s^{-1}}$, effective temperature of $\sim 5\times 10^3$ K, and timescale of 10-100 days) and shows a characteristic power-law decay, which may comprise a fraction of intermediate luminosity red transients (ILRTs) including AT 2017be. In addition to searching for the progenitor star in the archival data, we encourage X-ray follow-up observations of such ILRTs $\sim$ 1-10 yr after the collapse, targeting the fallback accretion disk.

Some massive stars experience episodic and intense mass loss phases with fluctuations in the luminosity. Ejected material forms circumstellar matter (CSM) around the star, and the subsequent core-collapse results in a type IIn supernova which is characterized by interaction between supernova ejecta and the CSM. The energy source which triggers these mass eruptions and dynamics of the outflow have not been clearly explained. Moreover, the mass eruption itself can alter the density structure of the envelope and affect the dynamics of the subsequent mass eruption if these events repeat. In fact, a large amount of observational evidence suggests multiple mass eruptions prior to core-collapse. We investigate the density structure of the envelope altered by the first mass eruption and the nature of the subsequent second mass eruption event in comparison with the first one. We deposited an extra energy twice at the bottom of the hydrogen envelope and calculate the time evolution by our radiation hydrodynamical simulation code. We do not deal with the origin of the energy source but focus on the dynamics of repeated mass eruptions from a single massive star. There are significant differences between the first and second mass eruptions in terms of the luminosity, color, amount of produced CSM. The second eruption leads to a redder burst event with the associated brightening phase lasting longer than the first one. The amount of ejected matter is different even with the same deposited energy in the first and second event, but the difference depends on the density structure of the progenitor star. Upcoming high cadence and deep transient surveys will provide us a lot of detailed pre-supernova activities, and some of them would show multi-peaked light curves. They should be interpreted taking the effect of density structure altered by the preceding outburst events into consideration.

Observations suggest that some massive stars experience violent and eruptive mass loss associated with significant brightening that cannot be explained by hydrostatic stellar models. This event seemingly forms dense circumstellar matter (CSM). The mechanism of eruptive mass loss has not been fully explained. We focus on the fact that the timescale of nuclear burning gets shorter than the dynamical timescale of the envelope a few years before core collapse for some massive stars. To reveal the properties of the eruptive mass loss, we investigate its relation to the energy injection at the bottom of the envelope supplied by nuclear burning taking place inside the core. In this study, we do not specify the actual mechanism for transporting energy from the site of nuclear burning to the bottom of the envelope. Instead, we parameterize the amount of injected energy and the injection time and try to extract information on these parameters from comparisons with observations. We carried out 1-D radiation hydrodynamical simulations for progenitors of red, yellow, and blue supergiants, and Wolf-Rayet stars. We calculated the evolution of the progenitors with a public stellar evolution code. We obtain the light curve associated with the eruption, the amount of ejected mass, and the CSM distribution at the time of core-collapse. The energy injection at the bottom of the envelope of a massive star within a period shorter than the dynamical timescale of the envelope could reproduce some observed optical outbursts prior to the core-collapse and form the CSM, which can power an interaction supernova (SN) classified as type IIn.

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

To determine the epoch of reionization precisely and to reveal the property of inhomogeneous reionization are some of the most important topics of modern cosmology. Existing methods to investigate reionization which use cosmic microwave background, Ly$\alpha$ emitters, quasars, or gamma ray bursts, have difficulties in terms of accuracy or event rate. We propose that recently discovered fast luminous blue transients (FLBTs) have a potential as a novel probe of reionization. We study the detectability of FLBTs at the epoch of reionization with upcoming WFIRST Wide-Field Instruments (WFI), using a star formation rate derived from galaxy observations and an event rate of FLBTs proportional to the star formation rate. We find that if FLBTs occur at a rate of 1\% of the core-collapse supernova rate, 2 (0.3) FLBTs per year per deg$^2$ at $z>6$ ($z>8$) can be detected by a survey with a limiting magnitude of 26.5 mag in the near-infrared band and a cadence of 10 days. We conclude that the WFIRST supernova deep survey can detect $\sim20$ FLBTs at the epoch of reionization in the near future.

We construct a numerical light curve model for interaction-powered supernovae that arise from an interaction between the ejecta and the circumstellar matter (CSM). In order to resolve the shocked region of an interaction-powered supernova, we solve the fluid equations and radiative transfer equation assuming the steady states in the rest frames of the reverse and forward shocks at each time step. Then we numerically solve the radiative transfer equation and the energy equation in the CSM with the thus obtained radiative flux from the forward shock as a radiation source. We also compare results of our models with observational data of two supernovae 2005kj and 2005ip classified as type IIn and discuss the validity of our assumptions. We conclude that our model can predict physical parameters associated with supernova ejecta and the CSM from the observed features of the light curve as long as the CSM is sufficiently dense. Furthermore, we found that the absorption of radiation in the CSM is an important factor to calculate the luminosity.

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

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

WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly constructed wind solution, where the optically thick outflow is launched from the carbon burning shell on an oxygen-neon core and accelerated by the rotating magnetic field to become supersonic and unbound well below the photosphere. Our model implies that WD J005311 has a mass of $M_* \sim 1.1\mbox{-}1.3\,M_\odot$, a magnetic field of $B_* \sim (2\mbox{-}5)\times 10^7\,\rm G$, and a spin angular frequency of $\Omega \sim 0.2\mbox{-}0.5 \,\rm s^{-1}$. The large magnetic field and fast spin support the carbon-oxygen WD merger origin. WD J005311 will neither explode as a type Ia supernova nor collapse into a neutron star. If the wind continues to blow another few kyr, WD J005311 will spin down significantly and join to the known sequence of slowly-rotating magnetic WDs. Otherwise it may appear as a fast-spinning magnetic WD and could be a new high energy source.

We present a bolometric light curve model of Type IIn supernovae powered by supernova ejecta colliding with a circumstellar medium. We estimate the conversion efficiency of the ejecta's kinetic energy to radiation at the reverse and forward shocks and find that a large density contrast makes a difference in the efficiency. The emission from the reverse shock can maintain high efficiency for a long time, and becomes important at the late phase of the light curve. We first construct a semi-analytical model that is applicable to the late phase of the light curve when the diffusion time of photons in the shocked region becomes negligible. We further develop radiation transfer simulations that incorporate these physical processes into the light curve. The numerical calculations predict light curves at early phases, which are testable by present and future short-cadence surveys. We compare our model with the bolometric light curve constructed from observations for a type IIn supernova 2005ip. Due to the reduced efficiency at the forward shock, we find from our model that the mass-loss rate of the progenitor star was $\approx 1\times 10^{-2}\ {\rm M_\odot \ yr^{-1}}$ for a wind velocity of $100\ {\rm km \ s^{-1}}$, an order of magnitude higher compared to previous work that used simple assumptions of the efficiency. This highlights the importance of taking these two components into account when extracting the physical parameters from observations.

We report on the current status of the radial velocity monitoring of nearby OB stars to look for binaries with small mass ratios. The combined data of radial velocities using the domestic 1-2 m-class telescopes seems to confirm the variations of radial velocities in a few weeks for four out of ten target single-lined spectroscopic binaries. More data are needed to estimate the exact periods and mass distributions.

We perform multi-wavelength light curve modeling of the recently discovered low-luminosity gamma-ray burst (GRB) 171205A. The emission model is based on the relativistic ejecta-circumstellar medium (CSM) interaction scenario. The collision of freely expanding spherical ejecta traveling at mildly relativistic velocities with the CSM produces the reverse and forward shocks, which dissipate a part of the kinetic energy of the mildly relativistic ejecta. We show that the early gamma-ray emission followed by an X-ray tail can be well explained by the radiation diffusing out from the shocked gas. Mildly relativistic ejecta with a kinetic energy of $5\times10^{50}$ erg and a wind-like CSM with a mass-loss rate of a few $10^{-4}\ M_\odot$ yr$^{-1}$ for a wind velocity of $10^3$ km s$^{-1}$, which extends up to $\sim 3\times 10^{13}$ cm, are required to account for the gamma-ray luminosity and duration of GRB 171205A. We also calculate the photospheric and non-thermal emission after the optically thick stage, which can fit the late-time X-ray, optical, and radio light curves. Our results suggest that the relativistic ejecta-CSM interaction can be a potential power source for low-luminosity GRBs and other X-ray bright transients.

Early-phase Type Ia supernovae (SNe Ia), especially those with luminosity enhancement within the first few days of explosions ("early-excess SNe Ia"), play an irreplaceable role in addressing the long-standing progenitor and explosion issue of SNe Ia. In this paper, we systematically investigate 11 early-excess SNe Ia from subluminous to luminous subclasses. Eight of them are selected from 23 SNe Ia with extremely early-phase optical light curves ("golden" early-phase SNe Ia), and three of them are selected from 40 SNe Ia (including 14 golden samples) with early-phase UV/NUV light curves. We found that previously discovered early-excess SNe Ia show a clear preference for specific SN Ia subclasses. In particular, the early-excess feature shown in all six luminous (91T- and 99aa-like) SNe Ia is in conflict with the viewing angle dependence predicted by the companion-ejecta interaction scenario. Instead, such a high early-excess fraction is likely related to the explosion physics of luminous SNe Ia; i.e. a more efficient detonation happening in the progenitor of luminous SNe Ia may consequently account for the early-excess feature powered by the radiation from a $^{56}$Ni-abundant outer layer. The diversity of early-excess features shown in different SN Ia subclasses suggests multiple origins of the discovered early-excess SNe Ia, challenging their applicability as a robust progenitor indicator. Further understanding of the early-excess diversity relies not only on multiband photometry and prompt-response spectroscopy of individual early-excess SNe Ia but also on investigations of the general early-phase light-curve behavior of each SN Ia subclass, which can be realized through ongoing/forthcoming transient survey projects in the near future.

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

The flow of fallback matter being shocked and repelled back by an energy deposition from a central object is discussed by using newly found self-similar solutions. We show that there exists a maximum mass accretion rate if the adiabatic index of the flow is less than or equal to 4/3. Otherwise we can find a solution with an arbitrarily large accretion rate by appropriately shrinking the energy deposition region. Applying the self-similar solution to supernova fallback, we discuss how the fate of newborn pulsars or magnetars depends on the fallback accretion and their spin-down power. Combining the condition for the fallback accretion to bury the surface magnetic field into the crust, we argue that supernova fallback with a rate of dM_{fb}/dt ~ 10^{-(4-6)} M_sun /s could be the main origin of the diversity of Galactic young neutron stars, i.e., rotation-powered pulsars, magnetars, and central compact objects.

Merging neutron stars generate shock waves that disintegrate heavy nuclei into nucleons especially in the outer envelope. It is expected that some of such neutrons having avoided capturing positrons remain as free neutrons even after disappearance of electron-positron pairs. To investigate how many free neutrons can be ejected from merging neutron stars, we performed special-relativistic Lagrangian hydrodynamics computation with simplified models of this phenomenon in which a spherically symmetric shock wave propagates in the hydrostatic envelope and emerge from the surface. We systematically study a wide parameter space of the size of the merging neutron stars and the energy involved in the shock waves. As a result, it is found that the mass of remaining free neutron is $10^{-7}$ to $10^{-6}\ M_{\odot}$, which is smaller than the previously expected by more than two orders of magnitude. There is a preferred energy of the order of $10^{48}$ erg that yields the maximum amount of free neutrons for large size of the envelope. We briefly discuss the emission from the free neutron layer and estimate the luminosity in the optical band to be about $7\times 10^{41}$ erg s $^{-1}(M_{\rm n}/10^{-6}\ M_\odot)$ in $\sim 30$ min after the merger.

The mechanism for the early-phase blue and excessive emission within the first few days, reported so far for a few type Ia supernovae (SNe Ia), has been suggested to be interaction of the SN ejecta either with a non-degenerate companion star or circumstellar media (CSM). Recently, another mechanism has been suggested within the context of the He detonation-triggered SN scenario (i.e., double detonation scenario or He-ignited violent merger), in which the radioactive isotopes in the outermost layer of the SN ejecta produce the early emission. In this paper, we investigate properties of the early-phase excessive emission predicted by these different scenarios. The He detonation scenario shows different behaviors in the early flash than the companion/CSM interaction scenarios. Especially clear diagnostics is provided once the behaviors in the UV and in the optical are combined. The spectra synthesized for the He detonation scenario are characterized by the absorptions due to the He detonation products, which especially develop in the decay phase. We further expect a relation between the properties of the early-phase flash and those of the maximum SN emission, in a way the brighter and slower initial flash is accompanied by a more substantial effect of the additional absorptions (and reddening). This relation, however, should be considered together with the maximum luminosity of the SN, since the larger luminosity suppresses the effect of the additional absorption. With these expected features, we address the possible origins of the observed excessive early-phase emission for a few SNe.

Recent observations found that electrons are accelerated to $\sim$10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected inferring that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Galactic disk by a factor $\sim$2.5. The cosmic rays (CR) produced by local CVs and novae will accumulate in the local bubble for $10^6$ - $10^7$ yrs. On these findings, we search for CRs from historic CVs and novae in the observed CR spectra. We model the CR spectra at the heliopause as sums of Galactic and local components based on observational data as much as possible. The initial Galactic CR electron and proton spectra are deduced from the gamma-ray emissivity, the local electron spectrum from the hard X-ray spectra at the CVs, and the local proton spectrum inferred by gamma-ray spectrum at novae. These spectral shapes are then expressed in a simple set of polynomial functions of CR energy and regressively fitted until the high-energy ($>$100 GeV) CR spectra near Earth and the Voyager-1 spectra at the heliopause are reproduced. We then extend the modeling to nuclear CR spectra and find that one spectral shape fits all local nuclear CRs and the apparent hardening of the nuclear CR spectra is caused by the roll-down of local nuclear spectra around 100 - 200 GeV. All local CR spectra populate in a limited energy band below 100 - 200 GeV and enhance gamma-ray emissivity below $\sim$10 GeV. Such an enhancement is observed in the inner Galaxy, suggesting the CR fluxes from CVs and novae are substantially higher there.

We investigate the properties of X-ray emission from shock breakout of a supernova in a stellar wind. We consider a simple model describing aspherical explosions, in which the shock front with an ellipsoidal shape propagates into the dense circumstellar matter. For this model, both X-ray light curves and spectra are simultaneously calculated using a Monte Carlo method. We show that the shock breakout occurs simultaneously in all directions in a steady and spherically symmetric wind. As a result, even for the aspherical explosion, the rise and decay timescales of the light curve do not significantly depend on the viewing angles. This fact suggests that the light curve of the shock breakout may be used as a probe of the wind mass loss rate. We compare our results with the observed spectrum and light curve of XRO 080109/SN 2008D. The observation can be reproduced by an explosion with a shock velocity of 60% of the speed of light and a circumstellar matter with a mass loss rate of 5.e-4 Msun/yr.

The stellar record of elemental abundances in satellite galaxies is important to identify the origin of r-process because such a small stellar system could have hosted a single r-process event, which would distinguish member stars that are formed before and after the event through the evidence of a considerable difference in the abundances of r-process elements, as found in the ultra-faint dwarf galaxy Reticulum II (Ret II). However, the limited mass of these systems prevents us from collecting information from a sufficient number of stars in individual satellites. Hence, it remains unclear whether the discovery of a remarkable r-process enrichment event in Ret II explains the nature of r-process abundances or is an exception. We perform high-resolution spectroscopic measurements of r-process abundances for twelve metal-poor stars in the Draco dwarf galaxy in the metallicity range of -2.5<[Fe/H]<-2. We found that these stars are separated into two groups with r-process abundances differing by one order of magnitude. A group of stars with high abundances of r-process elements was formed by a single r-process event that corresponds to the event evidenced in Ret II. On the other hand, the low r-process abundance group was formed by another sporadic enrichment channel producing a far fewer r-process elements, which is clearly identified for the first time. Accordingly, we identified two populations of stars with different r-process abundances, which are built by two r-process events that enriched gases at levels that differ by more than one order of magnitude.

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

Neutron star mergers (NSMs) are one of the most plausible sources of r-process elements in the universe. Therefore NSMs can also be a major source of ultra-heavy elements in cosmic rays. In this paper, we first estimate the contribution of r-process elements synthesized in NSMs to the ultra-heavy element cosmic rays (UHCRs) by calculating transport equations that take into account energy loss processes and spallations. We show that the flux of UHCRs accelerated by NSMs themselves fluctuates by many orders of magnitude on the timescale of several million years and can overwhelm UHCRs accelerated by supernova remnants (SNRs) after an NSM takes place within a few kilo-parsec from the solar system. Experiments with very long exposure times using meteorites as UHCR detectors can detect this fluctuation. As a consequence, we show that if NSMs are the primary source of UHCRs, future experiments using meteorites have a possibility to reveal the event history of NSMs in the solar vicinity. We also describe a possible difference in the abundance pattern and energy spectrum of UHCRs between NSM and SNR accelerations.

Hydrodynamical interaction of spherical ejecta freely expanding at mildly relativistic speeds into an ambient cold medium is studied in semi-analytical and numerical ways to investigate how ejecta produced in energetic stellar explosions dissipate their kinetic energy through the interaction with the surrounding medium. We especially focus on the case in which the circumstellar medium is well represented by a steady wind at a constant mass-loss rate having been ejected from the stellar surface prior to the explosion. As a result of the hydrodynamical interaction, the ejecta and circumstellar medium are swept by the reverse and forward shocks, leading to the formation of a geometrically thin shell. We present a semi-analytical model describing the dynamical evolution of the shell and compare the results with numerical simulations. The shell can give rise to bright emission as it gradually becomes transparent to photons. while it is optically thick. We develop an emission model for the expected emission from the optically thick shell, in which photons in the shell gradually diffuse out to the interstellar space. Then, we investigate the possibility that radiation powered by the hydrodynamical interaction is the origin of an underluminous class of gamma-ray bursts.

The main astronomical source of r-process elements has not yet been identified. One plausible site is neutron star mergers (NSMs), but from perspective of the Galactic chemical evolution, it has been pointed out that NSMs cannot reproduce the observed r-process abundance distribution of metal-poor stars at [Fe/H] $< -3$. Recently, Tsujimoto & Shigeyama (2014) pointed out that NSM ejecta can spread into much larger volume than ejecta from a supernova. We re-examine the enrichment of r-process elements by NSMs considering this difference in propagation using the chemical evolution model under the hierarchical galaxy formation. The observed r-process enhanced stars around [Fe/H] $\sim -3$ are reproduced if the star formation efficiency is lower for low-mass galaxies under a realistic delay time distribution for NSMs. We show that a significant fraction of NSM ejecta escape from its host proto-galaxy to pollute intergalactic matter and other proto-galaxies. The propagation of r-process elements over proto-galaxies changes the abundance distribution at [Fe/H] $< -3$ and obtains distribution compatible with observations of the Milky Way halo stars. In particular, the pre-enrichment of intergalactic medium explains the observed scarcity of EMP stars without Ba and abundance distribution of r-process elements at [Fe/H] $\lesssim -3.5$.

A two-dimensional special relativistic radiation-hydrodynamics code is developed and applied to numerical simulations of supernova shock breakout in bipolar explosions of a blue supergiant. Our calculations successfully simulate the dynamical evolution of a blast wave in the star and its emergence from the surface. Results of the model with spherical energy deposition show a good agreement with previous simulations. Furthermore, we calculate several models with bipolar energy deposition and compare their results with the spherically symmetric model. The bolometric light curves of the shock breakout emission are calculated by a ray-tracing method. Our radiation-hydrodynamic models indicate that the early part of the shock breakout emission can be used to probe the geometry of the blast wave produced as a result of the gravitational collapse of the iron core.

We present a novel semi-analytic solution describing the propagation of a spherical blast wave driven by a central energy source. The initial density profile has a power-law function of the distance from the center and the energy is injected only into the central region at a rate given by a power-law function of time. This solution is composed of three regions separated by the contact surface and the shock front. The innermost region is assumed to be uniform and the outside of the contact surface includes the shocked matter described by self-similar solutions. We analytically derive the applicable range of parameters of this solution from requirements to satisfy boundary conditions. We perform numerical simulations for flows with various values of parameters, some of which reside out of the thus derived applicable range, and compare with the semi-analytic solutions.

The progenitor systems for type Ia supernovae are still controversial. One of the methods to test the proposed scenario for the progenitor systems is to identify companions that are supposed to survive according to the so-called single degenerate scenario. These companions might be affected by supernova ejecta. We present several numerical simulations of surviving red-giant companions whose envelopes were stripped and heated. We find that red-giants with less-massive helium cores ($\lesssim0.30\,M_{\odot}$) can be so faint after the supernovae that we cannot detect them. In addition, we apply the results to the case of SNR 0509-67.5, and put constraints on the helium core mass, envelope stripping, and energy injection under the single degenerate scenario for type Ia supernovae.

In the single degenerate (SD) scenario of type Ia supernovae (SNe Ia), the collision of the ejecta with its companion results in stripping hydrogen rich matter from the companion star. This hydrogen rich matter might leave its trace in the light curves and/or spectra. In this paper, we perform radiation hydrodynamical simulations of this collision for three binary systems. As a result, we find that the emission from the shock-heated region is not as strong as in the previous study. This weak emission, however, may be a result of our underestimate of the coupling between the gas and radiation in the shock interaction. Therefore, though our results suggest that the observed early light curves of SNe Ia can not rule out binary systems with a short separation as the progenitor system, more elaborate numerical studies will be needed to reach a fair conclusion. Alternatively, our results indicate that the feature observed in the early phase of a recent type Ia SN 2014J might result from interaction of the ejecta with a red giant companion star. We also discuss the dependence of spectral features of H alpha and Si II absorption lines on viewing angles and investigate whether they can constrain the event rate of the SD progenitor.

Chemical abundance of r-process elements in nearby dwarf spheroidal (dSph) galaxies is a powerful tool to probe the site of r-process since their small-mass scale can sort out individual events producing r-process elements. A merger of binary neutron stars is a promising candidate of this site. In faint, or less massive dSph galaxies such as the Draco, a few binary neutron star mergers are expected to have occurred at most over the whole past. We have measured chemical abundances including Eu and Ba of three red giants in the Draco dSph by Subaru/HDS observation. The Eu detection for one star with [Fe/H]=-1.45 confirms a broadly constant [Eu/H] of ~-1.3 for stars with [Fe/H]>-2. This feature is shared by other dSphs with similar masses, i.e., the Sculptor and the Carina, and suggests that neutron star merger is the origin of r-process elements in terms of its rarity. In addition, two very metal-poor stars with [Fe/H]=-2.12 and -2.51 are found to exhibit very low Eu abundances such as [Eu/H]<-2 with an implication of a sudden increase of Eu abundance by more than 0.7 dex at [Fe/H] ~ -2.2 in the Draco dSph. In addition, the detection of Ba abundances for these stars suggests that the r-process enrichment initiated no later than the time when only a few % of stars in the present-day Draco dSph was formed. Though an identification of the origin of an early Eu production inside the Draco dSph should be awaited until more abundance data of stars with [Fe/H]<-2 in the Draco as well as other faint dSphs become available, the implied early emergence of Eu production event might be reconciled with the presence of extremely metal-poor stars enriched by r-process elements in the Galactic halo.

Growing interests in neutron star (NS) mergers as the origin of r-process elements have sprouted since the discovery of evidence for the ejection of these elements from a short-duration gamma-ray burst. The hypothesis of a NS merger origin is reinforced by a theoretical update of nucleosynthesis in NS mergers successful in yielding r-process nuclides with A>130. On the other hand, whether the origin of light r-process elements are associated with nucleosynthesis in NS merger events remains unclear. We find a signature of nucleosynthesis in NS mergers from peculiar chemical abundances of stars belonging to the Galactic globular cluster M15. This finding combined with the recent nucleosynthesis results implies a potential diversity of nucleosynthesis in NS mergers. Based on these considerations, we are successful in the interpretation of an observed correlation between [light r-process/Eu] and [Eu/Fe] among Galactic halo stars and accordingly narrow down the role of supernova nucleosynthesis in the r-process production site. We conclude that the tight correlation by a large fraction of halo stars is attributable to the fact that core-collapse supernovae produce light r-process elements while heavy r-process elements such as Eu and Ba are produced by NS mergers. On the other hand, stars in the outlier, composed of r-enhanced stars ([Eu/Fe]>+1) such as CS22892-052, were exclusively enriched by matter ejected by a subclass of NS mergers that is inclined to be massive and consist of both light and heavy r-process nuclides.

We investigate the hydrodynamical evolution of an ultra-relativistic fireball colliding with a freely expanding gas. The hydrodynamical interaction of the fireball and the gas results in the formation of a geometrically thin shell. We study the dynamical evolution of the shell by an analytical way and perform a numerical simulation equipped with an adaptive mesh refinement to investigate the internal structure of the shell. The shocked gas can give rise to bright emission in the X-ray and gamma-ray energy range. We propose that the breakout emission from the forward shock and the photospheric emission from the reverse-shocked fireball contribute to early gamma-ray emission from gamma-ray bursts.

While type Ia Supernovae (SNe Ia) have been used as precise cosmological distance indicators, their progenitor systems remain unresolved. One of the key questions is if there is a non-degenerate companion star at the time of a thermonuclear explosion of a white dwarf (WD). In this paper, we investigate if an interaction between the SN ejecta and the companion star may result in observable footprints around the maximum brightness and thereafter, by performing multi-dimensional radiation transfer simulations based on hydrodynamic simulations of the interaction. We find that such systems result in variations in various observational characteristics due to different viewing directions, while the predicted behaviors (redder and fainter for the companion direction) are opposite to what were suggested by the previous study. The variations are generally modest and within observed scatters. However, the model predicts trends between some observables different from observationally derived, thus a large sample of SNe Ia with small calibration errors may be used to constrain the existence of such a companion star. The variations in different colors in optical band passes can be mimicked by external extinctions, thus such an effect could be a source of a scatter in the peak luminosity and derived distance. After the peak, hydrogen-rich materials expelled from the companion will manifest themselves in hydrogen lines. H$\alpha$ is however extremely difficult to identify. Alternatively, we find that P$_{\beta}$ in post-maximum near-infrared spectra can potentially provide powerful diagnostics.

The origin of r-process elements remains unidentified and still puzzles us. The recent discovery of evidence for the ejection of r-process elements from a short-duration gamma-ray burst singled out neutron star mergers (NSMs) as their origin. In contrast, core-collapse supernovae are ruled out as the main origin of heavy r-process elements (A>110) by recent numerical simulations. However, the properties characterizing NSM events - their rarity and high yield of r-process elements per event - have been claimed to be incompatible with the observed stellar records on r-process elements in the Galaxy. We add to this picture with our results, which show that the observed constant [r-process/H] ratio in faint dwarf galaxies and one star unusually rich in r-process in the Sculptor galaxy agree well with this rarity of NSM events. Furthermore, we found that a large scatter in the abundance ratios of r-process elements to iron in the Galactic halo can be reproduced by a scheme that incorporates an assembly of various protogalactic fragments, in each of which r-process elements supplied by NSMs pervade the whole fragment while supernovae distribute heavy elements only inside the regions swept up by the blast waves. Our results demonstrate that NSMs occurring at Galactic rate of 12-23 per Myr are the main site of r-process elements, and we predict the detection of gravitational waves from NSMs at a high rate with upcoming advanced detectors.

We investigate the mechanism to reproduce notable spectral features at the ignition phase of nova explosion observed for a super-Eddington X-ray transient source MAXI J0158$-$744 in the Small Magellanic Cloud. These are a strong Ne IX emission line at 0.92 keV with a large equivalent width of $0.32^{+0.21}_{-0.11}$ keV and the absence of Ne X line at 1.02 keV. In this paper, we calculate the radiative transfer using a Monte Carlo code, taking into account the line blanketing effect due to transitions of N, O, Ne, Mg and Al ions in an accelerating wind emanating from a white dwarf with a structure based on a spherically symmetric stationary model. We found that the strong Ne IX line can be reproduced if the mass fraction of Ne is enhanced to $10^{-3}$ or more and that of O is reduced to $\sim5\times10^{-9}$ or less and that the absence of other lines including Ne X ions at higher energies can be also reproduced by the line blanketing effect. This enhancement of the Ne mass fraction indicates that the ejecta are enriched by Ne dredged up from the surface of the white dwarf composed of O, Ne, and Mg rather than C and O, as already pointed out in the previous work. We argue that the CNO cycle driving this nova explosion converted most of C and O into N and thus reduced the O mass fraction.

Bearing in mind the application to high-magnetic-field (high-B) radio pulsars, we investigate two-dimensional (2D) thermal evolutions of neutron stars (NSs). We pay particular attention to the influence of different equilibrium configurations on the surface temperature distributions. The equilibrium configurations are constructed in a systematic manner, in which both toroidal and poloidal magnetic fields are determined self-consistently with the inclusion of general relativistic effects. To solve the 2D heat transfer inside the NS interior out to the crust, we have developed an implicit code based on a finite-difference scheme that deals with anisotropic thermal conductivity and relevant cooling processes in the context of a standard cooling scenario. In agreement with previous studies, the surface temperatures near the pole become higher than those in the vicinity of the equator as a result of anisotropic heat transfer. Our results show that the ratio of the highest to the lowest surface temperatures changes maximally by one order of magnitude, depending on the equilibrium configurations. Despite such difference, we find that the area of such hot and cold spots is so small that the simulated X-ray spectrum could be well reproduced by a single temperature blackbody fitting.

Analyzing light curves of a complete sample of bright Swift long gamma-ray bursts (LGRBs) of which the peak photon fluxes constructed with the bin width of 1 second in the Swift 15-350 keV energy band exceed 2.6 photons cm$^{-2}$s$^{-1}$, we confirm that there do exist the third class in GRBs in addition to short and long GRBs. Being different from previous works based on the duration, fluence, etc. our classification method is based on two properties both quantified with light curve shapes of the prompt emission: the Absolute Deviation from the Constant Luminosity of their cumulative light curve $ADCL$, and the ratio of the mean counts to the maximum counts $\bar C/C_{\rm max}$. These are independent of the distance and the jet opening angle. A cluster analysis via the Gaussian mixture model detects three subclasses: one consisting of LGRBs with small $ADCL$ and large $\bar C/C_{\rm max}$ values referred to as Type I, one with large $ADCL$ and large $\bar C/C_{\rm max}$ referred to as Type II, and one with intermediate $ADCL$ and small $\bar C/C_{\rm max}$, which is composed of contaminating short GRBs with the extended emission. This result is reinforced by different temporal and spectral indices of their X-ray afterglows. The difference is prominent in the temporal index of the steep decay phase in particular: the indices for Type I LGRBs distribute between $-6$ and $-3$ while those for Type II LGRBs between $-3$ and $-2$. From these properties, we propose a possible scenario with different central engines: an accreting black hole and a magnetar.

We present theoretical expectations for non-thermal emission due to the bulk Comptonization at the ultra-relativistic shock breakout. We calculate the transfer of photons emitted from the shocked matter with a Monte Carlo code fully taking into account special relativity. As a hydrodynamical model, we use a self-similar solution of Nakayama & Shigeyama (2005). Our calculations reveal that the spectral shape exhibits a double peak or a single peak depending on the shock temperature at the shock breakout. If it is significantly smaller than the rest energy of an electron, the spectrum has a double peak. We also display a few example of light curves, and estimate the total radiation energy. In comparison with observations of gamma-ray bursts, a part of the higher energy component in the spectra and the total energy can be reproduced by some parameter sets. Meanwhile, the lower energy counterpart in the Band function is not reproduced by our results and the duration time seems too short to represent the entire event of a gamma-ray burst. Therefore the subsequent phase will constitute the lower energy part in the spectrum.

Overwhelming diversity of long gamma-ray bursts (LGRBs), discovered after the launch of {\it Swift} satellite, is a major obstacle to LGRB studies. Recently, it is shown that the prompt emission of LGRBs is classified into three subclasses: Type I, Type II LGRBs populating separate fundamental planes in a 3D space defined by the peak luminosity, the duration, and the spectral peak energy, and outliers not belonging to either of the planes. Here we show that Type I LGRBs (LGRBs I) exhibit different shapes of light curves from Type II LGRBs (LGRBs II). Furthermore, we demonstrate that this classification has uncovered a new scaling law in the light curves of LGRBs II spanning 8 orders of magnitude from the prompt to late X-ray afterglow emission. The scaled light curve has four distinct phases. The first phase has a characteristic time scale while the subsequent three phases exhibit power law behaviors with different exponents. We discuss its possible interpretation in terms of the emission from an optically thick fireball propagating in the cricumstellar matter at relativistic speeds and argue that the observed four phases correspond to its hydrodynamical phases. Our classification scheme can pin down the intrinsic luminosities of LGRBs II through the scaling law from a sample of polymorphic GRBs. Further refinement of this scheme and scaling law will make a subclass of LGRBs a standard candle as reliable and accurate as Type Ia supernovae in the more distant universe than supernovae can reach.

We performed a series of hydrodynamical calculations of an ultra-relativistic jet propagating through a massive star and the circumstellar matter to investigate the interaction between the ejecta and the circumstellar matter. We succeed in distinguishing two qualitatively different cases in which the ejecta are shocked and adiabatically cool. To examine whether the cocoon expanding at subrelativistic speeds emits any observable signal, we calculate expected photospheric emission from the cocoon. It is found that the emission can explain early thermal X-ray emission recently found in some long gamma-ray bursts. The result implies that the difference of the circumstellar environment of long gamma-ray bursts can be probed by observing their early thermal X-ray emission.

A time delay of Type Ia supernova (SN Ia) explosions hinders the imprint of their nucleosynthesis on stellar abundances. However, some occasional cases give birth to stars that avoid enrichment of their chemical compositions by massive stars and thereby exhibit a SN Ia-like elemental feature including a very low [Mg/Fe] (~-1). We highlight the elemental feature of Fe-group elements for two low-Mg/Fe objects detected in nearby galaxies, and propose the presence of a class of SNe Ia that yield the low abundance ratios of [Cr,Mn,Ni/Fe]. Our novel models of chemical evolution reveal that our proposed class of SNe Ia (slow SNe Ia) is associated with ones exploding on a long timescale after their stellar birth, and gives a significant impact on the chemical enrichment in the Large Magellanic Cloud (LMC). In the Galaxy, on the other hand, this effect is unseen due to the overwhelming enrichment by the major class of SNe Ia that explode promptly (prompt SNe Ia) and eject a large amount of Fe-group elements. This nicely explains the different [Cr,Mn,Ni/Fe] features between the two galaxies as well as the puzzling feature seen in the LMC stars exhibiting very low Ca but normal Mg abundances. Furthermore, the corresponding channel of slow SN Ia is exemplified by performing detailed nucleosynthesis calculations in the scheme of SNe Ia resulting from a 0.8+0.6 solar mass white dwarf merger.

We investigate how explosions of aspherical supernovae (A-SNe) can influence star formation histories and chemical evolution of dwarf galaxies by using a new chemodynamical model. We mainly present the numerical results of two comparative models so that the A-SN feedback effects on galaxies can be more clearly seen. SNe originating from stars with masses larger than 30M_sun are A-SNe in the "ASN" model whereas all SNe are spherical ones (S-SNe) in the "SSN" model. Each S-SN and A-SN are assumed to release feedback energy of 10^{51} erg and 10^{52} erg, respectively, and chemical yields and feedback energy of A-SN ejecta depend on angles between the axis of symmetry and the ejection directions. We find that star formation can become at least by a factor of ~3 lower in the ASN model in comparison with the SSN one owing to the more energetic feedback of A-SNe. As a result of this, chemical evolution can proceed very slowly in the ASN model. A-SN feedback effects can play a significant role in the formation of giant gaseous holes and energetic gaseous outflow and unique chemical abundances (e.g., high [Mg/Ca]). Based on these results, we provide a number of implications of the A-SN feedback effects on galaxy formation and evolution.

Using the Subaru 8.2m Telescope with an IRCS Echelle spectrograph, we obtained high-resolution (R=10,000) near-infrared (1.01-1.38 \mu m) spectra of images A and B of the gravitationally lensed QSO B1422+231 (z=3.628) consisting of four known lensed images. We detected MgII absorption lines at z=3.54, which show a large variance of column densities (~ 0.3 dex) and velocities (~ 10 km/s) between the sightlines A and B with a projected separation of only 8.4h_{70}^{-1} pc at the redshift. This is the smallest spatial structure of the high-z gas clouds ever detected after Rauch et al. found a 20-pc scale structure for the same z=3.54 absorption system using optical spectra of images A and C. The observed systematic variances imply that the system is an expanding shell as originally suggested by Rauch et al. By combining the data for three sightlines, we managed to constrain the radius and expansion velocity of the shell (~ 50-100 pc, 130 km/s), concluding that the shell is truly a supernova remnant (SNR) rather than other types of shell objects, such as a giant HII region. We also detected strong FeII absorption lines for this system, but with much broader Doppler width than that of \alpha-element lines. We suggest that this FeII absorption line originates in a localized FeII-rich gas cloud that is not completely mixed with plowed ambient interstellar gas clouds showing other \alpha-element low-ion absorption lines. Along with the Fe richness, we conclude that the SNR is produced by an SNIa explosion.

We present a new self-similar solution describing early evolution of an ultra-relativistic flow resulting from a collision of homologously expanding spherical ejecta with the circumstellar matter, in which a shock wave propagates in the circumstellar matter while a weak discontinuity propagates in the ejecta at the sound speed

We investigate effects of the magnetic field on the propagation of laminar flames of nuclear reactions taking place in white dwarfs (WDs) with the mass close to the Chandrasekhar limit. We calculate the velocities of laminar flames parallel and perpendicular to uniform magnetic fields as eigenvalues of steady solutions for magnetic hydrodynamical equations. As a result, we find that even when the magnetic pressure does not dominate the entire pressure it is possible for the magnetic field to suppress the flame propagation through the thermal conduction. Above the critical magnetic field, the flame velocity decreases with increasing magnetic field strength as $v \sim B^{-1}$. In media with densities of $10^{7}, 10^{8}, \,\mathrm{and}\,10^{9} \unit{g\,cm^{-3}}$, the critical magnetic fields are orders of $\sim 10^{10}, 10^{11}, \,\mathrm{and}\,10^{12} \unit{G}$, respectively.

To account for the non-thermal emission from the classical nova V2491 Cygni, we perform a series of numerical calculations of radiative transfer of \gamma-ray photons from the radioactive isotope $^{22}$Na in the matter ejected from a white dwarf. Using a simple wind model for the dynamical evolution of the ejecta and a monte-carlo code, we calculate radiative transfer of the \gamma-ray photons in the ejecta. Repeated scattering of the \gamma-ray photons by electrons in the ejecta, i.e., Compton degradation, results in an extremely flat spectrum in the hard X-ray range, which successfully reproduces the observed spectrum of the X-ray emission from V2491 Cygni. The amount of the isotope $^{22}$Na synthesized in the ejecta is required to be 3\times 10^{-5} M_\odot to account for the flux of the hard X-ray emission. Our model indicates that the ejecta become transparent to the \gamma-ray photons within several tens days. Using the results, we briefly discuss the detectability of the \gamma-ray line emission by the {\it INTEGRAL} gamma-ray observatory and the {\it Fermi} Gamma-ray Space Telescope.

We investigate the neutrino-process in an energetic Type Ic supernova (SN Ic) and resultant productions of the light elements including boron and its stable isotopes. SN Ic is a very unique boron source because it can produce boron through not only spallation reactions as discussed in Nakamura & Shigeyama (2004) but also the neutrino-process. The neutrino-process is considered to occur in core-collapse supernovae and previous studies were limited to Type II supernovae (SNe II). Although the progenitor star of an SN Ic does not posses a He envelope so that 7Li production via the neutrino-process is unlikely, 11B can be produced in the C-rich layers. We demonstrate a hydrodynamic simulation of SN Ic explosion and estimate the amounts of the light elements produced via the neutrino-process for the first time, and also the subsequent spallation reactions between the outermost layers of compact SN Ic progenitor and the ambient medium. We find that the neutrino-process in the current SN Ic model produces a significant amount of 11B, which is diluted by 10B from spallation reactions to get closer to B isotopic ratios observed in meteorites. We also confirm that high-temperature mu- and tau-neutrinos and their anti-neutrinos, reasonably suggested from the compact structure of SN Ic progenitors, enhance the light element production through the neutral-current reactions, which may imply an important role of SNe Ic in the Galactic chemical evolution.

We investigate roles of the bulk comptonization process in the propagation of thermal photons emitted at the shock breakout of a supernova explosion. We use a hydrodynamical model based on a self-similar solution for the shock breakout. The propagation of photons is treated by a Monte-Carlo technique. Results of the simulations successfully reproduce the power-law spectrum observed by X-ray observations for XRF 080109/SN 2008D, if a mildly relativistic shock velocity is assumed. Effects of some radiative processes, photoionization, radiative recombination, and free-free absorption on the propagation of emitted photons are also investigated. It is found that none of them hardly changes the spectrum regarding the progenitor stars of type Ib or Ic supernovae. Light curves calculated under the assumption of a spherical explosion indicate that the progenitor radius is required to be $/sim 10^{13}$ cm.

We investigate effects of aspherical energy deposition in core-collapse supernovae on the light curve of the supernova shock breakout. We performed two-dimensional hydrodynamical calculations of an aspherical supernova explosion to obtain the time when a shock wave generated in the stellar interior reaches the stellar surface in each radial direction. Using results of the calculations, light curves during the shock breakout are derived in an approximate way. We show that the light curve during the shock breakout can be a strong indicator of aspherical properties of core-collapse supernovae.

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

A new scheme for numerical integration of the 1D2V relativistic Vlasov-Maxwell system is proposed. Assuming that all particles in a cell of the phase space move with the same velocity as that of the particle located at the center of the cell at the beginning of each time step, we successfully integrate the system with no artificial loss of particles. Furthermore, splitting the equations into advection and interaction parts, the method conserves the sum of the kinetic energy of particles and the electromagnetic energy. Three test problems, the gyration of particles, the Weibel instability, and the wakefield acceleration, are solved by using our scheme. We confirm that our scheme can reproduce analytical results of the problems. Though we deal with the 1D2V relativistic Vlasov-Maxwell system, our method can be applied to the 2D3V and 3D3V cases.

We present results of a 2D3V kinetic Vlasov simulation of the Weibel instability. The kinetic Vlasov simulation allows us to investigate the velocity distribution of dilute plasmas, in which the effect of collisions between particles is negligible, and has the advantage that the accuracy of the calculated velocity distribution does not depend on the density of plasmas at each point in the physical space. We succeed in reproducing some features of the Weibel instability shown by other simulations, for example, the exponentially growing phase, the saturation of the magnetic field strength, the formation of filamentary structure, and the coalescence of the filaments. Especially, we concentrate on the behavior of the filaments after the saturation of the magnetic field strength and find that there is a kind of quasi-equilibrium states before the coalescence occurs. Furthermore, it is found that an analytical solution for stationary states of the 2D3V Vlasov-Maxwell system can reproduce some dominant features of the quasi-equilibrium, e.g, the configuration of the magnetic field and the velocity distribution at each point. The analytical expression could give a plausible model for the transition layer of a collisionless shock where a strong magnetic field generated by the Weibel instability provides an effective dissipation process instead of collisions between particles.

A novel method to derive stationary solutions of the Vlasov-Maxwell system is established. This method is based on the assumption that the deviation of the velocity distribution from the Maxwell-Boltzmann distribution can be expanded by the Hermite polynomials. By applying our method, a new two-dimensional equilibrium is derived, which may provide an initial setup for investigations of three-dimensional collisionless reconnection of magnetic fields.

The unveiled main-sequence splitting in omega Centauri as well as NGC 2808 suggests that matter highly-enriched in He (in terms of its mass fraction Y~0.4) was produced and made the color of some main-sequence stars bluer in these globular clusters (GCs). The potential production site for the He-rich matter is generally considered to be massive AGB stars that experience the second dredge-up. However, it is found that massive AGB stars provide the matter with Y~0.35 at most, while the observed blue-shift requires the presence of Y~0.4 matter. Here, we show that extra mixing, which operates in the red giant phase of stars less massive than ~2Msun, could be a mechanism that enhances He content in their envelopes up to Y~0.4. The extra mixing is supposed to be induced by red giant encounters with other stars in a collisional system like GCs. The Y~0.4 matter released in the AGB phase has alternative fates to (i) escape from a GC or (ii) be captured by kinematically cool stars through encounters. The AGB ejecta in omega Cen, which follows the latter case, can supply sufficient He to cause the observed blue-shift. Simultaneously, this scheme generates the extreme horizontal branch, as observed in omega Cen in response to the higher mass loss rates, which is also caused by stellar encounters.

We report our first results of photometric and spectroscopic observations for Tycho's supernova remnant (SNR Tycho) to search for the companion star of a type Ia supernova (SN Ia). From photometric observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag (in $V$-band) for spectroscopy, which are located within a circular region with the radius of 30" around the center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient amount of Fe I, we should be able to detect absorption lines at 3720 \AA and 3860 \AA associated with transitions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star of a SN Ia using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect significant absorption lines from two stars at 3720 \AA. Since widths of their absorption lines are broad, it is likely that the detected absorptions are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, we suggest another interpretation that the detected absorption lines might be caused by the peculiarity of stars. A star named Tycho(E) has the absorption line at 3720 \AA and its projected position is close to the center of SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's supernova.

We present a self-similar solution to describe the propagation of a shock wave whose energy is deposited or lost at the front. Both of the propagation of the shock wave in a medium having a power-law density profile and the expansion of the medium to a vacuum after the shock breakout are described with a Lagrangian coordinate. The Chapman-Jouguet detonation is found to accelerate the medium most effectively. The results are compared with some numerical simulations in the literature. We derive the fractions of the deposited/lost energy at the shock front in some specific cases, which will be useful when applying this solution to actual phenomena.

We investigate relativistic flows after a shock wave generated in a star arrives at the surface. First, the sphericity effect is involved through a successive approximation procedure by adding correction terms to an already known self-similar solution in the ultra-relativistic limit assuming the plane parallel geometry. We also find that the involvement of sphericity increases the acceleration in the early phase as compared with that of the original plane-parallel flow. Second, we obtain semi-analytic solutions for a mildly relativistic flow in which the rest mass energy density is not negligible in the equation of state. To take into account this, we use the enthalpy and the pressure instead of using the density and the pressure as thermodynamic variables. These solutions assume self-similar evolutions except for the initial conditions. Third, we have carried out numerical calculations with a special-relativistic hydrodynamical code based on the Godunov method in order to check the applicability of the above sphericity corrections and the semi-analytic solutions. The equation of state used in our calculations includes the rest mass energy density. Comparisons with results of numerical calculations support the validity of the sphericity correction terms. The evolutions of the pressure and the Lorentz factor of each fluid element of the semi-analytic solution for mildly relativistic flows match the numerical results at least in early phases. We also investigate the final free expansion phases by this code. We have derived the final energy distributions for these flows and compare them with previous works.

The origin of a double main-sequence (MS) in omega Centauri is explored. We have shown from theoretical calculations on the stellar evolution that the colors of MS stars are shifted to those of the observed blue MS if the surface layers are polluted by He-rich materials with the mass of ~ 0.1 solar mass. Stars are supposed to be polluted through numerous encounters with the ejecta descended from massive asymptotic giant-branch (AGB) stars. Two populations of stars with different kinematics exceptionally observed in omega Cen indicate that kinematically cooler stars are more polluted through encounters with AGB ejecta than kinematically hotter ones because the accretion rate is inversely proportional to the cube of the relative velocity. We propose that both of these factors split the MS in omega Cen. This theoretical scheme explains why only omega Cen exhibit a double MS and matches the amount of He necessary to produce the blue MS with that supplied from massive AGB stars. Furthermore, we predict that even if globular clusters (GCs) possess only one generation of stars, the velocity dispersion of stars broaden the MS in the color-magnitude diagram as long as the GCs are massive enough to keep the AGB ejecta after the burst of star formation. This view explains the broad MS recently found in the GC NGC 2808 which exhibits no scatter in [Fe/H] and thus is likely to consist of a single generation of stars unlike the case of omega Cen.

It is known that the carbon-enhanced, extremely metal-poor (CEMP) stars constitute a substantial proportion in the extremely metal-poor (EMP) stars of the Galactic Halo, by far larger than CH stars in Population II stars. We investigate their origin with taking into account an additional evolutionary path to the surface carbon-enrichment, triggered by hydrogen engulfment by the helium flash convection, in EMP stars of $[Fe/H] \lesssim -2.5$. This process is distinct from the third dredge-up operating in more metal-rich stars and also in EMP stars. In binary systems of EMP stars, the secondary stars become CEMP stars through mass transfer from the primary stars of low and intermediate masses, which have developed the surface carbon-enhancement. Our binary scenario can predict the variations in the abundances not only for carbon but also for nitrogen and s-process elements and reasonably explain the observed properties such as the stellar distributions with respect to the carbon abundances, the binary periods, and the evolutionary stages. Furthermore, from the observed frequencies of CEMP stars with and without s-process element enhancement, we demonstrate that the initial mass function of EMP stars need to give the mean mass $~10\msun$ under the reasonable assumptions on the distributions of orbital separations and mass ratio of binary components. This also indicates that the currently observed EMP stars were exclusively born as the secondary members of binaries, making up $\sim 10%$ remnants of EMP binary systems of mass $~10^8\msun$ in total; in addition to CEMP stars with white dwarf companions, a significant fraction of them have experienced supernova explosions of their companions. We discuss the implications of the present results in relation to the formation of Galactic halo.

We investigate energetic type Ic supernovae as production sites for Li6 and Be in the early stages of the Milky Way. Recent observations have revealed that some very metal-poor stars with [Fe/H]<-2.5 possess unexpectedly high abundances of Li6. Some also exbihit enhanced abundances of Be as well as N. From a theoretical point of view, recent studies of the evolution of metal-poor massive stars show that rotation-induced mixing can enrich the outer H and He layers with C, N, and O (CNO) elements, particularly N, and at the same time cause intense mass loss of these layers. Here we consider energetic supernova explosions occurring after the progeniter star has lost all but a small fraction of the He layer. The fastest portion of the supernova ejecta can interact directly with the circumstellar matter (CSM), both composed of He and CNO, and induce light element production through spallation and He-He fusion reactions. The CSM should be sufficiently thick to energetic particles so that the interactions terminate within its innermost regions. We calculate the resulting Li6/O and Be9/O ratios in the ejecta+CSM material out of which the very metal-poor stars may form. We find that they are consistent with the observed values if the mass of the He layer remaining on the pre-explosion core is 0.01-0.1 solar mass, and the mass fraction of N mixed in the He layer is about 0.01. Further observations of Li6, Be and N at low metallicity should provide critical tests of this production scenario.

We find self-similar solutions to describe the interaction of spherically symmetric ejecta expanding at relativistic speeds with an ambient medium having a power law density distribution. Using this solution, the time evolution of the Lorentz factor of the outer shock is derived as a function of the explosion energy, the mass of the ejecta, and parameters for the ambient medium. These solutions are an ultra-relativistic version of the solutions for the circumstellar interaction of supernova ejecta obtained by Chevalier and extensions of the relativistic blast wave solutions of Blandford & Mckee.

We propose a method to identify the companion stars of type Ia supernovae (SNe Ia) in young supernova remnants (SNRs) by recognizing distinct features of absorption lines due to Fe I appearing in the spectrum. If a sufficient amount of Fe I remains in the ejecta, Fe I atoms moving toward us absorb photons by transitions from the ground state to imprint broad absorption lines exclusively with the blue-shifted components in the spectrum of the companion star. To investigate the time evolution of column depth of Fe I in the ejecta, we have performed hydrodynamical calculations for SNRs expanding into the uniform ambient media, taking into account collisional ionizations, excitations, and photo-ionizations of heavy elements. As a result, it is found that the companion star in Tycho's SNR will exhibit observable features in absorption lines due to Fe I at $\lambda\lambda = 371.9935$ nm and 385.9911 nm if a carbon deflagration SN model is taken. However, these features may disappear by taking another model that emits a few times more intense ionizing photons from the shocked outer layers. To further explore the ionization states in the freely expanding ejecta, we need a reliable model to describe the structure of the outer layers.

Renewed interest in the first stars that were formed in the universe has led to the discovery of extremely iron-poor stars. Since several competing scenarios exist, our understanding of the mass range that determines the observed elemental abundances remains unclear. In this study, we consider three well-studied metal-poor stars in terms of the theoretical supernovae (SNe) model. Our results suggest that the observed abundance patterns in the metal-poor star BD +80 245 and the pair of stars HD 134439/40 agree strongly with the theoretical possibility that these stars inherited their heavy element abundance patterns from SNe initiated by thermonuclear runaways in the degenerate carbon-oxygen cores of primordial asymptotic giant branch stars with \~3.5-5 solar masses. Recent theoretical calculations have predicted that such SNe could be originated from metal-free stars in the intermediate mass range. On the other hand, intermediate mass stars containing some metals would end their lives as white dwarfs after expelling their envelopes in the wind due to intense momentum transport from outgoing photons to heavy elements. This new pathway for the formation of SNe requires that stars are formed from the primordial gas. Thus, we suggest that stars of a few solar masses were formed from the primordial gas and that some of them caused thermonuclear explosions when the mass of their degenerate carbon-oxygen cores increased to the Chandrasekhar limit without experiencing efficient mass loss.

We study the evolution of the ultra-relativistic shock wave in a plane-parallel atmosphere adjacent to a vacuum and the subsequent breakout phenomenon. When the density distribution has a power law with the distance from the surface, there is a self-similar motion of the fluid before and after the shock emergence. The time evolution of the Lorentz factor of the shock front is assumed to follow a power law when the time is measured from the moment at which the shock front reaches the surface. The power index is found to be determined by the condition for the flow to extend through a critical point. The energy spectrum of the ejected matter as a result of the shock breakout is derived and its dependence on the strength of the explosion is also deduced. The results are compared with the self-similar solution for the same problem with non-relativistic treatment.

Explosions of type Ic supernovae (SNe Ic) are investigated using a relativistic hydrodynamic code to study roles of their outermost layers of the ejecta in light element nucleosynthesis through spallation reactions as a possible mechanism of the "primary" process. We have confirmed that the energy distribution of the outermost layers with a mass fraction of only 0.001 % follows the empirical formula proposed by previous work when the explosion is furious. In such explosions, a significant fraction of the ejecta ($>$0.1 % in mass) have the energy greater than the threshold energy for spallation reactions. On the other hand, it is found that the outermost layers of ejecta become more energetic than the empirical formula would predict when the explosion energy per unit ejecta mass is smaller than $\sim 1.3\times 10^{51}{ergs/}\Msun$. As a consequence, it is necessary to numerically calculate explosions to estimate light element yields from SNe Ic. The usage of the empirical formula would overestimate the yields by a factor of $\gtsim 3$ for energetic explosions such as SN 1998bw and underestimate the yields by a similar factor for less energetic explosions like SN 1994I. The yields of light elements Li, Be, and B (LiBeB) from SNe Ic are estimated by solving the transfer equation of cosmic rays originated from ejecta of SNe Ic and compared with observations.

The origin of low alpha/Fe ratios in some metal-poor stars, so called low-alpha stars, is discussed. It is found that most of low-alpha stars in the Galaxy are on the main-sequence. This strongly suggests that these stars suffered from external pollution. It is also found that the abundance ratios Zn/Fe of low-alpha stars both in the Galaxy and in dwarf spheroidal galaxies are lower than the average value of Galactic halo stars whereas damped Ly alpha absorbers have higher ratios. This implies that some low-alpha stars accreted matter depleted from gas onto dust grains. To explain the features in these low-alpha stars, we have proposed that metal-poor stars harboring planetary systems are the origin of these low-alpha stars. Stars engulfing a small fraction of planetesimals enhance the surface content of Fe to exhibit low alpha/Fe ratios on their surfaces while they are on the main-sequence, because dwarfs have shallow surface convection zones where the engulfed matter is mixed. After the stars leave the main-sequence, the surface convection zones become deeper to reduce the enhancement of Fe. Eventually, when the stars ascend to the tip of the red giant branch, they engulf giant planets to become low-alpha stars again as observed in dwarf spheroidal galaxies. We predict that low-alpha stars with low Mn/Fe ratios harbor planetary systems.

The external pollution of the first stars in the Galaxy is investigated. The first stars were born in clouds composed of the pristine gas without heavy elements. These stars accreted gas polluted with heavy elements while they still remained in the cloud. As a result, it is found that they exhibit a distribution with respect to the surface metallicity. We have derived the actual form of this distribution function. This metallicity distribution function strongly suggests that the recently discovered most metal-deficient star HE0107-5240 with [Fe/H]=-5.3 was born as a metal-free star and accreted gas polluted with heavy elements. Thus the heavy elements such as Fe in HE0107-5240 must have been supplied from supernovae of later generations exploding inside the cloud in which the star had been formed. The elemental abundance pattern on the surface of stars suffering from such an external pollution should not be diverse but exhibit the average pattern of numerous supernovae. Future observations for a number of metal-deficient stars with [Fe/H]<-5 will be able to prove or disprove this external pollution scenario. Other possibilities to produce a star with this metallicity are also discussed.

The star formation history of the globular cluster Omega Centauri is investigated in the context of an inhomogeneous chemical evolution model in which supernovae induce star formation. The proposed model explains recent observations for Omega Cen stars, and divides star formation into three epochs. At the end of the first epoch, ~ 70% of the gas was expelled by supernovae. AGB stars then supplied s-process elements to the remaining gas during the first interval of ~300 Myr. This explains the observed sudden increase in Ba/Fe ratios in Omega Cen stars at [Fe/H] ~ -1.6. Supernovae at the end of the second epoch were unable to expel the gas. Eventually, Type Ia supernovae initiated supernova-induced star formation, and remaining gas was stripped when the cluster passed through the newly formed disk of the Milky Way. The formation of Omega Cen is also discussed in the framework of globular cluster formation triggered by cloud-cloud collisions. In this scenario, the relative velocity of clouds in the collision determines the later chemical evolution in the clusters. A head-on collision of proto-cluster clouds with a low relative velocity would have converted less than 1% of gas into stars and promoted the subsequent chemical evolution by supernova-driven star formation. This is consistent with present observed form of Omega Cen. In contrast, the other Galactic globular clusters are expected to have formed from more intense head-on collisions, and the resultant clouds would have been too thin for supernovae to accumulate enough gas to form the next generation of stars. This explains the absence of chemical evolution in these other globular clusters.

The unique elemental abundance pattern of the carbon-rich stars CS29498-043 and CS22949-037 is characterized by a large excess of magnesium and silicon in comparison with iron. This excess is investigated in the context of a supernova-induced star formation scenario, and it is concluded that these stars were born from the matter swept up by supernova remnants containing little iron and that such supernovae are similar to the least-luminous SNe ever observed, SNe 1997D and 1999br. Comparison of the observed abundance pattern in iron-group elements of subluminous supernovae with those of other supernovae leads to an intriguing implication for explosion, nucleosynthesis, and mixing in supernovae. The observed invariance of these ratios can not be accounted for by a spherically symmetric supernova model.

Stellar abundance pattern of neutron-capture elements such as barium is used as a powerful tool to infer how star formation proceeded in dwarf spheroidal (dSph) galaxies. It is found that the abundance correlation of barium with iron in stars belonging to dSph galaxies orbiting the Milky Way, i.e., Draco, Sextans, and Ursa Minor have a feature similar to the barium-iron correlation in Galactic metal-poor stars. The common feature of these two correlations can be realized by our inhomogeneous chemical evolution model based on the supernova-driven star formation scenario if dSph stars formed from gas with a velocity dispersion of ~26 km/s. This velocity dispersion together with the stellar luminosities strongly suggest that dark matter dominated dSph galaxies. The tidal force of the Milky Way links this velocity dispersion with the currently observed value <10 km/s by stripping the dark matter in dSph galaxies. As a result, the total mass of each dSph galaxy is found to have been originally ~25 times larger than at present. Our inhomogeneous chemical evolution model succeeds in reproducing the stellar [Fe/H] distribution function observed in Sextans. In this model, supernovae immediately after the end of the star formation epoch can expel the remaining gas over the gravitational potential of the dSph galaxy.

The origin of nucleosynthesis products of rapid neutron capture reactions (the r-process) is a longstanding astrophysical problem. Recent analyses of elemental abundances for extremely metal-poor stars shed light on the elemental abundances of individual supernovae. Comparison of the abundance distributions of some extremely metal-poor stars with those of the best-observed supernova SN 1987A clearly indicates that the overabundances of barium and strontium found in SN 1987A that have been ascribed to the slow neutron capture process must be results of r-process nucleosynthesis. The mass of freshly synthesized barium in SN 1987A is estimated to be 6x10^-6 solar mass based on the observed surface abundance and detailed hydrodynamical models for this supernova. These new findings lead to the conclusion that 20 solar mass stars, one of which is the progenitor star of SN 1987A, are the predominant production sites for r-process elements in the Galaxy and the r-process element donors for notable neutron-capture-rich giant stars, CS22892-052 and CS31082-001.

We investigate the chemical evolution of odd-numbered elements such as sodium (Na) and aluminum (Al) during the early epochs of the Galactic halo with the use of a model that reproduces the observed box-shaped distribution of extremely metal-poor stars in the [Na, Al/Mg] versus [Mg/H] plane. Our model is constructed under the assumptions that those stars retain the elemental abundance patterns produced by individual Type II supernovae (SNe), and that the yields of the odd elements depend on the initial metallicity, z, of their SN progenitors. As a result, recent abundance determinations that clarify how the [Na, Al/Mg] ratios of field stars have evolved to the solar values enable us to deduce how the yields of these odd elements depend on z. The observed trends in these abundances, in particular the very large scatter (over 1 dex in [Al/Mg]) requires that the Al yield scales as m_Al proportional to z^0.6 for [Mg/H]<-1.8, while the observed [Na/Mg] trend requires that the Na yield scales as m_Na proportional to z^-0.4 for [Mg/H]<-1.8 and m_Na proportional to z^0.4 for [Mg/H]>-1.8. It is found that the predicted frequency distribution of stars in the [Na/Mg] versus [Mg/H] diagram is very sensitive to the assumed form of the primordial IMF, and that its slope is steeper than the Salpeter IMF. The necessity to match the observed abundance patterns of odd elements and the frequency distribution of extremely metal-poor stars should provide useful constraints on nucleosynthesis calculations of metal-free massive stars as well as on theories of their formation.

The long-term evolution of supernova remnants (SNRs) in the primordial interstellar medium (ISM) with an inhomogeneous structure is calculated to investigate metal enrichment of the primordial gas. For this purpose, we have constructed a parallel 3-D hydrodynamics code incorporating the radiative cooling and self-gravity. The self-gravity and radiative cooling develop the inhomogeneous structure of the ISM from a small perturbation with a power-law spectrum. The resultant density ranges from 0.5 cm$^{-3}$ to 10$^6$ cm $^{-3}$. Calculations for a supernova (SN) with the progenitor mass of 20 solar mass are performed as the first step of a series of our study. It is found from the results that a single SN distributes some of newly synthesized heavy elements into a dense filament of the ISM with densities ranging from 100 to 10$^4$ cm$^{-3}$ depending on where the SN explodes. Thus, the metallicity [Mg/H] of the dense filaments polluted by the SN ejecta becomes $-2.7\pm0.5$. From these filaments, the first Population II stars will form. This value is in accordance with the previous analytical work (Shigeyama & Tsujimoto 1998) with accuracy of $\sim 0.3$ dex.

We suggest that if the astrophysical site for r-process nucleosynthesis in the early Galaxy is confined to a narrow mass range of Type II supernova (SN II) progenitors, with a lower mass limit of Mms = 20 Msun, a unique feature in the observed distribution of [Ba/Mg] vs.[Mg/H] for extremely metal-poor stars can be adequately reproduced. We associate this feature, a bifurcation of the observed elemental ratios into two branches in the Mg abundance interval -2.7 < [Mg/H] < -2.3, with two distinct processes. The first branch, which we call the ``y''-branch, is associated with the production of Ba and Mg from individual massive supernovae. We conclude that SNe II with Mms = 20 Msun are the dominant source of r-process nucleosynthesis in the early Galaxy. An SN-induced chemical evolution model with this Mms-dependent Ba yield creates the y-branch, reflecting the different nucleosynthesis yields of [Ba/Mg] for each SN II with Mms > 20 Msun. The second branch, which we call the ``i''-branch, is associated with the elemental abundance ratios of stars which were formed in the dense shells of the interstellar medium swept up by SNe II with Mms < 20 Msun that do not synthesize r-process elements, and applies to stars with observed Mg abundances in the range [Mg/H] < -2.7. The Ba abundances in these stars reflect those of the interstellar gas at the (later) time of their formation. The existence of a [Ba/Mg] i-branch strongly suggests that SNe II which are associated with stars of progenitor mass Mms < 20 Msun are infertile sources for the production of r-process elements. We predict the existence of this i-branch for other r-process elements, such as europium (Eu), to the extent that their production site is in common with Ba.

A model for Galactic chemical evolution, driven by supernova-induced star formation, is formulated and used to examine the nature of the Galactic halo at early epochs. In this model, new stars are formed following each supernova event, thus their abundance pattern is determined by the combination of heavy elements ejected from the supernova itself and those elements which are already present in the interstellar gas swept up by the supernova remnant. The end result is a prediction of large scatter in the abundance ratios among low-metallicity stars, reflecting a different nucleosynthesis yield for each Type II supernova with a different progenitor mass. Formation of new stars is terminated when supernova remnants sweep up too little gas to form shells. We show from calculations based on the above scenario that (i) the observed [Fe/H] distribution for the Galactic halo field stars can be reproduced without effectively decreasing the heavy-element yields from Type II supernovae by some manipulation required by previous models (e.g., via mass loss from the early Galaxy, or later mixing with ``pristine'' hydrogen clouds), (ii) the large observed scatter in the abundance ratio [Eu/Fe] for the most metal-poor stars can also be reproduced, and (iii) the frequency distribution of stars in the [Eu/Fe]-[Fe/H] plane can be predicted. Our model suggests that the probability of identifying essentially metal-free stars (Population III) in the local halo is around one in 10^{3-4}, provided that star formation in the halo is confined to individual gas clouds with mass of 10^{6-7} Msun and that the initial mass function of metal-free stars is not significantly different from the Salpeter mass function.

The supernova yields of several heavy elements including alpha-, iron-group, and r-process elements are obtained as a function of the mass of their progenitor main-sequence stars M_ms from the abundance patterns of extremely metal-poor stars with a procedure recently proposed by Shigeyama and Tsujimoto (1998). The ejected masses of alpha- and iron-group elements increase with M_ms, whereas more Eu is ejected from supernovae with lower M_ms. For these several heavy elements, it is shown that the average abundance ratios weighted by the Salpeter initial mass function coincide with the ratios observed in stars with -2<[Fe/H]<-1 within 0.1 dex. It follows that the correlations of stellar abundance ratios with the metallicity are twofold. One is the abundance ratios for [Fe/H] <-2.5 imprinted by the nucleosynthesis in individual supernovae on the timescale 10^7yr and the other for [Fe/H]>-2 results from the mixing of the products from a whole site of the nucleosynthesis, taking place on the timescale longer than 10^9yr.

Using results of nucleosynthesis calculations for theoretical core-collapse supernova models with various progenitor's masses, it is shown that abundance patterns of C, Mg, Si, Ca, and H seen in extremely metal-deficient stars with [Fe/H] < -2.5 follow those seen in the individual first generation supernova remnants (SNRs). This suggests that most of the stars with [Fe/H] < -2.5 were made from individual supernova (SN) events. To obtain the ratio of heavy elements to hydrogen, a formula is derived to estimate the mass of hydrogen swept up by a SNR when it occurs in the interstellar matter with the primordial abundances. We use [Mg/H] to indicate the metallicities instead of [Fe/H]. The metallicities [Mg/H] predicted from these SNRs range from ~-4 to ~-1.5 and the mass of Mg in a SN is well correlated with its progenitor's mass. Thus the observed [Mg/H] in an extremely metal deficient star has a correspondence to the progenitor's mass. A larger [Mg/H] corresponds to a larger progenitor's mass. Therefore, so called `age-metallicity relation' does not hold for stars with [Fe/H] < -2.5. In contrast, the [Mg/Fe] ratios in the theoretical SNRs have a different trend from those in extremely metal-deficient stars. It is also shown that the observed trend of [Mg/Fe] can predict the Fe yield of each SN given the correspondence of [Mg/H] to the progenitor's mass. The Fe yields thus obtained are consistent with those derived from SN light curve analyses. This indicates that there is still a problem in modelling a core-collapse supernova at its beginning of explosion or mass cut.

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.

Rayleigh-Taylor (R-T) instabilities in the explosion of SN 1993J are investigated by means of two-dimensional hydrodynamical simulations. It is found that the extent of mixing is sensitive to the progenitor's core mass and the envelope mass. Because the helium core mass (3 - 4 \ms) is smaller than that of SN 1987A, R-T instabilities at the He/C+O interfaces develop to induce a large scale mixing in the helium core, while the instability is relatively weak at the H/He interface due to the small envelope mass. The predicted abundance distribution, in particular the amount of the \ni~ mixing, is compared with those required in the theoretical light curves and the late time optical spectra. This enables us to specify the progenitor of SN 1993J in some detail.