Search

We perform multi-dimensional core-collapse supernova (CCSN) simulations in a massive scalar-tensor theory for the first time with a realistic equation of state and multi-energy neutrino radiation. Among the set of our models varying the scalar mass and the coupling strength between the scalar and gravitational fields, a particular model allows for recurrent spontaneous scalarizations (SSs) in the proto-neutron star (PNS). Each SS induces the PNS collapse and subsequent bounce, from which devastating shock waves emanate and eject the PNS envelope. The explosion energy can easily exceed $\mathcal O(10^{51})$ erg. This study reveals new aspects of SS as the explosion mechanism of CCSNe. We also discuss its characteristic multi-messenger signals: neutrinos and gravitational waves., 11 pages, 4 figures, accepted for publication in PRD

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., 6 pages, 5 figures, accepted for publication in MNRAS Letters

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\gamma} \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., Comment: 26 pages, 17 figures; v2 changes: inclusion of gravitational redshift on CC SN ALP emission spectra, impact of the creation of alpha particles during a CC SN on ALP emission spectra

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., 24 pages, 16 figures, Accepted for publication in ApJ

We present results of full general relativistic (GR), three-dimensional (3D) core-collapse simulation of a massive star with multi-energy neutrino transport. Using a 70Mȯ zero-metallicity star, we show that the black-hole (BH) formation occurs at ∼ 300 ms after bounce. At a few ∼ 10 ms before the BH formation, we find that the stalled bounce shock is revived by neutrino heating from the forming hot proto-neutron star (PNS), which is aided by vigorous convection behind the shock. Our numerical results present the first evidence to validate the BH formation by the so-called fallback scenario. Furthermore we present results from a rapidly rotating core-collapse model of a 27Mȯ star that is trending towards an explosion. We point out that the correlated neutrino and gravitational-wave signatures, if detected, could provide a smoking-gun evidence of rapid rotation of the newly-born PNS.

It was recently pointed out that very energetic subclasses of supernovae (SNe), like hypernovae and superluminous SNe, might host ultra-strong magnetic fields in their core. Such fields may catalyze the production of feebly interacting particles, changing the predicted emission rates. Here we consider the case of axion-like particles (ALPs) and show that the predicted large scale magnetic fields in the core contribute significantly to the ALP production, via a coherent conversion of thermal photons. Using recent state-of-the-art SN simulations including magnetohydrodynamics, we find that if ALPs have masses $m_a \sim {\mathcal O}(10)\, \rm MeV$, their emissivity via magnetic conversions is over two orders of magnitude larger than previously estimated. Moreover, the radiative decay of these massive ALPs would lead to a peculiar delay in the arrival times of the daughter photons. Therefore, high-statistics gamma-ray satellites can potentially discover MeV ALPs in an unprobed region of the parameter space and shed light on the magnetohydrodinamical nature of the SN explosion., 6 pages, 3 Figures

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., Comment: 20 pages, 19 figures

Through numerical simulations, it is predicted that the gravitational waves (GWs) reflect the characteristics of the core-collapse supernova (CCSN) explosion mechanism. There are multiple GW excitation processes that occur inside a star before its explosion, and it is suggested that the GWs originating from the CCSN have a mode for each excitation process in terms of time-frequency representation. Therefore, we propose an application of the Hilbert-Huang Transform (HHT), which is a high-resolution time-frequency analysis method, to analyze these GW modes for theoretically probing and increasing our understanding of the explosion mechanism. The HHT defines frequency as a function of time, and is not bound by the trade-off between time and frequency resolutions. In this study, we analyze a gravitational waveform obtained from a three-dimensional general-relativistic CCSN model that showed a vigorous activity of the standing-accretion-shock-instability (SASI). We succeed in extracting the SASI induced GWs with high resolution on a time-frequency representation using the HHT and we examine their instantaneous frequencies., Comment: 12 pages, 9 figures

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., 19 pages, 14 figures, accepted for publication in the ApJ

