We study thermal evolution of isolated neutron stars in scalar-tensor theories for the first time. Whether the rapid cooling due to the direct Urca process occurs or not is an interesting question in the viewpoint of the temperature observation of isolated neutron stars. The cooling effect is typically influenced by the proton fraction and the central density. If a fifth force is mediated due to modification of gravity, the relation between the central density and mass of neutron stars differs from one in general relativity, and the cooling curve is also naively expected to be varied. We find that an unscreened fifth force near the surface of neutron stars changes mass-central density relation, and the direct Urca process can be triggered even for neutron stars with smaller mass. This might enable us to test gravitational theories with neutron star cooling.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

We investigate the significance of isomers on a kilonova associated with neutron star mergers (NSMs) for the first time. We calculate the evolution of abundances and energy generation rates ($\dot{\epsilon}_{\rm nuc}$) of ejecta from NSMs, taking into account $\beta^-$ decay through isomers. We find that for ejecta with electron fraction ($Y_e$) of $0.2-0.3$, $\dot{\epsilon}_{\rm nuc}$ is appreciably changed from those without isomers, due to the large change in timing of $\beta^-$ decay through the isomeric states of ${}^{123, 125, 127} \rm Sn$ and ${}^{128} \rm Sb$. In particular, the effects of the isomers on $\dot{\epsilon}_{\rm nuc}$ are prominent for ejecta of $Y_e \sim 0.25$, which could emit a fraction of an early, blue component of a kilonova observed in GW170817. When the excitation from a ground state to its isomeric state is important, the isomers of ${}^{129} \rm Sb$ and ${}^{129, 131} \rm Te$ also cause additional and appreciable change in $\dot{\epsilon}_{\rm nuc}$. Furthermore, we demonstrate that larger amounts of lanthanide-free ejecta result in a better fit of an observed light curve of the kilonova in GW170817 if the isomers are taken into account.

We investigate X-ray bursts during the thermal evolution of an accreting neutron star which corresponds to the X-ray burster GS\ 1826-24. Physical quantities of the neutron star are included using an equation of state below and above the nuclear matter density. We adopt an equation of state and construct an approximate network that saves the computational time and calculates nuclear energy generation rates accompanying the abundance evolutions. The mass and radius of the neutron star are got by solving the stellar evolution equations from the center to the surface which involve necessary information such as the nuclear energy generation in accreting layers, heating from the crust, and neutrino emissions inside the stellar core. We reproduce the light curve and recurrence time of the X-ray burst from GS 1826-24 within the standard deviation of 1$\sigma$ for the assumed accretion rate, metallicity, and equation of state. It is concluded that the observed recurrence time is consistent with the theoretical model having metallicity of the initial CNO elements $Z_{\rm CNO}$ = 0.01. We suggest that the nuclear reaction rates responsible for the $rp$-process should be examined in detail, because the rates may change the shape of the light curve and our conclusion.

Whether fast cooling processes occur or not is crucial for the thermal evolution of neutron stars. In particular, the threshold of the direct Urca process, which is one of the fast cooling processes, is determined by the interior proton fraction $Y_p$, or the nuclear symmetry energy. Since recent observations indicate the small radius of neutron stars, a low value is preferred for the symmetry energy. In this study, simulations of neutron star cooling are performed adopting three models for equation of state (EoS): Togashi, Shen, and LS220 EoSs. The Togashi EoS has been recently constructed with realistic nuclear potentials under finite temperature, and found to account for the small radius of neutron stars. As a result, we find that, since the direct Urca process is forbidden, the neutron star cooling is slow with use of the Togashi EoS. This is because symmetry energy of Togashi EoS is lower than those of other EoSs. Hence, in order to account for observed age and surface temperature of isolated neutron stars (INS) with use of the Togashi EoS, other fast cooling processes are needed regardless of the surface composition.

