Researcher profile

Kei Kotake

Kei Kotake contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2026arXiv

Gravitational-Wave Signatures of Nonstandard Neutrino Properties in Collapsing Stellar Cores

We present a novel multi-messenger approach for probing nonstandard neutrino properties through the detection of gravitational waves (GWs) from collapsing stellar cores and associated supernova explosions. We show that neutrino flavor conversion inside the proto-neutron star (PNS), motivated by physics Beyond the Standard Model (BSM), can significantly boost PNS convection. This effect leads to large-amplitude GW emission over a wide frequency range during an otherwise relatively quiescent GW phase shortly after core bounce. Such a signal provides a promising new avenue for exploring nonstandard neutrino phenomena and other BSM physics impacting PNS convection.

preprint2022arXiv

3D template-based $Fermi$-LAT constraints on the diffuse supernova axion-like particle background

Axion-like particles (ALPs) may be abundantly produced in core-collapse (CC) supernovae (SNe), hence the cumulative signal from all past SN events can create a diffuse flux peaked at energies of about 25~MeV. We improve upon the modeling of the ALPs flux by including a set of CC SN models with different progenitor masses, as well as the effects of failed CC SNe -- which yield the formation of black holes instead of explosions. Relying on the coupling strength of ALPs to photons and the related Primakoff process, the diffuse SN ALP flux is converted into gamma rays while traversing the magnetic field of the Milky Way. The spatial morphology of this signal is expected to follow the shape of the Galactic magnetic field lines. We make use of this via a template-based analysis that utilizes 12 years of $Fermi$-LAT data in the energy range from 50 MeV to 500 GeV. In our benchmark case of the realization of astrophysical and cosmological parameters, we find an upper limit of $g_{aγ} \lesssim 3.76\times10^{-11}\;\mathrm{GeV}^{-1}$ at 95$\%$ confidence level for $m_a \ll 10^{-11}$ eV, while we find that systematic deviations from this benchmark scenario induce an uncertainty as large as about a factor of two. Our result slightly improves the CAST bound, while still being a factor of six (baseline scenario) weaker than the SN1987A gamma-ray burst limit.

preprint2022arXiv

Magnetic support for neutrino-driven explosion of 3D non-rotating core-collapse supernova models

The impact of the magnetic field on postbounce supernova dynamics of non-rotating stellar cores is studied by performing three-dimensional magnetohydrodynamics simulations with spectral neutrino transport. The explodability of strongly and weakly magnetized models of $20$ and $27$ $M_{\odot}$ pre-supernova progenitors are compared. We find that although the efficiency for the conversion of the neutrino heating into turbulent energy including magnetic fields in the gain region is not significantly different between the strong and weak field models, the amplified magnetic field due to the neutrino-driven convection on large hot bubbles just behind stalled shock results in a faster and more energetic explosion in the strongly magnetized models. In addition, by comparing the difference between the 2nd- and 5th-order spatial accuracy of the simulation in the strong field model for $27$ $M_{\odot}$ progenitor, we also find that the higher order accuracy in space is beneficial to the explosion because it enhances the growth of neutrino-driven convection in the gain region. Based on our results of core-collapse supernova simulations for the non-rotating model, a new possibility for the origin of the magnetic field of the protoneutron star (PNS) is proposed. The magnetic field is accumulated and amplified to magnetar level, that is, $\mathcal{O}(10^{14})$ G, in the convectively stable shell near the PNS surface.

preprint2022arXiv

Shock Revival in Core-collapse Supernovae Assisted by Heavy Axion-like Particles

Axion-like particles (ALPs) are a class of hypothetical pseudoscalar particles which feebly interact with ordinary matter. The hot plasma of core-collapse supernovae is a possible laboratory to explore physics beyond the standard model including ALPs. Once produced, some of the ALPs can be absorbed by the supernova matter and affect energy transfer. In this study, we calculate the ALP emission in core-collapse supernovae and the backreaction on supernova dynamics consistently. It is found that the stalled bounce shock can be revived if the coupling between ALPs and photons is as high as g_{ag}~10^{-9} GeV^{-1} and the ALP mass is 40-400 MeV. Most of the models result in more energetic explosions than the average observed supernova. While this can be used to place constraints on those ALPs, long-term simulations across multiple progenitors need to be further investigated to place robust limits.

preprint2022arXiv

Three-dimensional simulation of a core-collapse supernova for a binary star progenitor of SN 1987A

