Monte Carlo simulation with Geant4 for verification of rotational total skin electron therapy (TSET)

Detta arbete handlar om en behandlingsteknik för hudterapi där lågenergetiska elektroner utnyttjas för strålbehandling där patienten står eller sitter på en roterande plattform. För att behandlingen skall få ett lyckat resultat krävs en viss stråldos (den absorberade energin per massenhet) till patientens hud. För att kontrollera att dosen är rätt används normalt mätningar med olika detektorer. Dessa detektorer fungerar optimalt när de är riktade mot strålkällan men har en obestämd känslighet när de roteras. För att uppskatta den totala dosen under en rotation användes i stället ett datorsimuleringsprogram (Geant4). Programmet Geant4, för simulering av transport av partiklar i materia, utnyttjas av ett brett spektrum av användare inom högenergifysiken, rymdfysiken och för medicinska ändamål. Programmet har nyligen uppdaterats med nya rutiner för elektromagnetisk strålning (jonisation, bromsstrålning och multipelspridning) vilket har ökat precisionen av simuleringar för rymd- och medicinska applikationer. Ett stort antal användare i Europa utnyttjar nu programmet för beräkningar av lågenergetiska elektromagnetiska processer och i denna studie har vi använt programmet för att simulera behandlingstekniken som nämndes ovan. Resultaten av denna undersökning visar att Geant4 programmet ger god överensstämmelse i en stillastående referensgeometri när man gämnför med mätningar och med ett annat simuleringsprogram (VMC++). Resultaten visar också att Geant4 kan optimeras för beräkningar även vid mycket låga energier, och att programmet kan användas för att bestämma dosfördelningen när patienten roterar.Background and purpose: This study refers to a rotational radiotherapy technique for the treatment of mycosis fungoides. Patients are irradiated by a high dose-rate low-energy electron beam while they stand on a rotating platform. This report presents a first effort to use the Geant4 simulation program to calculate the relation between the stationary and rotational depth dose distributions in an ellipsoid phantom along the central beam axis. The absorbed dose delivered in one rotation can then be inferred from dose rate measurements performed in a stationary geometry. Method: Percentage depth dose and ionisation chamber measurements along the central beam axis were done at SSD = 100 cm and 250 cm (treatment distance). The measurement at SSD = 100 cm was used 1. to extract the beam characterization using VMC++ calculated percentage depth dose distributions and 2. to select the physics model of Geant4. A set of simulations were subsequently performed at the treatment distance in order to determine the following parameters in the simulation: 1. the number of primary electrons that are needed for the calculation of the ratio between the stationary and rotational depth dose distributions, R, to have a given level of statistical uncertainty (2%), 2. the distance off the central beam axis within which the energy should be scored, r, and 3. the step in beam incidence angle, . The rotational depth dose and R distributions were then calculated for one elliptical and one cylindrical phantom. The cylindrical phantom calculations were done for two production thresholds for the electrons, 10 keV and 100 keV. Results: This simulation study has contributed to the theoretical background of a rotational technique for total skin electron therapy (TSET). It was found that the low production threshold simulations give results that agree better with the measurements and the VMC++ calculation at SSD = 100 cm and also with the expected shape of the depth dose distribution for one rotation. Conclusion: The Monte Carlo calculations presented in this report were used to determine the the relation between the stationary and rotational depth dose distributions and thus contributed to a more accurate evaluation of the absorbed dose delivered at the prescription depth during treatment