Source author record

Kei Kotake

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

ResearcherUnclaimed source record

astro-ph.HE astro-ph astro-ph.SR gr-qc hep-ph nucl-th astro-ph.CO physics.flu-dyn

Catalog footprint

What is connected

62works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2026arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

Convection and Dynamo in Newly-born Neutron Stars

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

preprint2021arXiv

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

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

preprint2021arXiv

Impact of binary interactions on the diffuse supernova neutrino background

Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%-20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. The enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.

preprint2021arXiv

Pre-supernova Ultra-light Axion-like Particles

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

preprint2021arXiv

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

We present 3D hydrodynamics simulations of shell burning in two progenitors with zero-age main-sequence masses of 22 and 27 $M_{\odot}$ for $\sim$65 and 200 s up to the onset of gravitational collapse, respectively. The 22 and 27 $M_{\odot}$ stars are selected from a suite of 1D progenitors. The former and the latter have an extended Si- and O-rich layer with a width of $\sim$10$^9$ cm and $\sim$5$\times 10^9$ cm, respectively. Our 3D results show that turbulent mixing occurs in both of the progenitors with the angle-averaged turbulent Mach number exceeding $\sim$0.1 at the maximum. We observe that an episodic burning of O and Ne, which takes place underneath the convection bases, enhances the turbulent mixing in the 22 and 27 $M_\odot$ models, respectively. The distribution of nucleosynthetic yields is significantly different from that in 1D simulations, namely, in 3D more homogeneous and inhomogeneous in the radial and angular direction, respectively. By performing a spectrum analysis, we investigate the growth of turbulence and its role of material mixing in the convective layers. We also present a scalar spherical harmonics mode analysis of the turbulent Mach number. This analytical formula would be helpful for supernova modelers to implement the precollapse perturbations in core-collapse supernova simulations. Based on the results, we discuss implications for the possible onset of the perturbation-aided neutrino-driven supernova explosion.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Detectability of SASI activity in supernova neutrino signals

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

preprint2020arXiv

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

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

preprint2020arXiv

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

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the $9.6\, \rm M_\odot$ progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the $9.6\, \rm M_\odot$ progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

preprint2020arXiv

Time of Flight and Supernova Progenitor Effects on the Neutrino Halo

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

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

The criterion of supernova explosion revisited: the mass accretion history

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

preprint2014arXiv

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

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

preprint2014arXiv

An Investigation into Surface Temperature Distributions of High-B Pulsars

Bearing in mind the application to high-magnetic-field (high-B) radio pulsars, we investigate two-dimensional (2D) thermal evolutions of neutron stars (NSs). We pay particular attention to the influence of different equilibrium configurations on the surface temperature distributions. The equilibrium configurations are constructed in a systematic manner, in which both toroidal and poloidal magnetic fields are determined self-consistently with the inclusion of general relativistic effects. To solve the 2D heat transfer inside the NS interior out to the crust, we have developed an implicit code based on a finite-difference scheme that deals with anisotropic thermal conductivity and relevant cooling processes in the context of a standard cooling scenario. In agreement with previous studies, the surface temperatures near the pole become higher than those in the vicinity of the equator as a result of anisotropic heat transfer. Our results show that the ratio of the highest to the lowest surface temperatures changes maximally by one order of magnitude, depending on the equilibrium configurations. Despite such difference, we find that the area of such hot and cold spots is so small that the simulated X-ray spectrum could be well reproduced by a single temperature blackbody fitting.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2013arXiv

Effects of Resistivity on Magnetized Core-Collapse Supernovae

We studied roles of a turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out 2D-resistive-MHD simulations. The three cases with different initial strengths of magnetic field and rotation are investigated; 1. strongly magnetized rotating core; 2.moderately magnetized rotating core; 3. very strongly magnetized non-rotating core. In each case, both an ideal-MHD model and resistive-MHD models are computed. As a result of computations, each model shows a matter eruption helped by a magnetic acceleration (and also by a centrifugal acceleration in the rotating cases). We found that a resistivity attenuates the explosion in case~1 and 2, while it enhances the explosion in case~3. We also found that in the rotating cases, main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of magnetic field and a winding of poloidal magnetic field-lines by differential rotation, which are somewhat dampened down with the presence of a resistivity. Although the magnetorotational instability seems to occur in the rotating models, it will play only a minor role in a magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matters, $> 2.5$, attained in a ideal-MHD model is reduced to some extent in a resistive model. These results indicate that a resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

