Source author record

Takashi Yoshida

Takashi Yoshida appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2023arXiv

A Novel Estimation Method for Temperature of Magnetic Nanoparticles Dominated by Brownian Relaxation Based on Magnetic Particle Spectroscopy

This paper presents a novel method for estimating the temperature of magnetic nanoparticles (MNPs) based on AC magnetization harmonics of MNPs dominated by Brownian relaxation. The difference in the AC magnetization response and magnetization harmonic between the Fokker-Planck equation and the Langevin function was analyzed, and we studied the relationship between the magnetization harmonic and the key factors, such as Brownian relaxation time, temperature, magnetic field strength, core size and hydrodynamic size of MNPs, excitation frequency, and so on. We proposed a compensation function for AC magnetization harmonic with consideration of the key factors and the difference between the Fokker-Planck equation and the Langevin function. Then a temperature estimation model based on the compensation function and the Langevin function was established. By employing the least squares algorithm, the temperature was successfully calculated. The experimental results show that the temperature error is less than 0.035 K in the temperature range from 310 K to 320 K. The temperature estimation model is expected to improve the performance of the magnetic nanoparticle thermometer and be applied to magnetic nanoparticle-mediated hyperthermia.

preprint2022arXiv

Merger rate density of binary black holes through isolated Population I, II, III and extremely metal-poor binary star evolution

We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitational wave (GW) observations. Notably, Pop III and EMP ($< 10^{-2}$ $Z_\odot$) binary stars yield most of the pair instability (PI) mass gap events with $m_1 = 65$--$130$ $M_\odot$. Pop III binary stars contribute more to the PI mass gap events with increasing redshift, and all the PI mass gap events have the Pop III origin at redshifts $\gtrsim 8$. Our result can be assessed by future GW observations in the following two points. First, there are no binary BHs with $m_1=100$--$130$ $M_\odot$ in our result, and thus the $m_1$ distribution should suddenly drop in the range of $m_1=100$--$130$ $M_\odot$. Second, the PI mass gap event rate should increase toward higher redshift up to $\sim 11$, since those events mainly originate from the Pop III binary stars. We find that the following three assumptions are needed to reproduce the current GW observations: a top-heavy stellar initial mass function and the presence of close binary stars for Pop III and EMP binary stars, and inefficient convective overshoot in the main-sequence phase of stellar evolution. Without any of the above, the number of PI mass gap events becomes too low to reproduce current GW observations.

preprint2022arXiv

Neutrino emission from the collapse of $\sim 10^4$ $M_\odot$ population III supermassive stars

We calculate the neutrino signal from Population III supermassive star collapse using a neutrino transfer code originally developed for core collapse supernovae and massive star collapse. Using this code, we are able to investigate the supermassive star mass range thought to undergo neutrino trapping ($\sim 10^4$ M$_\odot$), a mass range which has been neglected by previous works because of the difficulty of neutrino transfer. For models in this mass range, we observe a neutrino-sphere with a large radius and low density compared to typical massive star neutrino-spheres. We calculate the neutrino light-curve emitted from this neutrino-sphere. The resulting neutrino luminosity is significantly lower than the results of a previous analytical model. We briefly discuss the possibility of detecting a neutrino burst from a supermassive star or the neutrino background from many supermassive stars and conclude that the former is unlikely with current technology, unless the SMS collapse is located as close as 1 Mpc, while the latter is also unlikely even under very generous assumptions. However, the supermassive star neutrino background is still of interest as it may serve as a source of noise in proposed dark matter direct detection experiments.

preprint2022arXiv

Stability analysis of supermassive primordial stars: a new mass range for general relativistic instability supernovae

Observed supermassive black holes in the early universe have several proposed formation channels, in part because most of these channels are difficult to probe. One of the more promising channels, the direct collapse of a supermassive star, has several possible probes including the explosion of a helium-core supermassive star triggered by a general relativistic instability. We develop a straightforward method for evaluating the general relativistic radial instability without simplifying assumptions and apply it to population III supermassive stars taken from a post Newtonian stellar evolution code. This method is more accurate than previous determinations and it finds that the instability occurs earlier in the evolutionary life of the star. Using the results of the stability analysis, we perform 1D general relativistic hydrodynamical simulations and we find two general relativistic instability supernovae fueled by alpha capture reactions as well as several lower mass pulsations, analogous to the puslational pair instability process. The mass range for the events (2.6-3.0 $\times 10^4$ ${\rm M}_\odot$) is lower than had been suggested by previous works (5.5 $\times 10^4$ ${\rm M}_\odot$) because the instability occurs earlier in the star's evolution. The explosion may be visible to, among others, JWST, while the discovery of the pulsations opens up additional possibilities for observation.

