BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Shoichi Yamada

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

ResearcherUnclaimed source record
astro-ph.HEastro-phnucl-thastro-ph.SRhep-phgr-qcastro-ph.EPmath.NANumerical Analysisastro-ph.COastro-ph.IMMachine Learningmath-phmath.MPphysics.comp-phphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

61works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

61 published item(s)

preprint2024arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A novel formulation for the evolution of relativistic rotating stars

We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational wave signals in the deci-Hz range from neutrinos during the proto-neutron star cooling phase

We investigate the gravitational waves (GWs) at low frequencies produced by neutrinos that are emitted anisotropically from the proto-neutron star (PNS) during its cooling phase that lasts about a minute. We first give a formulation based on the spherical-harmonic expansion of the neutrino luminosity to obtain the gravitational waveform as well as the characteristic strain. In the absence of multi-dimensional simulations of PNS cooling, from which we can extract reliable data on the neutrino luminosities as a function of solid angle, we construct them by hand. In the first model, the time evolution is approximated by piece-wise exponential functions (PEFs); in the second model we employ the time evolution obtained in a 1D cooling simulation for all harmonic components for simplicity; In both cases, we consider not only axisymmetric components but also non-axisymmetric ones; for the latter, in particular, we consider as the third model axisymetric neutrino emissions that are misaligned with the rotation axis and, as a result, rotate with the PNS. We find from the first model that the decay times in PEF at late phases can be inferred from the positions of bumps and dips in the characteristic strain of the GW in the case of a slow cooling, whereas it may be obtained by identifying the positions of slope change in the case of rapid cooling, which may be induced by convections in PNS. We also find from the first model that the rotation of PNS can be reflected in GW waveforms. Finally, we compared the characteristic strains with sensitivity curves GW detectors, finding that above phenomena are expected to be detected by DECIGO if the PNS is located at 10kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Singularity-Avoiding Multi-Dimensional Root-Finder

We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Deep Learning of the Eddington Tensor in the Core-collapse Supernova Simulation

We trained deep neural networks (DNNs) as a function of the neutrino energy density, flux, and the fluid velocity to reproduce the Eddington tensor for neutrinos obtained in our first-principles core-collapse supernova (CCSN) simulations. Although the moment method, which is one of the most popular approximations for neutrino transport, requires a closure relation, none of the analytical closure relations commonly employed in the literature captures all aspects of the neutrino angular distribution in momentum space. In this paper, we developed a closure relation by using the DNN that takes the neutrino energy density, flux, and the fluid velocity as the input and the Eddington tensor as the output. We consider two kinds of DNNs: a conventional DNN named a component-wise neural network (CWNN) and a tensor-basis neural network (TBNN). We found that the diagonal component of the Eddington tensor is reproduced better by the DNNs than the M1-closure relation especially for low to intermediate energies. For the off-diagonal component, the DNNs agree better with the Boltzmann solver than the M1 closure at large radii. In the comparison between the two DNNs, the TBNN has slightly better performance than the CWNN. With the new closure relations at hand based on the DNNs that well reproduce the Eddington tensor with much smaller costs, we opened up a new possibility for the moment method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

General Relativistic Effects on Hill Stability of Multi-Body Systems I: Stability of Three-Body Systems Containing a Massive Black Hole