preprint2012arXiv

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

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

preprint2012arXiv

Explosive Nucleosynthesis in Magnetohydrodynamical Jets from Collapsars II. Heavy-Element Nucleosynthesis of s, r, p-Processes

We investigate the nucleosynthesis in a massive star of 70 M_solar with solar metallicity in the main sequence stage. The helium core mass after hydrogen burning corresponds to 32 M_solar. Nucleosynthesis calculations have been performed during the stellar evolution and the jetlike supernova explosion of a collapsar model, where the weak s-, p-, and r-processes are taken into account. We confirm that s-elements of 60 < A < 90 are highly overproduced relative to the solar abundances in the hydrostatic nucleosynthesis. During oxygen burning, p-elements of A > 90 are produced via photodisintegrations of seed s-elements. However, the produced p-elements are disintegrated in later stages except for ^{180}Ta. In the explosive nucleosynthesis, elements of 90 < A < 160 are significantly overproduced relative to the solar values owing to the r-process. Only heavy p-elements (N > 50) are overproduced via the p-process. Compared with the previous study of r-process nucleosynthesis calculations in the collapsar model of 40 M_solar by Fujimoto et al. 2007, 2008, our jet model cannot contribute to the third peak of the solar r-elements and intermediate p-elements. Averaging the overproduction factors over the progenitor masses with the use of Salpeter's IMF, we suggest that the 70 M_solar star could contribute to the solar weak s-elements of 60 < A < 90 and neutron-rich elements of 90 < A < 160. We confirm the primary synthesis of light p-elements in the ejected matter of high peak temperature. The ejected matter has [Sr/Eu] \sim -0.4, which is different from that of a typical r-process-enriched star CS22892-052 ([Sr/Eu] \sim -1). We find that Sr-Y-Zr isotopes are primarily synthesized in the explosive nucleosynthesis in a similar process of the primary production of light p-elements, which has been considered as one of the sites of a lighter element primary process (LEPP).

preprint2012arXiv

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

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

preprint2012arXiv

Gravitational Wave Signatures of Hyperaccreting Collapsar Disks

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

preprint2012arXiv

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

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

preprint2012arXiv

Multiple physical elements to determine the gravitational-wave signatures of core-collapse supernovae

We review recent progress in the theoretical predictions of gravitational waves (GWs) of core-collapse supernovae. Following a brief summary of the methods in the numerical modeling, we summarize multiple physical elements that determine the GW signatures which have been considered to be important in extracting the information of the long-veiled explosion mechanism from the observation of the GWs. We conclude with a summary of the most urgent tasks to make the dream come true.

preprint2012arXiv

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

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

preprint2012arXiv

R-process Nucleosynthesis during the Magnetohydrodynamics Explosions of a Massive Star

We investigate the possibility of the r-process during the magnetohydrohynamical explosion of supernova in a massive star of 13 solar mass with the effects of neutrinos induced. We adopt five kinds of initial models which include properties of rotation and the toroidal component of the magnetic field . The simulations which succeed the explosions are limitted to a concentrated magnetic field and strong differential rotation. Low $Y_{e}$ ejecta produce heavy elements and the third peak can be reprocuced. However, the second peak is low because $Y_{e}$ distribution as a function of radius is steep and ejecta corresponding to middle $Y_{e}$ is very few.

preprint2012arXiv

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

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

preprint2011arXiv

Effects of Rotation on Stochasticity of Gravitational Waves in Nonlinear Phase of Core-Collapse Supernovae

By performing three-dimensional (3D) simulations that demonstrate the neutrino-driven core-collapse supernovae aided by the standing accretion shock instability (SASI), we study how the spiral modes of the SASI can have impacts on the properties of the gravitational-wave (GW) emission. To see the effects of rotation in the non-linear postbounce phase, we give a uniform rotation on the flow advecting from the outer boundary of the iron core, whose specific angular momentum is assumed to agree with recent stellar evolution models. We compute fifteen 3D models in which the initial angular momentum as well as the input neutrino luminosities from the protoneutron star are changed in a systematic manner. By performing a ray-tracing analysis, we accurately estimate the GW amplitudes generated by anisotropic neutrino emission. Our results show that the gravitational waveforms from neutrinos in models that include rotation exhibit a common feature otherwise they vary much more stochastically in the absence of rotation. The breaking of the stochasticity stems from the excess of the neutrino emission parallel to the spin axis. This is because the compression of matter is more enhanced in the vicinity of the equatorial plane due to the growth of the spiral SASI modes, leading to the formation of spiral flows circulating around the spin axis with higher temperatures. We point out that a recently proposed future space interferometers like Fabry-Perot type DECIGO would permit detection of these signals for a Galactic supernova.

