BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hiroki Nagakura

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

ResearcherUnclaimed source record
astro-ph.HEhep-phgr-qcnucl-thastro-ph.SRastro-ph.IMhep-exnucl-exphysics.data-an

Catalog footprint

What is connected

36works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

36 published item(s)

preprint2024arXiv

Collisional and Fast Neutrino Flavor Instabilities in Two-dimensional Core-collapse Supernova Simulation with Boltzmann Neutrino Transport

We present a comprehensive study on the occurrences of the collisional flavor instability (CFI) and the fast flavor instability (FFI) of neutrinos based on a two-dimensional (2D) core-collapse supernova (CCSN) simulation performed with a Boltzmann radiation hydrodynamics code. We find that CFI occurs in a region with the baryon-mass density of $10^{10}\lesssim ρ\lesssim 10^{12}\,\mathrm{g}\,\mathrm{cm}^{-3}$, which is similar to the previous results in one-dimensional (1D) CCSN models. In contrast to 1D, however, the CFI region varies with time vigorously in the 2D model, whereas it had a quiescent structure in 1D. This is attributed to the fact that the turbulent flows advected from a gain region account for the temporal variations. Another noticeable difference from the 1D models is the appearance of resonance-like CFI where number densities of $ν_e$, $\barν_e$ nearly coincide each other. The CFI growth rate there is enhanced and can reach $\sim10^{-3}\,\mathrm{cm}^{-1}$. As for FFI, on the other hand, it appears in three different regions; (1) the region overlapped with the resonance-like CFI, (2) neutrino decoupling regions where $\barν_e$'s are strongly emitted, and (3) optically thin regions where neutral-current scatterings dominate over charged-current reactions. Although overall properties for FFI are consistent with previous studies, we find that the number of electron-neutrinos lepton number crossing (ELN crossing) temporary becomes multiple, which can be assessed accurately only by multi-angle treatments in neutrino transport. We find that the growth rate of FFI is always higher than CFI if both of them occur, which suggests that the former is dominant for the linear evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Detectability of Late-time Supernova Neutrinos with Fallback Accretion onto Protoneutron star

We investigate the late-time neutrino emission powered by fallback mass accretion onto proto-neutron star (PNS), using neutrino radiation-hydrodynamic simulations with full Boltzmann neutrino transport. We follow the time evolution of accretion flow onto PNS until the system reaches a steady state. A standing shock wave is commonly formed in the accretion flow, whereas the shock radius varies depending on mass accretion rate and PNS mass. A sharp increase in temperature emerges in the vicinity of PNS ($\sim 10$ km), which characterizes neutrino emission. Both neutrino luminosity and average energy become higher with increasing mass accretion rate and PNS mass. The mean energy of emitted neutrinos is in the range of $10\lesssimε\lesssim20\,\mathrm{MeV}$, which is higher than that estimated from PNS cooling models ($\lesssim10\,\mathrm{MeV}$). Assuming a distance to core-collapse supernova of $10\,\mathrm{kpc}$, we quantify neutrino event rates for Super-Kamiokande (Super-K) and DUNE. The estimated detection rates are well above the background, and their energy-dependent features are qualitatively different from those expected from PNS cooling models. Another notable feature is that the neutrino emission is strongly flavor dependent, exhibiting that the neutrino event rate hinges on the neutrino oscillation model. We estimate them in the case with adiabatic Mikheev-Smirnov-Wolfenstein model, and show that the normal- and inverted mass hierarchy offer the large number of neutrino detection in Super-K and DUNE, respectively. Hence the simultaneous observation with Super-K and DUNE of the fallback neutrinos will provide a strong constraint on neutrino mass hierarchy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effects of nuclear matter and composition in core-collapse supernovae and long-term proto-neutron star cooling

We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and nucleons by using the same EOS by the variational method but employing two different treatments in computations of nuclear abundances. We perform numerical simulations of core-collapse supernovae adopting the three EOSs. We also perform numerical simulations of the long-term evolution over 70 s of the proto-neutron star cooling. We show that impacts by different modeling of composition are remarkable as in those by different treatments of uniform matter in the gravitational collapse, bounce, and shock propagation. The cooling of proto-neutron star and the resulting neutrino emission are also affected by the compositional difference even if the same treatment in computing uniform matter of EOS.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient estimation method for time evolution of proto-neutron star mass and radius from supernova neutrino signal