In this contribution, we present results from fully general-relativistic three-dimensional (3D) simulations of a non-rotating 15M⊙ star using different nuclear equations of state (EOSs). We show that the SASI (standing-accretion-shock-instability) activity occurs much more vigorously in models with softer EOS. By performing detailed analysis of the gravitational-wave (GW) emission, we find a new GW signature that is produced predominantly by the SASI-induced downflows to the proto-neutron star. We discuss the detectability of the GW signal by performing a coherent network analysis where multiple detectors including LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA are considered. We point out that the GW signal, whose typical frequency is in the best sensitivity range of the laser-interferometers, could potentially provide the live broadcast that pictures how the supernova shock is dancing in the core. The detection horizon of the signal is estimated as 2~3 kpc for the current generation detectors, which can extend up to ~100 kpc for the third generation detectors like Cosmic Explorer. We furthermore perform a correlation analysis between the SASI-modulated GW and neutrino signals. Our results show that the time correlation of the two signals becomes highest when we take into account the travel timescale of adverting material from the (average) neutrino-sphere to the proto-neutron star surface.

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

We present an analysis of the circular polarization of gravitational-waves (GWs) using results from three-dimensional (3D), general relativistic (GR) core-collapse simulations of a non-rotating $15 M_{\odot}$ star. For the signal detection, we perform a coherent network analysis taking into account the four interferometers of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA. We focus on the Stokes $V$ parameter, which directly characterizes the asymmetry of the GW circular polarization. We find that the amplitude of the GW polarization becomes bigger for our 3D-GR model that exhibits strong activity of the standing accretion shock instability (SASI). Our results suggest that the SASI-induced accretion flows to the proto-neutron star (PNS) lead to a characteristic, low-frequency modulation (100 $\sim$ 200 Hz) in both the waveform and the GW circular polarization. During the vigorous SASI phase, we observe that the GW polarization switches from the right- to left-handed mode, which is clearly visible in the spectrogram. By estimating the signal-to-noise ratio of the GW polarization, we demonstrate that the detection horizon of the circular polarization extends by more than a factor of several times farther comparing to that of the GW amplitude. Our results suggest that the GW circular polarization, if detected, could provide a new probe into the pre-explosion hydrodynamics such as the SASI activity and the $g$-mode oscillation of the PNS., 5 pages, 4 figures

We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star., Comment: 5 pages, 4 figures, accepted for publication in MNRAS letter

Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios., 20 pages, 11 figures, 4 Tables. Published in A&A

We present results from fully relativistic three-dimensional core-collapse supernova (CCSN) simulations of a non-rotating 15M⊙ star using three different nuclear equations of state (EoSs). From our simulations, we show that the development of the standing accretion shock instability (SASI) differs significantly depending on the stiffness of nuclear EoS. By evaluating the gravitational-wave (GW) emission, we find a new quasi-periodic GW signature on top of the previously identified high frequency (∼500–1000 Hz) one, which is originated from the g-mode oscillation of the photo-neutron star (PNS) surface. The newly found signal appears in relatively low frequency range from ∼100 to 200 Hz. 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 newly found GW emission frequency. This is because the SASI-induced temporary perturbed mass accretion strikes the PNS surface leading to the quasi-periodic GW emission.

We perform detailed nucleosynthesis calculations for two long-term, 2D simulations of core-collapse supernovae. We find that elements are produced up to Ru (Z = 44) and observe abundance patterns that are characteristic of a νp-process. One important characteristic of the long-term simulation is that there is still accretion of matter onto the proto-neutron star and unbinding of matter in some other regions at the time when the simulations stop (around 7s). Dividing the tracer particles into different bins according to their peak temperatures enables us to study and compare the nuclear compositions of these bins for the different simulations.

