Source author record

Yudai Suwa

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

39 published item(s)

preprint2023arXiv

Updating the $^{56}$Ni Problem in Core-collapse Supernova Explosion

Details of the core-collapse supernova (CCSN) explosion mechanism still need to be fully understood. There is an increasing number of successful examples of reproducing explosions in multidimensional hydrodynamic simulations, but subsequent studies pointed out that the growth rates of the explosion energy $\dot{E}_\mathrm{expl}$ of these simulations are insufficient to produce enough $^{56}$Ni to match observations. This issue is known as the `$^{56}$Ni problem' in CCSNe. Recently, however, some studies have suggested that this $^{56}$Ni problem is derived from the simplicity of the explosion model. In response, we investigate the effect of the explosion energy growth rate $\dot{E}_\mathrm{expl}$ on the behavior of nucleosynthesis in CCSNe in a more realistic model. We employ the 1D Lagrangian hydrodynamic code, in which we take neutrino heating and cooling terms into account with the light-bulb approximation. We reiterate that, consistent with previous rebuttal studies, there is the $^{56}$Ni problem: Although $^{56}$Ni is synthesized to almost the same mass coordinate independent of $\dot{E}_\mathrm{expl}$, some of the innermost material in the low-$\dot{E}_\mathrm{expl}$ model failed to escape, leading to a shift in the innermost mass coordinate of the ejecta to the outer positions. Comparing our results with observations, we find that while modern slow explosions can, in principle, reproduce observations of standard Type II SNe, this is not possible with stripped-envelope SNe. Our finding places a strong constraint on the explosion mechanism. There are significant differences in the progenitor structures and the explosion mechanism between Type II and stripped-envelope SNe.

preprint2022arXiv

Can Population III stars be major origins of both merging binary black holes and extremely metal poor stars?

Population (Pop) III stars, first stars, or metal-free stars are made of primordial gas. We have examined if they can be dominant origins of merging binary black holes (BHs) and extremely metal-poor stars. The abundance pattern of EMP stars is helpful to trace back the properties of Pop III stars. We have confirmed previous arguments that the observed BH merger rate needs Pop III star formation efficiency 10 times larger than theoretically predicted values, while the cosmic reionization history still permits such a high Pop III star formation efficiency. On the other hand, we have newly found that the elemental abundance pattern of EMP stars only allows the Pop III initial mass function with the minimum mass of $\sim 15 - 27$ $M_\odot$. In other words, the minimum mass must not deviate largely from the critical mass below and above which Pop III stars leave behind neutron stars and BHs, respectively. Pop III stars may be still a dominant origin of merging binary BHs but our study has reduced the allowed parameter space under a hypothesis that EMP stars are formed from primordial gas mixed with Pop III supernova ejecta.

preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. II. Impact of the Nuclear Equation of State

The late-time evolution of the neutrino event rate from supernovae is evaluated for Super-Kamiokande using simulated results of proto-neutron star (PNS) cooling. In the present work we extend the result of Suwa et al. (2019) [arXiv:1904.09996], which studied the dependence on the PNS mass, but focus on the impact of the nuclear equation of state (EOS). We find that the neutrino event rate depends on both the high-density and low-density EOS, where the former determines the radius of the PNS and the latter affects its surface temperature. Based on the present evaluation of the neutrino event rate we propose a new analysis method to extract the time variability of the neutrino average energy taking into account the statistical error in the observation.

preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. III. Extraction of Mass and Radius of Neutron Stars from Synthetic Data

Neutrinos are guaranteed to be observable from the next galactic supernova (SN). Optical light and gravitational waves are also observable, but may be difficult to observe if the location of the SN in the galaxy or the details of the explosion are unsuitable. The key to observing the next supernova is to first use neutrinos to understand various physical quantities and then link them to other signals. In this paper, we present Monte Carlo sampling calculations of neutrino events from galactic supernova explosions observed with Super-Kamiokande. The analytical solution of neutrino emission, which represents the long-term evolution of neutrino-light curve from supernovae, is used as a theoretical template. It gives the event rate and event spectrum through inverse beta decay interactions with explicit model parameter dependence. Parameter estimation is performed on these simulated sample data by fitting least squares using the analytical solution. The results show that the mass, radius and total energy of a remnant neutron star produced by a SN can be determined with an accuracy of $\sim 0.1M_\odot$, $\sim 1$ km, and $\sim 10^{51}$ erg, respectively, for a galactic SN at 8 kpc.

preprint2022arXiv

On the energy source of ultra-stripped supernovae

Ultra-stripped supernovae (USSN) with a relatively low ejecta mass of $\sim0.1M_\odot$ (e.g., iPTF14gqr and SN2019dge) are considered to originate from ultra-stripped carbon-oxygen (CO) cores in close binary systems and are likely to be progenitors of binary neutron stars (BNSs). Here we conduct long-term simulations of USSNe from ultra-stripped progenitors with various masses ($1.45\,M_\odot \leq M_\mathrm{CO} \leq 2.0\,M_\odot$) based on results of neutrino-radiation hydrodynamics simulations, and consistently calculate the nucleosynthesis and the SN light curves. We find that a USSN from a more massive progenitor has a larger ejecta mass but a smaller $^{56}$Ni mass mainly due to the fallback, which leads to the light curve being dimmer and slower. By comparing the synthetic light curves with the observed ones, we show that SN2019dge can be solely powered by $^{56}$Ni synthesized during the explosion of a progenitor with $M_\mathrm{CO} \lesssim 1.6\,M_\odot$ while iPTF14gqr cannot be explained by the $^{56}$Ni powered model; $\sim 0.05M_\odot$ of $^{56}$Ni inferred from the light curve fitting is argued to be difficult to synthesize for ultra-stripped progenitors. We consider fallback accretion onto and rotation-powered relativistic wind from the newborn NS as alternative energy sources and show that iPTF14gqr could be powered by a newborn NS with a magnetic field of $B_p \sim 10^{15}\,\mathrm{G}$ and an initial rotation period of $P_i \sim 0.1\,\mathrm{s}$.

preprint2022arXiv

X-raying the Birth of Binary Neutron Stars and Neutron Star-Black Hole Binaries

We consider fallback accretion after an ultra-stripped supernova (USSN) that accompanies formation of a binary neutron star (BNS) or a neutron star-black hole binary (NS-BH). The fallback matter initially accretes directly to the nascent NS, while it starts to accrete to the circumbinary disk, typically $0.1\mbox{-}1\, \mathrm{day}$ after the onset of the USSN explosion. The circumbinary disk mass further accretes, forming mini disks around each compact object, with a super-Eddington rate up to a few years. We show that such a system constitutes a binary ultraluminous X-ray source (ULX), and a fraction of the X rays can emerge through the USSN ejecta. We encourage follow-up observations of USSNe within $\lesssim 100\,\rm Mpc$ and $\sim 100\mbox{-}1,000\,\mathrm{day}$ after the explosion using Chandra, XMM Newton and NuSTAR, which could detect the X-ray counterpart with time variations representing the properties of the nascent compact binary, e.g., the orbital motion of the binary, the spin of the NS, and/or the quasiperiodic oscillation of the mini disks.

preprint2021arXiv

Analytic solutions for neutrino-light curves of core-collapse supernovae

Neutrinos are a guaranteed signal from supernova explosions in the Milky Way, and a most valuable messenger that can provide us with information about the deepest parts of supernovae. In particular, neutrinos will provide us with physical quantities, such as the radius and mass of protoneutron stars (PNS), which are the central engine of supernovae. This requires a theoretical model that connects observables such as neutrino luminosity and average energy with physical quantities. Here, we show analytic solutions for the neutrino-light curve derived from the neutrino radiation transport equation by employing the diffusion approximation and the analytic density solution of the hydrostatic equation for a PNS. The neutrino luminosity and the average energy as functions of time are explicitly presented, with dependence on PNS mass, radius, the total energy of neutrinos, surface density, and opacity. The analytic solutions provide good representations of the numerical models from a few seconds after the explosion and allow a rough estimate of these physical quantities to be made from observational data.

preprint2020arXiv

A Consistent Modeling of Neutrino-driven Wind with Accretion Flow onto a Protoneutron Star and its Implications for $^{56}$Ni Production

Details of the explosion mechanism of core-collapse supernovae (CCSNe) are not yet fully understood. There is now an increasing number of successful examples of reproducing explosions in the first-principles simulations, which have shown a slow increase of explosion energy. However, it was recently pointed out that the growth rates of the explosion energy of these simulations are insufficient to produce enough $^{56}$Ni mass to account for observations. We refer to this issue as the `nickel mass problem' (Ni problem, hereafter) in this paper. The neutrino-driven wind is suggested as one of the most promising candidates for the solution to the Ni problem in previous literature, but a multi-dimensional simulation for this is computationally too expensive to allow long-term investigations. In this paper, we first built a consistent model of the neutrino-driven wind with an accretion flow onto a protoneutron star (PNS), by connecting a steady-state solution of the neutrino-driven wind and a phenomenological mass accretion model. Comparing the results of our model with the results of first-principles simulations, we find that the total ejectable amount of the neutrino-driven wind is roughly determined within $\sim$ 1 sec from the onset of the explosion and the supplementable amount at a late phase ($t_e \gtrsim 1$ sec) remains $M_\mathrm{ej} \lesssim 0.01M_\odot$ at most. Our conclusion is that it is difficult to solve the Ni problem, by continuous injection of $^{56}$Ni by the neutrino-driven wind. We suggest that the total amount of synthesized $^{56}$Ni can be estimated robustly if simulations are followed up to $\sim 2$ seconds.

