The advent of general-purpose programmable GPUs has prompted the development of MC tools that can deliver a plan calculation within minutes [1]. The new hardware requires to substantially rewrite the core kernels of a MC calculation adopting the parallel execution model of a typical GPU. In the last three years a fast-MC tool called Fred has been developed at University of Rome "La Sapienza". The code can track proton pencil beams through a voxel 3D patient reconstruction grid based on high resolution CT scans. The kernel implements energy loss based on tabulated stopping power [2].
Different models (single Gaussian, double Gaussian, Gauss-Rutherford) can be used to reproduce Multiple Coulomb Scattering.  Secondary protons and deuterons are produced in the beam nuclear interaction following the multiplicity, energy and angular distributions provided by ICRU63 [3], and interpolated for other materials.

The voxel HU of the CT is converted in elemental composition following the method by Schneider[4]. Alpha and heavier fragments produced are treated as local dose deposition and (in this version) neutral particles are not produced. The proton RBE value is kept constant for direct comparison with commercial TPS results, but it can be also loaded from an external library of RBE values. The dose matrix can be processed by an optimizer based on a least-squares optimization algorithm as in [5].
The output/input system has been adapted to the CNAO environment[6]. All these solutions have been checked for efficacy and accuracy against a full-MC code[7], and the resulting kernel can deliver more than 1 million complete histories per second on a single GPU. The achieved tracking rate allows to perform a typical plan recalculation within minutes. 




[1] X. Jia et al., Phys. Med. Biol. (2012) 57:7783–7797; X. Jia et al., Phys. Med. Biol. (2014) 59:R151–R182; D. Giantsoudi et al., Phys. Med. Biol. (2015) 60:2257–2269 

[2] M.J.Berger et al, PSTAR, version 1.2.3 (2005). Available online: http://physics.nist.gov/Star

[3] Schneider U et al., Phys. Med. Biol. (1996) 41:111–24

[4] ICRU Report 63 (2000)

[5] A. Mairani et al., 2013 Phys Med Bio . 58, 2471-2490; Lomax et al 1999  Phys. Med. Bio. 44 185-205

[6] S. Molinelli et al., Phys. Med. Biol. (2013) 58:3837-47

[7] A. Ferrari et al., CERN-2005-10 (2005), INFN/TC_05/11, SLAC-R-773