Showing 1–50 of 105 results for author: Takiwaki, T

Universality in supernova gravitational waves with proto-neutron star properties

Authors: Hajime Sotani, Bernhard Müller, Tomoya Takiwaki

Abstract: Gravitational wave signals from core-collapse supernovae are one of the important observables for extracting the information of dense matter. To extract the properties of proto-neutron stars produced via core-collapse supernovae by asteroseismology, we perform a linear perturbation analysis using data obtained from two-dimensional numerical simulations. We employ 12 and 20 solar-mass progenitors a… ▽ More Gravitational wave signals from core-collapse supernovae are one of the important observables for extracting the information of dense matter. To extract the properties of proto-neutron stars produced via core-collapse supernovae by asteroseismology, we perform a linear perturbation analysis using data obtained from two-dimensional numerical simulations. We employ 12 and 20 solar-mass progenitors and compare two different treatments of gravity. One is a general relativistic one with a conformal flatness condition and the other is an effective gravitational potential mimicking the Tolman-Oppenheimer-Volkoff solution. We discuss how the frequencies of the proto-neutron star oscillations corresponding to the gravitational wave signals in the simulations depend on the proto-neutron star properties. In our models, we find that the gravitational wave frequencies of the proto-neutron stars determined with the Cowling approximation can be expressed to very good approximation as a function of the proto-neutron star average density almost independently of the progenitor mass, treatment of gravity in the simulations, and the interpolations in the simulations. On the other hand, if one considers the gravitational wave frequencies as a function of the surface gravity of proto-neutron stars, such a relation appears sensitive to the treatment of gravity and other numerical details in the simulations. Thus, the average density of proto-neutron stars seems more suitable for universally expressing the supernova gravitational wave frequencies, instead of the surface gravity. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: Accepted for publication in RPD

Report number: RIKEN-iTHEMS-Report-24

Three-dimensional Magneto-hydrodynamic Simulations of Core-collapse Supernovae: I. Hydrodynamic evolution and protoneutron star properties

Authors: Ko Nakamura, Tomoya Takiwaki, Jin Matsumoto, Kei Kotake

Abstract: We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of sixteen progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 $M_{\odot}$. The examined progenitors cover a wide range of the compactness parameter including a peak around $23 M_{\odot}$. We find that neutri… ▽ More We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of sixteen progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 $M_{\odot}$. The examined progenitors cover a wide range of the compactness parameter including a peak around $23 M_{\odot}$. We find that neutrino-driven explosions occur for all models within 0.3 s after bounce. We also find that the properties of the explosions and the central remnants are well correlated with the compactness. Early shock evolution is sensitive to the mass accretion rate onto the central core, reflecting the density profile of the progenitor stars. The most powerful explosions with diagnostic explosion energy $E_{\rm exp} \sim 0.75 \times 10^{51}$ erg are obtained by 23 and 24 $M_{\odot}$ models, which have the highest compactness among the examined models. These two models exhibit spiral SASI motions during 150-230 ms after bounce preceding a runaway shock expansion and leave a rapidly rotating neutron star with spin periods $\sim 50$ ms. Our models predict the gravitational masses of the neutron star ranging between $1.22 M_{\odot}$ and $1.67 M_{\odot}$ and their spin periods 0.04-4 s. The number distribution of these values roughly matches observation. On the other hand, our models predict small hydrodynamic kick velocity (15-260 km/s), although they are still growing at the end of our simulations. Further systematic studies, including rotation and binary effects, as well as long-term simulations up to several seconds, will enable us to explore the origin of various core-collapse supernova explosions. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 15 pages, 13+5 figures, 1 appendix

Two-dimensional models of core-collapse supernova explosions assisted by heavy sterile neutrinos

Authors: Kanji Mori, Tomoya Takiwaki, Kei Kotake, Shunsaku Horiuchi

Abstract: Core-collapse supernovae can be a copious source of sterile neutrinos, hypothetical particles that mix with active neutrinos. We develop two-dimensional stellar core-collapse models that incorporate the mixing between tau neutrinos and heavy sterile neutrinos -- those with the mass of 150--200 MeV -- to investigate signatures of sterile neutrinos in supernova observables. We find that the decay ch… ▽ More Core-collapse supernovae can be a copious source of sterile neutrinos, hypothetical particles that mix with active neutrinos. We develop two-dimensional stellar core-collapse models that incorporate the mixing between tau neutrinos and heavy sterile neutrinos -- those with the mass of 150--200 MeV -- to investigate signatures of sterile neutrinos in supernova observables. We find that the decay channel of a sterile neutrino into a pion and a tau neutrino can enhance the explosion energy and the synthesized nickel mass. Although the inclusion of sterile neutrinos considered in this study slightly reduce the neutrino and gravitational-wave signals, we find that they are still detectable for a Galactic event. Furthermore, we point out that if sterile neutrinos are as massive as ~200 MeV, they produce high-energy tau antineutrinos with energies of ~80 MeV, the detection of which can be a smoking signature of the sterile neutrinos and where Hyper-Kamiokande should play a pivotal role. △ Less

Submitted 14 June, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 18 pages, 14 figures, 3 tables, accepted for publication in PRD

Systematic local simulations of fast neutrino flavor conversions with scattering effects

Authors: Milad Delfan Azari, Hirokazu Sasaki, Tomoya Takiwaki, Hirotada Okawa

Abstract: We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) inhomogeneous and the homogeneous models as post processes by employing snapshots obtained by our self-consistent, realistic Boltzmann simulations in two spatial dimensions (2D). We show that the FFC growth rate is considerably larger in the inhomogeneous model than in the homogeneous model, as expec… ▽ More We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) inhomogeneous and the homogeneous models as post processes by employing snapshots obtained by our self-consistent, realistic Boltzmann simulations in two spatial dimensions (2D). We show that the FFC growth rate is considerably larger in the inhomogeneous model than in the homogeneous model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability under inhomogeneous conditions. On the other hand, in the homogeneous model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collision-induced enhancement occurs on a long time scale in the inhomogeneous model. Despite that collision-induced enhancement does not appear on the short time scale, that should be taken into account to predict the final conversion probability. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 15 pages, 16 figures

Report number: LA-UR-24-20443