preprint2022arXiv

Stellar models of Betelgeuse constrained using observed surface conditions

We study stellar models for Betelgeuse using the HR diagram and surface abundances as observational constraints. Previous studies on Betelgeuse have not systematically investigated the surface abundances, but we believe they can be impacted by, and thus be used as an observational constraint for various parameters such as initial mass, rotation, and overshoot scheme. We investigate stellar models with varying initial mass as they evolve past the main sequence, and we examine the red supergiant (RSG) properties in detail. For each mass, we vary the initial rotation up to $\sim$ 300 $km\:s^{-1}$, and test two different overshoot parameters. Overall, the acceptable initial mass range is 12 to 25 $M_\odot$, but for non-rotating models only, the range is decreased to 15 to 24 $M_\odot$. Also for rotating models, we find that $v/v_{\rm K} = 0.3$ is the upper limit for initial rotation, as more rapidly rotating models are unable to fit to Betelgeuse's surface abundances as an RSG. In addition, we report two possibilities for the current stage of evolution, core helium burning or core carbon burning and beyond. We find that certain 17 $M_\odot$ models could fit to both stages. Finally, we discuss the implications of our results in the context of merger scenarios which have been suggested as a mechanism to attain the observed surface velocity of Betelgeuse.

preprint2021arXiv

Possible Detection of X-Ray Emitting Circumstellar Material in the Synchrotron-Dominated Supernova Remnant RX J 1713.7-3946

We report on a discovery of an X-ray emitting circumstellar material knot inside the synchrotron dominant supernova remnant (SNR) RX J1713.7-3946. This knot was previously thought to be a Wolf-Rayet star (WR 85), but we realized that it is in fact $\sim$40$^{\prime\prime}$ away from WR 85, indicating no relation to WR 85. We performed high-resolution X-ray spectroscopy with the Reflection Grating Spectrometer (RGS) on board XMM-Newton. The RGS spectrum clearly resolves a number of emission lines, such as N Ly$α$, O Ly$α$, Fe XVIII, Ne X, Mg XI, and Si XIII. The spectrum can be well represented by an absorbed thermal emission model with a temperature of $k_{\rm B}T_{\rm e} = 0.65\pm 0.02$ keV. The elemental abundances are obtained to be ${\rm N/H} = 3.5\pm 0.8{\rm \left(N/H\right)_{\odot}}$, ${\rm O/H} = 0.5\pm0.1{\rm \left(O/H\right)_{\odot}}$, ${\rm Ne/H} = 0.9\pm0.1{\rm \left(Ne/H\right)_{\odot}}$, ${\rm Mg/H} = 1.0\pm0.1{\rm \left(Mg/H\right)_{\odot}}$, ${\rm Si/H} = 1.0\pm0.2{\rm \left(Si/H\right)_{\odot}}$, and ${\rm Fe/H} = 1.3\pm0.1{\rm \left(Fe/H\right)_{\odot}}$. The enhanced N abundance with others being about the solar values allows us to infer that this knot is circumstellar material ejected when the progenitor star evolved into a red supergiant. The abundance ratio of N to O is obtained to be $\rm N/O = 6.8_{-2.1}^{+2.5}\left(N/O\right)_{\odot}$. By comparing this to those in outer layers of red supergiant stars expected from stellar evolution simulations, we estimate the initial mass of the progenitor star to be $15\, \rm M_{\odot} \lesssim \rm M \lesssim 20\, \rm M_{\odot}$.

preprint2021arXiv

Three-dimensional Hydrodynamics Simulations of Precollapse Shell Burning in the Si- and O-rich Layers

