Medical Applications of Monte Carlo Methods with SuperComputing

1
Genova, 8 Marzo 2004
Radiotherapy treatment planning with Monte Carlo
on a distributed system
Stéphane Chauvie,
IRCC Mauriziano Hospital INFN S Croce e
Carle Hospital Turin, Italy
2
Contents

• Radiotherapy Treatment Planning
• Analitical algorithms for dose calculation
• Monte Carlo methods
• Cluster set-up
• Monte Carlo parallelization
• Data analisys and experimental measurements
  comparison open field and IM field
• Head and neck tumor with IMRT

Grant 2002-03/645
3
Radiotherapy Oncology
spare the surrounding healthy tissues.
deliver high dose to the target volume
allow local control of tumor
avoid side-effects
4
3D-CRT vs IMRT
Critical points - high dose gradients - strongly
unhomogeneous areas
3DCRT IMRT used in complex anatomical regions
How much is accurate the dose calculation ?
5
Dose calculation algorithms

• Pencil beam
• Convolution/Superposition

• Monte Carlo

Accurate but very slow Cheap (free)
Expensive Quick but inaccurate
Meas in ref pomint, neam stability flatness, CT
data, setup
6
Cluster Beowulf parallelisation


PC Ethernet
Th High performance networks of PCs are now
realistic alternative since offer parallel
processing of MC at a lower cost showing
competitive performances.
7
Cluster set-up
Hardware installation
Monte Carlo simulation
Software configuration
Benchmarking
Monte Carlo parallelisation
RUN
8
Installation, configuration benchmarking
Bios OS Disk conf Partition RAID Memory CPU Compil
ators Linking models
Parallelization LAM-MPI
Security SSH
9
Installation, configuration benchmarking
Efficiency Sup/ Nprocessors 0.997
10
Simulation geometry
?V 6 MV
e-
Varian 600C/D Millenium 120-leaf MLC
11
simulation physics

• Geant4 has only production thresholds, no
  tracking cuts
• all particles are tracked down to zero range
• energy, TOF ... cuts can be defined by the user

NO TUNING, NO CUT
12
Patient model
DICOM interface
Soft tissue - CT-tissue relationship ICRU
Bone - CT-?el linearity - cortical bone - bone
marrow diluition
Lung - CT-? linearity
13
Monte Carlo Parallelization
Take care of PRNG
IM simple field in homogeneous phantom
Phase Space Data
Water measurements
IM patient field in homogeneous phantom
Anthropomorphic phantom measurements
Simulation inside patient
IMRT treatment
14
Phase Space Data
PSD
(x,y,z) (px,py,pz)
15
Water measurements
PDD and dose profile in water
Scanner IC15 ionization chamber SSDSAD
16
Anthropomorphic phantom measurements
Microchamber A14SL SSDSAD
17
Patient simulation
TAC
X10 Y10 SSDSAD Gantry 0
18
IMRT treatment simulation
E0.9925
10X10 isocentric technique 7 field!
Every field segments no. 165,4?15,3 events
no. (15,5?0,5)107 hits no. (4,02 ? 0,39) 105
time (hours) 0,51?0,03
IMRT plan evaluation in 3,5 hours with 280000
hits and 3 nodes
19
Current Geant4 activities in Cuneo...
10 MeV Cyclotron
CT-PET