We study the effects of general relativistic gravity on the Hill stability, that is, the stability of a multi-body system against a close approach of one orbit to another, which has been hitherto studied mainly in Newtonian mechanics and applied to planetary systems. We focus in this paper on the three-body problem and extend the Newtonian analyses to the general relativistic regime in the post-Newtonian approximation. The approximate sufficient condition for the relativistic Hill stability of three-body systems is derived analytically and its validity and usefulness are confirmed numerically. In fact, relativity makes the system more unstable than Newtonian mechanics in the sense of the Hill stability as expected by our theoretical prediction. The criterion will be useful to analyze the results of large-scale N-body simulations of dense environments, in which the stability of three-body sub-systems is important.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Hydrodynamical Study on the Conversion of Hadronic Matter to Quark Matter: II. Diffusion-Induced Conversion

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM), regarding the conversion processes as combustion and describing them hydrodynamically. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider in this paper the conversion induced by diffusions of seed 3QM. This is a sequel to our previous paper, in which the shock-induced conversion was studied in the same frame work. We not only pay attention to the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions. We employ for HM the Shen's EOS, which is based on the relativistic mean field theory, and the bag model-based EOS for QM just as in the previous paper. We demonstrated in that paper that in this combination of EOS's the combustion will occur for a wide range of the bag constant and strong coupling constant in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the initial state. We find that weak deflagration nearly always occurs and that weak detonation is possible only when the diffusion constant is (unrealistically) large and the critical strange fraction is small. The velocities of the conversion front are ~ $10^3-10^7$ cm/s depending on the initial temperature and density as well as the parameters in the QM EOS and become particularly small when the final state is in the mixed phase. Finally we study linear stability of the laminar weak-deflagration front and find that it is unstable in the exothermic regime (Darrius-Landau instability) but stable in the endothermic regime, which is quite contrary to the ordinary combustions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Circular polarizations of gravitational waves from core-collapse supernovae: a clear indication of rapid rotation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rotational equilibria by Lagrangian variational principle: toward multi-dimensional stellar evolutions

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme for the barotropic case, and those by Fujisawa's self-consistent field scheme for the baroclinic case. Included in these rotational configurations are those with shellular-type rotations, which are commonly assumed in the evolution calculation of rotating stars. Although radiation processes, convections and meridional flows have not been taken into account in this study, we have in mind the application of this method to the two-dimensional evolution calculations of rotating stars, for which the Lagrangian formulation is best suited.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Evolution and Impacts of Magnetorotational Instability in Magnetized Core-Collapse Supernovae