The gravitational-wave (GW) asteroseismology is a powerful technique for extracting interior information of compact objects. In this work, we focus on spacetime modes, the so-called $w$-modes, of GWs emitted from a proto-neutron star (PNS) in the postbounce phase of core-collapse supernovae. Using results from recent three-dimensional supernova models, we study how to infer the properties of the PNS based on a quasi-normal mode analysis in the context of the GW asteroseismology. We find that the $w_1$-mode frequency multiplied by the PNS radius is expressed as a linear function with respect to the ratio of the PNS mass to the PNS radius. This relation is insensitive to the nuclear equation of state (EOS) employed in this work. Combining with another universal relation of the $f$-mode oscillations, we point out that the time dependent mass-radius relation of the PNS can be obtained by observing both the $f$- and $w_1$-mode GWs simultaneously. Our results suggest that the simultaneous detection of the two modes could provide a new probe into finite-temperature nuclear EOS that predominantly determines the PNS evolution., Comment: accepted for publication in PRD

We present our latest results of two-dimensional core-collapse supernova simulations for about 400 progenitors. Our self-consistent supernova models reveal the systematic features of core-collapse supernova properties such as neutrino luminosity and energy spectrum, explosion energy, remnant mass, and yield of radioactive 56Ni. We find that these explosion characteristics tend to show a monotonic increase as a function of mass accretion rate onto a shock. The accretion rate depends on the structure of the progenitor core and its envelope, which is well described by the compactness parameter.

Multi-dimensionality in the inner working of core-collapse supernovae has long been considered one of the most important ingredients to understand the explosion mechanism. We perform a series of numerical experiments to explore how rotation impacts the 3-dimensional hydrodynamics of core-collapse supernova. We employ a light-bulb scheme to trigger explosions and a three-species neutrino leakage scheme to treat deleptonization eects and neutrino losses from the neutron star interior. We nd that the rotation can help the onset of neutrino-driven explosions for models in which the initial angular momentum is matched to that obtained from recent stellar evolutionary calculations ( 0:3 - 3 rad s 1 at the center). For models with larger initial angular momenta, a shock surface deforms to be oblate due to larger centrifugal force. This makes a gain region, in which matter gains energy from neutrinos, more concentrated around the equatorial plane. As a result, the preferred direction of the explosion in 3-dimensional rotating models is perpendicular to the spin axis, which is in sharp contrast to the polar explosions around the axis that are often obtained from 2-dimensional simulations.

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

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., Comment: 7 pages, 5 figures, Accepted for publication in ApJL

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., 61 pages, 20 figures, accepted for publication in ApJS

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

We present results from general-relativistic (GR) three-dimensional (3D) core-collapse simulations with approximate neutrino transport for three non-rotating progenitors (11.2, 15, and 40 Msun) using different nuclear equations of state (EOSs). We find that the combination of progenitor's higher compactness at bounce and the use of softer EOS leads to stronger activity of the standing accretion shock instability (SASI). We confirm previous predications that the SASI produces characteristic time modulations both in neutrino and gravitational-wave (GW) signals. By performing a correlation analysis of the SASI-modulated neutrino and GW signals, we find that the correlation becomes highest when we take into account the time-delay effect due to the advection of material from the neutrino sphere to the proto-neutron star core surface. Our results suggest that the correlation of the neutrino and GW signals, if detected, would provide a new signature of the vigorous SASI activity in the supernova core, which can be hardly seen if neutrino-convection dominates over the SASI., 24 pages, 10 figures, Accepted for publication in ApJ

