Researcher profile

Tomoya Takiwaki

Tomoya Takiwaki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
51works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

51 published item(s)

preprint2026arXiv

Verifying the failing supernova constraint on dark photons with two-dimensional hydrodynamic simulations

Recent studies on the dark photon (DP) production in collapsing stars argue that the cooling effect induced by DPs can hinder supernova explosions and lead to a ``failing supernova" constraint on the photon-DP mixing parameter $ε$. In order to verify the idea, we perform two-dimensional neutrino-radiation hydrodynamic simulations coupled with the DP production with the masses of 0.3 and 0.45\,MeV. We find that the shock revival does not happen until the end of the simulations when $ε\gtrsim3\times10^{-9}$. The photon-DP mixing parameter above this value can be excluded by the failing supernova argument. Interestingly, our constraint roughly coincides with the one reported by the previous studies which adopted the post-processing framework. This result motivates one to investigate a wider parameter range of DPs with self-consistent simulations and evaluate uncertainties in the constraint.

preprint2022arXiv

Dynamics of fast neutrino flavor conversions with scattering effects: a detailed analysis

We calculate fast conversions of two flavor neutrinos by considering Boltzmann collisions of neutrino scatterings. In an idealized angular distribution of neutrinos with electron-lepton number crossing, we find that the collision terms of the neutrino scattering enhance the transition probability of fast flavor conversions as in the previous study. We analyze the dynamics of fast flavor conversions with collisions in detail based on the motion of polarization vectors in cylindrical coordinate analogous to a pendulum motion. The phase of the all polarization vector synchronizes in the linear evolution, and the phase deviation from the Hamiltonian governs the conversion of neutrino flavor. In the non-linear evolution, a closed orbit in the phase space is observed. The collision terms break the closed orbit and gradually make the phase space smaller. The flavor conversions are enhanced during this limit cycle. After the significant flavor conversion, all of the neutrino polarization vectors start to align with the z-axis owing to the collision effect within the time scale of the collision term irrespective of neutrino scattering angles. We also show the enhancement or suppression of the flavor conversions in various setups of the collision terms and verify consistency with previous studies. Though our analysis does not fully understand the dynamics of fast flavor conversion, the framework gives a new insight into this complicated phenomenon in further study.

preprint2022arXiv

Long-term evolution of a supernova remnant hosting a double neutron star binary

An ultra-stripped supernova (USSN) is a type of core-collapse SN explosion proposed to be a candidate formation site of a double neutron star (DNS) binary. We investigate the dynamical evolution of an ultra-stripped supernova remnant (USSNR), which should host a DNS at its center. By accounting for the mass-loss history of the progenitor binary using a model developed by a previous study, we construct the large-scale structure of the {circumstellar medium (CSM)} up to a radius $\sim 100\,{\rm pc}$, and simulate the explosion and subsequent evolution of a USSN surrounded by such a CSM environment. We find that the CSM encompasses an extended region characterized by a hot plasma with a temperature $\sim 10^8\,$K located around the termination shock of the wind from the progenitor binary ($\sim 10\,$pc), and the USSNR blastwave is drastically weakened while penetrating through this hot plasma. Radio continuum emission from a young USSNR is sufficiently bright to be detectable if it inhabits our Galaxy but faint compared to the observed Galactic SNRs, and thereafter declines in luminosity through adiabatic cooling. Within our parameter space, USSNRs typically exhibit a low radio luminosity and surface brightness compared to the known Galactic SNRs. Due to the small event rate of USSNe and their relatively short observable lifespan, we calculate that USSNRs account for only $\sim 0.1$-$1$ % of the total SNR population. This is consistent with the fact that no SNR hosting a DNS binary has been discovered in the Milky Way so far.

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

A systematic study on the rise time-peak luminosity relation for bright optical transients powered by wind shock breakout