preprint2020arXiv

Developing an end-to-end simulation framework of supernova neutrino detection

Massive stars can explode as supernovae at the end of their life cycle, releasing neutrinos whose total energy reaches $10^{53}$ erg. Moreover, neutrinos play key roles in supernovae, heating and reviving the shock wave as well as cooling the resulting protoneutron star. Therefore, neutrino detectors are waiting to observe the next galactic supernova and several theoretical simulations of supernova neutrinos are underway. While these simulation concentrate mainly on only the first one second after the supernova bounce, the only observation of a supernova with neutrinos, SN 1987A, revealed that neutrino emission lasts for more than 10 seconds. For this reason, long-time simulation and analysis tools are needed to compare theories with the next observation. Our study is to develop an integrated supernova analysis framework to prepare an analysis pipeline for treating galactic supernovae observations in the near future. This framework deals with the core-collapse, bounce and proto-neutron star cooling processes, as well as with neutrino detection on earth in a consistent manner. We have developed a new long-time supernova simulation in one dimension that explodes successfully and computes the neutrino emission for up to 20 seconds. Using this model we estimate the resulting neutrino signal in the Super-Kamiokande detector to be about 1,800 events for an explosion at 10 kpc and discuss its implications in this paper. We compare this result with the SN 1987A observation to test its reliability.