Fate of supernova progenitors in massive binary systems

Authors: Tomoya Kinugawa, Shunsaku Horiuchi, Tomoya Takiwaki, Kei Kotake

Abstract: How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate… ▽ More How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core-collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion induced collapse (AIC), rapidly rotating supernova (RSN), black hole formation, and gamma ray burst (GRB). We then conduct Binary Population Synthesis calculations for 12 models, varying the initial parameter distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.14$--$1.57) \times10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, RSN, and GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated GRB rates are higher than the observed long GRB rate, but very close to the low luminosity GRB rate. Furthering binary modeling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 21 pages, comments welcome

Spin-flavor precession of Dirac neutrinos in dense matter and its potential in core-collapse supernovae

Authors: Hirokazu Sasaki, Tomoya Takiwaki, A. Baha Balantekin

Abstract: We calculate the spin-flavor precession (SFP) of Dirac neutrinos induced by strong magnetic fields and finite neutrino magnetic moments in dense matter. As found in the case of Majorana neutrinos, the SFP of Dirac neutrinos is enhanced by the large magnetic field potential and suppressed by large matter potentials composed of the baryon density and the electron fraction. The SFP is possible irresp… ▽ More We calculate the spin-flavor precession (SFP) of Dirac neutrinos induced by strong magnetic fields and finite neutrino magnetic moments in dense matter. As found in the case of Majorana neutrinos, the SFP of Dirac neutrinos is enhanced by the large magnetic field potential and suppressed by large matter potentials composed of the baryon density and the electron fraction. The SFP is possible irrespective of the large baryon density when the electron fraction is close to 1/3. The diagonal neutrino magnetic moments that are prohibited for Majorana neutrinos enable the spin precession of Dirac neutrinos without any flavor mixing. With supernova hydrodynamics simulation data, we discuss the possibility of the SFP of both Dirac and Majorana neutrinos in core-collapse supernovae. The SFP of Dirac neutrinos occurs at a radius where the electron fraction is 1/3. The required magnetic field of the proto-neutron star for the SFP is a few $10^{14}$G at any explosion time. For the Majorana neutrinos, the required magnetic field fluctuates from $10^{13}$G to $10^{15}$G. Such a fluctuation of the magnetic field is more sensitive to the numerical scheme of the neutrino transport in the supernova simulation. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 14 pages, 10 figures

Report number: LA-UR-23-30314

Three-dimensional GRMHD Simulations of Rapidly Rotating Stellar Core-Collapse

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki, Tobias Fischer

Abstract: We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the stren… ▽ More We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Our results show that the MHD outflows are produced only for models (two out of four), to which magnetic field strengths of 10$^{12}$ G and rotation rates of 1 or 2 rad s$^{-1}$ are initially imposed in the core. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same order-of-magnitudes. Due to the memory effect, the frequency of the neutrino GW from our full-fledged 3D-MHD models is in the range less than $\sim$10 Hz. Toward the future GW detection for a Galactic core-collapse supernova, if driven by the MR mechanism, the planned next-generation detector as DECIGO is urgently needed to catch the low-frequency signals. △ Less

Submitted 17 November, 2023; v1 submitted 10 September, 2023; originally announced September 2023.

Comments: 12 pages, 10 figures

Neutrino-driven massive stellar explosions in 3D fostered by magnetic fields via turbulent $α$-effect

Authors: Jin Matsumoto, Tomoya Takiwaki, Kei Kotake

Abstract: We investigate the influence of magnetic field amplification on the core-collapse supernovae in highly magnetized progenitors through three-dimensional simulations. By considering rotating models, we observe a strong correlation between the exponential growth of the magnetic field in the gain region and the initiation of shock revival, with a faster onset compared to the non-rotating model. We hig… ▽ More We investigate the influence of magnetic field amplification on the core-collapse supernovae in highly magnetized progenitors through three-dimensional simulations. By considering rotating models, we observe a strong correlation between the exponential growth of the magnetic field in the gain region and the initiation of shock revival, with a faster onset compared to the non-rotating model. We highlight that the mean magnetic field experiences exponential amplification as a result of $α$-effect in the dynamo process, which works efficiently with the increasing kinetic helicity of the turbulence within the gain region. Our findings indicate that the significant amplification of the mean magnetic fields leads to the development of locally intense turbulent magnetic fields, particularly in the vicinity of the poles, thereby promoting the revival of the shock by neutrino heating. △ Less

Submitted 14 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

Multi-messenger signals of heavy axionlike particles in core-collapse supernovae: two-dimensional simulations

Authors: Kanji Mori, Tomoya Takiwaki, Kei Kotake, Shunsaku Horiuchi

Abstract: Core-collapse supernovae are a useful laboratory to probe the nature of exotic particles. If axionlike particles (ALPs) are produced in supernovae, they can affect the transfer of energy and leave traces in observational signatures. In this work, we present results from two-dimensional supernova models including the effects of the production and the absorption of ALPs that couple with photons. It… ▽ More Core-collapse supernovae are a useful laboratory to probe the nature of exotic particles. If axionlike particles (ALPs) are produced in supernovae, they can affect the transfer of energy and leave traces in observational signatures. In this work, we present results from two-dimensional supernova models including the effects of the production and the absorption of ALPs that couple with photons. It is found that the additional heating induced by ALPs can enhance the diagnostic energy of explosion, E_diag. For example, for moderate ALP-photon coupling, we find explosion energies ~0.6*10^51 erg compared to our reference model without ALPs of ~0.4*10^51 erg in the first ~0.5 s postbounce explored in this work. Our findings indicate that when the coupling constant is sufficiently high, the neutrino luminosities and mean energies are decreased because of the additional cooling of the proto-neutron star via ALPs. The gravitational wave amplitude is also reduced because the mass accretion on the proto-neutron star is suppressed. Although the ALP-photon coupling can foster explodability, including enhancing the explosion energy closer to recent observations, more long-term simulations in spatially three-dimension are needed to draw robust conclusions △ Less

Submitted 1 September, 2023; v1 submitted 22 April, 2023; originally announced April 2023.

Comments: 12 pages, 8 figures, accepted for publication in PRD

Monotonicity of the cores of massive stars