A number of astrophyical transients originating from stellar explosions are powered by the collision of the ejected material with the circumstellar medium, which efficiently produces thermal radiation via shock dissipation. We investigate how such interaction-powered transients are distributed in the peak bolometric luminosity vs the rise time phase space. Taking the advantage of less time-consuming one-dimensional simulations with spherical symmetry, we calculated more than 500 models with different circumstellar mass and radius, ejecta mass and energy, and chemical compositions. The peak bolometric luminosity, the total radiated energy, and the rise time of the interaction-powered emission are measured for each simulated light curve. We consider how these characteristic quantities are determined as a function of the model parameters and discuss possible implications for the observed populations of (potential) interaction-powered transients, such as type IIn supernovae and fast blue optical transients.

preprint2020arXiv

Detectability of Collective Neutrino Oscillation Signatures in the Supernova Explosion of a 8.8 $M_\odot$ star

In order to investigate the impact of collective neutrino oscillations (CNO) on the neutrino signal from a nearby supernova, we perform 3-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 \Msun progenitor star. We find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\barν_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to $\sim 10$ kpc. On the other hand, for the normal mass hierarchy, the spectrum of $ν_e$ becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to $\sim 10$ kpc. More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\barν_e$ in HK becomes softer due to CNO, the spectrum of $ν_e$ in DUNE becomes harder, and vice versa. This synergistic observations in $\barν_e$ and $ν_e$, by HK and DUNE respectively, will be an intriguing opportunity to test the occurrence of CNO.

preprint2020arXiv

Dimension dependence of numerical simulations on gravitational waves from protoneutron stars

We examine the eigenfrequencies of gravitational waves from the protoneutron stars (PNSs) provided via the core-collapse supernovae, focusing on how the frequencies depend on the dimension of numerical simulations especially for the early phase after core bounce. As expected, we find that the time evolution of gravitational wave frequencies depends strongly on the dimension of numerical simulation as well as the equation of state for high density matter. Even so, we find that the fundamental frequency as a function of PNS average density and the ratio of the specific eigenfrequencies to the fundamental frequency as a function of PNS properties are independent of the dimension of the numerical simulations, where the dependence of the equation of state is also very weak in this early postbounce phase. Thus, one can safely discuss the gravitational wave frequencies of the PNSs as a function of the PNS average density or compactness even with the frequencies obtained from the one dimensional simulations. We also provide phenomenological relation between the compactness and average density as well as the relation among the $f$-, $p$-, and $g$- mode frequencies.

preprint2020arXiv

Extremely energetic supernova explosions embedded in a massive circumstellar medium: the case of SN 2016aps

We perform one-dimensional radiation-hydrodynamic simulations of energetic supernova ejecta colliding with a massive circumstellar medium (CSM) aiming at explaining SN 2016aps, likely the brightest supernova observed to date. SN 2016aps was a superluminous Type-IIn SN, which released as much as $\gtrsim 5\times 10^{51}$ erg of thermal radiation. Our results suggest that the multi-band light curve of SN 2016aps is well explained by the collision of a $30\ M_\odot$ SN ejecta with the explosion energy of $10^{52}$ erg and a $\simeq 8\ M_\odot$ wind-like CSM with the outer radius of $10^{16}$ cm, i.e., a hypernova explosion embedded in a massive CSM. This finding indicates that very massive stars with initial masses larger than $40\ M_\odot$, which supposedly produce highly energetic SNe, occasionally eject their hydrogen-rich envelopes shortly before the core-collapse. We suggest that the pulsational pair-instability SNe may provide a natural explanation for the massive CSM and the energetic explosion. We also provide the relations between the peak luminosity, the radiated energy, and the rise time for interacting SNe with the kinetic energy of $10^{52}$ erg, which can be used for interpreting SN 2016aps-like objects in future surveys.

preprint2020arXiv

Fast Neutrino Flavor Conversion Modes in Multidimensional Core-collapse Supernova Models: the Role of the Asymmetric Neutrino Distributions

