Source author record

Takami Kuroda

Takami Kuroda 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

17works
5topics
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

17 published item(s)

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.

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

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.

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

Impact of magnetic field on neutrino-matter interactions in core-collapse supernova

We explore the impact of magnetic field on neutrino-matter interactions in core-collapse supernova. We first derive the modified source terms for neutrino-nucleon scattering and neutrino absorption and emission processes in the moment formalism. Then we perform full relativistic three-dimensional, magnetorotational core-collapse supernova simulations of a 20 $M_\odot$ star with spectral neutrino transport. Our simulations treat self-consistently the parity violation effects of weak interaction in the presence of external magnetic field. The result shows a significant global asymmetry, mostly confined in the proto-neutron star, with clearly reflecting the magnetic field structure. The asymmetric property arises from two factors: the angle between the neutrino flux and magnetic field, and the term, which is parallel to the magnetic field and is also proportional to the deviation of distribution function of neutrinos from thermal equilibrium. The typical correction value amounts to $\sim1$ % relative to the total neutrino-matter interaction rate for the magnetic field strength of $\sim 10^{15-16}$~G. Although these asymmetric properties do not immediately affect the explosion dynamics, our results imply that they would be significant once the neutrinos diffuse out the proto-neutron star core carrying those asymmetries away. We also show that, during our simulation time of $\sim370$ ms after bounce, our results indicate that the correction value due to the modified inelastic scattering process dominates over that of the modified neutrino absorption and emission process.

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.

preprint2016arXiv

A New Gravitational-Wave Signature from Standing Accretion Shock Instabilities in Supernovae

We present results from fully relativistic three-dimensional core-collapse supernova (CCSN) simulations of a non-rotating 15 Msun star using three different nuclear equations of state (EoSs). From our simulations covering up to ~350 ms after bounce, we show that the development of the standing accretion shock instability (SASI) differs significantly depending on the stiffness of nuclear EoS. Generally, the SASI activity occurs more vigorously in models with softer EoS. By evaluating the gravitational-wave (GW) emission, we find a new GW signature on top of the previously identified one, in which the typical GW frequency increases with time due to an accumulating accretion to the proto-neutron star (PNS). The newly observed quasi-periodic signal appears in the frequency range from ~100 to 200 Hz and persists for ~150 ms before neutrino-driven convection dominates over the SASI. By analyzing the cycle frequency of the SASI sloshing and spiral modes as well as the mass accretion rate to the emission region, we show that the SASI frequency is correlated with the GW frequency. This is because the SASI-induced temporary perturbed mass accretion strike the PNS surface, leading to the quasi-periodic GW emission. Our results show that the GW signal, which could be a smoking-gun signature of the SASI, is within the detection limits of LIGO, advanced Virgo, and KAGRA for Galactic events.

preprint2016arXiv

A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars

We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrinoelectron scattering, thermal neutrino production via pair annihilation and nucleonnucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino-matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock-stall of a 15Msun star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full-Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops-class supercomputers are at least needed.

preprint2016arXiv

Circular polarizations of gravitational waves from core-collapse supernovae: a clear indication of rapid rotation

We propose to employ the circular polarization of gravitational waves emitted by core-collapse supernovae as an unequivocal indication of rapid rotation deep in their cores. It has been demonstrated by three dimensional simulations that non-axisymmetric accretion flows may develop spontaneously via hydrodynamical instabilities in the post-bounce cores. It is not surprising then that the gravitational waves emitted by such fluid motions are circularly polarized. We show in this letter that a network of the second generation detectors of gravitational waves worldwide may be able to detect such polarizations up to the opposite side of Galaxy as long as the rotation period is shorter than a few seconds prior to collapse.

preprint2015arXiv

Coherent Network Analysis of Gravitational Waves from Three-Dimensional Core-Collapse Supernova Models

