JREKT':>@1;BLNSUBOT/+かに星雲 (SN1054)

 5*3hubblesite.org•1!&;HNUFQI2 /+•) $/+<APTGDMSTC468?0=.91!&;%7
( ", #$ , Vol. 21, 10 -

HPC-Phys@YITP, 2018/08/22

1 • – EPSON HC-40, CP/M
– Macintosh SE/30, SUN workstation• !# %– VPP500/28 @• !# %– KEK SR/BGQ -> YITP XC, RCNP SX– , $!"


1995 ~ 1D
2005 ~ 1D
2015
2007 ~ 2D, 3D~ 2Dn

2*ある日突然、非常に明るい星が現れる実は、重たい星の進化の最期$)('%%&

やがて暗くなり消える# $!$%"
A.K.Mann "Shadow of a Star" (W.H. Freeman and Company, 1997)

34

1#(>3=CBA@Crab Nebula (SN1054)– WPTKGRXV@-&• OVHX (9/5@1.4; 10km)• :YEFIJKG?$/Z– 5@MSXLUNW+DSN1987A• NXQV$%.\[[\• $/W7 @2 D<Prof. Koshiba
Kamiokanden– ,"'>@4*@0!
• #()* • 68@*@

http://www-sk.icrr.u-tokyo.ac.jp/

http://nobelprize.org/


http://periodictable.com/4

6 2](+QINMPD#*6<→8KHD6

→=EC3G→/%+)→(+907D10DD !

\ZZJUAdam Burrows (1995)

FeSiO, Ne, Mg
C, OHeHKH^5

D$4]B8KH
[ -1,FAD".E&B?G
X1: '<;:] @<YX>DLORTWLSVDY55,2/4'.17+
-*Fe core

重力崩壊

νn-閉込め!
'#(%268073 

""

νν ν

1000 kmν"ν

コアバウンス爆発ν!
$&)

超新星ν