preprint2011arXiv

Explosive nucleosynthesis in the neutrino-driven aspherical supernova explosion of a non-rotating 15$M_{\odot}$ star with solar metallicity

We investigate explosive nucleosynthesis in a non-rotating 15$M_\odot$ star with solar metallicity that explodes by a neutrino-heating supernova (SN) mechanism aided by both standing accretion shock instability (SASI) and convection. To trigger explosions in our two-dimensional hydrodynamic simulations, we approximate the neutrino transport with a simple light-bulb scheme and systematically change the neutrino fluxes emitted from the protoneutron star. By a post-processing calculation, we evaluate abundances and masses of the SN ejecta for nuclei with the mass number $\le 70$ employing a large nuclear reaction network. Aspherical abundance distributions, which are observed in nearby core-collapse SN remnants, are obtained for the non-rotating spherically-symmetric progenitor, due to the growth of low-mode SASI. Abundance pattern of the supernova ejecta is similar to that of the solar system for models whose masses ranges $(0.4-0.5) \Ms$ of the ejecta from the inner region ($\le 10,000\km$) of the precollapse core. For the models, the explosion energies and the \nuc{Ni}{56} masses are $ \simeq 10^{51} \rm erg$ and $(0.05-0.06) \Ms$, respectively; their estimated baryonic masses of the neutron star are comparable to the ones observed in neutron-star binaries. These findings may have little uncertainty because most of the ejecta is composed by matter that is heated via the shock wave and has relatively definite abundances. The abundance ratios for Ne, Mg, Si and Fe observed in Cygnus loop are well reproduced with the SN ejecta from an inner region of the $15\Ms$ progenitor.

preprint2011arXiv

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

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

preprint2011arXiv

Impacts of Collective Neutrino Oscillations on Supernova Explosions

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

preprint2010arXiv

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

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

preprint2010arXiv

Magnetic Field Decay Due to the Wave-Particle Resonances in the Outer Crust of the Neutron Star

Bearing in mind the application to the outer crust of the neutron stars (NSs), we investigate the magnetic field decay by means of the fully relativistic Particle-In-Cell simulations. Numerical computations are carried out in 2-dimensions, in which the initial magnetic fields are set to be composed both of the uniform magnetic fields that model the global fields penetrating the NS and of the turbulent magnetic fields that would be originated from the Hall cascade of the large-scale turbulence. Our results show that the whistler cascade of the turbulence transports the magnetic energy preferentially in the direction perpendicular to the uniform magnetic fields. It is also found that the distribution function of electrons becomes anisotropic because electrons with lower energies are predominantly heated in the direction parallel to the uniform magnetic fields due to the Landau resonance, while electrons with higher energies are heated mainly by the cyclotron resonance that makes the distribution function isotropic for the high energy tails. Furthermore we point out that the degree of anisotropy takes maximum as a function of the initial turbulent magnetic energy. As an alternative to the conventional ohmic dissipation, we propose that the magnetic fields in the outer crust of NSs, cascading down to the electron inertial scale via the whistler turbulence, would decay predominantly by the dissipation processes through the Landau damping and the cyclotron resonance.

preprint2010arXiv

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

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

preprint2009arXiv

General relativistic compact stars with exotic matter

We study structures of general relativistic compact stars with exotic matter. Our study is based on axisymetric and stationary formalism including purely toroidal magnetic field. We also study the finite size effects of quark-hadron mixed phase on structures of magnetars. For hybrid stars, we find a characteristic distribution of magnetic field, which has a discontinuity originated in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical phenomena, such as cooling processes.

preprint2009arXiv

Ray-Tracing Analysis of Anisotropic Neutrino Radiation for Estimating Gravitational Waves in Core-Collapse Supernovae