We carried out 2D-axisymmetric MHD simulations of core-collapse supernovae for rapidly-rotating magnetized progenitors. By changing both the strength of the magnetic field and the spatial resolution, the evolution of the magnetorotational instability (MRI) and its impacts upon the dynamics are investigated. We found that the MRI greatly amplifies the seed magnetic fields in the regime where not the Alfvén mode but the buoyant mode plays a primary role in the exponential growth phase. The MRI indeed has a powerful impact on the supernova dynamics. It makes the shock expansion faster and the explosion more energetic, with some models being accompanied by the collimated-jet formations. These effects, however, are not made by the magnetic pressure except for the collimated-jet formations. The angular momentum transfer induced by the MRI causes the expansion of the heating region, by which the accreting matter gain an additional time to be heated by neutrinos. The MRI also drifts low-$Y_p$ matter from the deep inside of the core to the heating region, which makes the net neutrino heating rate larger by the reduction of the cooling due to the electron capture. These two effects enhance the efficiency of the neutrino heating, which is found to be the key to boost the explosion. Indeed we found that our models explode far more weakly when the net neutrino heating is switched off. The contribution of the neutrino heating to the explosion energy could reach 60\% even in the case of strongest magnetic field in the current simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hydrodynamical Study on the Conversion of Hadronic Matter to Quark Matter: I. Shock-Induced Conversion

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM) locally, regarding the conversion processes as combustion and describing them hydrodynamically. Not only the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions as well as equations of state (EOS's) in the mixed phase. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider the transition via 2QM triggered by a rapid density rise in a shock wave. Based on the results, we discuss which combustion modes (strong/weak detonation) may be realized. HM is described by an EOS based on the relativistic mean field theory and 2, 3QM's are approximated by the MIT bag model. We demonstrate for a wide range of bag constant and strong coupling constant in this combination of EOS's that the combustion may occur in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the one for the shock wave in P-V plane, and which has no terrestrial counter part. We find that strong detonation always occurs. Depending on the EOS of quark matter (QM) as well as the density of HM and the Mach number of the detonation front, deconfinement from HM to 2QM is either completed or not completed in the shock wave. In the latter case, which is more likely if the EOS of QM ensures that deconfinement occurs above the nuclear saturation density and that the maximum mass of cold quark stars is larger than two solar mass, the conversion continues further via the mixing state of HM and 3QM on the time scale of weak interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exact and approximate expressions of energy generation rates and their impact on the explosion properties of Pair Instability Supernovae

Energetics of nuclear reaction is fundamentally important to understand the mechanism of pair instability supernovae (PISNe). Based on the hydrodynamic equations and thermodynamic relations, we derive exact expressions for energy conservation suitable to be solved in simulation. We also show that some formulae commonly used in the literature are obtained as approximations of the exact expressions. We simulate the evolution of very massive stars of ~100-320 Msun with zero- and 1/10 Zsun, and calculate further explosions as PISNe, applying each of the exact and approximate formulae. The calculations demonstrate that the explosion properties of PISN, such as the mass range, the 56Ni yield, and the explosion energy, are significantly affected by applying the different energy generation rates. We discuss how these results affect the estimate of the PISN detection rate, which depends on the theoretical predictions of such explosion properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Possible Signatures of Ejecta-Companion Interaction in iPTF 13bvn

We investigate the possible effects of the supernova ejecta hitting the companion star in iPTF 13bvn, focusing on the observable features when it becomes visible. iPTF 13bvn is a type Ib supernova that may become the first case that its progenitor is identified as a binary by near future observations. According to calculations by Bersten et al. (2014), the progenitor should have a mass $\approx3.5M_\odot$ to reproduce the supernova light curve, and such compact stars could only be produced via binary evolution. This is one of the reasons that we expect the progenitor to be a binary, but it should be confirmed by observing the remaining companion after the supernova. Their evolutionary calculations suggest that the companion star will be an overluminous OB star at the moment of supernova. With a combination of hydrodynamical and evolutionary simulations, we find that the secondary star will be heated by the supernova ejecta and expand to have larger luminosities and lower surface effective temperatures. The star will look rather like a red super giant, and this should be taken into account when searching for the companion star in the supernova ejecta in future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pre-supernova neutrino emissions from ONe cores in the progenitors of core-collapse supernovae: are they distinguishable from those of Fe cores?

Aiming to distinguish two types of progenitors of core collapse supernovae, i.e., one with a core composed mainly of oxygen and neon (abbreviated as ONe core) and the other with an iron core (or Fe core), we calculated the luminosities and spectra of neutrinos emitted from these cores prior to gravitational collapse, taking neutrino oscillation into account. We found that the total energies emitted as $\barν_e$ from the ONe core are $\lesssim 10^{46}\ {\rm erg}$, which is much smaller than $\sim 10^{47}\ {\rm erg}$ for Fe cores. The average energy, on the other hand, is twice as large for the ONe core as those for the Fe cores. The neutrinos produced by the plasmon decays in the ONe core are more numerous than those from the electron-positron annihilation in both cores but they have much lower average energies $\lesssim 1\ {\rm MeV}$. Although it is difficult to detect the pre-supernova neutrinos from the ONe core even if it is located within 200$\ $pc from the earth, we expect $ \sim 9 - 43$ and $\sim 7 - 61$ events for Fe cores at KamLAND and Super-Kamiokande, respectively, depending on the progenitor mass and neutrino-mass hierarchy. These numbers might be increased by an order of magnitude if we envisage next-generation detectors such as JUNO. We will hence be able to distinguish the two types of progenitors by the detection or non-detection of the pre-supernova neutrinos if they are close enough ($\lesssim1\ {\rm kpc}$).

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical Surface for Explosions of Rotational Core-Collapse Supernovae

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Influence of Magnetorotational Instability on Neutrino Heating: A New Mechanism for Weakly Magnetized Core-Collapse Supernovae

We investigated the impacts of magnetorotational instability (MRI) on the dynamics of weakly magnetized, rapidly rotating core-collapse by conducting high resolution MHD simulations in axisymmetry with simplified neutrino transfer. We found that an initially sub-magnetar class magnetic field is drastically amplified by MRI and substantially affects the dynamics thereafter. Although the magnetic pressure is not strong enough to eject matter, the amplified magnetic field efficiently transfers angular momentum from higher to lower latitudes, which causes the expansion of the heating region at low latitudes due to the extra centrifugal force. This then enhance the efficiency of neutrino heating and eventually leads to neutrino-driven explosion. This is a new scenario of core-collapse supernovae that has never been demonstrated by numerical simulations so far.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear analysis on the growth of non-spherical perturbations in supersonic accretion flows

We analyzed the growth of non-spherical perturbations in supersonic accretion flows. We have in mind the application to the post-bounce phase of core-collapse supernovae (CCSNe). Such non-spherical perturbations have been suggested by a series of papers by Arnett, who has numerically investigated violent convections in the outer layers of pre-collapse stars. Moreover, Couch & Ott (2013) demonstrated in their numerical simulations that such perturbations may lead to a successful supernova even for a progenitor that fails to explode without the fluctuations. This study investigated the linear growth of perturbations during the infall onto a stalled shock wave. The linearized equations are solved as an initial and boundary value problem with the use of Laplace transform. The background is a Bondi accretion flow whose parameters are chosen to mimic the 15 $\mathrm{M_\odot}$ progenitor model by Woosley & Heger (2007), which is supposed to be a typical progenitor of CCSNe. We found that the perturbations that are given at a large radius grow as they flow down to the shock radius; the density perturbations can be amplified by a factor of 30, for example. We analytically showed that the growth rate is proportional to $l$, the index of the spherical harmonics. We also found that the perturbations oscillate in time with frequencies that are similar to those of the standing accretion shock instability. This may have an implication for shock revival in CCSNe, which will be investigated in our forthcoming paper in more detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Outcome of Supernovae in Massive Binaries; Removed Mass, and its Separation Dependence

The majority of massive stars are formed in binary systems. It is hence reasonable to expect that most core-collapse supernovae (CCSNe) take place in binaries and the existence of a companion star may leave some imprints in observed features. Having this in mind, we have conducted two-dimensional hydrodynamical simulations of the collisions of CCSNe ejecta with the companion star in an almost-equal-mass ($\sim 10M_\odot$) binary to find out possible consequences of such events. In particular we pay attention to the amount of mass removed and its dependence on the binary separation. In contrast to the previous surmise, we find that the companion mass is stripped not by momentum transfer but by shock heating. Up to $25\%$ of the original mass can be removed for the closest separations and the removed mass decreases as $M_{ub} \propto a^{-4.3}$ with the binary separation $a$. By performing some experimental computations with artificially-modified densities of incident ejecta, we show that if the velocity of ejecta is fixed, the density of incident ejecta is the single important parameter that actually determines the removed mass as $M_{ub} \propto ρ_{ej} ^{1.4} $. On the other hand, another set of simulations with modified velocities of incident ejecta demonstrate that the strength of the forward shock, which heats up the stellar material and causes the mass loss of the companion star, is actually the key parameter for the removed mass.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of Resistivity on Magnetized Core-Collapse Supernovae

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global Simulations of Magnetorotational Instability in The Collapsed Core of A Massive Star

We performed the first numerical simulations of magnetorotational instability from a sub-magnetar-class seed magnetic field in core collapse supernovae. As a result of axisymmetric ideal MHD simulations, we found that the magnetic field is greatly amplified to magnetar-class strength. In saturation phase, a substantial part of the core is dominated by turbulence, and the magnetic field possesses dominant large scale components, comparable to the size of the proto-neutron star. A pattern of coherent chanel flows, which generally appears during exponential growth phase in previous local simulations, is not observed in our global simulations. While the approximate convergence in the exponential growth rate is attained by increasing spatial resolution, that of the saturation magnetic field is not achieved due to still large numerical diffusion. Although the effect of magnetic field on the dynamics is found to be mild, a simulation with a high-enough resolution might result in a larger impact.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New equations of state based on the liquid drop model of heavy nuclei and quantum approach to light nuclei for core-collapse supernova simulations

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an extension of the previous model, in which we adopted the liquid drop model to all nuclei under the nuclear statistical equilibrium. We reformulate the new liquid drop model so that the temperature dependences of bulk energies could be taken into account. Furthermore, we extend the region in the nuclear chart, in which shell affects are included, by using theoretical mass data in addition to experimental ones. We also adopt a quantum theoretical mass evaluation of light nuclei, which incorporates the Pauli- and self-energy shifts that are not included in the ordinary liquid drop model. The pasta phases for heavy nuclei are taken into account in the same way as in the previous model. We find that the abundances of heavy nuclei are modified by the shell effects of nuclei and temperature dependence of bulk energies. These changes may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. The abundances of light nuclei are also modified by the new mass evaluation, which may affect the heating and cooling rates of supernova cores and shocked envelopes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Photospheric emission from stratified jets

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stellar Core Collapse with Hadron-Quark Phase Transition

Hadronic matter undergoes a deconfinement transition to quark matter at high temperature and/or high density. It would be realized in collapsing cores of massive stars. In the framework of MIT bag model, the ambiguities of the interaction are encapsulated in the bag constant. Some progenitor stars that invoke the core collapses explode as supernovae, and other ones become black holes. The fates of core collapses are investigated for various cases. Equations of state including the hadron-quark phase transition are constructed for the cases of the bag constant B=90, 150 and 250 MeV fm^{-3}. To describe the mixed phase, the Gibbs condition is used. Adopting the equations of state with different bag constants, the core collapse simulations are performed for the progenitor models with 15 and 40Msolar. If the bag constant is small as B=90 MeV fm^{-3}, an interval between the bounce and black hole formation is shortened drastically for the model with 40Msolar and the second bounce revives the shock wave leading to explosion for the model with 15Msolar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Supernova Neutrino Light Curves and Spectra for Various Progenitor Stars: From Core Collapse to Proto-neutron Star Cooling

We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino Transfer in Three Dimensions for Core-Collapse Supernovae. I. Static Configurations

We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Regular and non-regular solutions of the Riemann problem in ideal magnetohydrodynamics

We have built a code to numerically solve the Riemann problem in ideal magnetohydrodynamics (MHD) for an arbitrary initial condition to investigate a variety of solutions more thoroughly. The code can handle not only regular solutions, in which no intermediate shocks are involved, but also all types of non-regular solutions if any. As a first application, we explored the neighborhood of the initial condition that was first picked up by Brio & Wu (1988) and has been frequently employed in the literature as a standard problem to validate numerical codes. Contrary to the conventional wisdom that there will always be a regular solution, we found an initial condition, for which there is no regular solution but a non-regular one. The latter solution has only regular solutions in its neighborhood and actually sits on the boundary of regular solutions. This implies that the regular solutions are not sufficient to solve the ideal MHD Riemann problem and suggests that at least some types of non-regular solutions are physical. We also demonstrate that the non-regular solutions are not unique. In fact, we found for the Brio & Wu initial condition that there are uncountably many non-regular solutions. This poses an intriguing question: why a particular non-regular solution is always obtained in numerical simulations? This has important ramifications to the discussion of which intermediate shocks are really admissible.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new baryonic equation of state at sub-nuclear densities for core-collapse simulations

We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with the atomic number up to ~ 1000. We have also taken into account the pasta phase, thanks to which the transition to uniform nuclear matter in our EOS occurs in the conventional manner: nuclei are not dissociated to nucleons but survive right up to the transition to uniform nuclear matter. We find that the free energy and other thermodynamical quantities are not very different from those given in the Shen's EOS, one of the standard EOS's that adopt the single nucleus approximation. The average mass is systematically different, on the other hand, which may have an important ramification to the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. It is also interesting that the root mean square of the mass number is not very different from the average mass number, since the former is important for the evaluation of coherent scattering rates on nuclei but has been unavailable so far. The EOS table is currently under construction, which will include the weak interaction rates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hyperon Matter and Black Hole Formation in Failed Supernovae

We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of Σ hyperons. Since hyperons soften the EOS, they shorten the time interval from the bounce to black hole formation, which corresponds to the duration of neutrino emission. This effect is larger for the attractive case than the repulsive case because Σ hyperons appear more easily. In addition, we investigate the impacts of pions to find that they also promotes the recollapse towards the black hole formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring Hadron Physics in Black Hole Formations: a New Promising Target of Neutrino Astronomy

The detection of neutrinos from massive stellar collapses can teach us a lot not only about source objects but also about microphysics working deep inside them. In this study we discuss quantitatively the possibility to extract information on the properties of dense and hot hadronic matter from neutrino signals coming out of black-hole-forming collapses of non-rotational massive stars. Based on our detailed numerical simulations we evaluate the event numbers for SuperKamiokande with neutrino oscillations being fully taken into account. We demonstrate that the event numbers from a Galactic event are large enough not only to detect it but also to distinguish one hadronic equation of state from another by our statistical method assuming the same progenitor model and non-rotation. This means that the massive stellar collapse can be a unique probe into hadron physics and will be a promising target of the nascent neutrino astronomy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of Quarks and Pions on Dynamics and Neutrino Signal of Black Hole Formation in Non-rotating Stellar Core Collapse

In the formation process of black holes, the density and temperature of matter become sufficiently high for quarks and pions to appear. In this study we numerically investigate stellar core collapse and black hole formation taking into account the equations of state involving quarks and/or pions. In our simulations, we utilize a code that solves the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. Initial models with three different masses, namely, 40, 100 and 375Msolar, are adopted. Our results show that quarks and pions shorten the duration of neutrino emission if the collapse bounces before black hole formation. In addition, pions increase the luminosity and average energy of neutrinos before black hole formation. We also find that the hadron-quark phase transition leads to an interesting evolution of temperature. Moreover, the neutrino event number is evaluated for the currently operating neutrino detector, SuperKamiokande, to confirm that it is not only detectable but also affected by the emergence of quarks and pions for Galactic events. While there are some issues, such as hyperons, beyond the scope of this study, this is the first serious attempt to assess the impact of quarks and pions in dynamical simulations of black hole formation and will serve as an important foundation for future studies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-layered configurations in differentially-rotational equilibrium

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Dynamics of Three-Planet Systems: an Approach from Dynamical System

We study in detail the motions of three planets interacting with each other under the influence of a central star. It is known that the system with more than two planets becomes unstable after remaining quasi-stable for long times, leading to highly eccentric orbital motions or ejections of some of the planets. In this paper, we are concerned with the underlying physics for this quasi-stability as well as the subsequent instability and advocate the so-called "stagnant motion" in the phase space, which has been explored in the field of dynamical system. We employ the Lyapunov exponent, the power spectra of orbital elements and the distribution of the durations of quasi-stable motions to analyze the phase space structure of the three-planet system, the simplest and hopefully representative one that shows the instability. We find from the Lyapunov exponent that the system is almost non-chaotic in the initial quasi-stable state whereas it becomes intermittently chaotic thereafter. The non-chaotic motions produce the horizontal dense band in the action-angle plot whereas the voids correspond to the chaotic motions. We obtain power laws for the power spectra of orbital eccentricities. Power-law distributions are also found for the durations of quasi-stable states. All these results combined together, we may reach the following picture: the phase space consists of the so-called KAM tori surrounded by satellite tori and imbedded in the chaotic sea. The satellite tori have a self-similar distribution and are responsible for the scale-free power-law distributions of the duration times. The system is trapped around one of the KAM torus and the satellites for a long time (the stagnant motion) and moves to another KAM torus with its own satellites from time to time, corresponding to the intermittent chaotic behaviors.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Oscillation and Future Detection of Failed Supernova Neutrinos from Black Hole Forming Collapse

Recently, stellar collapse involving black hole formation from massive stars is suggested to emit enormous fluxes of neutrinos on par with ordinary core-collapse supernovae. We investigate their detectability for the currently operating neutrino detector, SuperKamiokande. Neutrino oscillation is also taken into account for the evaluation. We find that the event number is larger than or comparable to that of supernova neutrinos and, hence, black hole formation is also a candidate for neutrino astronomy. Moreover, we find that the event number depends dominantly on the equation of state used in the computations of the black hole formation. This fact implies that the detection of neutrinos emitted from the black hole progenitors is very valuable to probe the properties of the equation of state for hot and/or dense matter.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Numerical Study on Stellar Core Collapse and Neutrino Emission: Probe into the Spherically Symmetric Black Hole Progenitors with 3 - 30Msun Iron Cores

The existence of various anomalous stars, such as the first stars in the universe or stars produced by stellar mergers, has been recently proposed. Some of these stars will result in black hole formation. In this study, we investigate iron core collapse and black hole formation systematically for the iron-core mass range of 3 - 30Msun, which has not been studied well so far. Models used here are mostly isentropic iron cores that may be produced in merged stars in the present universe but we also employ a model that is meant for a Population III star and is obtained by evolutionary calculation. We solve numerically the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. As a result, we find that massive iron cores with ~10Msun unexpectedly produce a bounce owing to the thermal pressure of nucleons before black hole formation. The features of neutrino signals emitted from such massive iron cores differ in time evolution and spectrum from those of ordinary supernovae. Firstly, the neutronization burst is less remarkable or disappears completely for more massive models because the density is lower at the bounce. Secondly, the spectra of neutrinos, except the electron type, are softer owing to the electron-positron pair creation before the bounce. We also study the effects of the initial density profile, finding that the larger the initial density gradient is, the more steeply the neutronization burst declines. Further more, we suggest a way to probe into the black hole progenitors from the neutrino emission and estimate the event number for the currently operating neutrino detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Effects of Magnetic Fields on Proto-Neutron Star Winds

We discuss effects of magnetic fields on proto-neutron star winds by performing numerical simulation. We assume that the atmosphere of proto-neutron star has a homogenous magnetic field (ranging from ~10^{12} G to ~10^{15} G) perpendicular to the radial direction and examine the dependence of the three key quantities (dynamical time scale, electron fraction, and entropy per baryon) for the successful r-process on the magnetic field strength. Our results show that even with a magneter-class field strength, ~10^{15} G, the feature of the wind dynamics varies only little from that of non-magnetic winds, and that the condition for successful r-process is not realized.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Gravitational Collapse and Neutrino Emission of Population III Massive Stars

Pop III stars are the first stars in the universe. They do not contain metals and their formation and evolution may be different from that of stars of later generations. In fact, according to the theory of star formation, Pop III stars might have very massive components ($\sim 100 - 10000M_\odot$). In this paper, we compute the spherically symmetric gravitational collapse of these Pop III massive stars. We solve the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail. Unlike supermassive stars ($\gtrsim 10^5 M_\odot$), the stars of concern in this paper become opaque to neutrinos. The collapse is simulated until after an apparent horizon is formed. We confirm that the neutrino transfer plays a crucial role in the dynamics of gravitational collapse, and find also that the $β$-equilibration leads to a somewhat unfamiliar evolution of electron fraction. Contrary to the naive expectation, the neutrino spectrum does not become harder for more massive stars. This is mainly because the neutrino cooling is more efficient and the outer core is more massive as the stellar mass increases. Here the outer core is the outer part of the iron core falling supersonically. We also evaluate the flux of relic neutrino from Pop III massive stars. As expected, the detection of these neutrinos is difficult for the currently operating detectors. However, if ever observed, the spectrum will enable us to obtain the information on the formation history of Pop III stars. We investigate 18 models covering the mass range of $300 - 10^4 M_\odot$, making this study the most detailed numerical exploration of spherical gravitational collapse of Pop III massive stars. This will also serve as an important foundation for multi-dimensional investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Gravitational Radiation from Axisymmetric Rotational Core Collapse

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

0 citations0 reviewsNo code linkedOpen access