NSν"(ν νShockwave
衝撃波の発生元素合成~10 km 649>5?1/&(Hans Bethe• 63- ) (Mcore~1.4Msolar)– RFe~5x103 km→ RNS~20 km• 49>5?'&,. -0/$ GM 2 GM 2 '53
ΔE Grav = −&−~10erg)
RNS (% RFe&(%49>5?&€• 49>5?:

http://www.nobelprize.org

1967:?;>

100[B]

/

Eexp ~ 1051 erg• 8<?7=:49>5?: En ~ 1053 erg

1[B]
99[B]

(*#2$8<?7=:& 0"!+
7

2F0_^]F8• &/)E (J E%B<@:`– LSVY\LWZF8>=•(UKQMTF7G:`– HACH6.DJ"?I8>=-+ : QV\PXRC!1F*';1)• LSVY\LWZ- [6- QV\PXR435• !&O\N- 99#,
- QV\PXR

$8•• %&•  $! ($#/"#)•' $! 9%!< 47:CFEK• # &IMPGNJ(%!2
5?=8→ QD

RD/3D2D

///IMPGNJ('*)



/.

+#

&→ 3D
GR

IMPGO,

• IMPGNJ('*9S)@>-3– $DAHLBC:";T• 2D/3D81

1:

@>6010

$A/R<?6C• 1D"':– ($:?7#%?

!#)A

-A,15Msolar

100 kmShen-EOS
LS-EOS0*3E →&$;F

1D%
LGOHA-9:>:D8• 2D, 3D:

Sumiyoshi et al. (2005)-A,

Takiwaki et al. (2014)

11Msolar
LS-EOS– 40*• QJMPINK • 2D vs 3DS →+2BD=– -A5BS

2D
3D
1D

100 km→ 1@B 0*9.

LGOHA
11

FKODLG'&+4.>=B?EJ@APFKODLGn• 8FKODLG-n– ν/6(;22• 1 3<!:7%*– 01

ProtoNS shockwave

nheating4,"→ )$3<$9
拡散ν加熱

十分に反応/散乱nn nnnn高温/高密度

50km

HMCIN>#5"
自由伝搬
1 ∂fν  
1 ' δfν *+ n ⋅ ∇fν = ) ,nc ∂tc ( δt + collision
衝撃波伝搬

100km

5000km$-4
12

6DCE@DF<'

Sumiyoshi & Yamada, ApJS (2012)fν (r,θ, φ; εν ,θν , φν ; t)Boltzmann eq.

€€

1 ∂fν  
1 ' δfν *+ n ⋅ ∇fν = ) ,c ∂tc ( δt + collision

2:?HA• 衝突項は複雑

+-K-L$-IνJ- =BE>HG& - /.,"- !
→%'• 近似法との比較が可能- 2D/3D: ), Ray-by-Ray0/9;513nnn- Ray-by-Ray68(<#743→ ν~20% 8*/

Sumiyoshi et al. ApJS (2015) 13

#+"CJOALD, )&<OB
Nagakura et alSApJS (2014, 2015)• 6DFM?HNR2D

R2D/G@N=IM– #+TB>EKO2(2moving mesh965CJOALD% P37'!Q;KOTYe, $ T.

;KOTν, $ Tν 6

50km

1 50 6CJOALD :*-548 cf. Ray-by-ray
14

T.
• 0

'LPFMQ/,)6*:;87, CAJBQD, +, "/, #&V

Furusawa EOS

vsLattimer-Swesty EOS

3!3 $-=-<>@2- ?54%(6#&<VCQKORE<>@USNGP9~400IRH1K-computer

Nagakura et al. ApJ (2018)15€

３次元でのニュートリノ輻射輸送 '?
Sumiyoshi & Yamada, ApJS (2012)

C>&LOGMR".3 2012, 10!1− µν2 cos φν ∂1− µν2 sin φν ∂fν
1 ∂fν µν ∂ 2+ 2 (r f ν ) +(sin θf ν ) +c ∂t r ∂rr sin θ∂θr sin θ∂φ
21−µ1 ∂1 ( δf ν +2ν cos θ ∂+[(1− µν ) f ν ] +(sin φν fν ) = * r ∂µνr sin θ∂φνc ) δt , collision
6D Boltzmann eq.

– T =>1; (Sn#)– 6A7=>81#

En,µnE´n,µn´nµν = cosθν€
1 $ δf ν '1= j emission (1− f ν ) −f ν + Cscattering [ ∫ f ν ( E ν- , µν- )dE ν- ]& )c % δt ( collision
λabsoption€– /+9T>BνRR5)•  BEJQFSR0<@A*2,• $)AD-RINSHPK %:(416• Linear equationAf ν = d• Neutrino distribution€Nn

Nn x Nspace~ 109– Nspace=nr x nq x nf– Nn=ne x nqn x nfn

Nvector~106 x 103• Memory size
ν-distribution >10GBmatrix: >1TB

• Iterative method– Pre-conditioner
Imakura et al. JSIAM (2012)

Kotake et al. PTEP (2012)17 
• Matrix: memory Nspace Ne x M2
operation Nspace Ne x M3
Space: Nspace, Neutrino: Nn=NeM
Parallel (MPI+OpenMP) by space r, q, f

MxM
M: block size
Angle variables

Sumiyoshi et al. ApJS (2011), Kotake et al. PTEP (2012)• 2D : 5D Boltzmann + 2D Hydrodynamics– Nspace=384x128, Nν=20x6x10 @ K-computer– 5M node-hour for ~0.5s evolution to find explosion– 1536 node, ~0.5 year
http://www.aics.riken.jp- Need Exa-flops machines for full 6D+3D simulations
18

!$( %#"&(' fn"&(• 1.0Core bouncef(Eν)ρc=100.613ρc=100.4ne12ρc=100.2113ρc=10 g/cm

Nene=1410100Eν [MeV]1.0 0.8"&(0.8nenµfν

0.60.4ne0.2
0.0-1.0-0.5

Nang=60.0µ = cosθ

0.5

1.021<µ>, <µ >0.0

– –   –
2LEQνi ~ Q0 ν ν
<µ>

0.8
14



0.60.4

<µ2>

0.2

1.0

0.0
5
10

<µ>
610

Nang=1210

7108

radius [cm]19:  • : Bi-CGSTAB
parameter optimized damped Jacobi-type preconditioner
Imakura et al. JSIAM (2012)

Boltzmann eq.
Solver in 3DTime step<10-7 sec→>10-5 sec

Optimized
Point-Jacobi

Improved forLong time evolution20)# •  &!$' %" • :  • ( • Super-Kamiokande• 爆発, 原始中性子星 n• 不発, ブラックホールnSN1987A

Disappearance of star!$' %"

Neutron starBlack holeAdams

N6946-BH1

https://www.spacetelescope.org

10Msun20Msun30Msun40Msun50Msun  21SNT• < *-fYa\eü'Nagk_hb)PUW*-• agk_hb?>BKE8
üci]dj1XJm@&N&• l%ci]dj1PUW;3K6
ü=/MQagk_hb?';3nü2 +;3R• FI"P(W;3Zk`QG-N4ü 9$7:&, ;3.NQC
üO^[_,H0#N!D Q5VAL22

Project in collaboration with• Numerical simulations• Supercomputing– H. Matsufuru, A. Imakura– H. Nagakura• EOS tables & neutrino rates– W. Iwakami– S. Furusawa,– H. Okawa– H. Shen, K. Oyamatsu, H. Toki
– A. Harada– S. X. Nakamura, T. Sato
– S. Yamada• Neutron star merger• Supernova research– Y. Sekiguchi, S. Fujibayashi– T. Takiwaki
Supported by– K. Nakazato- MEXT and JICFuS– K. Kotake- for K-computer and Post-K machine– K. Takahashi-K-Computer, Japan
http://www.aics.riken.jp

K-computer: hp180111, hp180179, hp180239
HPC resources at KEK, YITP, UT, RCNP

Grant-in-Aid for Scientific Research(15K05093, 17H06357, 17H06365)

23