Authors: Koh Takahashi, Tomoya Takiwaki, Takashi Yoshida

Abstract: Massive stars are linked with diverse astronomical processes and objects including star formation, supernovae and their remnants, cosmic rays, interstellar media, and galaxy evolution. Understanding their properties is of primary importance for modern astronomy, and finding simple rules that characterize them is especially useful. However, theoretical simulations have not yet realized such relatio… ▽ More Massive stars are linked with diverse astronomical processes and objects including star formation, supernovae and their remnants, cosmic rays, interstellar media, and galaxy evolution. Understanding their properties is of primary importance for modern astronomy, and finding simple rules that characterize them is especially useful. However, theoretical simulations have not yet realized such relations, instead finding that the late evolutionary phases are significantly affected by a complicated interplay between nuclear reactions, chemical mixing, and neutrino radiation, leading to non-monotonic initial mass dependencies of the iron core mass and the compactness parameter. We conduct a set of stellar evolution simulations, in which evolutions of He star models are followed until their central densities uniformly reach 10$^{10}$ g cm$^{-3}$, and analyze their final structures as well as their evolutionary properties including the lifetime, surface radius change, and presumable fates after core collapse. Based on the homogeneous data set, we have found that monotonicity is inherent in the cores of massive stars. We show that not only the density, entropy, and chemical distributions, but also their lifetimes and explosion properties such as the proto-neutron-star mass and the explosion energy can be simultaneously ordered into a monotonic sequence. This monotonicity can be regarded as an empirical principle that characterizes the cores of massive stars. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 31 pages. Accepted for publication in ApJ

Long-term gravitational wave asteroseismology of supernova: from core collapse to 20 seconds postbounce

Authors: Masamitsu Mori, Yudai Suwa, Tomoya Takiwaki

Abstract: We use an asteroseismology method to calculate the frequencies of gravitational waves in a long-term core-collapse supernova simulation, with a mass of 9.6 $M_\odot$. The simulation, which includes neutrino transport in general relativity is performed from core-collapse, bounce, explosion and cooling of protoneutron stars (PNSs) up to 20 s after the bounce self-consistently. Based on the hydrodyna… ▽ More We use an asteroseismology method to calculate the frequencies of gravitational waves in a long-term core-collapse supernova simulation, with a mass of 9.6 $M_\odot$. The simulation, which includes neutrino transport in general relativity is performed from core-collapse, bounce, explosion and cooling of protoneutron stars (PNSs) up to 20 s after the bounce self-consistently. Based on the hydrodynamics background, we calculate eigenmodes of the PNS oscillation through a perturbation analysis on fluid and metric. We classify the modes by the number of nodes and find that there are several eigenmodes. In the early phase before 1 s, there are a low-frequency g-mode around 0.5 kHz, a mid-frequency f-modes around 1 kHz, and high-frequency p-modes above them. Beyond 1 s, the g-modes drop too low in frequency and the p-modes become too high to be detected by ground-based interferometers. However, the f-mode persists at 1 kHz. We present a novel fitting formula for the ramp-up mode, comprising a mixture of g-mode and f-mode, using postbounce time as a fitting parameter. Our approach yields improved results for the long-term simulation compared to prior quadratic formulas. We also fit frequencies using combinations of gravitational mass and radius of the PNS. We test three types of fitting variables: compactness, surface gravity, and average density. We present results of the time evolution of each mode and the fitting for three different ranges, from 0.2 to 1, 4, and 20 s for each formula. We compare the deviation of the formulas from the eigenmodes to determine which fitting formula is the best. In conclusion, any variable fits the eigenmodes well to a similar degree. Comparing 3 variables in detail, the fitting with compactness is slightly the best among them. We also find that the fitting using less than 1 s of simulation data cannot be extrapolated to the long-term frequency prediction. △ Less

Submitted 5 February, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: submitted to PRD

Light Curves and Event Rates of Axion Instability Supernovae

Authors: Kanji Mori, Takashi J. Moriya, Tomoya Takiwaki, Kei Kotake, Shunsaku Horiuchi, Sergei I. Blinnikov

Abstract: It was recently proposed that exotic particles can trigger a new stellar instability which is analogous to the e-e+ pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses… ▽ More It was recently proposed that exotic particles can trigger a new stellar instability which is analogous to the e-e+ pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ~400 keV--2 MeV and coupling with photons of g_{ag}~10^{-5} GeV^{-1}. It is found that the 56Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10--20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7--2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function. △ Less

Submitted 15 December, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures, 1 table, accepted for publication in ApJ

The quest for new correlations in the realm of the Gamma-Ray Burst -- Supernova connection

Authors: M. G. Dainotti, B. De Simone, M. I. Khadir, K. Kawaguchi, T. J. Moriya, T. Takiwaki, N. Tominaga, A. Gangopadhyay

Abstract: Gamma-Ray Bursts (GRBs) are very energetic cosmological transients. Long GRBs are usually associated with Type Ib/c Supernovae (SNe), and we refer to them as GRB-SNe. Since the associated SN for a given GRB is observed only at low redshift, a possible selection effect exists when we consider intrinsically faint sources which cannot be observed at high redshift. Thus, it is important to explore the… ▽ More Gamma-Ray Bursts (GRBs) are very energetic cosmological transients. Long GRBs are usually associated with Type Ib/c Supernovae (SNe), and we refer to them as GRB-SNe. Since the associated SN for a given GRB is observed only at low redshift, a possible selection effect exists when we consider intrinsically faint sources which cannot be observed at high redshift. Thus, it is important to explore the possible relationships between GRB and SN parameters after these have been corrected for astrophysical biases due to the instrumental selection effects and redshift evolution of the variables involved. So far, only GRB prompt emission properties have been checked against the SNe Ib/c properties without considering the afterglow (AG). This work investigates the existence of relationships among GRB's prompt and AG and associated SN properties. We investigate 91 bidimensional correlations among the SN and GRB observables before and after their correction for selection biases and evolutionary effects. As a result of this investigation, we find hints of a new correlation with a Pearson correlation coefficient > 0.50 and a probability of being drawn by chance < 0.05. This correlation is between the luminosity at the end of the GRB optical plateau emission and the rest-frame peak time of the SN. According to this relation, the brightest optical plateaus are accompanied by the largest peak times. This correlation is corrected for selection biases and redshift evolution and may provide new constraints for the astrophysical models associated with the GRB-SNe connection. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 40 pages, 6 figures, 7 tables; accepted for publication in ApJ