We present 3D hydrodynamics simulations of shell burning in two progenitors with zero-age main-sequence masses of 22 and 27 $M_{\odot}$ for $\sim$65 and 200 s up to the onset of gravitational collapse, respectively. The 22 and 27 $M_{\odot}$ stars are selected from a suite of 1D progenitors. The former and the latter have an extended Si- and O-rich layer with a width of $\sim$10$^9$ cm and $\sim$5$\times 10^9$ cm, respectively. Our 3D results show that turbulent mixing occurs in both of the progenitors with the angle-averaged turbulent Mach number exceeding $\sim$0.1 at the maximum. We observe that an episodic burning of O and Ne, which takes place underneath the convection bases, enhances the turbulent mixing in the 22 and 27 $M_\odot$ models, respectively. The distribution of nucleosynthetic yields is significantly different from that in 1D simulations, namely, in 3D more homogeneous and inhomogeneous in the radial and angular direction, respectively. By performing a spectrum analysis, we investigate the growth of turbulence and its role of material mixing in the convective layers. We also present a scalar spherical harmonics mode analysis of the turbulent Mach number. This analytical formula would be helpful for supernova modelers to implement the precollapse perturbations in core-collapse supernova simulations. Based on the results, we discuss implications for the possible onset of the perturbation-aided neutrino-driven supernova explosion.

preprint2020arXiv

A Sub-Solar Metallicity Progenitor for Cassiopeia A, the remnant of a Type IIb Supernova

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

preprint2020arXiv

Fitting formulae for evolution tracks of massive stars under extreme metal poor environments for population synthesis calculations and star cluster simulations

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

preprint2020arXiv

Neutrino halo effect on collective neutrino oscillation in iron core-collapse supernova model of a 9.6 $M_{\odot}$ star

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

preprint2020arXiv

The Final Fate of Supermassive $M \sim 5 \times 10^4 \; M_\odot$ Pop III Stars: Explosion or Collapse?

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

preprint2020arXiv

Theoretical prediction of presupernova neutrinos and their detection

Almost 30 years have passed since the successful detection of supernova neutrinos from SN 1987A. In the last decades, remarkable progress has been made in neutrino detection technique, through which it may be possible to detect neutrinos from a new source, pre-supernova (pre-SN) neutrinos. They are emitted from a massive star prior to core bounce. Because neutrinos escape from the core freely, they carry information about the stellar physics directly. Pre-SN neutrinos may play an important role in verifying our understanding of stellar evolution for massive stars. Observations of pre-SN neutrinos, moreover, may serve as an alarm regarding a supernova explosion a few days in advance if the progenitor is located in our vicinity, enabling us to observe the next galactic supernova. In this review, we summarize the current status of pre-SN neutrino studies from both of the theoretical and observational points of view.

preprint2019arXiv

Matter Mixing in Aspherical Core-collapse Supernovae: Three-dimensional Simulations with Single Star and Binary Merger Progenitor Models for SN 1987A

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

preprint2016arXiv

Alpha-constrained QSE Nucleosynthesis in High-entropy and Fast-expanding Material

We investigate the nucleosynthesis process in high-entropy ($s/k_{\rm B}\gtrsim100$) and very fast-expanding ($τ_{\rm exp}\sim10^{-3}\ {\rm s}$) materials. In such a material with the electron fraction near 0.5, an interesting nucleosynthesis process occurs. In this process, the abundance distribution of heavy-nuclei of $A>100$ achieve quasi-statistical equilibrium (QSE) at high temperature and the abundances are frozen at the end of the nucleosynthesis. We explain this abundance distribution using the "alpha-constrained QSE" abundances formulated in this paper. We demonstrate that this nucleosynthesis would occur in neutrino-driven winds from massive proto-neutron stars in hypernovae, where $A\sim140$ $p$-nuclei are synthesized.

preprint2016arXiv

Mass ejection by pulsational pair instability in very massive stars and implications for luminous supernovae

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

preprint2016arXiv

Presupernova neutrino events relating to the final evolution of massive stars

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

preprint2015arXiv

Exact and approximate expressions of energy generation rates and their impact on the explosion properties of Pair Instability Supernovae

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

preprint2015arXiv

Neutrino-driven explosions of ultra-stripped type Ic supernovae generating binary neutron stars

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

preprint2015arXiv

Nucleosynthesis in neutrino-driven winds in hypernovae