We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(\nu)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $\nu$ self-interaction ($\nu$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li and heavy isotopes $^{92}$Nb, $^{98}$Tc and $^{138}$La are reduced by a factor of 1.5-2.0 by the $\nu$-SI. In contrast, $^{11}$B is relatively insensitive to the $\nu$-SI. We find that the abundance ratio of heavy to light nucleus, $^{138}$La/$^{11}$B, is sensitive to the neutrino mass hierarchy, and the normal mass hierarchy is more likely to be consistent with the solar abundances.

The neutrino process ($\nu$-process) for the production of 7Li and 11B in core-collapse supernovae (SNe) is extensively investigated. Initial abundances of s-nuclei and other physical conditions are derived from an updated calculation of the SN 1987A progenitor. The nuclear reaction network including neutrino reactions is constructed with the variable order Bader-Deuflhard integration method. We find that yields of 7Li and 11B significantly depend on the stellar metallicity while they are independent of the weak s-process during the stellar evolution. When the metallicity is high, there are more neutron absorbers, i.e., 56Fe, 14N (from initial CNO nuclei), and 54Fe, and the neutron abundance is small during the $\nu$-process. Since 7Be is predominantly destroyed via 7Be(n,p)7Li, a change in the neutron abundance results in different 7Be yields. Then, the calculated yield ratio 7Li/11B=0.93 for the solar metallicity is larger than that for the SN 1987A 7Li/11B=0.80 by 16 % in the inverted mass hierarchy case. We analyze contributions of respective reactions as well as abundance evolution, and clarify the $\nu$-process of 7Li and 11B.

We reinvestigate effects of neutrino oscillations on the production of 7Li and 11B in core-collapse supernovae (SNe). During the propagation of neutrinos from the proto-neutron star, their flavors change and the neutrino reaction rates for spallation of 12C and 4He are affected. In this work corrected neutrino spallation cross sections for 4He and 12C are adopted. Initial abundances involving heavy s-nuclei and other physical conditions are derived in a new calculation of the SN 1987A progenitor in which effects of the progenitor metallicity are included. A dependence of the SN nucleosynthesis and final yields of 7Li and 11B on the neutrino mass hierarchy are shown in several stellar locations. In the normal hierarchy case, the charged current reaction rates of electron neutrinos are enhanced, and yields of proton-rich nuclei, along with 7Be and 11C, are increased. In the inverted hierarchy case, the charged current reaction rates of electron antineutrinos are enhanced, and yields of neutron-rich nuclei, along with 7Li and 11B, are increased. We find that variation of the metallicity modifies the yields of 7Li, 7Be, 11B, and 11C. This effect is caused by changes in the neutron abundance during SN nucleosynthesis. Therefore, accurate calculations of Li and B production in SNe should take into account the metallicity of progenitor stars.

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

MAXI J0556-332 is the hottest transient accreting neutron star at the beginning of its quiescence. A theoretical model with crustal heating indicates that an additional shallow heat source of Qshallow > 6 MeV per accreted nucleon is required in the shallow outer crust with respect to the deeper star crust by considering the observed decline in accretion rate at the end of outburst. However, the physical source of this shallow heating is still unclear. In the present investigation, we performed stellar evolutionary calculations, adopting the effects of outburst behavior of the accretion rate. As a consequence, we find that the quiescent light curve of MAXI J0556-332 can be well explained by the nuclear energy generation due to the hot CNO cycle.

There are many kinds of models which describe the dynamics of dark energy (DE). Among all we adopt an equation of state (EoS) which varies as a function of time. We adopt Markov Chain Monte Carlo method to constrain the five parameters of our models. As a consequence, we can show the characteristic behavior of DE during the evolution of the universe. We constrain the EoS of DE with use of the avairable data of gamma-ray bursts and type Ia supernovae (SNe Ia) concerning the redshift-luminosity distance relations. As a result, we find that DE is quintessence-like in the early time and phantom-like in the present epoch or near the future, where the change occurs rather rapidly at $z\sim0.3$.

Quiescent X-ray luminosities are presented in low mass X-ray binaries with use of evolutionary calculations. The calculated luminosities are compared with observed ones in terms of timeaveraged mass accretion rate. It is shown that neutrino emission by strong pion condensation can explain quiescent X-ray luminosity of SAX J1808.4-3658 and we do not need direct Urca processes concerning nucleons and/or hyperons.

We show a scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star with high effective temperature, and it is consistent with the cooling without exotic phases. The Cas A observation also gives the mass range of $M \geq 1.5 M_\odot$. It may conflict with the current cooling scenarios of compact stars that heavy stars show rapid cooling. We include the effect of the color superconducting (CSC) quark matter phase on the thermal evolution of compact stars. We assume the gap energy of CSC quark phase is large ($\Delta \gtrsim \mathrm{10 MeV}$), and we simulate the cooling of compact stars. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency.

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

We investigate effects of triple-$\alpha$ and $^{12}\rm C(\alpha,\gamma) ^{16}O$ reaction rates on the production of supernova yields for a massive star of 25 $M_{\odot}$. We combine the reaction rates to examine the rate dependence, where the rates are considered to cover the possible variation of the rates based on experiments on the earth and theories. We adopt four combinations of the reaction rates from two triple-$\alpha$ reaction rates and two $^{12}\rm C(\alpha,\gamma)^{16}O$ ones. First, we examine the evolution of massive stars of 20 and 25 $M_{\odot}$ whose helium cores correspond to helium stars of 6 and 8 $M_{\odot}$, respectively. While the 25 $M_{\odot}$ stars evolve to the presupernova stages for all combinations of the reaction rates, evolutionary paths of the 20 $M_{\odot}$ stars proceed significantly different way for some combinations, which are unacceptable for progenitors of supernovae. Second, we perform calculations of supernova explosions within the limitation of spherical symmetry and compare the calculated abundance ratios with the solar system abundances. We can deduce some constraints to the reaction rates. As the results, a conventional rate is adequate for a triple-$\alpha$ reaction rate and a rather higher value of the reaction rate within the upper limit for the experimental uncertainties is favorable for a $^{12}\rm C(\alpha,\gamma)^{16}O$ rate.

Comparing the latest observed abundances of 4He and D, we make a ?2 analysis to see whether it is possible to reconcile primordial nucleosynthesis using up-to-date nuclear data of NACRE II and the mean-life of neutrons. If we adopt the observational data of ${}^4_{}$He by Izotov et al., we find that it is impossible to get reasonable concordance against the standard Big-Bang nucleosynthesis. However, including degenerate neutrinos, we have succeeded in obtaining consistent constraints between the neutrino degeneracy and the baryon-to-photon ratio from detailed comparison of calculated abundances with the observational data of ${}^4_{}$He and D: the baryon-to-photon ratio in units of $10^{-10}$ is found to be in the range 6.02 < $\eta^{}_{10}$ < 6:54 for the specified parameters of neutrino degeneracy.

Based on a scenario of the inhomogeneous big-bang nucleosynthesis (IBBN), we investigate the detailed nucleosynthesis that includes the production of heavy elements beyond Li-7. From the observational constraints on light elements of He4 and D for the baryon-to-photon ratio given by WMAP, possible regions found on the plane of the volume fraction of the high density region against the ratio between high- and low-density regions. In these allowed regions, we have confirmed that the heavy elements beyond Fe can be produced appreciably, where p- and/or r-process elements are produced well simultaneously compared to the solar system abundances. We suggest that recent observational signals such as He4 overabundance in globular clusters and high metallicity abundances in quasars could be partly due to the results of IBBN. Possible implications are given for the formation of the first generation stars

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.

We present a new scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star that has high effective temperature, and it is consistent with the cooling without exotic phases. The observation also gives the mass range of $M \geqslant 1.5 M_\odot$, which may conflict with the current plausible cooling scenario of compact stars. There are some cooled compact stars such as Vela or 3C58, which can be barely explained by the minimal cooling scenario, which includes the neutrino emission by nucleon superfluidity (PBF). Therefore, we invoke the exotic cooling processes, where a heavier star cools faster than lighter one. However, the scenario seems to be inconsistent with the observation of Cas A. Therefore, we present a new cooling scenario to explain the observation of Cas A by constructing models that include a quark color superconducting (CSC) phase with a large energy gap; this phase appears at ultrahigh density region and reduces neutrino emissivity. In our model, a compact star has CSC quark core with a low neutrino emissivity surrounded by high emissivity region made by normal quarks. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency without considering nucleon superfluidity. Furthermore, we show that our scenario is consistent with the recent observations of the effective temperature of Cas A during the last 10 years, including nucleon superfluidity.

We investigate the nucleosynthesis in a massive star of 70 M_solar with solar metallicity in the main sequence stage. The helium core mass after hydrogen burning corresponds to 32 M_solar. Nucleosynthesis calculations have been performed during the stellar evolution and the jetlike supernova explosion of a collapsar model, where the weak s-, p-, and r-processes are taken into account. We confirm that s-elements of 60 < A < 90 are highly overproduced relative to the solar abundances in the hydrostatic nucleosynthesis. During oxygen burning, p-elements of A > 90 are produced via photodisintegrations of seed s-elements. However, the produced p-elements are disintegrated in later stages except for ^{180}Ta. In the explosive nucleosynthesis, elements of 90 < A < 160 are significantly overproduced relative to the solar values owing to the r-process. Only heavy p-elements (N > 50) are overproduced via the p-process. Compared with the previous study of r-process nucleosynthesis calculations in the collapsar model of 40 M_solar by Fujimoto et al. 2007, 2008, our jet model cannot contribute to the third peak of the solar r-elements and intermediate p-elements. Averaging the overproduction factors over the progenitor masses with the use of Salpeter's IMF, we suggest that the 70 M_solar star could contribute to the solar weak s-elements of 60 < A < 90 and neutron-rich elements of 90 < A < 160. We confirm the primary synthesis of light p-elements in the ejected matter of high peak temperature. The ejected matter has [Sr/Eu] \sim -0.4, which is different from that of a typical r-process-enriched star CS22892-052 ([Sr/Eu] \sim -1). We find that Sr-Y-Zr isotopes are primarily synthesized in the explosive nucleosynthesis in a similar process of the primary production of light p-elements, which has been considered as one of the sites of a lighter element primary process (LEPP).

We investigate the cosmic thermal evolution with a vacuum energy which decays into photon at the low-redshift. We assume that the vacuum energy is a function of the scale factor that increases toward the early universe. We put on the constraints using recent observations of both type Ia supernovae (SNIa) by Union-2 compilation and the cosmic microwave background (CMB) temperature at the range of the redshift 0.01 < z < 3. From SNIa, we find that the effects of a decaying vacuum energy on the cosmic expansion rate should be very small but could be possible for z < 1.5. On the other hand, we obtain the severe constraints for parameters from the CMB temperature observations. Although the temperature can be still lower than the case of the standard cosmological model, it should only affect the thermal evolution at the early epoch.

We investigate the possibility of the r-process during the magnetohydrohynamical explosion of supernova in a massive star of 13 solar mass with the effects of neutrinos induced. We adopt five kinds of initial models which include properties of rotation and the toroidal component of the magnetic field . The simulations which succeed the explosions are limitted to a concentrated magnetic field and strong differential rotation. Low $Y_{e}$ ejecta produce heavy elements and the third peak can be reprocuced. However, the second peak is low because $Y_{e}$ distribution as a function of radius is steep and ejecta corresponding to middle $Y_{e}$ is very few.

The effects of a new triple-$\alpha$ reaction rate (OKK rate) on the helium flash of a helium accreting neutron star in a binary system have been investigated. Since the ignition points determine the properties of a thermonuclear flash of type I X-ray bursts, we examine the cases of different accretion rates, $dM/dt (\dot{M})$, of helium from $3\times10^{-10} M_{\odot} \rm yr^{-1}$ to $3\times10^{-8} M_{\odot} \rm yr^{-1}$, which could cover the observed accretion rates. We find that for the cases of low accretion rates, nuclear burnings are ignited at the helium layers of rather low densities. As a consequence, helium deflagration would be triggered for all cases of lower accretion rate than $\dot{M}\simeq 3\times10^{-8} M_{\odot} \rm yr^{-1}$. We find that OKK rate could be barely consistent with the available observations of the X-ray bursts on the helium accreting neutron star. However this coincidence is found to depend on the properties of crustal heating and the neutron star model.We suggest that OKK rate would be reduced by a factor of $10^{2-3}$ for $10^8$ K in the range of the observational errors.

Effects of a new triple-{\alpha} reaction rate on the s-process during the evolution of a massive star of 25 M\odot are investigated for the first time. Although the s-process in massive stars has been believed to be established with only minor change, we find that the s-process with use of the new rate during the core He-burning is very inefficient compared to the case with the previous triple-{\alpha} rate. However, the difference of the overproduction is found to be largely compensated by the subsequent C-burning. Since the s-process in massive stars has been attributed so far to the neutron irradiation during core He-burning, our finding reveals the importance of C-burning for the s-process during the evolution of massive stars.

We investigate explosive nucleosynthesis in a non-rotating 15$M_\odot$ star with solar metallicity that explodes by a neutrino-heating supernova (SN) mechanism aided by both standing accretion shock instability (SASI) and convection. To trigger explosions in our two-dimensional hydrodynamic simulations, we approximate the neutrino transport with a simple light-bulb scheme and systematically change the neutrino fluxes emitted from the protoneutron star. By a post-processing calculation, we evaluate abundances and masses of the SN ejecta for nuclei with the mass number $\le 70$ employing a large nuclear reaction network. Aspherical abundance distributions, which are observed in nearby core-collapse SN remnants, are obtained for the non-rotating spherically-symmetric progenitor, due to the growth of low-mode SASI. Abundance pattern of the supernova ejecta is similar to that of the solar system for models whose masses ranges $(0.4-0.5) \Ms$ of the ejecta from the inner region ($\le 10,000\km$) of the precollapse core. For the models, the explosion energies and the \nuc{Ni}{56} masses are $ \simeq 10^{51} \rm erg$ and $(0.05-0.06) \Ms$, respectively; their estimated baryonic masses of the neutron star are comparable to the ones observed in neutron-star binaries. These findings may have little uncertainty because most of the ejecta is composed by matter that is heated via the shock wave and has relatively definite abundances. The abundance ratios for Ne, Mg, Si and Fe observed in Cygnus loop are well reproduced with the SN ejecta from an inner region of the $15\Ms$ progenitor.

The Brans-Dicke model with a variable cosmological term ($BD\Lambda$) has been investigated with use of the coupling constant of $\omega=10^4$. Parameters inherent in this model are constrained from comparison between Big Bang nucleosynthesis and the observed abundances. Furthermore, the magnitude redshift ($m-z$) relations are studied for $BD\Lambda$ with and without another constant cosmological term in a flat universe. Observational data of Type Ia Supernovae are used in the redshift range of $0.01<z<2$. It is found that our model with energy density of the constant cosmological term with the value of 0.7 can explain the SNIa observations, though the model parameters are insensitive to the $m-z$ relation.

Effects of a new triple-alpha reaction rate on the ignition of carbon-oxygen white dwarfs accreting helium in a binary systems have been investigated. The ignition points determine the properties of a thermonuclear explosion of a Type Ia supernova. We examine the cases of different accretion rates of helium and different initial masses of the white dwarf, which was studied in detail by Nomoto. We find that for all cases from slow to intermediate accretion rates, nuclear burnings are ignited at the helium layers. As a consequence, carbon deflagration would be triggered for the lower accretion rate compared to that of $dM/dt\simeq 4\times10^{-8} M_{\odot} \rm yr^{-1}$ which has been believed to the lower limit of the accretion rate for the deflagration supernova. Furthermore, off-center helium detonation should result for intermediate and slow accretion rates and the region of carbon deflagration for slow accretion rate is disappeared.

This is a reply report to astro-ph/0604264. We studied heavy element production in high baryon density region in early universe astro-ph/0507439. However it is claimed in astro-ph/0604264 that small scale but high baryon density region contradicts observations for the light element abundance or in order not to contradict to observations high density region must be so small that it cannot affect the present heavy element abundance. In this paper we study big bang nucleosynthesis in high baryon density region and show that in certain parameter spaces it is possible to produce enough amount of heavy element without contradiction to CMB and light element observations.

We constrain the thermal evolution of the universe with a decaying cosmological term by using the method of the analysis for the Wilkinson Microwave Anisotropy Probe (WMAP) observation data. The cosmological term is assumed to be a function of the scale factor that increases toward the early universe, and the radiation energy density is lower compared to that in the model with the standard cosmological "constant" (LCDM). The decrease in the radiation density affects the thermal history of the universe; e.g. the photon decoupling occurs at higher-z compared to the case of the standard LCDM model. As a consequence, a decaying cosmological term affects the cosmic microwave background anisotropy. Thanks to the Markov-chain Monte Carlo method, we compare the angular power spectrum in the decaying LCDM model with the CMB data, and we get severe constraints on parameters of the model.

We have made detailed calculations of the composition of magnetically driven jets ejected from collapsars, or rapidly rotating massive stars, based on long-term magnetohydrodynamic simulations of their core collapse with various distributions of magnetic field and angular momentum before collapse. We follow the evolution of the abundances of about 4000 nuclides from the collapse phase to the ejection phase and through the jet generation phase using a large nuclear reaction network. We find that the r-process successfully operates only in energetic jets (> 1e51 ergs), such that U and Th are synthesized abundantly, even when the collapsar has a relatively weak magnetic field (1e10 G) and a moderately rotating core before the collapse. The abundance patterns inside the jets are similar to those of the r-elements in the solar system. About 0.01-0.06 Msun neutron-rich, heavy nuclei are ejected from a collapsar with energetic jets. The higher energy jets have larger amounts of Ni56, varying from 0.00037 to 0.06Msun. Less energetic jets, which eject small amounts of Ni56, could induce a gamma-ray burst (GRB) a supernova, such as GRB 060505 or GRB 060614. Considerable amounts of r-elements are likely to be ejected from GRBs with hypernovae, if both the GRB and hypernova are induced by jets that are driven near the black hole.

We investigate the evolution of a universe with a decaying cosmological term (vacuum energy) that is assumed to be a function of the scale factor. In this model, while the cosmological term increases to the early universe, the radiation energy density is lower than the model with the cosmological "constant". We find that the effects of the decaying cosmological term on the expansion rate at the redshift z<2 is negligible. However, the decrease in the radiation density affects on the thermal history of the universe; e.g. the photon decoupling occurs at higher $z$ compared to the case of the standard \Lambda CDM model. As a consequence, a decaying cosmological term affects on the cosmic microwave background anisotropy. We show the angular power spectrum in D\Lambda CDM model and compare with the Wilkinson Microwave Anisotropy Probe (WMAP) data.

We perform two-dimensional, magnetohydrodynamical core-collapse simulations of massive stars accompanying the QCD phase transition. We study how the phase transition affects the gravitational waveforms near the epoch of core-bounce. As for initial models, we change the strength of rotation and magnetic fields. Particularly, the degree of differential rotation in the iron core (Fe-core) is changed parametrically. As for the microphysics, we adopt a phenomenological equation of state above the nuclear density, including two parameters to change the hardness before the transition. We assume the first order phase transition, where the conversion of bulk nuclear matter to a chirally symmetric quark-gluon phase is described by the MIT bag model. Based on these computations, we find that the phase transition can make the maximum amplitudes larger up to $\sim$ 10 percents than the ones without the phase transition. On the other hand, the maximum amplitudes become smaller up to $\sim$ 10 percents owing to the phase transition, when the degree of the differential rotation becomes larger. We find that even extremely strong magnetic fields $\sim 10^{17}$ G in the protoneutron star do not affect these results.

We perform two-dimensional, axisymmetric, magnetohydrodynamic simulations of the collapse of a rotating star of 40 Msun and in the light of the collapsar model of gamma-ray burst. Considering two distributions of angular momentum, up to \sim 10^{17} cm^2/s, and the uniform vertical magnetic field, we investigate the formation of an accretion disk around a black hole and the jet production near the hole. After material reaches to the black hole with the high angular momentum, the disk is formed inside a surface of weak shock. The disk becomes in a quasi-steady state for stars whose magnetic field is less than 10^{10} G before the collapse. We find that the jet can be driven by the magnetic fields even if the central core does not rotate as rapidly as previously assumed and outer layers of the star has sufficiently high angular momentum. The magnetic fields are chiefly amplified inside the disk due to the compression and the wrapping of the field. The fields inside the disk propagate to the polar region along the inner boundary near the black hole through the Alfv{\'e}n wave, and eventually drive the jet. The quasi-steady disk is not an advection-dominated disk but a neutrino cooling-dominated one. Mass accretion rates in the disks are greater than 0.01 Msun/sec with large fluctuations. The disk is transparent for neutrinos. The dense part of the disk, which locates near the hole, emits neutrino efficiently at a constant rate of < 8 \times 10^{51} erg/s. The neutrino luminosity is much smaller than those from supernovae after the neutrino burst.

We investigate synthesis of heavy elements in a collapsar. We have calculated detailed composition of magnetically driven jets ejected from a collapsar, which is based on long-term, magneto-hydrodynamic simulations of a rapidly rotating massive star of 40Msun during core collapse. We follow evolution of abundances of about 4000 nuclides from the collapse phase to the ejection phase through the jet generation phase with use of two large nuclear reaction networks. We find that the r-process successfully operates in the jets, so that U and Th are synthesized abundantly when the progenitor has large magnetic field of 10^{12} G and rapidly rotating core. Abundance pattern inside the jets is similar compared to that of r-elements in the solar system. Heavy neutron-rich nuclei \sim 0.01Msun can be ejected from the collapsar. The detailed abundances depend on nuclear properties of mass model, beta-decay rate, and fission, for nuclei near the neutron drip line. Furthermore, we find that p-nuclei are produced without seed nuclei: not only light p-nuclei, such as Se74, Kr78, Sr84, and Mo92, but also heavy p-nuclei, In113, Sn115, and La138, can be abundantly synthesized in the jets. The amounts of p-nuclei in the ejecta are much greater than those in core-collapse supernovae (SNe). In particular, Mo92, In113, Sn115, and La138 deficient in the SNe, are significantly produced in the ejecta.

We present a new astrophysical site of the big bang nucleosynthesis (BBN) that are very peculiar compared with the standard BBN. Some models of the baryogenesis suggest that very high baryon density regions were formed in the early universe. On the other hand, recent observations suggest that heavy elements already exist in high red-shifts and the origin of these elements become a big puzzle. Motivated by these, we investigate BBN in very high baryon density regions. BBN proceeds in proton-rich environment, which is known to be like the p-process. However, by taking very heavy nuclei into account, we find that BBN proceeds through both the p-process and the r-process simultaneously. P-nuclei such as 92Mo, 94Mo, 96Ru, 98Ru whose origin is not well known are also synthesized.

We investigate the big-bang nucleosynthesis in a Brans-Dicke model with a varying $\Lambda$ term using the Monte-Carlo method and likelihood analysis. It is found that the cosmic expansion rate differs appreciably from that of the standard model. The produced abundances of $^4$He, D, and barely Li are consistent with the observed ones within the uncertainties in nuclear reaction rates when the baryon to photon ratio $\eta = (5.47 - 6.64)\times 10^{-10}$, which is in agreement with the value deduced from WMAP.

We present a scenario of formation of strange stars due to spin-down of {\it rapidly rotating} neutron stars left after supernova explosions . By assuming a process where the total baryon mass is conserved but the angular momentum is lost due to emission of gravitational waves and/or the magnetic braking, we find that the transition from rapidly rotating neutron stars to slowly rotating strange stars is possible; a large amount of energy $\sim 10^{53} ergs $ could be released. The liberated energy might become a new energy source for a delayed explosion of supernova. Furthermore, our scenario suggests that the supernova associated with gamma-ray bursts could become candidates for targets in the future observation of gravitational waves.

We investigate $r$-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to the case of the spherical explosion, jet-like explosion due to the combined effects of the rotation and magnetic field lowers the electron fraction significantly inside the layers above the iron core. We find that the ejected material of low electron fraction responsible for the $r$-process comes out from the silicon rich layer of the presupernova model. This leads to the production up to the third peak in the solar $r$-process elements. We examine whether the fission affects the $r$-process paths by using the full nuclear reaction network with both the spontaneous and $\beta$-delayed fission included. Moreover, we pay particular attention how the mass formula affects the $r$-process peaks with use of two mass formulae. It is found that both formulae can reproduce the global abundance pattern up to the third peak though detailed distributions are rather different. We point out that there are variations in the $r$-process nucleosynthesis if the MHD effects play an important role in the supernova explosion.

We investigate nucleosynthesis inside both a gamma-ray burst accretion disk and a wind launched from an inner region of the disk using one-dimensional models of the disk and wind and a nuclear reaction network. Far from a central black hole, the composition of accreting gas is taken to be that of an O-rich layer of a massive star before core collapse. We find that the disk consists of five layers characterized by dominant elements: O16, Si28, Fe54 (and Ni56), He4, and nucleons, and the individual layers shift inward with keeping the overall profiles of compositions as the accretion rate decreases. Ni56 are abundantly ejected through the wind from the inner region of the disk with the electron fraction \simeq 0.5. In addition to iron group, elements heavier than Cu, in particular Cu63 and Zn64, are massively produced through the wind. Various neutron-rich nuclei can be also produced in the wind from neutron-rich regions of the disk, though the estimated yields have large uncertainties.

We investigate p-process nucleosynthesis in a supercritical accretion disk around a compact object of 1.4 M_solar, using the self-similar solution of an optically thick advection dominated flow. Supercritical accretion is expected to occur in a supernova with fallback material accreting onto a new-born compact object. It is found that appreciable amounts of p-nuclei are synthesized via the p-process in supernova-driven supercritical accretion disks (SSADs) when the accretion rate m_dot = M_dot c^2/(16 L_Edd) >10^5, where L_Edd is the Eddington luminosity. Abundance profiles of p-nuclei ejected from SSADs have similar feature to those of the oxygen/neon layers in Type II supernovae when the abundance of the fallback gas far from the compact object is that of the oxygen/neon layers in the progenitor. The overall abundance profile is in agreement with that of the solar system. Some p-nuclei, such as Mo, Ru, Sn, and La, are underproduced in the SSADs as in Type II supernovae. If the fallback gas is mixed with a small fraction of proton through Rayleigh-Taylor instability during the explosion, significant amounts of Mo92 are produced inside the SSADs. Ru96 and La138 are also produced when the fallback gas contains abundant proton though the overall abundance profile of p-nuclei is rather different from that of the solar system. The p-process nucleosynthesis in SSADs contributes to chemical evolution of p-nuclei, in particular Mo92, if several percents of fallback matter are ejected via jets and/or winds.

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

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

Nucleosynthesis in ONeMg novae has been investigated with the wide ranges of three parameters, i.e., the white dwarf mass, the envelope mass at ignition, and the initial composition. A quasi-analytic one-zone approach is used with an up-to-date nuclear reaction network. The nucleosynthesis results show correlation with the peak temperatures or the cooling timescales during outbursts. Among the combinations of white dwarf and envelope masses which give the same peak temperature, the explosion is more violent for a lower white dwarf mass owing to its smaller gravitational potential. Comparison of the nucleosynthesis results with observations implies that at least two-third of the white dwarf masses for the observed ONeMg novae are $\simeq 1.1 M_\odot$, which are significantly lower than estimated by previous hydrodynamic studies but consistent with the observations of V1974 Cyg. Moreover, the envelope masses derived from the comparison are $\gtrsim 10^{-4} M_\odot$, which are in good agreement with the ejecta masses estimated from observations but significantly higher than in previous hydrodynamic studies. With such a low mass white dwarf and a high mass envelope, the nova can produce interesting amounts of $\gamma$-ray emitters $^7$Be, $^{22}$Na, and $^{26}$Al. We suggest that V1974 Cyg has produced $^{22}$Na as high as the upper limit derived from the COMPTEL survey. In addition, a non-negligible part of the Galactic $^{26}$Al may originate from ONeMg novae, if not the major contributors. Both the future INTEGRAL survey for these $\gamma$-ray emitters and abundance estimates derived from ultraviolet, optical, and near infrared spectroscopies will impose a severe constraint on the current nova models.

We have performed 1-dimensional calculations for explosive nucleosynthesis in collapse-driven supernova and investigated its sensitivity to the initial form of the shock wave. We have found the tendency that the peak temperature becomes higher around the mass cut if the input energy is injected more in the form of kinetic energy rather than internal energy. Then, the mass cut becomes larger, and, as a result, neutron-rich matter is less included in the ejecta; this is favorable for producing the observational data compared with a previous model. Our results imply that the standard method to treat various processes for stellar evolution, such as convection and electron capture during the silicon burning stage, are still compatible with the calculation of explosive nucleosynthesis.

The large abundance ratio of $^{44}Ti/^{56}Ni$ in Cas A is puzzling. In fact, the ratio seems to be larger than the theoretical constraint derived by Woosley & Hoffman (1991). However, this constraint is obtained on the assumption that the explosion is spherically symmetric, whereas Cas A is famous for the asymmetric form of the remnant. Recently, Nagataki et al. (1997) calculated the explosive nucleosynthesis of axisymmetrically deformed collapse-driven supernova. They reported that the ratio of $^{44}Ti/^{56}Ni$ was enhanced by the stronger alpha-rich freezeout in the polar region. In this paper, we apply these results to Cas A and examine whether this effect can explain the large amount of $^{44}Ti$ and the large ratio of $^{44}Ti/^{56}Ni$. We demonstrate that the conventional spherically symmetric explosion model can not explain the $^{44}$Ti mass produced in Cas A if its lifetime is shorter than $\sim$ 80 years and the intervening space is transparent to the gamma-ray line from the decay of $^{44}$Ti. On the other hand, we show the axisymmetric explosion models can solve the problem. We expect the same effect from a three dimensionally asymmetric explosion, since the stronger alpha-rich freezeout will also occur in that case in the region where the larger energy is deposited.

We have estimated the dispersion of the inferred relative frequencies of Type-Ia and Type- II supernovae ($N_{Ia}/ N_{II}$) in our Galaxy by fitting the numerical results of supernova nucleosynthesis to the solar-system abundances. The ratio $N_{Ia}/ N_{II}$ is estimated to be 0.056--0.14, which is consistent with the observation, if the model of Woosley and Weaver (1995, WW95) is adopted for Type-II supernovae (SNe II). On the other hand, the upper limit of $N_{Ia}/ N_{II}$ becomes too large for the model of Hashimoto (1995, Ha95). However, Ha95 can fit the solar values better than WW95 as far as the abundant nuclei are concerned. These results mean that Ha95 can reproduce the main nuclei of the solar-system well and that WW95 can reproduce the solar values over a wide mass-number range. We also note that $N_{Ia}/ N_{II}$ tends to become smaller if the delayed detonation model is adopted for Type-Ia supernovae (SNe Ia). The dispersion of $N_{Ia}/ N_{II}$ obtained in the present investigation is larger than that concluded by Tsujimoto et al. (1995), which will have an influence on the estimate of the average life time of SNe Ia's progenitors and/or the star-formation rate history in our galaxy.

This review concentrates on nucleosynthesis processes in general and their applications to massive stars and supernovae. A brief initial introduction is given to the physics in astrophysical plasmas which governs composition changes. We present the basic equations for thermonuclear reaction rates and nuclear reaction networks. The required nuclear physics input for reaction rates is discussed, i.e. cross sections for nuclear reactions, photodisintegrations, electron and positron captures, neutrino captures, inelastic neutrino scattering, and beta-decay half-lives. We examine especially the present state of uncertainties in predicting thermonuclear reaction rates, while the status of experiments is discussed by others in this volume (see M. Wiescher). It follows a brief review of hydrostatic burning stages in stellar evolution before discussing the fate of massive stars, i.e. the nucleosynthesis in type II supernova explosions (SNe II). Except for SNe Ia, which are explained by exploding white dwarfs in binary stellar systems (which will not be discussed here), all other supernova types seem to be linked to the gravitational collapse of massive stars (M$>$8M$_\odot$) at the end of their hydrostatic evolution. SN1987A, the first type II supernova for which the progenitor star was known, is used as an example for nucleosynthesis calculations. Finally, we discuss the production of heavy elements in the r-process up to Th and U and its possible connection to supernovae.

Explosive nucleosynthesis under the axisymmetric explosion in Type II supernova has been performed by means of two dimensional hydrodynamical calculations. We have compared the results with the observations of SN 1987A. Our chief findings are as follows: (1) $^{44}Ti$ is synthesized so much as to explain the tail of the bolometric light curve of SN 1987A. We think this is because the alpha-rich freezeout takes place more actively under the axisymmetric explosion. (2) $^{57}Ni$ and $^{58}Ni$ tend to be overproduced compared with the observations. However, this tendency relies strongly on the progenitor's model. We have also compared the abundance of each element in the mass number range $A= 16-73$ with the solar values. We have found three outstanding features. (1) For the nuclei in the range $A=16-40$, their abundances are insensitive to the initial form of the shock wave. This insensitivity is favored since the spherical calculations thus far can explain the solar system abundances in this mass range. (2) There is an enhancement around A=45 in the axisymmetric explosion compared with the spherical explosion fairly well. In particular, $^{44}Ca$, which is underproduced in the present spherical calculations, is enhanced significantly. (3) In addition, there is an enhancement around A=65. This tendency does not rely on the form of the mass cut but of the initial shock wave. This enhancement may be the problem of the overproduction in this mass range, although this effect would be relatively small since Type I supernovae are chiefly responsible for this mass number range.

Presupernova evolution and explosive nucleosynthesis in massive stars for main-sequence masses from 13 $M_\odot$ to 70 $M_\odot$ are calculated. We examine the dependence of the supernova yields on the stellar mass, $^{12}C(\alpha, \gamma) ^{16}O}$ rate, and explosion energy. The supernova yields integrated over the initial mass function are compared with the solar abundances.

Oxygen to iron abundance ratios of metal-poor stars provide information on nucleosynthesis yields from massive stars which end in Type II supernova explosions. Using a standard model of chemical evolution of the Galaxy we have reproduced the solar neighborhood abundance data and estimated the oxygen and iron yields of genuine SN II origin. The estimated yields are compared with the theoretical yields to derive the relation between the lower and upper mass limits in each generation of stars and the IMF slope. Independently of this relation, we furthermore derive the relation between the lower mass limit and the IMF slope from the stellar mass to light ratio in the solar neighborhood. These independent relations unambiguously determine the upper mass limit of $m_u=50 \pm 10 M_sun$ and the IMF slope index of 1.3 - 1.6 above 1 M_sun. This upper mass limit corresponds to the mass beyond which stars end as black holes without ejecting processed matter into the interstellar medium. We also find that the IMF slope index below 0.5 M_sun cannot be much shallower than 0.8.