TY - JOUR
T1 - Principal-axis Analysis of the Eddington Tensor for the Early Post-bounce Phase of Rotational Core-collapse Supernovae
AU - Iwakami, Wakana
AU - Harada, Akira
AU - Nagakura, Hiroki
AU - Akaho, Ryuichiro
AU - Okawa, Hirotada
AU - Furusawa, Shun
AU - Matsufuru, Hideo
AU - Sumiyoshi, Kohsuke
AU - Yamada, Shoichi
N1 - Funding Information:
This research used the K and Fugaku supercomputers and the high-performance computing resources of FX100 at Nagoya University ICTS, Grand Chariot at Hokkaido University, and Oakforest-PACS at JCAHPC through the HPCI System Research Project (Project ID: hp130025, 140211, 150225, 150262, 160071, 160211, 170031, 170230, 170304, 180111, 180179, 180239, 190100, 190160, 200102, 200124, 210050, 210051, 210164, 220047, 220173). The new system Supercomputer “Flow” at Nagoya University ICTS, NEC SX Aurora Tsubasa at KEK, the Research Center for Nuclear Physics (RCNP) at Osaka University, and XC50 and the general common-use computer system operated by CfCA at the National Astronomical Observatory of Japan (NAOJ) also contribute to this study. Large-scale storage of numerical data is provided by Japan Lattice Data Grid (JLDG). This work was supported in part by Grants-in-Aid for Scientific Research (26104006, 15K05093, and 19K03837) and the Grant-in-Aid for Scientific Research on Innovative areas “Gravitational Wave Physics and Astronomy: Genesis” (17H06357 and 17H06365) and “Unraveling the History of the Universe and Matter Evolution with Underground Physics” (19H05802 and 19H05811) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. This work was also partly supported by research programs at the K-computer of the RIKEN R-CCS, HPCI Strategic Program of Japanese MEXT, Priority Issue on Post-K-computer (Elucidation of the Fundamental Laws and Evolution of the universe), the Joint Institute for Computational Fundamental Sciences (JICFus), and the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Young Scientists (Start-up, JP19K23435). A.H. was supported in part by JSPS Grant-in-Aid for Research Activity Start-up (19K23435) and for Early-Career Scientists (21K13913). S.F. was supported by JSPS KAKENHI (19K14723). S.Y. is supported by the Institute for Advanced Theoretical and Experimental Physics, Waseda University, and the Waseda University Grant for Special Research Projects (project No. 2020-C273).
Publisher Copyright:
© 2022. The Author(s). Published by the American Astronomical Society.
PY - 2022/7/1
Y1 - 2022/7/1
N2 - Using full Boltzmann neutrino transport, we performed 2D core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2 and 15.0 M ⊙, both rotational and nonrotational. We employed the results obtained in the early post-bounce phase (t 2 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 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.
AB - Using full Boltzmann neutrino transport, we performed 2D core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2 and 15.0 M ⊙, both rotational and nonrotational. We employed the results obtained in the early post-bounce phase (t 2 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 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.
UR - http://www.scopus.com/inward/record.url?scp=85134812766&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85134812766&partnerID=8YFLogxK
U2 - 10.3847/1538-4357/ac714b
DO - 10.3847/1538-4357/ac714b
M3 - Article
AN - SCOPUS:85134812766
SN - 0004-637X
VL - 933
JO - Astrophysical Journal
JF - Astrophysical Journal
IS - 1
M1 - 91
ER -