Journal ref: ApJ, 938, 41 (2022)

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

Authors: Tomoki Matsuoka, Shiu-Hang Lee, Keiichi Maeda, Tomoya Takiwaki, Takashi J. Moriya

Abstract: 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 const… ▽ More 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. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 23 pages, 10 figures, 3 tables, accepted for publication in The Astrophysical Journal

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

Authors: Jin Matsumoto, Yuta Asahina, Tomoya Takiwaki, Kei Kotake, Hiroyuki R. Takahashi

Abstract: 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 neu… ▽ More 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. △ Less

Submitted 18 August, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 16 pages, 21 figures, accepted for publication in MNRAS

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

Authors: Ko Nakamura, Tomoya Takiwaki, Kei Kotake

Abstract: 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, tota… ▽ More 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. △ Less

Submitted 7 June, 2022; v1 submitted 13 February, 2022; originally announced February 2022.

Comments: 14 pages, 14 figures, accepted for publication in MNRAS

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

Authors: Kanji Mori, Tomoya Takiwaki, Kei Kotake, Shunsaku Horiuchi

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 10 pages, 6 figures, accepted for publication in PRD

Modelling stellar convective transport with plumes: I. Non-equilibrium turbulence effect in double-averaging formulation

Authors: Nobumitsu Yokoi, Youhei Masada, Tomoya Takiwaki

Abstract: Plumes in a convective flow are considered to be relevant to the turbulent transport in convection. The effective mass, momentum, and heat transports in the convective turbulence are investigated in the framework of time--space double averaging procedure, where a field quantity is decomposed into three parts: the spatiotemporal mean (spatial average of the time-averaged) field, the dispersion or c… ▽ More Plumes in a convective flow are considered to be relevant to the turbulent transport in convection. The effective mass, momentum, and heat transports in the convective turbulence are investigated in the framework of time--space double averaging procedure, where a field quantity is decomposed into three parts: the spatiotemporal mean (spatial average of the time-averaged) field, the dispersion or coherent fluctuation, and the random or incoherent fluctuation. With this framework, turbulent correlations in the mean-field equations are divided into the dispersion/coherent and random/incoherent correlation part. By reckoning the plume as the coherent fluctuation, a transport model for the convective turbulence is constructed with the aid of the non-equilibrium effect, in which the change of turbulence characteristics along the mean stream is taken into account for the modelling of the turbulent transport coefficients. In this work, for the first time, change of turbulence properties along plume motions is incorporated into the expression of the turbulent transport coefficients. This non-equilibrium model is applied to a stellar convective flow. One of the prominent characteristics of a surface cooling-driven convection, the enhanced and localised turbulent mass flux below the surface layer, which cannot be reproduced at all by the usual eddy-diffusivity model with mixing length theory (MLT), is well reproduced by the present model. Our results show that the incorporation of plume motion into turbulent transport model is an important and very relevant extension of mean-field theory beyond the heuristic gradient transport model with MLT. △ Less

Submitted 27 April, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 18 pages, 8 figures

Universal relation for supernova gravitational waves

Authors: Hajime Sotani, Tomoya Takiwaki, Hajime Togashi

Abstract: Using the numerical simulation data for two-dimensional core-collapse supernova, we examine the protoneutron star (PNS) asteroseismology with the relativistic Cowling approximation. As shown in the previous study, the gravitational wave signals appearing in the numerical simulation can be well identified with the gravity (fundamental) oscillation in the early (later) phase before (after) the avoid… ▽ More Using the numerical simulation data for two-dimensional core-collapse supernova, we examine the protoneutron star (PNS) asteroseismology with the relativistic Cowling approximation. As shown in the previous study, the gravitational wave signals appearing in the numerical simulation can be well identified with the gravity (fundamental) oscillation in the early (later) phase before (after) the avoided crossing between the gravity and fundamental oscillations. On the other hand, the time evolution of supernova gravitational waves strongly depends on the PNS models, such as the progenitor mass and the equation of state for dense matter. Nevertheless, we find that the fundamental and gravity mode frequencies according to the gravitational wave signals appearing in the numerical simulations can be expressed as a function of the protoneutron star average density independently of the PNS models. Using the average density, we derive the empirical formula for supernova gravitational wave frequency. In addition, we confirm that the dependence of the PNS surface density on the PNS average density is almost independent of the PNS models and also discuss how the different treatment of the non-uniform matter in the equation of state affects the observables. △ Less

Submitted 25 November, 2021; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: accepted for publication in PRD; For the waveform data of gravitational waves, see https://doi.org/10.5281/zenodo.5555497

Report number: RIKEN-iTHEMS-Report-21

On the treatment of phenomenological turbulent effects in one dimensional simulations of core-collapse supernovae

Authors: Shunsuke Sasaki, Tomoya Takiwaki

Abstract: We have developed a phenomenological turbulent model with one-dimensional (1D) simulation based on Reynolds decomposition. Using this method, we have systematically studied models with different effects of compression, mixing length parameters, and diffusion coefficient of internal energy, turbulence energy and electron fraction. With employed turbulent effects, supernova explosion can be achieved… ▽ More We have developed a phenomenological turbulent model with one-dimensional (1D) simulation based on Reynolds decomposition. Using this method, we have systematically studied models with different effects of compression, mixing length parameters, and diffusion coefficient of internal energy, turbulence energy and electron fraction. With employed turbulent effects, supernova explosion can be achieved in 1D geometry, which can mimic the evolution of shock in the 3D simulations. We found that enhancement of turbulent energy by compression affects the early shock evolution. The diffusion coefficients of internal energy and turbulent energy also affect the explodability. The smaller diffusion makes the shock revival faster. Our comparison between the two reveals that the diffusion coefficients of internal energy has a greater impact. These simulations would help understand the role of turbulence in core-collapse supernovae. △ Less

Submitted 23 January, 2024; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: 13 pages, 18 figures, accepted in MNRAS

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

Authors: Hirokazu Sasaki, Tomoya Takiwaki