preprint2016arXiv

Parametric initial conditions for core-collapse supernova simulations

We investigate a method to construct parametrized progenitor models for core-collapse supernova simulations. Different from all modern core-collapse supernova studies, which rely on progenitor models from stellar evolution calculations, we follow the methodology of Baron & Cooperstein (1990) to construct initial models. Choosing parametrized spatial distributions of entropy and electron fraction as a function of mass coordinate and solving the equation of hydrostatic equilibrium, we obtain the initial density structures of our progenitor models. First, we calculate structures with parameters fitting broadly the evolutionary model s11.2 of Woosley et al. (2002). We then demonstrate the reliability of our method by performing general relativistic hydrodynamic simulations in spherical symmetry with the isotropic diffusion source approximation to solve the neutrino transport. Our comprehensive parameter study shows that initial models with a small central entropy ($\lesssim 0.4\,k_B$ nucleon$^{-1}$) can explode even in spherically symmetric simulations. Models with a large entropy ($\gtrsim 6\,k_B$ nucleon$^{-1}$) in the Si/O layer have a rather large explosion energy ($\sim 4\times 10^{50}$ erg) at the end of the simulations, which is still rapidly increasing.

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

How Much Can $^{56}$Ni Be Synthesized by Magnetar Model for Long Gamma-ray Bursts and Hypernovae?

A rapidly rotating neutron star with strong magnetic fields, called magnetar, is a possible candidate for the central engine of long gamma-ray bursts and hypernovae (HNe). We solve the evolution of a shock wave driven by the wind from magnetar and evaluate the temperature evolution, by which we estimate the amount of $^{56}$Ni that produces a bright emission of HNe. We obtain a constraint on the magnetar parameters, namely the poloidal magnetic field strength ($B_p$) and initial angular velocity ($Ω_i$), for synthesizing enough $^{56}$Ni mass to explain HNe ($M_{^{56}\mathrm{Ni}}\gtrsim 0.2M_\odot$), i.e. $(B_p/10^{16}~\mathrm{G})^{1/2}(Ω_i/10^4~\mathrm{rad~s}^{-1})\gtrsim 0.7$.

