BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shoichi Yamada

Shoichi Yamada contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcastro-phmath.NANumerical Analysisastro-ph.IMhep-phMachine Learningnucl-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 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
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
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
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
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