We present results from a self-consistent, non-rotating core-collapse supernova simulation in three spatial dimensions using a binary evolution progenitor model of SN 1987A by Urushibata et al. (2018). This 18.3 solar-mass progenitor model is evolved from a slow-merger of 14 and 8 solar-mass stars, and it satisfies most of the observational constraints such as red-to-blue evolution, lifetime, total mass and position in the Hertzsprung-Russell diagram at collapse, and chemical anomalies. Our simulation is initiated from a spherically symmetric collapse and mapped to the three-dimensional coordinates at 10 ms after bounce to follow the non-spherical hydrodynamics evolution. We obtain the neutrino-driven shock revival for this progenitor at 350 ms after bounce, leading to the formation of a newly-born neutron star with average gravitational mass of 1.35 solar mass and spin period of 0.1 s. We also discuss the detectability of gravitational wave and neutrino signals for a Galactic event with the same characteristics as SN 1987A. At our final simulation time (660 ms postbounce), the diagnostic explosion energy, though still growing, is smaller (0.15 foe) compared to the observed value (1.5 foe). The 56Ni mass obtained from the simulation (0.01 solar mass) is also smaller than the reported mass from SN 1987A (0.07 solar mass). Long-term simulation including several missing physical ingredients in our 3D models such as rotation, magnetic fields, or more elaborate neutrino opacity should be done to bridge the gap between the theoretical predictions and the observed values.

preprint2021arXiv

Characteristic Time Variability of Gravitational-Wave and Neutrino Signals from Three-dimensional Simulations of Non-Rotating and Rapidly Rotating Stellar Core-Collapse

We present results from full general relativistic three-dimensional hydrodynamics simulations of stellar core collapse of a 70 M$_\odot$ star with spectral neutrino transport. To investigate the impact of rotation on non-axisymmetric instabilities, we compute three models by parametrically changing the initial strength of rotation. The most rapidly rotating model exhibits a transient development of the low-$T/|W|$ instability with one-armed spiral flow at the early postbounce phase. Subsequently, the two-armed spiral flow appears, which persists during the simulation time. The moderately rotating model also shows the growth of the low-$T/|W|$ instability, but only with the two-armed spiral flow. In the nonrotating model, a vigorous activity of the standing accretion-shock instability (SASI) is only observed. The SASI is first dominated by the sloshing mode, which is followed by the spiral SASI until the black hole formation. We present a spectrogram analysis of the gravitational waves (GWs) and neutrinos, focusing on the time correlation. Our results show that characteristic time modulations in the GW and neutrino signals can be linked to the growth of the non-axisymmetric instabilities. We find that the degree of the protoneutron star (PNS) deformation, depending upon which modes of the non-axisymmetric instabilities develop, predominantly affects the characteristic frequencies of the correlated GW and neutrino signals. We point out that these signals would be simultaneously detectable by the current-generation detectors up to $\sim10$ kpc. Our findings suggest that the joint observation of GWs and neutrinos is indispensable for extracting information on the PNS evolution preceding the black hole formation.

preprint2021arXiv

Convection and Dynamo in Newly-born Neutron Stars

To study properties of magneto-hydrodynamic (MHD) convection and resultant dynamo activities in proto-neutron stars (PNSs), we construct a "PNS in a box" simulation model with solving compressible MHD equation coupled with a nuclear equation of state (EOS) and a simplified leptonic transport. As a demonstration, we apply it to two types of PNS models with different internal structures: fully-convective model and spherical-shell convection model. By varying the spin rate of models, the rotational dependence of convection and dynamo that operate inside the PNS is investigated. We find that, as a consequence of turbulent transport by rotating stratified convection, large-scale structures of flow and thermodynamic fields are developed in all models. Depending on the spin rate and the convection zone depth, various profiles of the large-scale structures are obtained, which can be physically understood as steady-state solutions to the "mean-field" equation of motion. Additionally to those hydrodynamic structures, the large-scale magnetic component with $\mathcal{O}(10^{15})$ G is also spontaneously organized in disordered tangled magnetic fields in all models. The higher the spin rate, the stronger the large-scale magnetic component is built up. Intriguingly, as an overall trend, the fully-convective models have a stronger large-scale magnetic component than that in the spherical-shell convection models. The deeper the convection zone extends, the larger the size of the convection eddies becomes. As a result, the rotationally-constrained convection seems to be more easily achieved in the fully-convective model, resulting in the higher efficiency of the large-scale dynamo there. To gain a better understanding of the origin of the diversity of NS's magnetic field, we need to study the PNS dynamo in a wider parameter range.

preprint2021arXiv

Core-collapse supernova simulations and the formation of neutron stars, hybrid stars, and black holes