preprint2015arXiv

Neutrino-driven explosions of ultra-stripped type Ic supernovae generating binary neutron stars

We study explosion characteristics of ultra-stripped supernovae (SNe), which are candidates of SNe generating binary neutron stars (NSs). As a first step, we perform stellar evolutionary simulations of bare carbon-oxygen cores of mass from 1.45 to 2.0 $M_\odot$ until the iron cores become unstable and start collapsing. We then perform axisymmetric hydrodynamics simulations with spectral neutrino transport using these stellar evolution outcomes as initial conditions. All models exhibit successful explosions driven by neutrino heating. The diagnostic explosion energy, ejecta mass, Ni mass, and NS mass are typically $\sim 10^{50}$ erg, $\sim 0.1 M_\odot$, $\sim 0.01M_\odot$, and $\approx 1.3 M_\odot$, which are compatible with observations of rapidly-evolving and luminous transient such as SN 2005ek. We also find that the ultra-stripped SN is a candidate for producing the secondary low-mass NS in the observed compact binary NSs like PSR J0737-3039.

preprint2015arXiv

Optical Synchrotron Precursors of Radio Hypernovae

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

preprint2015arXiv

Probing Rotation of Core-collapse Supernova with Concurrent Analysis of Gravitational Waves and Neutrinos

The next time a core-collapse supernova (SN) explodes in our galaxy, vari- ous detectors will be ready and waiting to detect its emissions of gravitational waves (GWs) and neutrinos. Current numerical simulations have successfully introduced multi-dimensional effects to produce exploding SN models, but thus far the explosion mechanism is not well understood. In this paper, we focus on an investigation of progenitor core rotation via comparison of the start time of GW emission and that of the neutronization burst. The GW and neutrino de- tectors are assumed to be, respectively, the KAGRA detector and a co-located gadolinium-loaded water Cherenkov detector, either EGADS or GADZOOKS!. Our detection simulation studies show that for a nearby supernova (0.2 kpc) we can confirm the lack of core rotation close to 100% of the time, and the presence of core rotation about 90% of the time. Using this approach there is also po- tential to confirm rotation for considerably more distant Milky Way supernova explosions.

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

Anisotropic neutrino effect on magnetar spin: constraint on inner toroidal field

The ultra-strong magnetic field of magnetars modifies the neutrino cross section due to the parity violation of the weak interaction and can induce asymmetric propagation of neutrinos. Such an anisotropic neutrino radiation transfers not only the linear momentum of a neutron star but also the angular momentum, if a strong toroidal field is embedded inside the stellar interior. As such, the hidden toroidal field implied by recent observations potentially affects the rotational spin evolution of new-born magnetars. We analytically solve the transport equation for neutrinos and evaluate the degree of anisotropy that causes the magnetar to spin-up or spin-down during the early neutrino cooling phase. Supposing that after the neutrino cooling phase the dominant process causing the magnetar spin-down is the canonical magnetic dipole radiation, we compare the solution with the observed present rotational periods of anomalous X-ray pulsars 1E 1841-045 and 1E 2259+586, whose poloidal (dipole) fields are $\sim 10^{15}$ G and $10^{14}$ G, respectively. Combining with the supernova remnant age associated with these magnetars, the present evaluation implies a rough constraint of global (average) toroidal field strength at $B^ϕ\lesssim 10^{15}$ G.

preprint2014arXiv

Soft X-ray Extended Emissions of Short Gamma-Ray Bursts as Electromagnetic Counterparts of Compact Binary Mergers; Possible Origin and Detectability