Recent hydrodynamic studies of core-collapse supernovae imply that the neutrino-heated ejecta from a nascent neutron star develop to supersonic outflows. These supersonic winds are influenced by the reverse shock from the preceding supernova ejecta, forming the wind termination shock. We investigate the effects of the termination shock in neutrino-driven winds and its role in the r-process. Supersonic outflows are calculated with a semianalytic neutrino-driven wind model. Subsequent termination-shocked, subsonic outflows are obtained by applying the Rankine-Hugoniot relations. We find a couple of effects that can be relevant for the r-process. First is the sudden slowdown of the temperature decrease by the wind termination. Second is the entropy jump by termination-shock heating, up to several hundred NAk. Calculations of nucleosynthesis in the obtained winds are performed to examine these effects on the r-process. We find that the slowdown of the temperature decrease plays a decisive role in determining the r-process abundance curves. This is due to the strong dependences of the nucleosynthetic path on the temperature during the r-process freezeout phase. Our results suggest that only the termination-shocked winds with relatively small shock radii (~500 km) are relevant for the bulk of the solar r-process abundances (A ≈ 100-180). The heaviest part of the solar r-process curve (A ≈ 180-200), however, can be reproduced in both shocked and unshocked winds. These results may help to constrain the mass range of supernova progenitors relevant for the r-process. We find, on the other hand, a negligible role of the entropy jump in the r-process. This is because the sizable entropy increase takes place only at a large shock radius (10,000 km), where the r-process has already ceased.

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., Comment: PRD in press

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., Comment: 15pages, 13 figures, 1 table, revised and submitted to Astrophysical Journal after referee report

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 $\xi$ that characterizes the structure of the progenitors is also a key in 2D to diagnose the properties of neutrino-driven explosions. Models with high $\xi$ 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 $\xi$, 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 $\xi$., Comment: 15 pages, 22 figures, 277 progenitors added, accepted to PASJ

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

Tomoyuki Maruyama1, Jun Hidaka2, Toshitaka Kajino2, Nobutoshi Yasutake3, Takami Kuroda2, Myung-Ki Cheoun4, Chung-Yeol Ryu5 and Grant J. Mathews6 1College of Bioresource Sciences, Nihon University, Fujisawa 252-8510, Japan 2National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 3Department of Physics, Chiba Institute of Technology, Narashino 275-0023, Japan 4Department of Physics, Soongsil University, Seoul, 156-743, Korea 5General Education Curriculum Center, Hanyang University, Seoul, 133-791, Korea 6Center of Astrophysics, University of Notre Dame, Notre Dame, IN 46556, USA

We calculate the spin deceleration of proto-neutron stars by asymmetric neutrino absorption in the context of a fully relativistic mean field theory. We calculate for the first time the spin deceleration of neutron stars due to asymmetric neutrino absorption in a toroidal magnetic field configurations. We find a surprising effect that the deceleration can much larger for asymmetric neutrino absorption in a toroidal magnetic field than the usually presumed braking due to magnetic dipole radiation. This may explain the observation that magnetars appear to have had a more rapid deceleration in the past., Comment: 13 pages, 1 figure, 1 table

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., Comment: 47 pages, 20 figures, Accepted by PRD

We investigate the effects of rotation on the evolution of core-collapse supernova explosion using a 15 solar mass progenitor model with a variety of neutrino luminosity and rotational velocity. Stars should have some amount of angular momentum, which would affect stellar evolution and its final explosion. In this paper we focus on the effect of rotation on gravitational collapse of a core, on a core bounce of accreting matter, and on subsequent generation and evolution of a shock wave. We find that the rotation plays a positive role for supernova explosions. More rapidly rotating models present more rapid expansion of the shock front and more energetic explosions. When the rotational speed is moderate, the shock once stalls at about 200 km away from the center similarly to a non-rotating model. Then the rotating progenitor experiences effective neutrino heating especially around an equatorial plane and explodes even with somewhat low neutrino luminosity for which the non-rotating model cannot overcome accreting matter and finally collapses. When the rotational speed is fast, the shock expands to about 300 km immediately after the core bounce and then evolves to move outward without shock stalling. We conclude that this positive effect of rotation to explosions is dominant against some possible negative aspects, for example, lower neutrino luminosity caused by less contraction of the rotating core.

