Introduction a la simulation Monte Carlo avec PRIMO

Qu'est-ce que PRIMO ?

PRIMO est un logiciel gratuit de simulation Monte Carlo specialement concu pour la dosimetrie des accelerateurs lineaires medicaux (LINAC). Developpe par l'equipe du Dr. Lorenzo Brualla a l'Universite de Duisburg-Essen en Allemagne, PRIMO rend accessible la simulation Monte Carlo aux physiciens medicaux et aux etudiants sans necessiter de competences avancees en programmation.

Contrairement aux codes Monte Carlo generalistes comme GEANT4 ou EGSnrc qui exigent une maitrise de la programmation et de longues phases de configuration, PRIMO fournit une interface graphique intuitive permettant de configurer, d'executer et d'analyser une simulation complete en quelques etapes. Le logiciel est construit autour du moteur Monte Carlo PENELOPE (PENetration and Energy LOss of Positrons and Electrons), un code reconnu pour sa precision dans le transport des photons et des electrons a basse et moyenne energie.

Le moteur PENELOPE

PENELOPE est un code Monte Carlo developpee par l'Universitat de Barcelona sous l'egide de la Nuclear Energy Agency (NEA/OECD). Il simule le transport couple d'electrons, de positrons et de photons dans des materiaux arbitraires.

Physique simulee

PENELOPE modelise les interactions suivantes avec une grande rigueur :

Pour les photons :

• Diffusion Rayleigh (diffusion coherente)
• Diffusion Compton (diffusion incoherente)
• Effet photoelectrique
• Production de paires electron-positron

Pour les electrons et les positrons :

• Diffusion elastique (modele de section efficace partielle ondes)
• Collisions inelastiques (ionisation et excitation)
• Rayonnement de freinage (bremsstrahlung)
• Annihilation des positrons

Algorithme de transport mixte

PENELOPE utilise un algorithme de transport dit "mixte" pour les particules chargees. Les interactions sont classifiees en deux categories :

• Interactions dures (hard events) : les interactions avec un transfert d'energie ou une deviation angulaire superieur a un seuil defini par l'utilisateur. Ces interactions sont simulees individuellement, une par une.
• Interactions molles (soft events) : les interactions avec de faibles transferts, regroupees et traitees de maniere condensee entre deux interactions dures a l'aide de distributions de diffusion multiple.

$\text{Parcours total} = \sum_{i=1}^{N_{hard}} s_i$

ou $s_i$ est le libre parcours entre deux interactions dures consecutives, et les interactions molles sont traitees de facon condensee le long de chaque segment $s_i$.

Les parametres de transport $C_1$, $C_2$, $W_{cc}$ et $W_{cr}$ controlent la frontiere entre interactions dures et molles :

• $C_1$ : deflexion angulaire moyenne maximale entre deux interactions dures (typiquement 0,05 a 0,1)
• $C_2$ : perte d'energie fractionnelle moyenne maximale entre deux interactions dures (typiquement 0,05 a 0,1)
• $W_{cc}$ : energie de coupure pour les collisions inelastiques dures (en eV)
• $W_{cr}$ : energie de coupure pour les emissions de bremsstrahlung dures (en eV)

Architecture d'une simulation PRIMO

PRIMO decompose la simulation en trois segments independants, ce qui permet de reutiliser les resultats intermediaires et d'optimiser le temps de calcul.

Segment 1 : la tete de l'accelerateur (partie superieure)

Le premier segment simule le transport des particules depuis la cible jusqu'a la sortie du systeme de collimation primaire. Il inclut :

• La cible en tungstene ou en cuivre ou les photons sont produits par bremsstrahlung
• Le collimateur primaire qui delimite le cone utile
• Le filtre egalisateur (flattening filter) qui homogeneise la fluence
• La chambre moniteur qui mesure le nombre d'unites moniteur

PRIMO propose des modeles predefinis pour les accelerateurs les plus repandus (Varian Clinac, Elekta Synergy, Siemens Oncor, etc.). L'utilisateur selectionne le modele et l'energie souhaitee, et PRIMO genere automatiquement la geometrie correspondante.