We investigate the possible origin of extended emissions (EEs) of short gamma-ray bursts with an isotropic energy of ~ 10^(50-51) erg and a duration of a few 10 s to ~ 100 s, based on a compact binary (neutron star (NS)-NS or NS-black hole (BH)) merger scenario. We analyze the evolution of magnetized neutrino-dominated accretion disks of mass ~ 0.1 M_sun around BHs formed after the mergers, and estimate the power of relativistic outflows via the Blandford-Znajek (BZ) process. We show that a rotation energy of the BH up to > 10^52 erg can be extracted with an observed time scale of > 30 (1+z) s with a relatively small disk viscosity parameter of alpha < 0.01. Such a BZ power dissipates by clashing with non-relativistic pre-ejected matter of mass M ~ 10^-(2-4) M_sun, and forms a mildly relativistic fireball. We show that the dissipative photospheric emissions from such fireballs are likely in the soft X-ray band (1-10 keV) for M ~ 10^-2 M_sun possibly in NS-NS mergers, and in the BAT band (15-150 keV) for M ~ 10^-4 M_sun possibly in NS-BH mergers. In the former case, such soft EEs can provide a good chance of ~ 6 yr^-1 for simultaneous detections of the gravitational waves with a ~ 0.1 deg angular resolution by soft X-ray survey facilities like Wide-Field MAXI.

preprint2014arXiv

Symmetry energy impact in simulations of core-collapse supernovae

We present a review of a broad selection of nuclear matter equations of state (EOSs) applicable in core-collapse supernova studies. The large variety of nuclear matter properties, such as the symmetry energy, which are covered by these EOSs leads to distinct outcomes in supernova simulations. Many of the currently used EOS models can be ruled out by nuclear experiments, nuclear many-body calculations, and observations of neutron stars. In particular the two classical supernova EOS describe neutron matter poorly. Nevertheless, we explore their impact in supernova simulations since they are commonly used in astrophysics. They serve as extremely soft and stiff representative nuclear models. The corresponding supernova simulations represent two extreme cases, e.g., with respect to the protoneutron star (PNS) compactness and shock evolution. Moreover, in multi-dimensional supernova simulations EOS differences have a strong effect on the explosion dynamics. Because of the extreme behaviors of the classical supernova EOSs we also include DD2, a relativistic mean field EOS with density-dependent couplings, which is in satisfactory agreement with many current nuclear and observational constraints. This is the first time that DD2 is applied to supernova simulations and compared with the classical supernova EOS. We find that the overall behaviour of the latter EOS in supernova simulations lies in between the two extreme classical EOSs. As pointed out in previous studies, we confirm the impact of the symmetry energy on the electron fraction. Furthermore, we find that the symmetry energy becomes less important during the post bounce evolution, where conversely the symmetric part of the EOS becomes increasingly dominating, which is related to the high temperatures obtained. Moreover, we study the possible impact of quark matter at high densities and light nuclear clusters at low and intermediate dens

preprint2013arXiv

An Alternative Numerical Method for the Stationary Pulsar Magnetosphere

Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. In 1999, Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line. This current sheet is artificially added to make the current system closed. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetosphere models without a current sheet. In our method, the pulsar equation is decomposed into Ampere's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model and find that the force-free condition breaks down in a vicinity of the light cylinder due to dissipation that is included implicitly in the numerical method.

preprint2013arXiv

Blue Supergiant Model for Ultra-Long Gamma-Ray Burst with Superluminous-Supernova-Like Bump

Long GRBs (LGRBs) have typical duration of ~ 30 s and some of them are associated with hypernovae, like Type Ic SN 1998bw. Wolf-Rayet stars are the most plausible LGRB progenitors, since the free-fall time of the envelope is consistent with the duration, and the natural outcome of the progenitor is a Type Ic SN. While a new population of ultra-long GRBs (ULGRBs), GRB 111209A, GRB 101225A, and GRB 121027A, has a duration of ~ 10^4 s, two of them are accompanied by superluminous-supernova (SLSN) like bumps, which are <~ 10 times brighter than typical hypernovae. Wolf-Rayet progenitors cannot explain ULGRBs because of too long duration and too bright SN-like bump. A blue supergiant (BSG) progenitor model, however, can explain the duration of ULGRBs. Moreover, SLSN-like bump can be attributed to the so-called cocoon-fireball photospheric emissions (CFPEs). Since a large cocoon is inevitably produced during the relativistic jet piercing though the BSG envelope, this component can be a smoking-gun evidence of BSG model for ULGRBs. In this paper, we examine u, g, r, i, and J-band light curves of three ULGRBs and demonstrate that they can be fitted quite well by our BSG model with the appropriate choices of the jet opening angle and the number density of the ambient gas. In addition, we predict that for 121027A, SLSN-like bump could have been observed for ~ 20 - 80 days after the burst. We also propose that some SLSNe might be CFPEs of off-axis ULGRBs without visible prompt emission.