We calculate the neutrino production cross-section in the proto-neutron-star matter under a strong magnetic field in the relativistic mean-field approach. We introduce a new parameter-set which can reproduce the 1.96 solar mass neutron star. We find that the production process increases emitted neutrinos along the direction parallel to the magnetic field and decrease those along its opposite direction. It means that resultant asymmetry due to the neutrino absorption and scattering process in the magnetic field becomes larger by the addition of the neutrino production process.

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., Comment: 50pages, 20 figures, Accepted by ApJ. Latest version with following the referee's suggestions and comments

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, aswell as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are non-linearly 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 access to the most advanced numerical resources such as the “K computer”. To tackle this problem, we have embarked on efforts on multiple fronts. In particular, we report in this paper our recent progresses in the treatment of multidimensional (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 quantitative comparisons will be made between 2D and 3D models. 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 example test simulations. We will also address the perspectives of exascale computations on the next generation supercomputers.

Core‐collapse supernovae (CCSNe) are the final fate of the massive stars, but their explosion mechanisms are still uncertain. One of the clues to the solution of the explosion mechanism is to examine the asymmetric effects. This is because most of observed CCSNe are asymmetric explosions. One of the factors to the asymmetric explosions are the magnetorotational effects. The magnetic fields are amplified intensively along the rotational axsis during the collapse, and it leads to the bipolar outflows which may eject outer mantle. To understand the role of magnetorotational effects during CCSNe, we have developed a new multidimensional magnetohydrodynamic(MHD) code and calculate collapse of a 25 M⊙ star with various magnetic field and rotational velocity.

Recent hydrodynamic studies of core-collapse supernovae imply that the neutrino-heated ejecta from a nascent neutron star develops to supersonic outflows. These supersonic winds are influenced by the reverse shock from the preceding supernova ejecta, forming the wind termination shock. We investigate the effects of the termination shock in neutrino-driven winds and its roll on the r-process. Supersonic outflows are calculated with a semi-analytic neutrino-driven wind model. Subsequent termination-shocked, subsonic outflows are obtained by applying the Rankine-Hugoniot relations. We find a couple of effects that can be relevant for the r-process. First is the sudden slowdown of the temperature decrease by the wind termination. Second is the entropy jump by termination-shock heating, up to several 100NAk. Nucleosynthesis calculations in the obtained winds are performed to examine these effects on the r-process. We find that 1) the slowdown of the temperature decrease plays a decisive roll to determine the r-process abundance curves. This is due to the strong dependences of the nucleosynthetic path on the temperature during the r-process freezeout phase. Our results suggest that only the termination-shocked winds with relatively small shock radii (~500km) are relevant for the bulk of the solar r-process abundances (A~100-180). The heaviest part in the solar r-process curve (A~180-200), however, can be reproduced both in shocked and unshocked winds. These results may help to constrain the mass range of supernova progenitors relevant for the r-process. We find, on the other hand, 2) negligible roles of the entropy jump on the r-process. This is a consequence that the sizable entropy increase takes place only at a large shock radius (~10,000km) where the r-process has already ceased., 11 pages, 7 figures, submitted to ApJ, revised following referee's comments,Accepted for publication in ApJ

Optical spectroscopy and photometry of SN 2006aj have been performed with the Subaru telescope at t > 200 days after GRB060218, the X-ray Flash with which it was associated. Strong nebular emission-lines with an expansion velocity of v ~ 7,300 km/s were detected. The peaked but relatively broad [OI]6300,6363 suggests the existence of ~ 2 Msun of materials in which ~1.3 Msun is oxygen. The core might be produced by a mildly asymmetric explosion. The spectra are unique among SNe Ic in (1) the absence of [CaII]7291,7324 emission, and (2) a strong emission feature at ~ 7400A, which requires ~ 0.05 Msun of newly-synthesized 58Ni. Such a large amount of stable neutron-rich Ni strongly indicates the formation of a neutron star. The progenitor and the explosion energy are constrained to 18 Msun < Mms < 22 Msun and E ~ (1 - 3) 10^{51} erg, respectively., Accepted for Publication in the Astrophysical Journal Letters (2007, ApJ, 658, L5). 8 pages, including 1 table and 3 figures. Typos corrected