-
Protoneutron Star Convection Simulated with a New General Relativistic Boltzmann Neutrino Radiation-Hydrodynamics Code
Authors:
Ryuichiro Akaho,
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The origi…
▽ More
We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial condition is linearly unstable to convection. We find in our 2D simulation that convection is instigated there indeed and expands inward after being full-fledged. The convection is then settled to a quasi-steady state in $\sim100\,\mathrm{ms}$, being sustained by the negative $Y_e$ gradient, which is in turn maintained by neutrino emissions. It enhances the luminosities and mean energies of all species of neutrinos compared to 1D. Taking advantage of the Boltzmann solver, we analyze the possible occurrence of the neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in the regions, where $Y_e$ is lower, and that the growth rate can be as high as $\sim 10^{-1}\,{\mathrm{cm}^{-1}}$.
△ Less
Submitted 20 December, 2022; v1 submitted 3 June, 2022;
originally announced June 2022.
-
Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae
Authors:
Wakana Iwakami,
Akira Harada,
Hiroki Nagakura,
Ryuichiro Akaho,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
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 meth…
▽ More
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.
△ Less
Submitted 21 May, 2022; v1 submitted 13 September, 2021;
originally announced September 2021.
-
Gravitational Wave Physics and Astronomy in the nascent era
Authors:
Makoto Arimoto,
Hideki Asada,
Michael L. Cherry,
Michiko S. Fujii,
Yasushi Fukazawa,
Akira Harada,
Kazuhiro Hayama,
Takashi Hosokawa,
Kunihito Ioka,
Yoichi Itoh,
Nobuyuki Kanda,
Koji S. Kawabata,
Kyohei Kawaguchi,
Nobuyuki Kawai,
Tsutomu Kobayashi,
Kazunori Kohri,
Yusuke Koshio,
Kei Kotake,
Jun Kumamoto,
Masahiro N. Machida,
Hideo Matsufuru,
Tatehiro Mihara,
Masaki Mori,
Tomoki Morokuma,
Shinji Mukohyama
, et al. (28 additional authors not shown)
Abstract:
The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion,…
▽ More
The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics.
△ Less
Submitted 6 April, 2021;
originally announced April 2021.
-
Multidimensional Boltzmann Neutrino Transport Code in Full General Relativity for Core-collapse Simulations
Authors:
Ryuichiro Akaho,
Akira Harada,
Hiroki Nagakura,
Kohsuke Sumiyoshi,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse super…
▽ More
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our code's capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code.
△ Less
Submitted 19 March, 2021; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Simulations of the Early Post-Bounce Phase of Core-Collapse Supernovae in Three-Dimensional Space with Full Boltzmann Neutrino Transport
Authors:
Wakana Iwakami,
Hirotada Okawa,
Hiroki Nagakura,
Akira Harada,
Shun Furusawa,
Kosuke Sumiyoshi,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation…
▽ More
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there.
△ Less
Submitted 31 August, 2020; v1 submitted 5 April, 2020;
originally announced April 2020.
-
The Boltzmann-radiation-hydrodynamics Simulations of Core-collapse Supernovae with Different Equations of State: the Role of Nuclear Composition and the Behavior of Neutrinos
Authors:
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Kohsuke Sumiyoshi,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. On…
▽ More
Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The protoneutron star seen in the FS model is more compact than that in the LS model because the existence of multinuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the motion of matter determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semitransparent region.
△ Less
Submitted 17 November, 2020; v1 submitted 19 March, 2020;
originally announced March 2020.
-
On the Neutrino Distributions in Phase Space for the Rotating Core-collapse Supernova Simulated with a Boltzmann-neutrino-radiation-hydrodynamics Code
Authors:
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of…
▽ More
With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.
△ Less
Submitted 6 March, 2019; v1 submitted 29 October, 2018;
originally announced October 2018.
-
Simulations of core-collapse supernovae in spatial axisymmetry with full Boltzmann neutrino transport
Authors:
Hiroki Nagakura,
Wakana Iwakami,
Shun Furusawa,
Hirotada Okawa,
Akira Harada,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideo Matsufuru,
Akira Imakura
Abstract:
We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing differ…
▽ More
We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r θ-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the θθ-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.
△ Less
Submitted 12 February, 2018; v1 submitted 6 February, 2017;
originally announced February 2017.
-
Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars II. The Implementation of moving-mesh for neutron star kicks
Authors:
Hiroki Nagakura,
Wakana Iwakami,
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideo Matsufuru,
Akira Imakura
Abstract:
We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolu…
▽ More
We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.
△ Less
Submitted 2 May, 2016;
originally announced May 2016.
-
Multi-dimensional Features of Neutrino Transfer in Core-Collapse Supernovae
Authors:
K. Sumiyoshi,
T. Takiwaki,
H. Matsufuru,
S. Yamada
Abstract:
We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six dimensional (6D) phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolutions to stationary states of neutrino distribution funct…
▽ More
We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six dimensional (6D) phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolutions to stationary states of neutrino distribution functions by our Boltzmann solver. Taking advantage of the multi-angle and multi-energy feature realized by the S$_n$ method in our code, we reveal the genuine characteristics of spatially three dimensional (3D) neutrino transfer such as non-radial fluxes and non-diagonal Eddington tensors. In addition, we assess the ray-by-ray approximation, turning off the lateral-transport terms in our code. We demonstrate that the ray-by-ray approximation tends to propagate fluctuations in thermodynamical states around the neutrino-sphere along each radial ray and overestimate the variations between the neutrino distributions on different radial rays. We find that the difference in the densities and fluxes of neutrinos between the ray-by-ray approximation and the full Boltzmann transport becomes ~20%, which is also the case for the local heating rate, whereas the volume-integrated heating rate in the Boltzmann transport is found to be only slightly larger (~2%) than the counterpart in the ray-by-ray approximation due to cancellation among different rays. These results suggest that we had better assess carefully the possible influences of various approximations in the neutrino transfer employed in the current simulations on supernova dynamics. Detailed information on the angle and energy moments of neutrino distribution functions will be profitable for the future development of numerical methods in neutrino-radiation hydrodynamics.
△ Less
Submitted 12 November, 2014; v1 submitted 18 March, 2014;
originally announced March 2014.
-
Nucleon sigma term and strange quark content from lattice QCD with exact chiral symmetry
Authors:
H. Ohki,
H. Fukaya,
S. Hashimoto,
T. Kaneko,
H. Matsufuru,
J. Noaki,
T. Onogi,
E. Shintani,
N. Yamada,
for JLQCD Collaboration
Abstract:
We calculate the nucleon sigma term in two-flavor lattice QCD utilizing the Feynman-Hellman theorem. Both sea and valence quarks are described by the overlap fermion formulation, which preserves exact chiral and flavor symmetries on the lattice. We analyse the lattice data for the nucleon mass using the analytical formulae derived from the baryon chiral perturbation theory. From the data at vale…
▽ More
We calculate the nucleon sigma term in two-flavor lattice QCD utilizing the Feynman-Hellman theorem. Both sea and valence quarks are described by the overlap fermion formulation, which preserves exact chiral and flavor symmetries on the lattice. We analyse the lattice data for the nucleon mass using the analytical formulae derived from the baryon chiral perturbation theory. From the data at valence quark mass set different from sea quark mass, we may extract the sea quark contribution to the sigma term, which corresponds to the strange quark content. We find that the strange quark content is much smaller than the previous lattice calculations and phenomenological estimates.
△ Less
Submitted 21 August, 2008; v1 submitted 29 June, 2008;
originally announced June 2008.