Abstract: 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 conversion… ▽ More 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. △ Less

Submitted 24 February, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 16 pages, 13 figures

Report number: LA-UR-21-29563

Journal ref: Prog Theor Exp Phys (2022)

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

Authors: Takami Kuroda, Tobias Fischer, Tomoya Takiwaki, Kei Kotake

Abstract: 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… ▽ More 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. △ Less

Submitted 20 October, 2021; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 24 pages, 16 figures, Accepted for publication in ApJ

Pre-supernova Ultra-light Axion-like Particles

Authors: Kanji Mori, Tomoya Takiwaki, Kei Kotake

Abstract: 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 oxyge… ▽ More 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. △ Less

Submitted 15 December, 2021; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 8 pages, 8 figures, 1 table, accepted for publication in PRD

Insights into non-axisymmetric instabilities in three-dimensional rotating supernova models with neutrino and gravitational-wave signatures

Authors: Tomoya Takiwaki, Kei Kotake, Thierry Foglizzo

Abstract: We present a detailed analysis to clarify what determines the growth of the low-$T/|W|$ instability in the context of rapidly rotating core-collapse of massive stars. To this end, we perform three-dimensional core-collapse supernova (CCSN) simulations of a $27 M_{\odot}$ star including several updates in the general relativistic correction to gravity, the multi-energy treatment of heavy-lepton neu… ▽ More We present a detailed analysis to clarify what determines the growth of the low-$T/|W|$ instability in the context of rapidly rotating core-collapse of massive stars. To this end, we perform three-dimensional core-collapse supernova (CCSN) simulations of a $27 M_{\odot}$ star including several updates in the general relativistic correction to gravity, the multi-energy treatment of heavy-lepton neutrinos, and the nuclear equation of state. Non-axisymmetric deformations are analyzed from the point of view of the time evolution of the pattern frequency and the corotation radius. The corotation radius is found to coincide with the convective layer in the proto neutron star (PNS). We propose a new mechanism to account for the growth of the low-$T/|W|$ instability in the CCSN environment. Near the convective boundary where a small Brunt-Väisälä frequency is expected, Rossby waves propagating in the azimuthal direction at mid latitude induce non-axisymmetric unstable modes, in both hemispheres. They merge with each other and finally become the spiral arm in the equatorial plane. We also investigate how the growth of the low-$T/|W|$ instability impacts the neutrino and gravitational-wave signatures. △ Less

Submitted 8 September, 2021; v1 submitted 6 July, 2021; originally announced July 2021.

Comments: 20 pages, 22 figures, submitted to MNRAS. For associated movie files, see https://doi.org/10.5281/zenodo.5463660. For the waveform data of gravitational waves, see https://doi.org/10.5281/zenodo.5489428

A three-dimensional hydrodynamics simulation of oxygen-shell burning in the final evolution of a fast-rotating massive star

Authors: Takashi Yoshida, Tomoya Takiwaki, David R. Aguilera-Dena, Kei Kotake, Koh Takahashi, Ko Nakamura, Hideyuki Umeda, Norbert Langer

Abstract: We perform for the first time a 3D hydrodynamics simulation of the evolution of the last minutes pre-collapse of the oxygen shell of a fast-rotating massive star. This star has an initial mass of 38 M$_\odot$, a metallicity of $\sim$1/50 Z$_\odot$, an initial rotational velocity of 600 km s$^{-1}$, and experiences chemically homogeneous evolution. It has a silicon- and oxygen-rich (Si/O) convectiv… ▽ More We perform for the first time a 3D hydrodynamics simulation of the evolution of the last minutes pre-collapse of the oxygen shell of a fast-rotating massive star. This star has an initial mass of 38 M$_\odot$, a metallicity of $\sim$1/50 Z$_\odot$, an initial rotational velocity of 600 km s$^{-1}$, and experiences chemically homogeneous evolution. It has a silicon- and oxygen-rich (Si/O) convective layer at (4.7-17)$\times 10^{8}$ cm, where oxygen-shell burning takes place. The power spectrum analysis of the turbulent velocity indicates the dominance of the large-scale mode ($\ell \sim 3$), which has also been seen in non-rotating stars that have a wide Si/O layer. Spiral arm structures of density and silicon-enriched material produced by oxygen-shell burning appear in the equatorial plane of the Si/O shell. Non-axisymmetric, large-scale ($m \le 3$) modes are dominant in these structures. The spiral arm structures have not been identified in previous non-rotating 3D pre-supernova models. Governed by such a convection pattern, the angle-averaged specific angular momentum becomes constant in the Si/O convective layer, which is not considered in spherically symmetrical stellar evolution models. Such spiral arms and constant specific angular momentum might affect the ensuing explosion or implosion of the star. △ Less

Submitted 30 June, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 6 pages, 6 figures, accepted for publication in MNRAS Letters

Neutrino-antineutrino oscillations induced by strong magnetic fields in dense matter

Authors: Hirokazu Sasaki, Tomoya Takiwaki

Abstract: We simulate neutrino-antineutrino oscillations caused by strong magnetic fields in dense matter. With the strong magnetic fields and large neutrino magnetic moments, Majorana neutrinos can reach flavor equilibrium. We find that the flavor equilibration of neutrino-antineutrino oscillations is sensitive to the values of the baryon density and the electron fraction inside the matter. The neutrino-an… ▽ More We simulate neutrino-antineutrino oscillations caused by strong magnetic fields in dense matter. With the strong magnetic fields and large neutrino magnetic moments, Majorana neutrinos can reach flavor equilibrium. We find that the flavor equilibration of neutrino-antineutrino oscillations is sensitive to the values of the baryon density and the electron fraction inside the matter. The neutrino-antineutrino oscillations are suppressed in the case of the large baryon density in neutron (proton)-rich matter. On the other hand, the flavor equilibration occurs when the electron fraction is close to $0.5$ even in the large baryon density. From the simulations, we propose a necessary condition for the equilibration of neutrino-antineutrino oscillations in dense matter. We also study whether such necessary condition is satisfied near the proto-neutron star by using results of neutrino hydrodynamic simulations of core-collapse supernovae. In our explosion model, the flavor equilibration would be possible if the magnetic field on the surface of the proto-neutron star is larger than $10^{14}$ G which is the typical value of the magnetic fields of magnetars. △ Less