A dense neutrino gas, such as the one anticipated in the supernova environment, can experience fast neutrino flavor conversions on scales much shorter than those expected in vacuum probably provided that the angular distributions of $ν_e$ and $\barν_e$ cross each other. We perform a detailed investigation of the neutrino angular distributions obtained by solving the Boltzmann equations for fixed matter profiles of some representative snapshots during the post-bounce phase of core-collapse supernovae in multidimensional calculations of an $11.2\mathrm{M}_{\odot}$ and a $27\mathrm{M}_{\odot}$ progenitor models. Although the $11.2\mathrm{M}_{\odot}$ model features $ν_e - \barν_e$ angular crossings and the associated fast modes at different time snapshots, the $27\mathrm{M}_{\odot}$ model does not show any crossings within the decoupling region. We show that this can be understood by studying the multipole components of the neutrino distributions. In fact, there is a higher chance for the occurrence of $ν_e - \barν_e$ angular crossings for the zones where the multipole components of the neutrino distributions are strong enough. We also show that there can exist more than one crossings between the angular distributions of $ν_e$ and $\barν_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $ν_e-\barν_e$ angular crossings which appears very deep inside the proto-neutron star.

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

On the occurrence of fast neutrino flavor conversions in multidimensional supernova models

The dense neutrino medium in a core-collapse supernova or a neutron-star merger event can experience fast flavor conversions on time/distance scales that are much smaller than those of vacuum oscillations. It is believed that fast neutrino flavor transformation occurs in the region where the angular distributions of $ν_e$ and $\barν_e$ cross each other. We present the first study of this crossing phenomenon and the fast neutrino flavor conversions in multidimensional (multi-D) supernova models. We examine the neutrino distributions obtained by solving the Boltzmann transport equation for several fixed profiles which are representative snapshots taken from separate 2D and 3D supernova simulations with an $11.2 M_\odot$ progenitor model. Our research shows that the spherically asymmetric patterns of the $ν_e$ and $\barν_e$ fluxes in multi-D models can assist the appearance of the crossing between the $ν_e$ and $\barν_e$ angular distributions. In the models that we have studied, there exist unstable neutrino oscillation modes in and beyond the neutrino decoupling region which have amplitude growth rates as large as an $e$-fold per nanosecond. This finding can have important consequences for the explosion mechanism, nucleosynthesis, and neutrino signals of core-collapse supernovae.

preprint2020arXiv

Systematic investigation of the effect of 56Ni mixing in the early photospheric velocity evolution of stripped-envelope supernovae

Mixing of 56Ni, whose nuclear decay energy is a major luminosity source in stripped-envelope supernovae, is known to affect the observational properties of stripped-envelope supernovae such as light-curve and color evolution. Here we systematically investigate the effect of 56Ni mixing on the photospheric velocity evolution in stripped-envelope supernovae. We show that 56Ni mixing significantly affects the early photospheric velocity evolution. The photospheric velocity, which is often used to constrain the ejecta mass and explosion energy, significantly varies by just changing the degree of 56Ni mixing. In addition, the models with a small degree of 56Ni mixing show a flattening in the early photospheric velocity evolution, while the fully mixed models show a monotonic decrease. The velocity flattening appears in both helium and carbon+oxygen progenitor explosions with a variety of ejecta mass, explosion energy, and 56Ni mass. Some stripped-envelope supernovae with early photospheric velocity information do show such a flattening. We find that Type Ib SN 2007Y, which has early photospheric velocity information, has a signature of a moderate degree of 56Ni mixing in the photospheric velocity evolution and about a half of the ejecta is mixed in it. The immediate spectroscopic follow-up observations of stripped-envelope supernovae shortly after the explosion providing the early photospheric evolution give an important clue to constrain 56Ni mixing in the ejecta.

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.

preprint2019arXiv

Supernova ejecta interacting with a circumstellar disk. I. two-dimensional radiation-hydrodynamic simulations

We perform a series of two-dimensional radiation-hydrodynamic simulations of the collision between supernova ejecta and circumstellar media (CSM). The hydrodynamic interaction of a fast flow and the surrounding media efficiently dissipates the kinetic energy of the fast flow and considered as a dominant energy source for a specific class of core-collapse supernovae. Despite some observational evidence for aspherical ejecta and/or CSM structure, multi-dimensional effects in the ejecta-CSM interaction are relatively unexplored. Our numerical simulations equipped with an adaptive mesh refinement technique successfully reproduce hydrodynamic instabilities developing around the ejecta-CSM interface. We also investigate effects of disk-like CSM on the dynamical evolution of supernova ejecta and bolometric light curves. We find that emission powered by ejecta-disk interaction exhibits significant viewing angle dependence. For a line of sight close to the symmetry axis, the observer directly sees the supernova ejecta, leading to a short brightening timescale. For an observer seeing the emission through the CSM disk, thermal photons diffuse throughout the CSM and thus the light curve is severely smeared out.

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