We investigate the nucleosynthesis in the neutrino-driven winds blown off from a $3M_\odot$ massive proto-neutron star (mPNS) temporarily formed during the collapse of a $100M_\odot$ star. Such mPNSs would be formed in hypernovae. We construct steady and spherically symmetric wind solutions. We set large neutrino luminosities of $\sim 10^{53}\ {\rm erg\ s^{-1}}$ and average energies of electron neutrinos and antineutrinos in the ranges of $ε_{ν_e}=9-16\ {\rm MeV}$ and $ε_{\barν_e}=11-18\ {\rm MeV}$ based on a recent numerical relativity simulation. The wind solutions indicate much shorter temperature-decrease timescale than that of the winds from ordinary PNSs and, depending on $ε_ν$, the winds can be both neutron-rich and proton-rich. In the neutron-rich wind, the $r$-process occurs and the abundance distribution of a fiducial wind model of the mPNS gives an approximate agreement with the abundance pattern of metal-poor weak $r$ star HD~122563, although the third-peak elements are produced only when the $\barν_e$ energy is much larger than the $ν_e$ energy. In the proton-rich wind, the strong $νp$-process occurs and $A>100$ nuclides are synthesized. The synthesized nuclei can be neutron-rich in some cases because the large neutrino luminosity of the mPNS supplies a sufficient amount of neutrons.

preprint2015arXiv

Pre-supernova neutrino emissions from ONe cores in the progenitors of core-collapse supernovae: are they distinguishable from those of Fe cores?

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

preprint2014arXiv

Stellar Yields of Rotating First Stars. I. Yields of Weak Supernovae and Abundances of Carbon-enhanced Hyper Metal Poor Stars

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

preprint2013arXiv

Evolution of progenitors for electron capture supernovae

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

preprint2013arXiv

Type Ic core-collapse supernova explosions evolved from very massive stars

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

preprint2012arXiv

Massive Star Evolution and Nucleosynthesis -Lower End of Fe-Core Collapse Supernova Progenitors and Remnant Neutron Star Mass Distribution-

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

preprint2011arXiv

A Progenitor for the Extremely Luminous Type Ic Supernova 2007bi

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

preprint2011arXiv

Beta Decays of Isotones with Neutron Magic Number of N=126 and R-process Nucleosynthesis

Beta decays of the isotones with N=126 are studied by shell model calculations taking into account both the Gamow-Teller (GT) and first-forbidden (FF) transitions. The FF transitions are found to be important to reduce the half-lives, by nearly twice to several times, from those by the GT contributions only. Possible implications of the short half-lives of the waiting point nuclei on the r-process nucleosynthesis during the supernova explosions are discussed. A slight shift of the third peak of the element abundances in the r-process toward higher mass region is found.

preprint2011arXiv

Evolution of Fluorine Abundances with the ν-Process

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

preprint2011arXiv

Invariants of Collective Neutrino Oscillations

We consider the flavor evolution of a dense neutrino gas by taking into account both vacuum oscillations and self interactions of neutrinos. We examine the system from a many-body perspective as well as from the point of view of an effective one-body description formulated in terms of the neutrino polarization vectors. We show that, in the single angle approximation, both the many-body picture and the effective one-particle picture possess several constants of motion. We write down these constants of motion explicitly in terms of the neutrino isospin operators for the many-body case and in terms of the polarization vectors for the effective one-body case. The existence of these constants of motion is a direct consequence of the fact that the collective neutrino oscillation Hamiltonian belongs to the class of Gaudin Hamiltonians. This class of Hamiltonians also includes the (reduced) BCS pairing Hamiltonian describing superconductivity. We point out the similarity between the collective neutrino oscillation Hamiltonian and the BCS pairing Hamiltonian. The constants of motion manifest the exact solvability of the system. Borrowing the well established techniques of calculating the exact BCS spectrum, we present exact eigenstates and eigenvalues of both the many-body and the effective one-particle Hamiltonians describing the collective neutrino oscillations. For the effective one-body case, we show that spectral splits of neutrinos can be understood in terms of the adiabatic evolution of some quasi-particle degrees of freedom from a high density region where they coincide with flavor eigenstates to the vacuum where they coincide with mass eigenstates. We write down the most general consistency equations which should be satisfied by the effective one-body eigenstates and show that they reduce to the spectral split consistency equations for the appropriate initial conditions.

preprint2011arXiv

Supernova neutrino signals by liquid Argon detector and neutrino magnetic moment

We study electron-neutrino and electron-antineutrino signals from a supernova with strong magnetic field detected by a 100 kton liquid Ar detector. The change of neutrino flavors by resonant spin-flavor conversions, matter effects, and neutrino self-interactions are taken into account. Different neutrino signals, characterized by neutronization burst event and the total event numbers of electron-neutrinos and electron-antineutrinos, are expected with different neutrino oscillation parameters and neutrino magnetic moment. Observations of supernova neutrino signals by a 100 kton liquid Ar detector would constrain oscillation parameters as well as neutrino magnetic moment in either normal and inverted mass hierarchies.

