@article{oai:nifs-repository.repo.nii.ac.jp:00010620,
 author = {SATAKE,  Shinsuke and Nakata, Motoki and Pianpanit, Theerasarn and SUGAMA,  Hideo and NUNAMI,  Masanori and MATSUOKA,  Seikichi and Ishiguro,  Seiji and KANNO,  Ryutaro},
 journal = {Computer Physics Communications},
 month = {Mar},
 note = {0000-0002-9387-9987, A numerical method to implement a linearized Coulomb collision operator in the two-weight Monte Carlo method for multi-ion-species neoclassical transport simulation is developed. The conservation properties and the self-adjoint property of the operator in the collisions between two particle species with different temperatures are verified. The linearized operator in a Monte Carlo code is benchmarked with other two kinetic simulations, a continuum gyrokinetic code with the same linearized collision operator and a full-f PIC code with Nanbu collision operator. The benchmark simulations of the equilibration process of plasma flow and temperature fluctuation among several particle species show very good agreement between Monte Carlo code and the other two codes. An error in the H-theorem in the two-weight Monte Carlo method is found, which is caused by the weight spreading phenomenon inherent in the two-weight method. It is demonstrated that the weight averaging method serves to restoring the H-theorem without causing side effect.},
 title = {Benchmark of a new multi-ion-species collision operator for δf Monte Carlo neoclassical simulation},
 volume = {252},
 year = {2020}
}