Gravitational wave asteroseismology with protoneutron stars

We examine the time evolution of the frequencies of the gravitational wave after the bounce within the framework of relativistic linear perturbation theory using the results of one dimensional numerical simulations of core-collapse supernovae. Protoneutron star models are constructed in such a way that the mass and radius of protoneutron star become equivalent to the results obtained from the numerical simulations. Then, we find that the frequencies of gravitational waves radiating from protoneutron stars strongly depend on the mass and radius of protoneutron stars, but almost independently of the profiles of electron fraction and entropy per baryon inside the star. Additionally, we find that the frequencies of gravitational waves can be characterized by the square root of the average density of protoneutron star irrespectively the progenitor models, which are completely different from the empirical formula for cold neutron stars. The dependence of the spectra on the mass and radius is different from that of the $g$-mode: the oscillations around the surface of protoneutron stars due to the convection and the standing accretion-shock instability. Careful observations of the these modes of gravitational waves can determine the evolution of the mass and radius of protoneutron stars after core-bounce. Furthermore, the expected frequencies of gravitational waves are around a few hundred hertz in the early stage after bounce, which must be a good candidate for the ground-based gravitational wave detectors.

preprint2016arXiv

Three-dimensional simulations of rapidly rotating core-collapse supernovae: finding a neutrino-powered explosion aided by non-axisymmetric flows

We report results from a series of three-dimensional (3D) rotational core-collapse simulations for $11.2$ and 27 $M_{/odot}$ stars employing neutrino transport scheme by the isotropic diffusion source approximation. By changing the initial strength of rotation systematically, we find a rotation-assisted explosion for the 27$M_{/odot}$ progenitor, which fails in the absence of rotation. The unique feature was not captured in previous two-dimensional (2D) self-consistent rotating models because the growing non-axisymmetric instabilities play a key role. In the rapidly rotating case, strong spiral flows generated by the so-called low $T/|W|$ instability enhance the energy transport from the proto-neutron star (PNS) to the gain region, which makes the shock expansion more energetic. The explosion occurs more strongly in the direction perpendicular to the rotational axis, which is different from previous 2D predictions.

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 $ξ$.

preprint2015arXiv

The criterion of supernova explosion revisited: the mass accretion history

By performing neutrino-radiation hydrodynamic simulations in spherical symmetry (1D) and axial symmetry (2D) with different progenitor models by Woosley & Heger (2007) from 12 $M_{\odot}$ to 100 $M_{\odot}$, we find that all 1D runs fail to produce an explosion and several 2D runs succeed. The difference in the shock evolutions for different progenitors can be interpreted by the difference in their mass accretion histories, which are in turn determined by the density structures of progenitors. The mass accretion history has two phases in the majority of the models: the earlier phase in which the mass accretion rate is high and rapidly decreasing and the later phase with a low and almost constant accretion rate. They are separated by the so-called turning point, the origin of which is a change of the accreting layer. We argue that shock revival will most likely occur around the turning point and hence that its location in the $\dot M$-$L_ν$ plane will be a good measure for the possibility of shock revival: if the turning point lies above the critical curve and the system stays there for a long time, shock revival will obtain. In addition, we develop a phenomenological model to approximately evaluate the trajectories in the $\dot M$-$L_ν$ plane, which, after calibrating free parameters by a small number of 1D simulations, reproduces the location of the turning point reasonably well by using the initial density structure of progenitor alone. We suggest the application of the phenomenological model to a large collection of progenitors in order to infer without simulations which ones are more likely to explode.

preprint2015arXiv

The r-process nucleosynthesis in the various jet-like explosions of magnetorotational core-collapse supernovae

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.

preprint2014arXiv

