Electromagnetic wave-plasma interaction has drawn much attention recently due to numerous important technologies and applications, taking advantage of phenomena such as electromagnetic waves being reflected or absorbed in a plasma medium. The physics of wave-plasma interaction can be complicated, when non-uniform, non-equilibrium, or anisotropic plasmas are involved, in which numerical simulations can be used to fill the gaps between theoretical solutions and experimental measurements. Among many numerical methods, the particle-in-cell method, which can solve accurately both the electromagnetic fields and particle trajectories
self-consistently, would be the best choice to study wave-plasma interaction problems as long as the computational cost can be accepted. However, the applications of PIC on wave-plasma interaction remain rare, and the numerical effects of the PIC method on accurately evaluating the wave attenuation have not been studied in depth. In this paper, a number of numerical parameters and physical parameters are tested using a 1D electromagnetic PIC method plus Monte Carlo collision model. It is found that as long the as the basic PIC criterion is met, the PIC results can be trustable, and the numerical noise due to limited number of particles has a minor effect. The physical parameters of the EM wave frequency, amplitude, the plasma temperature, thickness, and collision type are studied, and their effects on the wave attenuation are presented. In addition, strategies on establishing simulation setup and evaluating the wave attenuation in terms of power or energy are discussed.