Submitted 22 June, 2021; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 14 pages, 8 figures, accepted for publication in PRD

Journal ref: Phys. Rev. D 104, 023018 (2021)

Gravitational Wave Physics and Astronomy in the nascent era

Authors: Makoto Arimoto, Hideki Asada, Michael L. Cherry, Michiko S. Fujii, Yasushi Fukazawa, Akira Harada, Kazuhiro Hayama, Takashi Hosokawa, Kunihito Ioka, Yoichi Itoh, Nobuyuki Kanda, Koji S. Kawabata, Kyohei Kawaguchi, Nobuyuki Kawai, Tsutomu Kobayashi, Kazunori Kohri, Yusuke Koshio, Kei Kotake, Jun Kumamoto, Masahiro N. Machida, Hideo Matsufuru, Tatehiro Mihara, Masaki Mori, Tomoki Morokuma, Shinji Mukohyama , et al. (28 additional authors not shown)

Abstract: The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion,… ▽ More The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 89pages, Prog. Theor. Exp. Phys. accepted

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

Authors: Akihiro Suzuki, Matt Nicholl, Takashi J. Moriya, Tomoya Takiwaki

Abstract: 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… ▽ More 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. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Comments: 17 pages, 9 figures, accepted for publication in ApJ

Three-dimensional Hydrodynamics Simulations of Precollapse Shell Burning in the Si- and O-rich Layers

Authors: Takashi Yoshida, Tomoya Takiwaki, Kei Kotake, Koh Takahashi, Ko Nakamura, Hideyuki Umeda

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2021; v1 submitted 24 December, 2020; originally announced December 2020.

Comments: 23 pages, 26 figures, 5 tables. Accepted for publication in Astrophys. J

Impact of binary interactions on the diffuse supernova neutrino background

Authors: Shunsaku Horiuchi, Tomoya Kinugawa, Tomoya Takiwaki, Koh Takahashi, Kei Kotake

Abstract: 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 de… ▽ More 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. △ Less

Submitted 27 January, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 12 pages, 8 figures, matches published version

Journal ref: Phys. Rev. D 103, 043003 (2021)

Three-flavor collective neutrino conversions with multi-azimuthal-angle instability in an electron-capture supernova model

Authors: Masamichi Zaizen, Shunsaku Horiuchi, Tomoya Takiwaki, Kei Kotake, Takashi Yoshida, Hideyuki Umeda, John F. Cherry

Abstract: We investigate the multi-azimuthal angle (MAA) effect on collective neutrino oscillation by considering the three-dimensional neutrino momentum distribution in a realistic electron-capture supernova model with an $8.8 M_{\odot}$ progenitor. We find that the MAA effect induces collective flavor conversions at epochs when it is completely suppressed under the axial-symmetric approximation. This nove… ▽ More We investigate the multi-azimuthal angle (MAA) effect on collective neutrino oscillation by considering the three-dimensional neutrino momentum distribution in a realistic electron-capture supernova model with an $8.8 M_{\odot}$ progenitor. We find that the MAA effect induces collective flavor conversions at epochs when it is completely suppressed under the axial-symmetric approximation. This novel activity is switched on/off by the growth of the MAA instability and imprints additional time evolution in the expected neutrino event rate. We validate our results by extending the linear stability analysis into the three-flavor scheme including mixing angles, and confirm that the onset of collective neutrino oscillation matches the steep growth of flavor instability. We discuss how the MAA effect alters neutrino detection at Super-Kamiokande and DUNE. △ Less

Submitted 18 February, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 15 pages with 28 figures, Version accepted in Phys. Rev. D

Journal ref: Phys. Rev. D 103, 063008 (2021)

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

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: 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 o… ▽ More 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. △ Less

Submitted 13 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 20 pages, 19 figures

Accuracy of relativistic Cowling approximation in protoneutron star asteroseismology

Authors: Hajime Sotani, Tomoya Takiwaki

Abstract: The relativistic Cowling approximation, where the metric perturbations are neglected during the fluid oscillations, is often adopted for considering the gravitational waves from the protoneutron stars (PNSs) provided via core-collapse supernova explosions. In this study, we evaluate how the Cowling approximation works well by comparing the frequencies with the Cowling approximation to those withou… ▽ More The relativistic Cowling approximation, where the metric perturbations are neglected during the fluid oscillations, is often adopted for considering the gravitational waves from the protoneutron stars (PNSs) provided via core-collapse supernova explosions. In this study, we evaluate how the Cowling approximation works well by comparing the frequencies with the Cowling approximation to those without the approximation. Then, we find that the behavior of the frequencies with the approximation is qualitatively the same way as that without the approximation, where the frequencies with the approximation can totally be determined within $\sim 20\%$ accuracy. In particular, the fundamental mode with the Cowling approximation is overestimated. In addition, we also discuss the damping time of various eigenmodes in gravitational waves from the PNSs, where the damping time for the PNSs before the avoided crossing between the $f$- and $g_1$-modes, is quite different from that for cold neutron stars, but it is more or less similar to that for cold neutron stars in the later phase. The damping time is long enough compared to the typical time interval of short-Fourier transformation that often used in the analysis, and that ideally guarantees the validity of the transformation. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Report number: RIKEN-iTHEMS-Report-20

Two-dimensional numerical study for magnetic field dependence of neutrino-driven core-collapse supernova models

Authors: Jin Matsumoto, Tomoya Takiwaki, Kei Kotake, Yuta Asahina, Hiroyuki R. Takahashi