Le resultat de ce segment est un fichier d'espace de phase (phase space file, PSF) qui enregistre les proprietes (position, direction, energie, type) de chaque particule traversant un plan situe au-dessus du systeme de collimation secondaire.

Segment 2 : la collimation secondaire

Le deuxieme segment transporte les particules du PSF a travers le systeme de collimation secondaire :

• Machoires (jaws) definissant la taille du champ rectangulaire
• Collimateur multilame (MLC) pour les champs conformes et modules

L'utilisateur definit les positions des machoires et des lames du MLC pour obtenir la forme de champ desiree. PRIMO supporte les principaux modeles de MLC (Varian Millennium 120, Varian HD120, Elekta MLCi, etc.).

Le resultat est un second PSF enregistre a la sortie du MLC, pret a etre utilise pour le transport dans le fantome.

Segment 3 : le fantome et le calcul de dose

Le troisieme segment transporte les particules dans un fantome d'eau (ou un autre materiau) et calcule la distribution de dose tridimensionnelle. L'utilisateur configure :

• Les dimensions du fantome
• La resolution des voxels (taille des elements de volume dans lesquels la dose est calculee)
• La distance source-surface (SSD)

La dose absorbee dans chaque voxel est calculee en sommant l'energie deposee par toutes les particules traversant le voxel, divisee par la masse du voxel :

$D_{voxel} = \frac{\sum_{i=1}^{N} \Delta E_i}{\rho \times V_{voxel}}$

ou $\Delta E_i$ est l'energie deposee par la particule $i$, $\rho$ la masse volumique du milieu, et $V_{voxel}$ le volume du voxel.

Configurer une simulation : guide pratique

Etape 1 : selection de l'accelerateur

Au lancement de PRIMO, l'utilisateur selectionne le modele d'accelerateur lineaire parmi la bibliotheque disponible. Chaque modele a ete valide par les developpeurs contre des mesures experimentales publiees.

Pour un resultat fidele, il est important de connaitre les specifications exactes de l'accelerateur clinique : modele, energie nominale, type de MLC, et si possible les parametres de la source d'electrons (taille focale, energie moyenne, divergence).

Etape 2 : ajustement de la source d'electrons

Les parametres de la source d'electrons primaires influencent fortement la distribution de dose resultante. Les principaux parametres a ajuster sont :

• L'energie moyenne : typiquement proche de l'energie nominale du faisceau (ex. 6,2 MeV pour un faisceau nominal de 6 MV)
• La largeur du spectre en energie : modele gaussien ou uniforme
• La taille focale : full width at half maximum (FWHM) de la distribution spatiale de la source, typiquement 1 a 3 mm

Etape 3 : definition du champ et du fantome

L'utilisateur definit ensuite :

• La taille du champ (ouverture des machoires et position des lames MLC)
• Les dimensions du fantome d'eau (typiquement 40 x 40 x 40 cm^3 pour la dosimetrie de reference)
• La resolution des voxels (typiquement 2 mm dans le plan transverse et 2 mm en profondeur)
• La SSD (generalement 100 cm)

Etape 4 : parametres de simulation

Les parametres critiques incluent :

• Le nombre de particules primaires : determine l'incertitude statistique du resultat. Pour une incertitude de 1 % dans les voxels du champ, il faut typiquement 10^8 a 10^9 particules.
• Les energies de coupure : en dessous desquelles les particules sont arretees et leur energie restante deposee localement. Des valeurs typiques sont 50 keV pour les electrons et 10 keV pour les photons.
• La technique de reduction de variance : PRIMO propose le splitting, qui consiste a multiplier les particules dans les regions d'interet pour ameliorer la statistique sans augmenter le nombre de particules primaires.

Etape 5 : execution et analyse

L'execution de la simulation peut durer de quelques heures a plusieurs jours selon le nombre de particules, la resolution et la complexite de la geometrie. PRIMO supporte l'execution parallele sur plusieurs coeurs de processeur.