We propose a ray-tracing method to estimate gravitational waves (GWs) generated by anisotropic neutrino emission in supernova cores. To calculate the gravitational waveforms, we derive analytic formulae in a useful form, which are applicable also for three-dimensional computations. Pushed by evidence of slow rotation prior to core-collapse, we focus on asphericities in neutrino emission and matter motions outside the protoneutron star. Based on the two-dimensional (2D) models, which mimic SASI-aided neutrino heating explosions, we compute the neutrino anisotropies via the ray-tracing method in a post-processing manner and calculate the resulting waveforms. With these computations, it is found that the waveforms exhibit more variety in contrast to the ones previously estimated by the ray-by-ray analysis (e.g., Kotake et al. (2007)). In addition to a positively growing feature, which was predicted to determine the total wave amplitudes predominantly, the waveforms are shown to exhibit large negative growth for some epochs during the growth of SASI. These features are found to stem from the excess of neutrino emission in lateral directions, which can be precisely captured by the ray-tracing calculation. Due to the negative contributions and the neutrino absorptions appropriately taken into account by the ray-tracing method, the wave amplitudes become more than one-order-of magnitude smaller than the previous estimation, thus making their detections very hard for a galactic source.On the other hand, it is pointed out that the GW spectrum from matter motions have its peak near $\sim 100$ Hz, which could be characteristic for the SASI-induced supernova explosions.(abridged)

preprint2009arXiv

Relativistic hybrid stars with super-strong toroidal magnetic fields: An evolutionary track with QCD phase transition

We investigate structures of hybrid stars, which feature quark core surrounded by a hadronic matter mantle, with super-strong toroidal magnetic fields in full general relativity. Modeling the equation of state (EOS) with a first order transition by bridging the MIT bag model for the description of quark matter and the nuclear EOS by Shen et al., we numerically construct thousands of the equilibrium configurations to study the effects of the phase transition. It is found that the appearance of the quark phase can affect distributions of the magnetic fields inside the hybrid stars, making the maximum field strength about up to 30 % larger than for the normal neutron stars. Using the equilibrium configurations, we explore the possible evolutionary paths to the formation of hybrid stars due to the spin-down of magnetized rotating neutron stars. We find that the energy release by the phase transition to the hybrid stars is quite large ($\la 10^{52} \rm erg$) even for super strongly magnetized compact stars. Our results suggest that the strong gravitational-wave emission and the sudden spin-up signature could be observable signals of the QCD phase transition, possibly for a source out to Megaparsec distances.

preprint2009arXiv

Relativistic stars with purely toroidal magnetic fields with realistic equations of state

We investigate equilibrium sequences of relativistic stars containing purely toroidal magnetic fields with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). We numerically construct thousands of equilibrium configurations. Particularly we pay attention to the equilibrium sequences of constant baryon mass and/or constant magnetic flux, which model evolutions of an isolated neutron star. Important properties obtained in this study are summarized as follows ; (1) The dependence of the mass-shedding angular velocity on the EOSs is determined from that of the non-magnetized case. The stars with Shen(FPS) EOS reach the mass-shedding limit at the smallest(largest) angular velocity, while the stars with SLy or Lattimer-Swesty EOSs take the moderate values. (2) For the supramassive sequences, the equilibrium configurations are found to be generally oblate for the realistic EOSs in sharp contrast to the polytropic stars. For FPS(LS) EOS, the parameter region which permits the prolately deformed stars is widest(narrowest). For SLy and Shen EOS, it is in medium. Furthermore, the angular velocities $Ω_{\rm up}$, above which the stars start to spin up as they lose angular momentum, are found to depend sharply on the realistic EOSs. Our analysis indicates that the hierarchy of this spin up angular velocity is $Ω_{\rm up,SLy} > Ω_{\rm up,FPS} > Ω_{\rm up,LS}>Ω_{\rm up,Shen}$ and this relation holds even if the sequences have strong magnetic fields. Our results suggest the EOSs within the relativistic stars containing purely toroidal magnetic fields can be constrained by observing the angular velocity, the gravitational wave, and the signature of the spin up.

preprint2009arXiv

Stochastic Nature of Gravitational Waves from Supernova Explosions with Standing Accretion Shock Instability

