Source author record

Koh Takahashi

Koh Takahashi 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

16works
4topics
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

16 published item(s)

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

Impact of binary interactions on the diffuse supernova neutrino background

Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%-20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. The enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.

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

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

Modeling of Magneto-Rotational Stellar Evolution I. Method and first applications

While magnetic fields have long been considered to be important for the evolution of magnetic non-degenerate stars and compact stars, it has become clear in recent years that actually all of the stars are deeply affected. This is particularly true regarding their internal angular momentum distribution, but magnetic fields may also influence internal mixing processes and even the fate of the star. We propose a new framework for stellar evolution simulations, in which the interplay between magnetic field, rotation, mass loss, and changes in the stellar density and temperature distributions are treated self-consistently. For average large-scale stellar magnetic fields which are symmetric to the axis of rotation of the star, we derive 1D evolution equations for the toroidal and poloidal components from the mean-field MHD equation by applying Alfven's theorem, and a conservative form of the angular momentum transfer due to the Lorentz force is formulated. We implement our formalism into a numerical stellar evolution code and simulate the magneto-rotational evolution of 1.5 M$_\odot$ stars. The Lorentz force aided by the $Ω$ effect imposes torsional Alfven waves propagating through the magnetized medium, leading to near-rigid rotation within the Alfven timescale. Our models with different initial spins and B-fields can reproduce the main observed properties of Ap/Bp stars. Calculations continued to the red-giant regime show a pronounced core-envelope coupling, which reproduces the core and surface rotation periods determined by asteroseismic observations.

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.

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

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

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.

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.