Une fois la simulation terminee, PRIMO affiche :

• Les courbes de PDD le long de l'axe central
• Les profils lateraux a differentes profondeurs
• Les distributions de dose 2D et 3D
• L'incertitude statistique dans chaque voxel

Comparaison avec les mesures experimentales

La validation d'une simulation Monte Carlo passe obligatoirement par la comparaison avec des mesures experimentales. Les criteres de validation typiques sont :

L'analyse gamma est l'outil standard pour comparer quantitativement deux distributions de dose. Elle combine un critere de difference de dose et un critere de distance :

$\gamma(\vec{r}) = \min \left\{ \sqrt{\frac{|\vec{r} - \vec{r}_m|^2}{\Delta d^2} + \frac{|D(\vec{r}) - D_m(\vec{r}_m)|^2}{\Delta D^2}} \right\}$

ou $\Delta d$ est le critere de distance (typiquement 2 ou 3 mm), $\Delta D$ le critere de dose (typiquement 2 ou 3 %), et l'indice $m$ designe la distribution mesuree. Un point est considere comme "passant" si $\gamma \leq 1$.

Les criteres d'acceptation standards en dosimetrie clinique sont :

Avantages de la simulation Monte Carlo

La simulation Monte Carlo avec PRIMO offre des avantages considerables par rapport aux seules mesures experimentales :

• Acces a la dose en tout point : la simulation fournit la dose dans chaque voxel du fantome, y compris dans des regions difficiles a mesurer (interfaces, petits champs, zones de buildup).
• Decomposition des contributions : possibilite de separer les contributions des photons primaires, des photons diffuses et des electrons de contamination.
• Geometries complexes : possibilite d'etudier des configurations non standard (heterogeneites, bolus, compensateurs).
• Cout nul en consommables : pas de films, pas de chambre d'ionisation, pas de fantome physique necessaire.
• Reproductibilite parfaite : les resultats sont entierement reproductibles pour un meme jeu de parametres et de graines aleatoires.

Lien avec DosiPlot

Les donnees de reference integrees dans DosiPlot ont ete generees avec PRIMO pour des configurations cliniques standard. Cette approche permet aux utilisateurs de DosiPlot d'acceder a des distributions de dose de haute precision, validees par simulation Monte Carlo, sans avoir a installer ni a configurer PRIMO.

Pour les utilisateurs souhaitant aller plus loin, la comparaison entre les donnees DosiPlot et les resultats de leurs propres simulations PRIMO constitue un excellent exercice d'apprentissage et de validation.

FAQ

PRIMO est-il vraiment gratuit ?

Oui, PRIMO est distribue gratuitement a des fins de recherche et d'enseignement. Il peut etre telecharge depuis le site officiel apres inscription. Le logiciel fonctionne sous Windows et ne necessite pas de licence payante.

Combien de temps dure une simulation typique ?

Le temps d'execution depend du nombre de particules, de la resolution des voxels et de la puissance de calcul disponible. Pour une simulation de reference (champ 10 x 10 cm^2, fantome d'eau, resolution 2 mm, 10^8 particules), il faut compter environ 2 a 8 heures sur un PC de bureau moderne equipe d'un processeur multi-coeurs. L'utilisation du splitting peut reduire ce temps de maniere significative.

Peut-on simuler des plans de traitement complets avec PRIMO ?

PRIMO est principalement concu pour la dosimetrie de reference et le commissioning des faisceaux. Il ne remplace pas un TPS pour la planification de traitement. Cependant, il peut simuler des champs individuels avec des configurations de MLC definies, ce qui le rend utile pour la verification de plans simples et la validation de l'algorithme du TPS.

Quelle precision peut-on atteindre avec PRIMO ?

Avec un modele d'accelerateur correctement ajuste et un nombre suffisant de particules, PRIMO permet d'atteindre un accord avec les mesures de l'ordre de 1 % pour les PDD et de 1 a 2 mm pour la position des isodoses dans les profils lateraux. Ces performances sont comparables a celles des autres codes Monte Carlo de reference (EGSnrc, GEANT4).