A Comparison of Two- and Three-dimensional Neutrino-hydrodynamics simulations of Core-collapse Supernovae

We present numerical results on two- (2D) and three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2$M_\odot$ star. By changing numerical resolutions and seed perturbations systematically, we study how the postbounce dynamics is different in 2D and 3D. The calculations were performed with an energy-dependent treatment of the neutrino transport based on the isotropic diffusion source approximation scheme, which we have updated to achieve a very high computational efficiency. All the computed models in this work including nine 3D models and fifteen 2D models exhibit the revival of the stalled bounce shock, leading to the possibility of explosion. All of them are driven by the neutrino-heating mechanism, which is fostered by neutrino-driven convection and the standing-accretion-shock instability (SASI). Reflecting the stochastic nature of multi-dimensional (multi-D) neutrino-driven explosions, the blast morphology changes from models to models. However, we find that the final fate of the multi-D models whether an explosion is obtained or not, is little affected by the explosion stochasticity. In agreement with some previous studies, higher numerical resolutions lead to slower onset of the shock revival in both 3D and 2D. Based on the self-consistent supernova models leading to the possibility of explosions, our results systematically show that the revived shock expands more energetically in 2D than in 3D.

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.

preprint2014arXiv

Magnetohydrodynamic Turbulence Powered by Magnetorotational Instability in Nascent Proto-Neutron Stars

Magnetorotational instability (MRI) in a convectively-stable layer around the neutrinosphere is simulated by a three-dimensional model of supernova core. To resolve MRI-unstable modes, a thin layer approximation considering only the radial global stratification is adopted. Our intriguing finding is that the convectively-stable layer around the neutrinosphere becomes fully-turbulent due to the MRI and its nonlinear penetration into the strongly-stratified MRI-stable region. The intensity of the MRI-driven turbulence increases with magnetic flux threading the core, but is limited by a free energy stored in the differential rotation. The turbulent neutrinosphere is a natural consequence of rotating core-collapse and could exert a positive impact on the supernova mechanism.

preprint2014arXiv

Revisiting Impacts of Nuclear Burning for Reviving Weak Shocks in Neutrino-Driven Supernovae

We revisit potential impacts of nuclear burning on the onset of the neutrino-driven explosions of core-collapse supernovae. By changing the neutrino luminosity and its decay time to obtain parametric explosions in one-(1D) and two-dimensional (2D) models with or without a 13-isotope alpha network, we study how the inclusion of nuclear burning could affect the postbounce dynamics for four progenitor models; three for 15.0 Msun stars, one for an 11.2 Msun star. We find that the energy supply due to nuclear burning of infalling material behind the shock can energize the shock expansion especially for models that produce only marginal explosions in the absence of nuclear burning. These models are energized by nuclear energy deposition when the shock front passes through the silicon-rich layer and/or later it touches the oxygen-rich layer. Depending on the neutrino luminosity and its decay time, a diagnostic energy of explosion increases up to a few times 10^50 erg for models with nuclear burning compared to the corresponding models without. We point out that these features are most remarkable for the Limongi-Chieffi progenitor in both 1D and 2D, because the progenitor model possesses a massive oxygen layer with its inner-edge radius being smallest among the employed progenitors, so that the shock can touch the rich fuel on a shorter timescale after bounce. The energy difference is generally smaller (~0.1-0.2 times 10^51 erg) in 2D than in 1D (at most ~0.6 times 10^51 erg). This is because neutrino-driven convection and the shock instability in 2D models enhance the neutrino heating efficiency, which makes the contribution of nuclear burning relatively smaller compared to 1D models. Considering uncertainties in progenitor models, our results indicate that nuclear burning should remain as one of the important ingredients to foster the onset of neutrino-driven explosions.

preprint2014arXiv

The red supergiant and supernova rate problems: implications for core-collapse supernova physics