preprint2013arXiv

From supernovae to neutron stars

The gravitational collapse, bounce, the explosion of an iron core of an 11.2 $M_{\odot}$ star is simulated by two-dimensional neutrino-radiation hydrodynamic code. The explosion is driven by the neutrino heating aided by multi-dimensional hydrodynamic effects such as the convection. Following the explosion phase, we continue the simulation focusing on the thermal evolution of the protoneutron star up to $\sim$70 s when the crust of the neutron star is formed using one-dimensional simulation. We find that the crust forms at high-density region ($ρ\sim10^{14}$ g cm$^{-3}$) and it would proceed from inside to outside. This is the first self-consistent simulation that successfully follows from the collapse phase to the protoneutron star cooling phase based on the multi-dimensional hydrodynamic simulation.

preprint2013arXiv

Luminous supernova-like UV/optical/infrared transients associated with ultra-long gamma-ray bursts from metal-poor blue supergiants

Metal-poor massive stars may typically end up their lives as blue supergiants (BSGs). Gamma-ray bursts (GRBs) from such progenitors could have ultra-long duration of relativistic jets. For example Population III (Pop III) GRBs at z ~ 10-20 might be observable as X-ray rich events with a typical duration of T_90 ~ 10^4(1+z) sec. Recent GRB111209A at z = 0.677 has an ultra long duration of T_90 ~ 2.5*10^4 sec so that it have been suggested that the progenitor might be a metal-poor BSGs in the local universe. Here, we suggest luminous UV/optical/infrared emissions associated with such a new class of GRB from metal poor BSGs. Before the jet head breaks out the progenitor envelope, the energy injected by the jet is stored in a hot-plasma cocoon, which finally emerges and expands as a baryon-loaded fireball. We show that the photospheric emissions from the cocoon fireball could be intrinsically very bright (L_peak ~ 10^(42-44) erg/sec) in UV/optical bands (E_peak ~ 10 eV) with a typical duration of ~ 100 days in the rest frame. Such cocoon emissions from Pop III GRB might be detectable in infrared bands at ~ years after Pop III GRBs at up to z ~ 15 by up-coming facilities like JWST. We also suggest that GRB111209A might have been rebrightening in UV/optical bands up to an AB magnitude of < 26. The cocoon emissions from local metal-poor BSGs might have been already observed as luminous supernovae without GRB since they can be seen from the off-axis direction of the jet.

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

Long Duration X-Ray Flash and X-Ray Rich Gamma Ray Burst from Low Mass Population III Star

Recent numerical simulations suggest that Population III (Pop III) stars were born with masses not larger than $\sim 100 M_{\odot}$ but typically $\sim 40M_{\odot}$. By self-consistently considering the jet generation and propagation in the envelope of these low mass Pop III stars, we find that a Pop III blue super giant star has the possibility to raise a gamma-ray burst (GRB) even though it keeps a massive hydrogen envelope. We evaluate observational characters of Pop III GRBs and predict that Pop III GRBs have the duration of $\sim 10^5$ sec in the observer frame and the peak luminosity of $\sim 5 \times 10^{50} {\rm erg} {\rm sec}^{-1}$. Assuming that the $E_p-L_p$ (or $E_p-E_{γ, \rm iso}$) correlation holds for Pop III GRBs, we find that the spectrum peak energy falls $\sim$ a few keV (or $\sim 100$ keV) in the observer frame. We discuss the detectability of Pop III GRBs by future satellite missions such as EXIST and Lobster. If the $E_p-E_{γ, \rm iso}$ correlation holds, we have the possibility to detect Pop III GRBs at $z \sim 9$ as long duration X-ray rich GRBs by EXIST. On the other hand, if the $E_p-L_p$ correlation holds, we have the possibility to detect Pop III GRBs up to $z \sim 19$ as long duration X-ray flashes by Lobster.

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

Neutrino Acceleration by Bulk Matter Motion and Explosion Mechanism of Gamma-Ray Bursts