We study properties of gravitational waves based on the three-dimensional simulations, which demonstrate the neutrino-driven explosions aided by the standing accretion shock instability (SASI). Pushed by evidence supporting slow rotation prior to core-collapse, we focus on the asphericities in neutrino emissions and matter motions outside the protoneutron star. By performing a ray-tracing calculation in 3D, we estimate accurately the gravitational waveforms from anisotropic neutrino emissions. In contrast to the previous work assuming axisymmetry, we find that the gravitational waveforms vary much more stochastically because the explosion anisotropies depend sensitively on the growth of the SASI which develops chaotically in all directions. Our results show that the gravitational-wave spectrum has its peak near $\sim 100$ Hz, reflecting the SASI-induced matter overturns of $\sim O(10)$ ms. We point out that the detection of such signals, possibly visible to the LIGO-class detectors for a Galactic supernova, could be an important probe into the long-veiled explosion mechanism.

preprint2007arXiv

Equilibrium Configurations of Strongly Magnetized Neutron Stars with Realistic Equations of State

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). Employing the Tomimura-Eriguchi scheme to construct the equilibrium configurations. we study the basic physical properties of the sequences in the framework of Newton gravity. In addition we newly take into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g., structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyze this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than $10^{15}\rm{g} \rm{cm}^{-3}$. The maximum baryon mass of the magnetized stars with axis ratio $q\sim 0.7$ increases about up to twenty percents for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly-born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.

preprint2007arXiv

Gravitational Wave Background from Population III Stars

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

preprint2007arXiv

Magnetorotational Collapse of Population III Stars

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

preprint2007arXiv

Three-Dimensional Simulations of Standing Accretion Shock Instability in Core-Collapse Supernovae

We have studied non-axisymmetric standing accretion shock instability, or SASI, by 3D hydrodynamical simulations. This is an extention of our previous study on axisymmetric SASI. We have prepared a spherically symmetric and steady accretion flow through a standing shock wave onto a proto-neutron star, taking into account a realistic equation of state and neutrino heating and cooling. This unperturbed model is supposed to represent approximately the typical post-bounce phase of core-collapse supernovae. We then have added a small perturbation (~1%) to the radial velocity and computed the ensuing evolutions. Not only axisymmetric but non-axisymmetric perturbations have been also imposed. We have applied mode analysis to the non-spherical deformation of the shock surface, using the spherical harmonics. We have found that (1) the growth rates of SASI are degenerate with respect to the azimuthal index m of the spherical harmonics Y_l^m, just as expected for a spherically symmetric background, (2) nonlinear mode couplings produce only m=0 modes for the axisymmetric perturbations, whereas m=!0 modes are also generated in the non-axisymmetric cases according to the selection rule for the quadratic couplings, (3) the nonlinear saturation level of each mode is lower in general for 3D than for 2D because a larger number of modes are contributing to turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the equi-partition is nearly established among different m modes in the nonlinear phase, (6) the spectra with respect to l obey power laws with a slope slightly steeper for 3D, and (7) although these features are common to the models with and without a shock revival at the end of simulation, the dominance of low l modes is more remarkable in the models with a shock revival.

preprint2006arXiv

Gravitational Radiation from Standing Accretion Shock Instability in Core-Collapse Supernovae

We present the results of numerical experiments, in which we study how the asphericities induced by the growth of the standing accretion shock instability (SASI) produce the gravitational waveforms in the postbounce phase of core-collapse supernovae. To obtain the neutrino-driven explosions, we parameterize the neutrino fluxes emitted from the central protoneutron star and approximate the neutrino transfer by a light-bulb scheme. We find that the waveforms due to the anisotropic neutrino emissions show the monotonic increase with time, whose amplitudes are up to two order-of-magnitudes larger than the ones from the convective matter motions outside the protoneutron stars. We point out that the amplitudes begin to become larger when the growth of the SASI enters the nonlinear phase, in which the deformation of the shocks and the neutrino anisotropy become large. From the spectrum analysis of the waveforms, we find that the amplitudes from the neutrinos are dominant over the ones from the matter motions at the frequency below $\sim 100$ Hz, which are suggested to be within the detection limits of the detectors in the next generation such as LCGT and the advanced LIGO for a supernova at 10 kpc. As a contribution to the gravitational wave background, we show that the amplitudes from this source could be larger at the frequency above $\sim$ 1 Hz than the primordial gravitational wave backgrounds, but unfortunately, invisible to the proposed space-based detectors.

preprint2006arXiv

