Particle-in-cell simulation of transport and energy deposition of intense proton beams in solid-state materials

A particle-in-cell (PIC) simulation code is used to investigate the transport and energy deposition of an intense proton beam in solid-state material. This code is able to simulate close particle interactions by using a Monte Carlo binary collision model. Such a model takes into account all related interactions between the incident protons and material particles, e.g., proton-nucleus, proton-bound-electron, and proton-free-electron collisions. This code also includes a Monte Carlo model for the collisional ionization and electron-ion recombination as well as the depression of the ionization potential by shielding of surrounding particles. Moreover, for intense proton beams, in order to include collective electromagnetic effects, significantly speed up the simulation, and simultaneously avoid numerical instabilities, an approach that combines the PIC method with a reduced model of high-density plasma based on Ohm's law is used. Simulation results indicate that the collective electromagnetic effects have a significant influence on the transport and energy deposition of proton beams. The Ohmic electric field would increase the stopping power and leads to a shortened range of proton beams in solid. The magnetic field would localize the energy deposition by collimating proton beams, which would otherwise be deflected by the collisions with nuclei.