The neutrino annihilation is one of the most promising candidates for the jet production process of gamma-ray bursts. Although neutrino interaction rates depend strongly on the neutrino spectrum, the estimations of annihilation rate have been done with an assumption of the neutrino thermal spectrum based on the presence of the neutrinospheres, in which neutrinos and matter couple strongly. We consider the spectral change of neutrinos caused by the scattering by infalling materials and amplification of the annihilation rate. We solve the kinetic equation of neutrinos in spherically symmetric background flow and find that neutrinos are successfully accelerated and partly form nonthermal spectrum. We find that the accelerated neutrinos can significantly enhance the annihilation rate by a factor of $\sim 10$, depending on the injection optical depth.

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

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

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

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

Afterglow of binary neutron star merger

The merger of two neutron stars results often in a rapidly and differentially rotating hypermassive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromagnetic radiation with a large luminosity ~ 0.1 B^2 R^3 Omega is emitted with baryon (B, R, and Omega are poloidal magnetic-field strength at stellar surface, stellar radius, and angular velocity of a HMNS). The predicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole direction, is ~ 10^{47} (B/10^{13} G)^2(R/10 km)^3(Omega/10^4 rad/s) ergs/s, that is comparable to the luminosity of quasars.

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

preprint2011arXiv

Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity

A truncated moment formalism for general relativistic radiation hydrodynamics, based on the Thorne's moment formalism, is derived. The fluid rest frame is chosen to be the fiducial frame for defining the radiation moments. Then, zeroth-, first-, and second-rank radiation moments are defined from the distribution function with a physically reasonable assumption for it in the optically thin and thick limits. The source terms are written, focusing specifically on the neutrino transfer and neglecting higher harmonic angular dependence of the reaction angle. Finally, basic equations for a truncated moment formalism for general relativistic radiation hydrodynamics in a closed covariant form are derived assuming a closure relation among the radiation stress tensor, energy density, and energy flux, and a variable Eddington factor, which works well.

preprint2010arXiv

Can Gamma-Ray Burst Jets Break Out the First Stars?

We show that a relativistic gamma-ray burst (GRB) jet can potentially pierce the envelope of very massive first generation star (Population III; Pop III) by using the stellar density profile to estimate both the jet luminosity (via accretion) and its penetrability. The jet breakout is possible even if the Pop III star has a supergiant hydrogen envelope without mass loss, thanks to the long-lived powerful accretion of the envelope itself. While the Pop III GRB is estimated to be energetic, E_{gamma,iso} ~ 10^{55} erg, the supergiant envelope hides the initial bright phase into the cocoon component, leading to a GRB with a long duration ~ 1000(1+z) sec and an ordinary isotropic luminosity ~ 10^{52} erg s^{-1} (~ 10^{-9} erg cm^{-2} s^{-1} at redshift z ~ 20). The neutrino-annihilation is not effective for Pop III GRBs because of a low central temperature, while the magnetic mechanism is viable. We also derive analytic estimates of the breakout conditions, which are applicable to various progenitor models. The GRB luminosity and duration are found to be very sensitive to the core and envelope mass, providing possible probes of the first luminous objects at the end of the high redshift dark ages.

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.

preprint2009arXiv

Probing the central engine of long gamma-ray bursts and hypernovae with gravitational waves and neutrinos

There are the two common candidates as the viable energy source for the central engine of long gamma-ray bursts (GRBs) and hypernovae (HNe), neutrino annihilation and magnetic fields. We investigate gravitational wave (GW) emission accompanied by these two mechanisms. Especially, we focus on GW signals produced by neutrinos from a hyper-accreting disk around a massive black hole. We show that neutrino-induced GWs are detectable for $\sim$1 Mpc events by LISA and $\sim$ 100 Mpc by DECIGO/BBO, if the central engine is powered by neutrinos. The GW signals depend on the viewing angle and they are anti-correlated with neutrino ones. But, simultaneous neutrino detections are also expected, and helpful for diagnosing the explosion mechanism when later electromagnetic observations enable us to identify the source. GW and neutrino observations are potentially useful for probing choked jets that do not produce prompt emission, as well as successful jets. Even in non-detection cases, observations of GWs and neutrinos could lead to profitable implications for the central engine of GRBs and HNe.

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.