Inelastic Neutrino-Helium Scatterings and Standing Accretion Shock Instability in Core-Collapse Supernovae

We present the results of numerical experiments, in which we have investigated the influence of the inelastic neutrino-helium interactions on the standing accretion shock instability supposed to occur in the post-bounce supernova core. The axisymmetric hydrodynamical simulations of accretion flows through the standing accretion shock wave onto the protoneutron star show that the interactions are relatively minor and the linear growth of the shock instability is hardly affected. The extra heating given by the inelastic reactions becomes important for the shock revival after the instability enters the non-linear regime, but only when the neutrino luminosity is very close to the critical value, at which the shock would be revived without the interactions. We have also studied the dependence of the results on the initial amplitudes of perturbation and the temperatures of mu and tau neutrinos.

preprint2005arXiv

Numerical Analysis on Standing Accretion Shock Instability with Neutrino Heating in the Supernova Cores

We have numerically studied the instability of the spherically symmetric standing accretion shock wave against non-spherical perturbations. We have in mind the application to the collapse-driven supernovae in the post bounce phase, where the prompt shock wave generated by core bounce is commonly stalled. We take an experimental stand point in this paper. Using spherically symmetric, completely steady, shocked accretion flows as unperturbed states, we have clearly observed both the linear growth and the subsequent nonlinear saturation of the instability. In so doing, we have employed a realistic equation of state together with heating and cooling via neutrino reactions with nucleons. We have done a mode analysis based on the spherical harmonics decomposition and found that the modes with l=1, 2 are dominant not only in the linear regime, but also after the nonlinear couplings generate various modes and the saturation occurs. Varying the neutrino luminosity, we have constructed the unperturbed states both with and without a negative entropy-gradient. We have found that in both cases the growth of the instability is similar, suggesting the convection does not play a dominant role, which also appears to be supported by the recent linear analysis of the convection in accretion flows by Foglizzo et al. The real part of the eigen frequency seems to be mainly determined by the advection time rather than by the sound-crossing time. Whatever the cause may be, the instability is favorable for the shock revival.

preprint2004arXiv

Gravitational Radiation from Rotational Core Collapse: Effects of Magnetic Fields and Realistic Equation of States

We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic simulations of the rotational collapse of a supernova core. In order to calculate the waveforms of the gravitational wave, we derive the quadrupole formula including the contributions from the electromagnetic fields. Recent stellar evolution calculations imply that the magnetic fields of the toroidal components are much stronger than those of the poloidal ones at the presupernova stage. Thus, we systematically investigate the effects of the toroidal magnetic fields on the amplitudes and waveforms. Furthermore, we employ the two kinds of the realistic equation of states, which are often used in the supernova simulations. Then, we investigate the effects of the equation of states on the gravitational wave signals. With these computations, we find that the peak amplitudes are lowered by an order of 10% for the models with the strongest toroidal magnetic fields. However, the peak amplitudes are mostly within sensitivity range of laser interferometers such as TAMA and the first LIGO for a source at a distance of 10 kpc. Furthermore, we point out that the amplitudes of second peaks are still within the detection limit of the first LIGO for the source, although the characteristics of second peaks are reduced by the magnetic fields. We stress the importance of the detection, since it will give us information about the angular momentum distribution of massive evolved stars. When we compare the gravitational waves from the two realistic equation of states, significant differences are not found, except that the typical frequencies of the gravitational wave become slightly higher for the softer equation of state.

preprint2004arXiv

North-South Neutrino Heating Asymmetry in Strongly Magnetized and Rotating Stellar Cores

We perform a series of two-dimensional magnetohydrodynamic simulations of supernova cores. Since the distributions of the angular momentum and the magnetic fields of strongly magnetized stars are quite uncertain, we systematically change the combinations of the strength of the angular momentum, the rotations law, the degree of differential rotation, and the profiles of the magnetic fields to construct the initial conditions. By so doing, we estimate how the rotation-induced anisotropic neutrino heating are affected by the strong magnetic fields through parity-violating effects and first investigate how the north-south asymmetry of the neutrino heating in a strongly magnetized supernova core could be. As for the microphysics, we employ a realistic equation of state based on the relativistic mean field theory and take into account electron captures and the neutrino transport via the neutrino leakage scheme. With these computations, we find that the parity-violating corrections reduce $ \lesssim 0.5 %$ of the neutrino heating rate than that without the magnetic fields in the vicinity of the north pole of a star, on the other hand, enhance about $ \lesssim 0.5 %$ in the vicinity of the south pole. If the global asymmetry of the neutrino heating in the both of the poles develops in the later phases, the newly born neutron star might be kicked toward the north pole in the subsequent time.