In this paper we present a novel method to estimate the time evolution of proto-neutron star (PNS) structure from the neutrino signal in core-collapse supernovae (CCSN). Employing recent results of multi-dimensional CCSN simulations, we delve into a relation between total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting TONE to PNS mass and radius as a function of time. By combining another fitting function representing the correlation between TONE and cumulative number of event at each neutrino observatory, PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of CCSN neutrino signal.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Explosive nucleosynthesis with fast neutrino-flavor conversion in core-collapse supernovae

Fast neutrino-flavor conversion (FFC) is a possible game-changing ingredient in core-collapse supernova (CCSN) theory. In this paper, we examine the impact of FFC on explosive nucleosynthesis by including the effects of FFC in conjunction with asymmetric neutrino emission into nucleosynthetic computations in a parametric way. We find that the ejecta compositions are not appreciably affected by FFC for elements lighter than Co, while the compositions are influenced by FFC for the heavier elements. We also find that the role of FFC varies depending on the asymmetric degree of neutrino emission ($m_{\rm asy}$) and the degree of neutrino-flavor mixing. The impact of FFC is not monotonic to $m_{\rm asy}$; The change in the ejecta composition increases for higher $m_{\rm asy}$ up to $\sim 10\%$ compared with that without FFC, whereas FFC has little effect on the nucleosynthesis in very large asymmetric neutrino emission ($\gtrsim 30 \%$). Our results suggest that FFC facilitates the production of neutron-rich ejecta in most cases, although it makes the ejecta more proton-rich if anti-neutrino conversion is more vigorous than that of neutrino. The key ingredient accounting for this trend is neutrino absorption, whose effects on nucleosynthesis can be quantified by simple diagnostics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Self-consistency in models of neutrino scattering and fast flavor conversion

Several recent numerical studies have examined the effects of neutrino neutral-current scattering on fast flavor conversion (FFC). Those studies are in apparent conflict, with some finding enhancement of flavor conversion and others finding suppression. The ones that report enhancement all use homogeneous models, and we discuss in detail the self-consistency issues that they face as a result. We reproduce scattering-enhanced FFC in our own homogeneous calculations, showing that it occurs for both neutral- and charged-current scattering and that it is due to isotropization of the angular distributions over time. Because this is the very feature of the calculations that is not self-consistent, we conclude that the enhancement effect is of unclear astrophysical relevance and may not occur in natural environments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-thermal neutrinos created by shock acceleration in successful and failed core-collapse supernova

We present a comprehensive study of neutrino shock acceleration in core-collapse supernova (CCSN). The leading players are heavy leptonic neutrinos, $ν_μ$ and $ν_τ$; the former and latter potentially gain the energy up to $\sim 100$ MeV and $\sim 200$ MeV, respectively, through the shock acceleration. Demonstrating the neutrino shock acceleration by Monte Carlo neutrino transport, we make a statement that it commonly occurs in the early post bounce phase ($\lesssim 50$ ms after bounce) for all massive stellar collapse experiencing nuclear bounce and would reoccur in the late phase ($\gtrsim 100$ ms) for failed CCSNe. This opens up a new possibility to detect high energy neutrinos by terrestrial detectors from Galactic CCSNe; hence, we estimate the event counts for Hyper(Super)-Kamiokande, DUNE, and JUNO. We find that the event count with the energy of $\gtrsim 80$ MeV is a few orders of magnitude higher than that of the thermal neutrinos regardless of the detectors, and muon production may also happen in these detectors by $ν_μ$ with the energy of $\gtrsim 100$ MeV. The neutrino signals provide a precious information on deciphering the inner dynamics of CCSN and placing a constraint on the physics of neutrino oscillation; indeed, the detection of the high energy neutrinos through charged current reaction channels will be a smoking gun evidence of neutrino flavor conversion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Prospects of fast flavor neutrino conversion in rotating core-collapse supernovae