preprint2010arXiv

Big Bang Nucleosynthesis with long-lived strongly interacting relic particles

We study effects of relic long-lived strongly interacting massive particles (X particles) on big bang nucleosynthesis (BBN). The X particle is assumed to have existed during the BBN epoch, but decayed long before detected. The interaction strength between an X and a nucleon is assumed to be similar to that between nucleons. Rates of nuclear reactions and beta decay of X-nuclei are calculated, and the BBN in the presence of neutral charged X^0 particles is calculated taking account of captures of X^0 by nuclei. As a result, the X^0 particles form bound states with normal nuclei during a relatively early epoch of BBN leading to the production of heavy elements. Constraints on the abundance of X^0 are derived from observations of primordial light element abundances. Particle models which predict long-lived colored particles with lifetimes longer than about 200 s are rejected. This scenario prefers the production of 9Be and 10B. There might, therefore, remain a signature of the X particle on primordial abundances of those elements. Possible signatures left on light element abundances expected in four different models are summarized.

preprint2010arXiv

Boron Synthesis in Type Ic Supernovae

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.

preprint2010arXiv

New results on catalyzed BBN with a long-lived negatively-charged massive particle

It has been proposed that the apparent discrepancies between the inferred primordial abundances of 6Li and 7Li and the predictions of big bang nucleosynthesis (BBN) can be resolved by the existence of a negatively-charged massive unstable supersymmetric particle (X-) during the BBN epoch. Here, we present new BBN calculations with an X- particle utilizing an improved nuclear reaction network including captures of nuclei by the particle, nuclear reactions and beta-decays of normal nuclei and nuclei bound to the X- particles (X-nuclei), and new reaction rates derived from recent rigorous quantum many-body dynamical calculations. We find that this is still a viable model to explain the observed 6Li and 7Li abundances. However, contrary to previous results, neutral X-nuclei cannot significantly affect the BBN light-element abundances. We also show that with the new rates the production of heavier nuclei is suppressed and there is no signature on abundances of nuclei heavier than Be in the X--particle catalyzed BBN model as has been previously proposed. We also consider the version of this model whereby the X- particle decays into the present cold dark matter. We analyze the this paradigm in light of the recent constraints on the dark-matter mass deduced from the possible detected events in the CDMS-II experiment. We conclude that based upon the inferred range for the dark-matter mass, only X- decay via the weak interaction can achieve the desired 7Li destruction while also reproducing the observed 6Li abundance.

preprint2008arXiv

The X^- Solution to the ^6Li and ^7Li Big Bang Nucleosynthesis Problems

The $^6$Li abundance observed in metal poor halo stars appears to exhibit a plateau as a function of metallicity similar to that for $^7$Li, suggesting a big bang origin. However, the inferred primordial abundance of $^6$Li is $\sim$1000 times larger than that predicted by standard big bang nucleosynthesis for the baryon-to-photon ratio inferred from the WMAP data. Also, the inferred $^7$Li primordial abundance is 3 times smaller than the big bang prediction. We here describe in detail a possible simultaneous solution to both the problems of underproduction of $^6$Li and overproduction of $^7$Li in big bang nucleosynthesis. This solution involves a hypothetical massive, negatively-charged leptonic particle that would bind to the light nuclei produced in big bang nucleosynthesis, but would decay long before it could be detected. We consider only the $X$-nuclear reactions and assume that the effect of decay products is negligible, as would be the case if lifetime were large or the mass difference between the charged particle and its daughter were small. An interesting feature of this paradigm is that, because the particle remains bound to the existing nuclei after the cessation of the usual big bang nuclear reactions, a second longer epoch of nucleosynthesis can occur among $X$-nuclei. We confirm that reactions in which the hypothetical particle is transferred can occur that greatly enhance the production of $^6$Li while depleting $^7$Li. We also identify a new reaction that destroys large amounts of $^7$Be, and hence reduces the ultimate $^7$Li abundance. Thus, big-bang nucleosynthesis in the presence of these hypothetical particles, together with or without an event of stellar processing, can simultaneously solve the two Li abundance problems.