Using predictions from three-dimensional (3D) hydrodynamics simulations of core-collapse supernovae (CCSNe), we present a coherent network analysis to detection, reconstruction, and the source localization of the gravitational-wave (GW) signals. We use the {\tt RIDGE} pipeline for the analysis, in which the network of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA is considered. By combining with a GW spectrogram analysis, we show that several important hydrodynamics features in the original waveforms persist in the waveforms of the reconstructed signals. The characteristic excess in the spectrograms originates not only from rotating core-collapse, bounce and the subsequent ring down of the proto-neutron star (PNS) as previously identified, but also from the formation of magnetohydrodynamics jets and non-axisymmetric instabilities in the vicinity of the PNS. Regarding the GW signals emitted near at the rotating core bounce, the horizon distance extends up to $\sim$ 18 kpc for the most rapidly rotating 3D model in this work. Following the rotating core bounce, the dominant source of the GW emission shifts to the non-axisymmetric instabilities. The horizon distances extend maximally up to $\sim$ 40 kpc seen from the spin axis. With an increasing number of 3D models trending towards explosion recently, our results suggest that in addition to the best studied GW signals due to rotating core-collapse and bounce, the time is ripe to consider how we can do science from GWs of CCSNe much more seriously than before. Particularly the quasi-periodic signals due to the non-axisymmetric instabilities and the detectability should deserve further investigation to elucidate the inner-working of the rapidly rotating CCSNe.

preprint2015arXiv

Systematic Features of Axisymmetric Neutrino-Driven Core-Collapse Supernova Models in Multiple Progenitors

We present an overview of two-dimensional (2D) core-collapse supernova simulations employing neutrino transport scheme by the isotropic diffusion source approximation. We study 101 solar-metallicity, 247 ultra metal-poor, and 30 zero-metal progenitors covering zero-age main sequence mass from $10.8 M_{\odot} $ to $75.0 M_{\odot} $. Using the 378 progenitors in total, we systematically investigate how the differences in the structures of these multiple progenitors impact the hydrodynamics evolution. By following a long-term evolution over 1.0 s after bounce, most of the computed models exhibit neutrino-driven revival of the stalled bounce shock at $\sim$ 200 - 800 ms postbounce, leading to the possibility of explosion. Pushing the boundaries of expectations in previous one-dimensional (1D) studies, our results confirm that the compactness parameter $ξ$ that characterizes the structure of the progenitors is also a key in 2D to diagnose the properties of neutrino-driven explosions. Models with high $ξ$ undergo high ram pressure from the accreting matter onto the stalled shock, which affects the subsequent evolution of the shock expansion and the mass of the protoneutron star under the influence of neutrino-driven convection and the standing accretion-shock instability. We show that the accretion luminosity becomes higher for models with high $ξ$, which makes the growth rate of the diagnostic explosion energy higher and the synthesized nickel mass bigger. We find that these explosion characteristics tend to show a monotonic increase as a function of the compactness parameter $ξ$.

preprint2014arXiv

Asymmetric Neutrino Production in Strongly Magnetized Proto-Neutron Stars

We calculate the neutrino production cross-section through the direct URCA process in proto-neutron star matter in the presence of a strong magnetic field. We assume isoentropic conditions and introduce a new equation of state parameter-set in the relativistic mean-field approach that can reproduce neutron stars with $M > 1.96$ M$_\odot$ as required by observations. We find that the production process increases the flux of emitted neutrinos along the direction parallel to the magnetic field and decreases the flux in the opposite direction. This means that the neutrino flux asymmetry due to the neutrino absorption and scattering processes in a magnetic field becomes larger by the inclusion of the neutrino production process.

preprint2014arXiv

Gravitational Wave Signatures from Low-mode Spiral Instabilities in Rapidly Rotating Supernova Cores

We study properties of gravitational waves (GWs) from rotating core-collapse of a 15M_odot star by performing three-dimensional general-relativistic hydrodynamic simulations with an approximate neutrino transport. By parametrically changing the precollapse angular momentum, we focus on the effects of rotation on the GW signatures in the early postbounce evolution. Regarding three-flavor neutrino transport, we solve the energy-averaged set of radiation energy and momentum based on the Thorne's momentum formalism. In addition to the gravitational quadrupole radiation from matter motions, we take into account GWs from anisotropic neutrino emission. With these computations, our results present several supporting evidences for the previous anticipation that non-axisymmetric instabilities play an essential role in determining the postbounce GW signatures. During prompt convection, we find that the waveforms show narrow-band and highly quasi-periodic signals which persist until the end of simulations. We point out that such feature reflects the growth of the one-armed spiral modes. The typical frequency of the quasi-periodic waveforms can be well explained by the propagating acoustic waves between the stalled shock and the rotating proto-neutron star surface, which suggests the appearance of the standing-accretion-shock instability. Although the GW signals exhibit strong variability between the two polarizations and different observer directions, they are within the detection limits of next generation detectors such as by KAGRA and Advanced LIGO, if the source with sufficient angular momentum is located in our Galaxy.