There is mounting evidence that neutrinos undergo fast flavor conversion (FFC) in core-collapse supernova (CCSN). In this letter, we investigate the roles of stellar rotation on the occurrence of FFC by carrying out axisymmetric CCSN simulations with full Boltzmann neutrino transport. Our result suggests that electron neutrino lepton number (ELN) angular crossings, which are the necessary and sufficient condition to trigger FFC, preferably occur in the equatorial region for rotating CCSNe. By scrutinizing the neutrino--matter interaction and neutrino radiation field, we find some pieces of evidence that the stellar rotation facilitates the occurrence of FFC. The low-electron-fraction region in the post-shock layer expands by centrifugal force, enhancing the disparity of neutrino absorption between electron-type neutrinos ($ν_{\rm e}$) and their anti-particles ($\barν_{\rm e}$). This has a significant impact on the angular distribution of neutrinos in momentum space, in which $ν_{\rm e}$ tends to be more isotropic than $\barν_{\rm e}$; consequently, ELN crossings emerge. The ELN crossing found in this study is clearly associated with rotation, which motivates further investigation on how the subsequent FFC influences explosion dynamics, nucleosynthesis, and neutrino signals in rotating CCSNe.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The impact of asymmetric neutrino emissions on nucleosynthesis in core-collapse supernovae II -- progenitor dependences --

We investigate the impact of asymmetric neutrino-emissions on explosive nucleosynthesis in core-collapse supernovae (CCSNe) of progenitors with a mass range of 9.5 to 25$M_{\odot}$. We perform axisymmetric, hydrodynamic simulations of the CCSN explosion with a simplified neutrino-transport, in which anti-correlated dipolar emissions of $ν_{\rm e}$ and ${\bar ν}_{\rm e}$ are imposed. We then evaluate abundances and masses of the CCSN ejecta in a post-processing manner. We find that the asymmetric $ν$-emission leads to the abundant ejection of $p$- and $n$-rich matter in the high-$ν_{\rm e}$ and -${\bar ν}_{\rm e}$ hemispheres, respectively. It substantially affects the abundances of the ejecta for elements heavier than Ni regardless of progenitors, although those elements lighter than Ca are less sensitive. Based on these results, we calculate the IMF-averaged abundances of the CCSN ejecta with taking into account the contribution from Type Ia SNe. For $m_{\rm asy} = 10/3\%$ and $10\%$, where $m_{\rm asy}$ denotes the asymmetric degree of the dipole components in the neutrino emissions, the averaged abundances for elements lighter than Y are comparable to those of the solar abundances, whereas those of elements heavier than Ge are overproduced in the case with $m_{\rm asy} \ge 30\%$. Our result also suggests that the effect of the asymmetric neutrino emissions is imprinted in the difference of abundance ratio of [Ni/Fe] and [Zn/Fe] between the high-$ν_{\rm e}$ and -${\bar ν}_{\rm e}$ hemispheres, indicating that the future spectroscopic X-ray observations of a CCSN remnant will bring evidence of the asymmetric neutrino emissions if exist.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations

This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical Fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behavior emerges. We find that the vigor of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS $-$ heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of $ν_μ$, $ν_τ$, $\barν_μ$, and $\barν_τ$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $ν_e$ and $\barν_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $ν_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutrino oscillations in supernovae: angular moments and fast instabilities

Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero-crossing in $n_{ν_e} - n_{\barν_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using Fornax simulation data, we show that fast collective oscillations qualitatively depend on how forward-peaked the neutrino angular distributions are.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutrino transport with Monte Carlo method: I. Towards fully consistent implementation of nucleon recoils in core-collapse supernova simulations

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The sensitivity of presupernova neutrinos to stellar evolution models

We examine the sensitivity of neutrino emissions to stellar evolution models for a 15$M_\odot$ progenitor, paying particular attention to a phase prior to the collapse. We demonstrate that the number luminosities in both electron-type neutrinos ($ν_e$) and their anti-partners ($\barν_e$) differ by more than an order of magnitude by changing spatial resolutions and nuclear network sizes on stellar evolution models. We also develop a phenomenological model to capture the essential trend of the diversity, in which neutrino luminosities are expressed as a function of central density, temperature and electron fraction. In the analysis, we show that neutrino luminosity can be well characterized by these central quantities. This analysis also reveals that the most influential quantity to the time evolution of $ν_e$ luminosity is matter density, while it is temperature for $\barν_e$. These qualitative trends will be useful and applicable to constrain the physical state of progenitors at the final stages of stellar evolution from future neutrino observations, although more detailed systematic studies including various mass progenitors are required to assess the applicability.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A new possibility of the fast neutrino-flavor conversion in the pre-shock region of core-collapse supernova

We make a strong case that the fast neutrino-flavor conversion, one of the collective flavor oscillation modes, commonly occurs in core-collapse supernovae (CCSNe). It is confirmed in the numerical data obtained in realistic simulations of CCSNe but the argument is much more generic and applicable universally: the coherent neutrino-nucleus scattering makes the electron lepton number (ELN) change signs at some inward direction and trigger the flavor conversion in the outward direction in the pre-shock region. Although the ELN crossing is tiny and that is why it has eluded recognition so far, it is still large enough to induce the flavor conversion. Our findings will have an important observational consequences for CCSNe neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fast-pairwise collective neutrino oscillations associated with asymmetric neutrino emissions in core-collapse supernova

We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos ($ν_{\rm e}$) and their anti-particles ($\barν_{\rm e}$), in which the asymmetry of $ν_{\rm e}$ and $\barν_{\rm e}$ is anti-correlated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter they appear only in the hemisphere of higher $\barν_{\rm e}$ emissions (the same hemisphere with stronger shock expansion). We analyze in depth the characteristics of electron-lepton-number (ELN) crossing by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, $\barν_{\rm e}$ ($ν_{\rm e}$) dominates over $ν_{\rm e}$ ($\barν_{\rm e}$) in the forward (backward) direction: at the larger radius the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-D transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modelling of CCSNe.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Overarching Framework of Core-Collapse Supernova Explosions as Revealed by 3D Fornax Simulations

We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular vis à vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y$_e = 0.5$, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hyperbolic Self-Gravity Solver for Large Scale Hydrodynamical Simulations

A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optical Synchrotron Precursors of Radio Hypernovae

We examine the bright radio synchrotron counterparts of low-luminosity gamma-ray bursts (llGRBs) and relativistic supernovae (SNe) and find that they can be powered by spherical hypernova (HN) explosions. Our results imply that radio-bright HNe are driven by relativistic jets that are choked deep inside the progenitor stars or quasi-spherical magnetized winds from fast-rotating magnetars. We also consider the optical synchrotron counterparts of radio-bright HNe and show that they can be observed as precursors several days before the SN peak with an r-band absolute magnitude of M_r ~ -14 mag. While previous studies suggested that additional trans-relativistic components are required to power the bright radio emission, we find that they overestimated the energy budget of the trans-relativistic component by overlooking some factors related to the minimum energy of non-thermal electrons. If an additional trans-relativistic component exists, then a much brighter optical precursor with M_r ~ -20 mag can be expected. Thus, the scenarios of radio-bright HNe can be distinguished by using optical precursors, which can be detectable from < 100 Mpc by current SN surveys like the Kiso SN Survey, Palomar Transient Factory, and Panoramic Survey Telescope & Rapid Response System.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical Surface for Explosions of Rotational Core-Collapse Supernovae

The effect of rotation on the explosion of core-collapse supernovae is investigated systematically in three-dimensional simulations. In order to obtain the critical conditions for explosion as a function of mass accretion rate, neutrino luminosity, and specific angular momentum, rigidly rotating matter was injected from the outer boundary with an angular momentum, which is increased every 500 ms. It is found that there is a critical value of the specific angular momentum, above which the standing shock wave revives, for a given combination of mass accretion rate and neutrino luminosity, i.e. an explosion can occur by rotation even if the neutrino luminosity is lower than the critical value for a given mass accretion rate in non-rotational models. The coupling of rotation and hydrodynamical instabilities plays an important role to characterize the dynamics of shock revival for the range of specific angular momentum that are supposed to be realistic. Contrary to expectations from past studies, the most rapidly expanding direction of the shock wave is not aligned with the rotation axis. Being perpendicular to the rotation axis on average, it can be oriented in various directions. Its dispersion is small when the spiral mode of the standing accretion shock instability (SASI) governs the dynamics, while it is large when neutrino-driven convection is dominant. As a result of the comparison between 2D and 3D rotational models, it is found that m=!0 modes of neutrino-driven convection or SASI are important for shock revival around the critical surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Jet collimation in the ejecta of double neutron star merger: new canonical picture of short gamma-ray bursts

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10°. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parametric Study of Flow Patterns behind the Standing Accretion Shock Wave for Core-Collapse Supernovae

In this study, we conduct three-dimensional hydrodynamic simulations systematically to investigate the flow patterns behind the accretion shock waves that are commonly formed in the post-bounce phase of core-collapse supernovae. Adding small perturbations to spherically symmetric, steady, shocked accretion flows, we compute the subsequent evolutions to find what flow pattern emerges as a consequence of hydrodynamical instabilities such as convection and standing accretion shock instability (SASI) for different neutrino luminosities and mass accretion rates. Depending on these two controlling parameters, various flow patterns are indeed realized. We classify them into three basic patterns and two intermediate ones; the former includes sloshing motion (SL), spiral motion (SP) and multiple buoyant bubble formation (BB); the latter consists of spiral motion with buoyant-bubble formation (SPB) and spiral motion with pulsationally changing rotational velocities (SPP). Although the post-shock flow is highly chaotic, there is a clear trend in the pattern realization. The sloshing and spiral motions tend to be dominant for high accretion rates and low neutrino luminosities, and multiple buoyant bubbles prevail for low accretion rates and high neutrino luminosities. It is interesting that the dominant pattern is not always identical between the semi-nonlinear and nonlinear phases near the critical luminosity; the intermediate cases are realized in the latter case. Running several simulations with different random perturbations, we confirm that the realization of flow pattern is robust in most cases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of advection and collision terms in the Boltzmann equations is the source of difficulties, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce and shock-stall for a 15M_{sun} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects in our method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Post-shock-revival evolutions in the neutrino-heating mechanism of core-collapse supernovae

We perform some experimental simulations in spherical symmetry and axisymmetry to understand the post-shock-revival evolution of core-collapse supernovae. Assuming that the stalled shock wave is relaunched by neutrino heating and employing the so-called light bulb approximation, we induce shock revival by raising the neutrino luminosity by hand up to the critical value, which is also de- termined by dynamical simulations. A 15M_{sun} progenitor model is employed. We incorporate nuclear network calculations with a consistent equation of state in the simulations to account for the energy release by nuclear reactions and their feedback to hydrodynamics. Varying the shock-relaunch time rather arbitrarily, we investigate the ensuing long-term evolutions systematically, paying particular attention to the explosion energy and nucleosynthetic yields as a function of this relaunch time, or equivalently the accretion rate at shock revival. We study in detail how the diagnostic explosion energy approaches the asymptotic value and which physical processes contribute to the explosion energy in what proportions as well as their dependence on the relaunch time and the dimension of dynamics. We find that the contribution of nuclear reactions to the explosion energy is comparable to or greater than that of neutrino heating. In particular, recombinations are dominant over burnings in the contributions of nuclear reactions. Interestingly 1D models studied in this paper cannot produce the appropriate explosion energy and nickel mass simultaneously, overproducing nickels, whereas this problem is resolved in 2D models if the shock is relaunched at 300-400ms after bounce.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Semi-dynamical approach to the shock revival in core-collapse supernovae

We develop a new semi-dynamical method to study shock revival by neutrino heating in core- collapse supernovae. Our new approach is an extension of the previous studies that employ spherically symmetric, steady, shocked accretion flows together with the light bulb approximation. The latter has been widely used in the supernova community for the phenomenological investigation of the criteria for successful supernova explosions. In the present approach, on the other hand, we get rid of the steady-state condition and take into account shock wave motions instead. We have in mind the scenario that not the critical luminosity but the critical fluctuation generated by hydrodynamical instabilities such as SASI and neutrino-driven convection in the post-shock region determines the onset of shock revival. After confirming that the new approach indeed captures the dynamics of revived shock wave qualitatively, we then apply the method to various initial conditions and find that there is a critical fluctuation for shock revival, which can be well fit by the following formula: f_crit ~ 0.8 * (M_in/1.4M_sun) * {1- (rsh/10^8cm)}, in which fcrit denotes the critical pressure fluctuation normalized by the unperturbed post-shock value. Min and rsh stand for the mass of the central compact object and the shock radius, respectively. The critical fluctuation decreases with the shock radius, whereas it increases with the mass of the central object. We discuss the possible implications of our results for 3D effects on shock revival, which is currently controversial in the supernova community.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability in core-collapse supernovae

We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Population III Gamma-Ray Bursts and Breakout Criteria for Accretion-Powered Jets

We investigate the propagation of accretion-powered jets in various types of massive stars such as Wolf-Rayet stars, light Population III (Pop III) stars, and massive Pop III stars, all of which are the progenitor candidates of Gamma-Ray Bursts (GRBs). We perform two dimensional axisymmetric simulations of relativistic hydrodynamics taking into account both the envelope collapse and the jet propagation (i.e., the negative feedback of the jet on the accretion). Based on our hydrodynamic simulations, we show for the first time that the accretion-powered jet can potentially break out relativistically from the outer layers of Pop III progenitors. In our simulations, the accretion rate is estimated by the mass flux going through the inner boundary, and the jet is injected with a fixed accretion-to-jet conversion efficiency $η$. By varying the efficiency $η$ and opening angle $θ_{op}$ for more than 40 models, we find that the jet can make a relativistic breakout from all types of progenitors for GRBs if a simple condition $η\gtrsim 10^{-4} (θ_{op}/8^{\circ})^2$ is satisfied, which is consistent with analytical estimates. Otherwise no explosion or some failed spherical explosions occur.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The propagation of neutrino-driven jets in Wolf-Rayet stars

We numerically investigate the jet propagation through a rotating collapsing Wolf-Rayet star with detailed central engine physics constructed based on the neutrino-driven collapsar model. The collapsing star determines the evolution of mass accretion rate, black hole mass and spin, all of which are important ingredients for determining the jet luminosity. We reveal that neutrino-driven jets in rapidly spinning Wolf-Rayet stars are capable of breaking out from the stellar envelope, while those propagating in slower rotating progenitors fail to jet breakout due to insufficient kinetic power. For progenitor models with successful jet breakouts, the kinetic energy accumulated in the cocoon could be as large as ~1051erg and might significantly contribute to the luminosity of the afterglow emission or to the kinetic energy of the accompanying supernova if nickel production takes place. We further analyze the post breakout phase using a simple analytical prescription and conclude that the relativistic jet component could produce events with an isotropic-luminosity Lp(iso) ~1052erg/s and isotropic-energy Ej(iso) ~1054erg. Our findings support the idea of rapidly rotating Wolf-Rayet stars as plausible progenitors of GRBs, while slowly rotational ones could be responsible for low luminosity or failed GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Jet propagations, breakouts and photospheric emissions in collapsing massive progenitors of long duration gamma ray bursts

We investigate by two-dimensional axisymmetric relativistic hydrodynamical simulations (1) jet propagations through an envelope of a rapidly rotating and collapsing massive star, which is supposed to be a progenitor of long duration gamma ray bursts (GRBs), (2) breakouts and subsequent expansions into stellar winds and (3) accompanying photospheric emissions. We find that if the envelope rotates uniformly almost at the mass shedding limit, its outer part stops contracting eventually when the centrifugal force becomes large enough. Then another shock wave is formed, propagates outwards and breaks out of the envelope into the stellar wind. Which breaks out earlier, the jet or the centrifugal bounce-induced shock, depends on the timing of jet injection. If the shock breakout occurs earlier owing to a later injection, the jet propagation and subsequent photospheric emissions are affected substantially. We pay particular attention to observational consequences of the difference in the timing of jet injection. We calculate optical depths to find the location of photospheres, extracting densities and temperatures at appropriate retarded times from the hydrodynamical data. We show that the luminosity and observed temperature of the photospheric emissions are both much lower than those reported in previous studies. Although luminosities are still high enough for GRBs, the observed temperature are lower than the energy at the spectral peak expected by the Yonetoku-relation. This may imply that energy exchanges between photons and matter are terminated deeper inside or some non-thermal processes are operating to boost photon energies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Direct time radio variability induced by non-axisymmetric standing accretion shock instability: implications for M87

We show for the first time the direct time variable radio images in the context of shocked accretion flows around a black hole under the general relativistic treatment of both hydrodynamics and radiation transfer. Time variability around a black hole can be induced by the non-axisymmetric standing accretion shock instability (namely Black Hole SASI). Since the spiral arm shock waves generate the density and temperature waves at the post shock region, they cause time variability in the black hole vicinity. Based on our dynamical simulations, we discuss a possibility of detection for the time variable radio images of M87 by the future space telescope VSOP2/ASTRO-G satellite. The most luminous part of the images is predicted to be near 15 Schwarzschild radii for some snapshots. We show that our results are consistent with existing observational data such as timeaveraged radio spectra, VLBA images and variability timescale for M87. We also discuss observations of M87 with millimeter and sub-millimeter interferometers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-layered configurations in differentially-rotational equilibrium

We present a new formula to numerically construct configurations in rotational equilibrium, which consist of multiple layers. Each layer rotates uniformly or differentially according to cylindrical rotation-laws that are different from layer to layer. Assuming a different barotropic equation of state (EOS) for each layer, we solve the Bernoulli equation in each layer separately and combine the solutions by imposing continuity of the pressure at each boundary of the layers. It is confirmed that a single continuous barotropic EOS is incompatible with the junction condition. Identifying appropriate variables to be solved, we construct a convergent iteration scheme. For demonstration, we obtain two-layered configurations, each layer of which rotates rapidly with either an "$Ω$-constant law" or a "$j$-constant law" or a "$v$-constant law". Other rotation laws and/or a larger number of layers can be treated similarly. We hope that this formula will be useful in studying the stellar evolution in multi-dimension with the non-spherical configuration induced by rotation being fully taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic Dissipative Accretion Flow onto Black Hole

Dissipations, e.g. heat flow and bulk and shear viscosities, cause the transport of energy, momentum and angular momentum, which is the essence of accretion of matters onto celestial objects. Dissipations are usually described by the Fourier and Navier-Stokes laws ("classic laws" of dissipations). However the classic laws result in an infinitely fast propagation of dissipations. In relativistic formulation, the classic laws of dissipations violate the causality, and hence no relativistic theory of accretion flow onto celestial object is formulated. In this short report, we summarize the causal dissipative hydrodynamics, so-called "Extended Irreversible Thermodynamics" (EIT), with a supplemental comment of which the original works of EIT are not aware, and then show two theorems about relativistic dissipative flows around a Schwarzschild black hole. By these theorems, a significant property of EIT in contrast with classic laws of dissipations is clarified, and a dissipative instability of an exact solution of relativistic perfect fluid flow is also obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hysteresis of Backflow Imprinted in Collimated Jets

We report two different types of backflow from jets by performing 2D special relativistic hydrodynamical simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), and the other is bent path of the backflow (bent backflow). We find that the former appears when the head advance speed is comparable to or higher than the local sound speed at the hotspot while the latter appears when the head advance speed is slower than the sound speed bat the hotspot. Bent backflow collides with the unshocked jet and laterally squeezes the jet. At the same time, a pair of new oblique shocks are formed at the tip of the jet and new bent fast backflows are generated via these oblique shocks. The hysteresis of backflow collisions is thus imprinted in the jet as a node and anti-node structure. This process also promotes broadening of the jet cross sectional area and it also causes a decrease in the head advance velocity. This hydrodynamic process may be tested by observations of compact young jets.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Standing Accretion Shock Instability in the Disk around the Kerr Black Hole

This paper is a sequel to our previous work for accretion onto a Schwarzschild black hole and the so-called standing accretion shock instability (SASI), in this paper we investigate non-axisymmetric perturbations for a Kerr black hole. The linear and non-linear phases for the shock evolution are analyzed in detail by both 2D general relativistic hydrodynamical simulations and linear analysis. Since the structure of steady axisymmetric accretion flows with a standing shock wave is very sensitive to the inner transonic flow, their properties such as Mach numbers, which are important for the stability, depend on the Kerr parameter very much. Although the essential features of the instability do not differ from the previous results for the Schwarzschild black hole, the frame dragging effects specific to the Kerr black hole is also evident. Interestingly, the oscillation periods of the fundamental unstable modes are dependent only on the shock radius irrespective of the injection parameters.

0 citations0 reviewsNo code linkedOpen access