Mapping supernovae to their progenitors is fundamental to understanding the collapse of massive stars. We investigate the red supergiant problem, which concerns why red supergiants with masses $\sim16$-$30 M_\odot$ have not been identified as progenitors of Type IIP supernovae, and the supernova rate problem, which concerns why the observed cosmic supernova rate is smaller than the observed cosmic star formation rate. We find key physics to solving these in the compactness parameter, which characterizes the density structure of the progenitor. If massive stars with compactness above $ξ_{2.5} \sim 0.2$ fail to produce canonical supernovae, (i) stars in the mass range $16$-$30 M_\odot$ populate an island of stars that have high $ξ_{2.5}$ and do not produce canonical supernovae, and (ii) the fraction of such stars is consistent with the missing fraction of supernovae relative to star formation. We support this scenario with a series of two- and three-dimensional radiation hydrodynamics core-collapse simulations. Using more than 300 progenitors covering initial masses $10.8$-$75 M_\odot$ and three initial metallicities, we show that high compactness is conducive to failed explosions. We then argue that a critical compactness of $\sim 0.2$ as the divide between successful and failed explosions is consistent with state-of-the-art three-dimensional core-collapse simulations. Our study implies that numerical simulations of core collapse need not produce robust explosions in a significant fraction of compact massive star initial conditions.

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.

preprint2012arXiv

Gravitational Wave Signatures of Hyperaccreting Collapsar Disks

By performing two-dimensional special relativistic (SR) magnetohydrodynamic simulations, we study possible signatures of gravitational waves (GWs) in the context of the collapsar model for long-duration gamma-ray bursts. In our SR simulations, the central black hole is treated as an absorbing boundary. By doing so, we focus on the GWs generated by asphericities in neutrino emission and matter motions in the vicinity of the hyperaccreting disks. We compute nine models by adding initial angular momenta and magnetic fields parametrically to a precollapse core of a $35 M_{\odot}$ progenitor star. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multiflavor neutrino leakage scheme. To accurately estimate GWs produced by anisotropic neutrino emission, we perform a ray-tracing analysis in general relativity by a post-processing procedure. By employing a stress formula that includes contributions both from magnetic fields and special relativistic corrections, we study also the effects of magnetic fields on the gravitational waveforms. We find that the GW amplitudes from anisotropic neutrino emission show a monotonic increase with time, whose amplitudes are much larger than those from matter motions of the accreting material. We show that the increasing trend of the neutrino GWs stems from the excess of neutrino emission in the direction near parallel to the spin axis illuminated from the hyperaccreting disks. We point out that a recently proposed future space-based interferometer like Fabry-Perot type DECIGO would permit the detection of these GW signals within $\approx$ 100 Mpc.

preprint2012arXiv

Multimessengers from core-collapse supernovae: multidimensionality as a key to bridge theory and observation

Core-collapse supernovae are dramatic explosions marking the catastrophic end of massive stars. The only means to get direct information about the supernova engine is from observations of neutrinos emitted by the forming neutron star, and through gravitational waves which are produced when the hydrodynamic flow or the neutrino flux is not perfectly spherically symmetric. The multidimensionality of the supernova engine, which breaks the sphericity of the central core such as convection, rotation, magnetic fields, and hydrodynamic instabilities of the supernova shock, is attracting great attention as the most important ingredient to understand the long-veiled explosion mechanism. Based on our recent work, we summarize properties of gravitational waves, neutrinos, and explosive nucleosynthesis obtained in a series of our multidimensional hydrodynamic simulations and discuss how the mystery of the central engines can be unraveled by deciphering these multimessengers produced under the thick veils of massive stars.

preprint2012arXiv

On the Importance of the Equation of State for the Neutrino-Driven Supernova Explosion Mechanism

By implementing widely-used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 1D simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models of a 15 M_sun progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, which is not the case using the SHEN EOS. Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 11.2 M_sun progenitor produce neutrino-driven explosions for all the EOS under investigation. Models using the LS EOS are slightly more energetic compared to those with the SHEN EOS. The more efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity and a larger mass that is enclosed within the gain region. The models based on the LS EOS also show a more vigorous and aspherical downflow of accreting matter to the surface of the protoneutron star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure waves, that can push in turn the shock to larger radii and provide more favorable conditions for the explosion. [abbreviated]

preprint2012arXiv

Three-dimensional Hydrodynamic Core-Collapse Supernova Simulations for an $11.2 M_{\odot}$ Star with Spectral Neutrino Transport

