Comparison of heavy-ion transport simulations; Collision integral with pions and resonances in a box

重イオン輸送シミュレーションの比較; パイオンとデルタ共鳴を含む箱内の衝突積算

小野 章*; Xu, J.*; Colonna, M.*; Danielewicz, P.*; Ko, C. M.*; Tsang, M. B.*; Wang, Y,-J.*; Wolter, H.*; Zhang, Y.-X.*; Chen, L.-W.*; Cozma, D.*; Elfner, H.*; Feng, Z.-Q.*; 池野 なつ美*; Li, B.-A.*; Mallik, S.*; 奈良 寧*; 小川 達彦   ; 大西 明*; Oliinychenko, D.*; Su, J.*; Song, T.*; Zhang, F.-S.*; Zhang, Z.*

Ono, Akira*; Xu, J.*; Colonna, M.*; Danielewicz, P.*; Ko, C. M.*; Tsang, M. B.*; Wang, Y,-J.*; Wolter, H.*; Zhang, Y.-X.*; Chen, L.-W.*; Cozma, D.*; Elfner, H.*; Feng, Z.-Q.*; Ikeno, Natsumi*; Li, B.-A.*; Mallik, S.*; Nara, Yasushi*; Ogawa, Tatsuhiko; Onishi, Akira*; Oliinychenko, D.*; Su, J.*; Song, T.*; Zhang, F.-S.*; Zhang, Z.*

2017年4月に開催された国際会議Transport2017において、重イオン核反応モデルの国際的な比較が議論された。重イオン加速器の安全評価や宇宙飛行士の被ばく評価等で重要な役割を果たすため、世界中で重イオン核反応の様々な理論モデルが開発されている。本研究では、辺の長さが20fmの直方体に320個の中性子と陽子をランダム配置し、それらが70fm/cの間に起こす散乱の回数やエネルギーを計算した。ここでは、特にパイオンやその前駆体であるデルタ共鳴の生成に注目して比較を行った。参加コードは、個々の粒子の時間発展を追うQMD型コードと、粒子の位置や運動量の確率分布を決めておき、散乱や崩壊が発生したときそれらを乱数サンプリングするBUU型コードがあり、発表者が用いたJQMDは前者に属する。本研究により、計算における時間刻みが各コードによる結果の差の主な原因であることが分かった。さらに、今後のJQMDの改良方針の策定に有益な知見を得ることができた。

International comparison of heavy-ion induced reaction models were discussed in the international conference "Transport2017" held in April 2017. Owing to their importance for safety assessment of heavy-ion accelerators and dosimetry of astronauts, various models to simulate heavy-ion induced reaction models are developed. This study is intended to clarify the difference among them to pinpoint their problems. In the comparison study, 320 protons and neutrons were packed in a 20-fm-large cube to calculate the number and energies of collisions during the time evolution. The author contributed to this study by running calculation using JQMD (JAERI Quantum Molecular Dynamics). This study showed that time step in the calculation is one of the biggest causes of the discrepancies. For example, the calculation by JQMD comprises 1-fm/c time steps, each of which is composed of transport, scattering and decay phases. Therefore a sequence of scattering, and decay followed by another scattering in 1 fm/c cannot be considered. Moreover, in JQMD particles are labeled by sequential numbers and scattering reactions are simulated by the order. Therefore scattering between low ID numbers, that between high ID numbers and that between the first (low ID) pair is overlooked in JQMD. Above indications obtained in this study must be kept in our mind for future JQMD upgrades.