Abstract: We study the effects of the magnetic field on the dynamics of non-rotating stellar cores by performing two-dimensional (2D), magnetohydrodynamics (MHD) simulations. To this end, we have updated our neutrino-radiation-hydrodynamics supernova code to include MHD employing a divergence cleaning method with both careful treatments of finite volume and area reconstructions. By changing the initial stre… ▽ More We study the effects of the magnetic field on the dynamics of non-rotating stellar cores by performing two-dimensional (2D), magnetohydrodynamics (MHD) simulations. To this end, we have updated our neutrino-radiation-hydrodynamics supernova code to include MHD employing a divergence cleaning method with both careful treatments of finite volume and area reconstructions. By changing the initial strength of the magnetic field, the evolution of $15.0$, $18.4$ and $27.0$ $M_\odot$ presupernova progenitors is investigated. An intriguing finding in our study is that the neutrino-driven explosion occurs regardless of the strength of the initial magnetic field. For the 2D models presented in this work, the neutrino heating is the main driver for the explosion, whereas the magnetic field secondary contributes to the pre-explosion dynamics. Our results show that the strong magnetic field weakens the growth of the neutrino-driven turbulence in the small scale compared to the weak magnetic field. This results in the slower increase of the turbulent kinetic energy in the postshock region, leading to the slightly delayed onset of the shock revival for models with the stronger initial magnetic field. △ Less

Submitted 6 October, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: 23 pages, 10 figures, accepted for publication in MNRAS

Avoided crossing in gravitational wave spectra from protoneutron star

Authors: Hajime Sotani, Tomoya Takiwaki

Abstract: The ramp up signals of gravitational waves appearing in the numerical simulations could be important signals to estimate parameters of the protoneutron star (PNS) at supernova explosions. To identify the signals with PNS oscillations, we make a linear perturbation analysis and compare the resultant eigenfrequencies with the ramp up signals obtained via the two-dimensional numerical simulations. Th… ▽ More The ramp up signals of gravitational waves appearing in the numerical simulations could be important signals to estimate parameters of the protoneutron star (PNS) at supernova explosions. To identify the signals with PNS oscillations, we make a linear perturbation analysis and compare the resultant eigenfrequencies with the ramp up signals obtained via the two-dimensional numerical simulations. Then, we find that the ramp up signals correspond well to the $g_1$-mode in the early phase and to the $f$-mode, to which the $g_1$-mode is exchanged via the avoided crossing. We also confirm that the $f$- and $g_1$-modes are almost independent of the selection of the PNS surface density in the later phase after core bounce. In addition, we successfully find that the fitting formula of $g_1$- and $f$-modes, which correspond to the ramp up signals in the numerical simulation, as a function of the PNS average density. That is, via the direct observation of the gravitational waves after supernova explosion, one could extract the time evolution of the PNS average density by using our fitting formula. △ Less

Submitted 8 October, 2021; v1 submitted 2 August, 2020; originally announced August 2020.

Comments: For the waveform data of gravitational waves, see https://doi.org/10.5281/zenodo.5481811

Report number: RIKEN-iTHEMS-Report-20

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

Authors: Takashi J. Moriya, Akihiro Suzuki, Tomoya Takiwaki, Yen-Chen Pan, Sergei I. Blinnikov

Abstract: 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 af… ▽ More 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. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: 9 pages, 13 figures, 1 table, accepted by Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 497, Issue 2, pp.1619-1626 (2020)

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

Authors: Akihiro Suzuki, Takashi J. Moriya, Tomoya Takiwaki

Abstract: 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-consumi… ▽ More 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. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 24 pages, 11 figures, accepted for publication in ApJ

Dimension dependence of numerical simulations on gravitational waves from protoneutron stars

Authors: Hajime Sotani, Tomoya Takiwaki

Abstract: 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… ▽ More 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. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Report number: RIKEN-iTHEMS-Report-20

Journal ref: Phys. Rev. D 102, 023028 (2020)

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

Authors: Takami Kuroda, Almudena Arcones, Tomoya Takiwaki, Kei Kotake

Abstract: 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.… ▽ More 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. △ Less

Submitted 12 May, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: 19 pages, 14 figures, accepted for publication in the ApJ

Convection and Dynamo in Newly-born Neutron Stars

Authors: Youhei Masada, Tomoya Takiwaki, Kei Kotake

Abstract: 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-conve… ▽ More 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. △ Less

Submitted 3 November, 2021; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: 18 pages, 14 figures, accepted for publication in ApJ

Time of Flight and Supernova Progenitor Effects on the Neutrino Halo

Authors: John F. Cherry, George M. Fuller, Shunsaku Horiuchi, Kei Kotake, Tomoya Takiwaki, Tobias Fischer

Abstract: 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 po… ▽ More 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. △ Less

Submitted 25 May, 2020; v1 submitted 24 December, 2019; originally announced December 2019.

Comments: 7 pages, 3 figures and 1 table

Journal ref: Phys. Rev. D 102, 023022 (2020)

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

Authors: Akihiro Suzuki, Takashi J. Moriya, Tomoya Takiwaki

Abstract: 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 evid… ▽ More 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. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Comments: 35 pages, 22 figures, accepted for publication in ApJ

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

Authors: Sajad Abbar, Huaiyu Duan, Kohsuke Sumiyoshi, Tomoya Takiwaki, Maria Cristina Volpe

Abstract: 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 m… ▽ More 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. △ Less

Submitted 10 February, 2020; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: Version accepted in PRD

Journal ref: Phys. Rev. D 101, 043016 (2020)

A New Gravitational Wave Signature of Low-$T/|W|$ Instability in Rapidly Rotating Stellar Core Collapse

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: 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… ▽ More 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. △ Less

Submitted 6 February, 2020; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: 6 pages, 5 figures, accepted for publication in MNRAS Letters

Neutrino halo effect on collective neutrino oscillation in iron core-collapse supernova model of a 9.6 $M_{\odot}$ star

Authors: Masamichi Zaizen, John F. Cherry, Tomoya Takiwaki, Shunsaku Horiuchi, Kei Kotake, Hideyuki Umeda, Takashi Yoshida

Abstract: 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.… ▽ More 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. △ Less

Submitted 17 May, 2020; v1 submitted 28 August, 2019; originally announced August 2019.

Comments: 29 pages with 32 figures, Accepted for publication in JCAP

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

Authors: Hirokazu Sasaki, Tomoya Takiwaki, Shio Kawagoe, Shunsaku Horiuchi, Koji Ishidoshiro

Abstract: 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… ▽ More 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. △ Less

Submitted 5 March, 2020; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 18 pages, 17 figures, accepted for publication in PRD

Journal ref: Phys. Rev. D 101, 063027 (2020)

Dependence of outer boundary condition on protoneutron star asteroseismology with gravitational-wave signatures

