I would like to simulate energy deposition as a result of the stopping of high energy particles in a dense plasma. For example the energy deposited when 10MeV protons are stopped in a hydrogen plasma of density 300g/cc. I think combining a Particle in Cell code with a monte carlo code would give the appropriate physics. However I have not found a code available to do this kind of simulation. Does anyone know of any appropriate simulation tools for this problem?