preprint2014arXiv

Impacts of Rotation on Three-dimensional Hydrodynamics of Core-collapse Supernovae

We perform a series of simplified numerical experiments to explore how rotation impacts on the three-dimensional (3D) hydrodynamics of core-collapse supernovae. For the sake of our systematic study, we employ a light-bulb scheme to trigger explosions and a three-flavor neutrino leakage scheme to treat deleptonization effects and neutrino losses from proto-neutron star interior. Using a 15 solar mass progenitor, we compute thirty models in 3D with a wide variety of initial angular momentum and light-bulb neutrino luminosity. We find that the rotation can help onset of neutrino-driven explosions for the models in which the initial angular momentum is matched to that obtained in recent stellar evolutionary calculations (0.3-3 rad/s at the center). For the models with larger initial angular momentum, the shock surface deforms to be more oblate due to larger centrifugal force. This makes not only a gain region more concentrated around the equatorial plane, but also the mass in the gain region bigger. As a result, buoyant bubbles tend to be coherently formed and rise in the equatorial region, which pushes the revived shock ever larger radii until a global explosion is triggered. We find that these are the main reasons that the preferred direction of explosion in 3D rotating models is often perpendicular to the spin axis, which is in sharp contrast to the polar explosions around the axis that was obtained in previous 2D simulations.

preprint2012arXiv

Core-Collapse Supernovae as Supercomputing Science: a status report toward 6D simulations with exact Boltzmann neutrino transport in full general relativity

This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers.

preprint2012arXiv

Fully General Relativistic Simulations of Core-Collapse Supernovae with An Approximate Neutrino Transport

We present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport. Using a M1 closure scheme with an analytic variable Eddington factor, we solve the energy-independent set of radiation energy and momentum based on the Thorne's momentum formalism. To simplify the source terms of the transport equations, a methodology of multiflavour neutrino leakage scheme is partly employed. Our newly developed code is designed to evolve the Einstein field equation together with the GR radiation hydrodynamic equations. We follow the dynamics starting from the onset of gravitational core-collapse of a 15 $M_{\odot}$ star, through bounce, up to about 100 ms postbounce in this study to study how the spacial multi-dimensionality and GR would affect the dynamics in the early postbounce phase. Our 3D results support the anticipation in previous 1D results that the neutrino luminosity and average neutrino energy of any neutrino flavor in the postbounce phase increase when switching from SR to GR hydrodynamics. This is because the deeper gravitational well of GR produces more compact core structures, and thus hotter neutrino spheres at smaller radii. By analyzing the residency timescale to the neutrino-heating timescale in the gain region, we show that the criterion to initiate neutrino-driven explosions can be most easily satisfied in 3D models, irrespective of SR or GR hydrodynamics. Our results suggest that the combination of GR and 3D hydrodynamics provides the most favorable condition to drive a robust neutrino-driven explosion.

preprint2010arXiv

Three Dimensional Magneto Hydrodynamical Simulations of Gravitational Collapse of a 15Msun Star

We introduce our newly developed two different, three dimensional magneto hydrodynamical codes in detail. One of our codes is written in the Newtonian limit (NMHD) and the other is in the fully general relativistic code (GRMHD). Both codes employ adaptive mesh refinement and, in GRMHD, the metric is evolved with the "Baumgarte-Shapiro-Shibata-Nakamura" formalism known as the most stable method at present. We did several test problems and as for the first practical test, we calculated gravitational collapse of a $15M_\odot$ star. Main features found from our calculations are; (1) High velocity bipolar outflow is driven from the proto-neutronstar and moves through along the rotational axis in strongly magnetized models; (2) A one-armed spiral structure appears which is originated from the low-$|T/W|$ instability; (3) By comparing GRMHD and NMHD models, the maximum density increases about $\sim30%$ in GRMHD models due to the stronger gravitational effect. These features agree very well with previous studies and our codes are thus reliable to numerical simulation of gravitational collapse of massive stars.