Authors: Hajime Sotani, Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: To obtain the eigenfrequencies of a protoneutron star (PNS) in the postbounce phase of core-collapse supernovae (CCSNe), we perform a linear perturbation analysis of the angle-averaged PNS profiles using results from a general relativistic CCSN simulation of a $15 M_{\odot}$ star. In this work, we investigate how the choice of the outer boundary condition could affect the PNS oscillation modes in… ▽ More To obtain the eigenfrequencies of a protoneutron star (PNS) in the postbounce phase of core-collapse supernovae (CCSNe), we perform a linear perturbation analysis of the angle-averaged PNS profiles using results from a general relativistic CCSN simulation of a $15 M_{\odot}$ star. In this work, we investigate how the choice of the outer boundary condition could affect the PNS oscillation modes in the linear analysis. By changing the density at the outer boundary of the PNS surface in a parametric manner, we show that the eigenfrequencies strongly depend on the surface density. By comparing with the gravitational wave (GW) signatures obtained in the hydrodynamics simulation, the so-called surface $g$-mode of the PNS can be well ascribed to the fundamental oscillations of the PNS. The frequency of the fundamental oscillations can be fitted by a function of the mass and radius of the PNS similar to the case of cold neutron stars. In the case that the position of the outer boundary is chosen to cover not only the PNS but also the surrounding postshock region, we obtain the eigenfrequencies close to the modulation frequencies of the standing accretion-shock instability (SASI). However, we point out that these oscillation modes are unlikely to have the same physical origin of the SASI modes seen in the hydrodynamics simulation. We discuss possible limitations of applying the angle-averaged, linear perturbation analysis to extract the full ingredients of the CCSN GW signatures. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: accepted for publication in PRD

Journal ref: Phys. Rev. D 99, 123024 (2019)

Long-term Simulations of Multi-Dimensional Core-collapse Supernovae: Implications for Neutron Star Kicks

Authors: Ko Nakamura, Tomoya Takiwaki, Kei Kotake

Abstract: Core-collapse supernovae (CCSNe) are the final stage of massive stars, marking the birth of neutron stars (NSs). The aspherical mass ejection drives a natal kick of the forming NS. In this work we study the properties of the NS kick based on our long-term hydrodynamics CCSN simulations. We perform two-dimensional (2D) simulations for ten progenitors from a 10.8 to 20 $\, M_{\odot}$ star covering a… ▽ More Core-collapse supernovae (CCSNe) are the final stage of massive stars, marking the birth of neutron stars (NSs). The aspherical mass ejection drives a natal kick of the forming NS. In this work we study the properties of the NS kick based on our long-term hydrodynamics CCSN simulations. We perform two-dimensional (2D) simulations for ten progenitors from a 10.8 to 20 $\, M_{\odot}$ star covering a wide range of the progenitor's compactness parameter, and two three-dimensional (3D) simulations for an 11.2 $\, M_{\odot}$ star. Our 2D models present a variety of explosion energies between $\sim 1.3 \times 10^{50}$ erg and $\sim 1.2 \times 10^{51}$ erg, and NS kick velocities between $\sim 100$ km s$^{-1}$ and $\sim 1500$ km s$^{-1}$. For the 2D exploding models, we find that the kick velocities tend to become higher with the progenitor's compactness. This is because the high progenitor compactness results in high neutrino luminosity from the proto-neutron star (PNS), leading to more energetic explosions. Since high-compactness progenitors produce massive PNSs, we point out that the NS masses and the kick velocities can be correlated, which is moderately supported by observation. Comparing 2D and 3D models of the 11.2 $\, M_{\odot}$ star, the diagnostic explosion energy in 3D is, as previously identified, higher than that in 2D, whereas the 3D model results in a smaller asymmetry in the ejecta distribution and a smaller kick velocity than in 2D. Our results confirm the importance of self-consistent CCSN modeling covering a long-term postbounce evolution in 3D for a quantitative prediction of the NS kicks. △ Less

Submitted 5 August, 2019; v1 submitted 17 April, 2019; originally announced April 2019.

Comments: 13 pages, 10 figures, accepted for publication in PASJ

One-, Two-, and Three-dimensional Simulations of Oxygen Shell Burning Just Before the Core-Collapse of Massive Stars

Authors: Takashi Yoshida, Tomoya Takiwaki, Kei Kotake, Koh Takahashi, Ko Nakamura, Hideyuki Umeda

Abstract: We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution code, we first calculate the evolution of massive stars with an initial mass of 9-40 $M_\odot$. Four different overshoot parameters are applied, and CO core mass… ▽ More We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution code, we first calculate the evolution of massive stars with an initial mass of 9-40 $M_\odot$. Four different overshoot parameters are applied, and CO core mass trend similar to previous works is obtained in the 1D models. Selecting eleven 1D models that have a silicon and oxygen coexisting layer, we perform 2D hydrodynamics simulations of the evolution $\sim$100 s until the onset of core-collapse. We find that convection with large-scale eddies and the turbulent Mach number $\sim$0.1 is obtained in the models having a Si/O layer with a scale of 10$^8$ cm, whereas most models that have an extended O/Si layer up to a few $\times 10^9$ cm exhibit lower turbulent velocity. Our results indicate that the supernova progenitors that possess a thick Si/O layer could provide a preferable condition for perturbation-aided explosions. We perform 3D simulation of a 25 $M_\odot$ model, which exhibits large-scale convection in the 2D models. The 3D model develops large ($\ell = 2$) convection similar to the 2D model, however, the turbulent velocity is lower. By estimating the neutrino emission properties of the 3D model, we point out that a time modulation of the event rates, if observed in KamLAND and Hyper-Kamiokande, would provide an important information about structural changes in the presupernova convective layer. △ Less

Submitted 19 June, 2019; v1 submitted 18 March, 2019; originally announced March 2019.

Comments: 25 pages, 17 figures, 2 tables. Accepted for publication in ApJ

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

Authors: Sajad Abbar, Huaiyu Duan, Kohsuke Sumiyoshi, Tomoya Takiwaki, Maria Cristina Volpe

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2020; v1 submitted 17 December, 2018; originally announced December 2018.

Journal ref: Phys. Rev. D 100, 043004 (2019)