We investigate observable signatures of a first-order quantum chromodynamics (QCD) phase transition in the context of core collapse supernovae. To this end, we conduct axially symmetric numerical relativity simulations with multi-energy neutrino transport, using a hadron-quark hybrid equation of state (EOS). We consider four non-rotating progenitor models, whose masses range from $9.6$ to $70$ M$_\odot$. We find that the two less massive progenitor stars (9.6 and 11.2 M$_\odot$) show a successful explosion, which is driven by the neutrino heating. They do not undergo the QCD phase transition and leave behind a neutron star (NS). As for the more massive progenitor stars (50 and 70 M$_\odot$), the proto-neutron star (PNS) core enters the phase transition region and experiences the second collapse. Because of a sudden stiffening of the EOS entering to the pure quark matter regime, a strong shock wave is formed and blows off the PNS envelope in the 50 M$_\odot$ model. Consequently, the remnant becomes a quark core surrounded by hadronic matters, leading to the formation of the hybrid star. However for the 70 M$_\odot$ model, the shock wave cannot overcome the continuous mass accretion and it readily becomes a black hole. We find that the neutrino and gravitational wave (GW) signals from supernova explosions driven by the hadron-quark phase transition are detectable for the present generation of neutrino and GW detectors. Furthermore, the analysis of the GW detector response reveals unique kHz signatures, which will allow us to distinguish this class of supernova explosions from failed and neutrino-driven explosions.

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

Pre-supernova Ultra-light Axion-like Particles

We calculate the production of ultra-light axion-like particles (ALPs) in a nearby supernova progenitor. Once produced, ALPs escape from the star and a part of them is converted into photons during propagation in the Galactic magnetic field. It is found that the MeV photon flux that reaches Earth may be detectable by gamma ray telescopes for ALPs lighter than ~1 neV when Betelgeuse undergoes oxygen and silicon burning. The dependence of the gamma ray flux on the stellar mass is much smaller than the uncertainty that originates from the Galactic magnetic field. If ALPs are lighter than ~0.1 neV and the supernova progenitor is close enough to the Solar System, the gamma ray flux is insensitive to the distance d because the ALP-photon conversion probability is proportional to d^2. (Non-)detection of gamma rays from a supernova progenitor with next-generation gamma ray telescopes just after pre-supernova neutrino alerts would lead to an independent constraint on ALP parameters as stringent as a SN 1987A limit.

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 New Gravitational Wave Signature of Low-$T/|W|$ Instability in Rapidly Rotating Stellar Core Collapse

We present results from a full general relativistic three-dimensional hydrodynamics simulation of rapidly rotating core-collapse of a 70 M$_{\odot}$ star with three-flavor spectral neutrino transport. We find a strong gravitational wave (GW) emission that originates from the growth of the one- and two-armed spiral waves extending from the nascent proto-neutron star (PNS). The GW spectrogram shows several unique features that are produced by the non-axisymmetric instabilities. After bounce, the spectrogram first shows a transient quasi-periodic time modulation at $\sim$ 450 Hz. In the second active phase, it again shows the quasi-periodic modulation but with the peak frequency increasing with time, which continues until the final simulation time. From our detailed analysis, such features can be well explained by a combination of the so-called low-$T/|W|$ instability and the PNS core contraction.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

Detectability of SASI activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of Standing Accretion Shock Instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda $\it {et\, al.}$ (2017) as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to $\sim 6$ kpc for IceCube and and up to $\sim 3$ kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

preprint2020arXiv

Magnetorotational Explosion of A Massive Star Supported by Neutrino Heating in General Relativistic Three Dimensional Simulations

We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of progenitor's rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M$_{\odot}$ star. While we find no shock revival in our two non-magnetized models during our simulation times ($\sim500$ ms after bounce), the magnetorotationally (MR) driven shock expansion immediately initiates after bounce in our rapidly rotating and strongly magnetized model. We show that the expansion of the MR-driven flows toward the polar directions is predominantly driven by the magnetic pressure, whereas the shock expansion toward the equatorial direction is supported by neutrino heating. Our detailed analysis indicates that the growth of the so-called kink instability may hinder the collimation of jets, resulting in the formation of broader outflows. Furthermore we find a dipole emission of lepton number, only in the MR explosion model, whose asymmetry is consistent with the explosion morphology. Although it is similar to the lepton-number emission self-sustained asymmetry (LESA), our analysis shows that the dipole emission occurs not from the protoneutron star convection zone but from above the neutrino sphere indicating that it is not associated with the LESA. We also report several unique neutrino signatures, which are significantly dependent on both the time and the viewing angle, if observed, possibly providing a rich information regarding the onset of the MR-driven explosion.

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

Time of Flight and Supernova Progenitor Effects on the Neutrino Halo

We argue that the neutrino halo, a population of neutrinos that have undergone direction-changing scattering in the stellar envelope of a core-collapse supernova (CCSNe), is sensitive to neutrino emission history through time of flight. We show that the constant time approximation, commonly used in calculating the neutrino halo, does not capture the spatiotemporal evolution of the halo neutrino population and that correcting for time of flight can produce conditions which may trigger fast neutrino flavor conversion. We also find that there exists a window of time early in all CCSNe where the neutrino halo population is sufficiently small that it may be negligible. This suggests that collective neutrino oscillation calculations which neglect the Halo may be well founded at sufficiently early times.