preprint2003arXiv

Gravitational Radiation from Axisymmetric Rotational Core Collapse

We have done a series of two-dimensional hydrodynamic simulations of the rotational collapse of a supernova core in axisymmetry. We have employed a realistic equation of state (EOS) and taken into account electron captures and neutrino transport by the so-called leakage scheme. It is an important progress to apply the realistic EOS coupled with the microphysics to 2-D simulations for computing gravitational radiation in rotational core collapse. We have used the quadrupole formula to calculate the amplitudes and the waveforms of gravitational wave assuming the Newtonian gravity. From these computations, we have extended the conventional category of the gravitational waveforms. Our results have shown that the peak amplitudes of gravitational wave are mostly within the sensitivity range of the laser interferometers such as TAMA and first LIGO for a source at a distance of 10 kpc. Furthermore we have found that the amplitudes of the second peaks are within the detection limit of first LIGO for the source and first pointed out the importance of the detections, since they will give us the information as to the angular momentum distribution of evolved massive stars.

Kei Kotake

What is connected

Connect this record

See the researcher in context

Building this map preview

62 published item(s)

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

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

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

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

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

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

Convection and Dynamo in Newly-born Neutron Stars

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

Impact of binary interactions on the diffuse supernova neutrino background

Pre-supernova Ultra-light Axion-like Particles

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

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

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

Detectability of SASI activity in supernova neutrino signals

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

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

Time of Flight and Supernova Progenitor Effects on the Neutrino Halo

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

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

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

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

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

The criterion of supernova explosion revisited: the mass accretion history

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

An Investigation into Surface Temperature Distributions of High-B Pulsars

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

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

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

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

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

Effects of Resistivity on Magnetized Core-Collapse Supernovae

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

Explosive Nucleosynthesis in Magnetohydrodynamical Jets from Collapsars II. Heavy-Element Nucleosynthesis of s, r, p-Processes

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

Gravitational Wave Signatures of Hyperaccreting Collapsar Disks

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

Multiple physical elements to determine the gravitational-wave signatures of core-collapse supernovae

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

R-process Nucleosynthesis during the Magnetohydrodynamics Explosions of a Massive Star

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

Effects of Rotation on Stochasticity of Gravitational Waves in Nonlinear Phase of Core-Collapse Supernovae

Explosive nucleosynthesis in the neutrino-driven aspherical supernova explosion of a non-rotating 15$M_{\odot}$ star with solar metallicity

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

Impacts of Collective Neutrino Oscillations on Supernova Explosions

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

Magnetic Field Decay Due to the Wave-Particle Resonances in the Outer Crust of the Neutron Star

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

General relativistic compact stars with exotic matter

Ray-Tracing Analysis of Anisotropic Neutrino Radiation for Estimating Gravitational Waves in Core-Collapse Supernovae

Relativistic hybrid stars with super-strong toroidal magnetic fields: An evolutionary track with QCD phase transition

Relativistic stars with purely toroidal magnetic fields with realistic equations of state

Stochastic Nature of Gravitational Waves from Supernova Explosions with Standing Accretion Shock Instability

Equilibrium Configurations of Strongly Magnetized Neutron Stars with Realistic Equations of State

Gravitational Wave Background from Population III Stars

Magnetorotational Collapse of Population III Stars

Three-Dimensional Simulations of Standing Accretion Shock Instability in Core-Collapse Supernovae

Gravitational Radiation from Standing Accretion Shock Instability in Core-Collapse Supernovae

Inelastic Neutrino-Helium Scatterings and Standing Accretion Shock Instability in Core-Collapse Supernovae

Numerical Analysis on Standing Accretion Shock Instability with Neutrino Heating in the Supernova Cores

Gravitational Radiation from Rotational Core Collapse: Effects of Magnetic Fields and Realistic Equation of States

North-South Neutrino Heating Asymmetry in Strongly Magnetized and Rotating Stellar Cores

Gravitational Radiation from Axisymmetric Rotational Core Collapse