We present numerical results on three-dimensional (3D) hydrodynamic core-collapse simulations of an $11.2 M_{\odot}$ star. By comparing one-(1D) and two-dimensional(2D) results with those of 3D, we study how the increasing spacial multi-dimensionality affects the postbounce supernova dynamics. The calculations were performed with an energy-dependent treatment of the neutrino transport that is solved by the isotropic diffusion source approximation scheme. By performing a tracer-particle analysis, we show that the maximum residency time of material in the gain region is shown to be longer for 3D due to non-axisymmetric flow motions than 2D, which is one of advantageous aspects of 3D models to obtain neutrino-driven explosions. Our results show that convective matter motions below the gain radius become much more violent in 3D than 2D, making the neutrino luminosity larger for 3D. Nevertheless the emitted neutrino energies are made smaller due to the enhanced cooling. Our results indicate whether these advantages for driving 3D explosions could or could not overwhelm the disadvantages is sensitive to the employed numerical resolutions. An encouraging finding is that the shock expansion tends to become more energetic for models with finer resolutions. To draw a robust conclusion, 3D simulations with much more higher numerical resolutions and also with more advanced treatment of neutrino transport as well as of gravity is needed, which could be hopefully practicable by utilizing forthcoming Petaflops-class supercomputers.

preprint2011arXiv

Gravitational Wave Signatures of Magnetohydrodynamically-Driven Core-Collapse Supernova Explosions

By performing a series of two-dimensional, special relativistic magnetohydrodynamic (MHD) simulations, we study signatures of gravitational waves (GWs) in the magnetohydrodynamically-driven core-collapse supernovae. In order to extract the gravitational waveforms, we present a stress formula including contributions both from magnetic fields and special relativistic corrections. By changing the precollapse magnetic fields and initial angular momentum distributions parametrically, we compute twelve models. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multiflavor neutrino leakage scheme. With these computations, we find that the total GW amplitudes show a monotonic increase after bounce for models with a strong precollapse magnetic field ($10^{12}$G) also with a rapid rotation imposed. We show that this trend stems both from the kinetic contribution of MHD outflows with large radial velocities and also from the magnetic contribution dominated by the toroidal magnetic fields that predominantly trigger MHD explosions. For models with weaker initial magnetic fields, the total GW amplitudes after bounce stay almost zero, because the contribution from the magnetic fields cancels with the one from the hydrodynamic counterpart. These features can be clearly understood with a careful analysis on the explosion dynamics. We point out that the GW signals with the increasing trend, possibly visible to the next-generation detectors for a Galactic supernova, would be associated with MHD explosions with the explosion energies exceeding $10^{51}$ erg.

preprint2011arXiv

Impacts of Collective Neutrino Oscillations on Supernova Explosions

By performing a series of one- and two-dimensional (1-, 2D) hydrodynamic simulations with spectral neutrino transport, we study possible impacts of collective neutrino oscillations on the dynamics of core-collapse supernovae. To model the spectral swapping which is one of the possible outcome of the collective neutrino oscillations, we parametrize the onset time when the spectral swap begins, the radius where the spectral swap occurs, and the threshold energy above which the spectral interchange between heavy-lepton neutrinos and electron/anti-electron neutrinos takes place, respectively. By doing so, we systematically study how the neutrino heating enhanced by the spectral swapping could affect the shock evolution as well as the matter ejection. We also investigate the progenitor dependence using a suite of progenitor models (13, 15, 20, and 25 $M_\odot$). We find that there is a critical heating rate induced by the spectral swapping to trigger explosions, which significantly differs between the progenitors. The critical heating rate is generally smaller for 2D than 1D due to the multidimensionality that enhances the neutrino heating efficiency. For the progenitors employed in this paper, the final remnant masses are estimated to range in 1.1-1.5$M_\odot$. For our 2D model of the $15M_\odot$ progenitor, we find a set of the oscillation parameters that could account for strong supernova explosions ($\sim 10^{51}$ erg), simultaneously leaving behind the remnant mass close to $\sim 1.4 M_\odot$.

preprint2010arXiv

Explosion geometry of a rotating 13 $M_{\odot}$ star driven by the SASI-aided neutrino-heating supernova mechanism

By performing axisymmetric hydrodynamic simulations of core-collapse supernovae with spectral neutrino transport based on the isotropic diffusion source approximation scheme, we support the assumption that the neutrino-heating mechanism aided by the standing accretion shock instability and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results show that bipolar explosions are more likely to be associated with models which include rotation. We point out that models, which form a north-south symmetric bipolar explosion, can lead to larger explosion energies than for the corresponding unipolar explosions.

preprint2010arXiv

Neutrino Pair Annihilation in Collapsars; Ray-Tracing Method in Special Relativity

We develop a numerical scheme and code for estimating the energy and momentum transfer via neutrino pair annihilation ($ν+ {\bar ν} \to e^{-}+ e^{+}$), bearing in mind the application to the collapsar models of gamma-ray bursts (GRBs). To calculate the neutrino flux illuminated from the accretion disk, we perform a ray-tracing calculation in the framework of special relativity. The numerical accuracy of the developed code is certificated by several tests, in which we show comparisons with the corresponding analytical solutions. Using hydrodynamical data in our collapsar simulation, we estimate the annihilation rates in a post-processing manner. We show that the neutrino energy deposition and momentum transfers are strongest near the inner edge of the accretion disk. The beaming effects of special relativity are found to change the annihilation rates by several factors in the polar funnel region. After the accretion disk settles into a stationary state (typically later than $\sim 9$ s from the onset of gravitational collapse), we find that the neutrino-heating timescale in the vicinity of the polar funnel ($\lesssim 80$ km) can become shorter than the hydrodynamical timescale, indicating that the neutrino-heated outflows can be launched there. We point out that the momentum transfer can play as important role as the energy deposition for the efficient acceleration of neutrino-driven outflows. Our results suggest that the neutrino pair annihilation has a potential importance equal to the conventional magnetohydrodynamic mechanism for igniting the GRB fireballs.

preprint2007arXiv

Gravitational Wave Background from Population III Stars

We estimate the gravitational wave (GW) background from Population III (Pop III) stars using the results from our hydrodynamic simulations (Suwa et al. 2007). We calculate gravitational waveforms from matter motions and anisotropic neutrino emissions for single Pop III stars. We find that the GW amplitudes from matter motions are dominant until just after bounce, but those from neutrinos dominate later on at frequencies below $\sim 10$ Hz in the GW spectrum. Computing the overall signal produced by the ensemble of such Pop III stars, we find that the resultant density parameter of the GW background peaks at the amplitude of $Ω_\mathrm{GW}\sim 10^{-10}$ in the frequency interval $\sim 1-10$ Hz. We show that such signals, depending on the formation rate of Pop III stars, can be within the detection limits of future planned interferometers such as DECIGO and BBO in the frequency interval of $\sim 0.1-1$ Hz. Our results suggest that the detection of the GW background from Pop III stars can be an important tool to supply the information about the star formation history in the early universe.

preprint2007arXiv

Magnetorotational Collapse of Population III Stars

We perform a series of two-dimensional magnetorotational core-collapse simulations of Pop III stars. Changing the initial distributions of rotation and magnetic fields prior to collapse in a parametric manner, we compute 19 models. By so doing, we systematically investigate how rotation and magnetic fields affect the collapse dynamics and explore how the properties of the black-hole formations and neutrino emissions could be affected. As for the microphysics, we employ a realistic equation of state and approximate the neutrino transfer by a multiflavour leakage scheme. With these computations, we find that the jet-like explosions are obtained by the magnetodriven shock waves if the initial magnetic field is as large as $10^{12}$ G. We point out that the black-hole masses at the formation decrease with the initial field strength, on the other hand, increase with the initial rotation rates. As for the neutrino properties, we point out that the degree of the differential rotation plays an important role to determine which species of the neutrino luminosity is more dominant than the others. Furthermore, we find that the stronger magnetic fields make the peak neutrino luminosities smaller, because the magnetic pressure acts to halt the collapse in the central regions, leading to the suppression of the releasable gravitational binding energies.