Solutions were developed to solve the special photon transport problems. A respective Monte Carlo code were implemented. The photon transport calculations were made for simulation of light distribution in tissues, for 3D estimation of dosis brachytherapy source, for determination of scattering coincidences in PET and for solving of one special problem for ambient dosimetry. The developed calculation methods are based on a purely statistical approach and can therefore universal applied. The efficiency in respect of precision, numerical effectiveness, as well as memory requirement were optimised and verified by the calculations done. The modification and enhancement can be easy realised thanks to modular, object orientated implementation of the program and enables development of new application fields in physics and medicine. / Zur Lösung speziellen Probleme des Photonentransports wurden grundlegende Ansätze gewonnen und in einem Transportprogramm umgesetzt. Die Photonen-Transportrechnungen zur Simulation des Lichttransports in trüben Medien, für die dreidimensionale Dosisberechnung für interstitielle Brachytherapiequellen, für die Ermittlung der Streustrahlungsverteilungen in PET-Scannern und für die Lösung eines speziellen Problems der Umgebungsdosimetrie wurden durchgeführt. Die entwickelten Berechnungsmethoden basieren auf einer rein probabilistischen Herangehensweise und lassen sich deshalb universell verwenden. Ihre Leistungsfähigkeit hinsichtlich der Genauigkeit, der numerischen Effektivität sowie des Bedarfs an Rechenresourcen wurde optimiert und konnte durch die durchgeführten Berechnungen bestätigt werden. Durch den modularen, objektorientierten Programmaufbau sind Modifikationen und Erweiterungen relativ einfach durchzuführen. Das ermöglicht eine Erschließung von neuen Anwendungsgebieten in der Physik und Medizin.

Ce travail de thèse répond à un besoin de qualification des outils permettant de calculer les échauffements photoniques dans les réacteurs nucléaires. En effet, la problématique des échauffements g dans les matériaux de structure a pris de l’importance ces dernières années, notamment pour la sûreté des réacteurs de 3ème génération dans lesquels un réflecteur lourd en acier est introduit.Les photons présents dans le coeur sont tous directement ou indirectement issus des interactions des neutrons avec la matière. Ils sont créés au moment de l’interaction ou en différé par l’intermédiaire de noyaux créés au moment de l’interaction. Par conséquent, le premier axe de travail a été d’effectuer une analyse critique des données de production photonique dans les bibliothèques de données nucléaires standards. La découverte d’omissions dans la bibliothèque JEFF-3.1.1 nous a amené à proposer une méthode de production de nouvelles évaluations contenant de nouveaux spectres d’émission de photong. Ces nouvelles évaluations ont ensuite été proposées et en partie acceptées pour la nouvelle version de la bibliothèque JEFF.Il existe deux codes de transport de particules développés au CEA : TRIPOLI4 etAPOLLO2. Le deuxième axe de travail a été de qualifier ces deux codes. Pour cela, nous avons interprété les mesures d’échauffement g effectuées dans le cadre du programme expérimental PERLE. Des détecteurs thermoluminescents (TLD) ont été introduits dans un réflecteur lourd en acier entourant un réseau de crayons combustibles. Nous avons dû proposer un schéma de calcul spécifique aux deux codes afin de calculer la réponse des TLD.Les comparaisons calcul-mesure ont montré que TRIPOLI4 permettait decorrectement estimer l’échauffement dans le réflecteur relativement à l’échauffement dans lazone fissile. En effet, les écarts calcul-mesure sont inférieurs à l’incertitude expérimentale à1s. Pour le calcul APOLLO2, nous avons tout d’abord commencé par une phase de validation par rapport à TRIPOLI4 afin d’estimer les biais liés aux approximations imposées par le traitement déterministe du transport des particules. Après cette phase de validation,nous avons pu montrer qu’APOLLO2, comme TRIPOLI4, permettait d’estimer correctement l’échauffement dans le réflecteur avec des écarts calcul-mesure comparables à l’incertitude expérimentale. / During the last few years, the g-heating issue has gained in stature, mainly for thesafety of the 3rd-generation reactors in which a stainless steel reflector is inserted. Thepurpose of this work is the qualification of the needed tools for calculation of the g-heating inthe nuclear reactors.In a nuclear reactor, all the photons are directly or indirectly produced by the neutronmatterinteractions. Thus, the first phase of this work is a critical analysis of the photonproduction data in the standard nuclear data library. New evaluations have been proposed tothe next version of the JEFF library after that some omissions have been found. They havepartly been accepted for JEFF-3.2.Two particle-transport codes are currently developed in the CEA: the deterministiccode APOLLO2 and the Monte Carlo code TRIPOLI4. The second part of this work is thequalification of both these codes by interpreting an integral experiment called PERLE. Theexperimental set-up is made by a LWR pin assembly surrounded by a stainless steelreflector in which the g heating is measured by Thermo-luminescent Detector (TLD). Acalculation scheme has been proposed for both APOLLO2 and TRIPOLI4 in order tocalculate the TLD’s responses.Comparisons between calculations and measurements have shown that TRIPOLI4gives a satisfactory estimation of the g heating in the reflector. These discrepancies arewithin the experimental 1s uncertainty. Before the qualification, APOLLO2 has beenpreviously validated against TRIPOLI4 reference calculation. This validation gives anestimation of the bias due to the deterministic approximations of the transport equationresolution. The qualification has shown that the discrepancies between APOLLO2predictions and TLD’s measurements are in the same range as experimental uncertainties.

The general context of this work is the use of low energy photon sealed sources in brachytherapy. We have worked in particular on two isotopes: I-125 (mean energy of 27 keV) and Pd-103 (mean energy of 21 keV). The sealed sources are prepared as cylindrical seeds 4.5 mm in length and 0.8 mm in radius. Even if the external dimensions are standard, the internal design can be extremely different from one model to the other as the manufacturers try to improve the dosimetric characteristics of their sources. These isotopes are used mainly as permanent implants for prostate tumours but can also be used in the treatment of eye tumours. Compared to the higher energy photon sources, they offer the physical advantages of a safer manipulation from a radioprotection point of view and of the reduction of the dose to the surrounding healthy tissues. When performing a clinical treatment, it is absolutely mandatory to be able to report very precisely various parameters that can have an impact on the patient treatment outcomes. These parameters are, for example, the prescribed dose, the doses at different organs, the degree of uniformity that has been achieved on the target or some dose-volume information. The brachytherapy treatment planning systems (TPS) also permit more and more to conform the treatment to the patient anatomy, like in external treatments. In the case of prostate tumours, it has been possible for a few years, using ultrasound imaging, to check the positioning of the seeds and to calculate the dose distribution in real time during the implantation procedure. It is clear that to achieve a good precision in the treatment itself, the sources dosimetric characteristics must be known with a maximum of precision. As at these low energies the dosimetric characteristics are very dependant on the internal design, this implies a thorough study of any new source design. This is the subject of the first part of this work with the determination of the dosimetric characteristics of two new models of seed produced by the IBt Company. That determination has been done using measurements with thermoluminescent detectors and Monte Carlo calculations with two codes: MCNP4C and EGSnrc. The comparison of these two codes with the measurements has brought into evidence the necessity to use up to date cross section libraries during the calculations and to take into account the presence of the detectors during the measurements. However, the dosimetry is only one part of the problem when dealing with these low energy sources in radiotherapy. The irradiated materials are complex living tissues, composed of many substructures on which radiations have not always the same impact. It is now largely accepted that the most radiosensitive part of a living cell is the DNA. When photons interact with matter, they produce electrons that deposit their energy as ionisations and excitations of the atoms. These ionisations, if they occur with or close to the DNA molecule, are responsible for damaging it. That structure can support a certain amount of damage and stay functional thanks to repair mechanisms, but these mechanisms have limits in the way they can handle multiple breaks in the DNA strands. If these breaks are too close in space and time, repair is not possible and the cell dies. The density of ionisations increases when the energy of the incoming photon decreases so it is expected that low energy photons will have a higher efficiency at killing cells and so a higher Relative Biological Effect than high energy photons. To study that subject, one must reduce the volume of matter considered during energy deposition events to reach the sizes relevant for cells, micron or even nanometre volumes. That is why the part of radiation physics dealing with that problem is called microdosimetry. The second part of this work is dedicated to a theoretical microdosimetric study of these low energy photon sources using the microdosimetric functions for volumes of 1µm and a cluster analysis study for nanometric volumes. The two types of studies show that the photons of 20-30 keV are effectively more biologically efficient than high energy photons. The microdosimetric results give a ratio of the relative to Co-60 of 2.6 for Pd-103 and 2.2 for I-125, which is concordant with the experimental values published by Wuu et al. (1996). The cluster analysis also shows that the electrons generated by photons of 20-30 keV produce more high order clusters than electrons of 300 keV. The mean cluster order for clusters of 10 nm is 3.0 for Pd-103 and is 3.3 for I-125 compared to 2.1 for electrons of 300 keV. In this case, I-125 shows a higher mean cluster order, which is related to a potentially higher biological efficiency. This is explained by the fact that I-125 photons interact with a probability of 51% by Compton Effect, when that probability decreases to 27% for Pd-103. Compton interactions generate a high number of very low electrons that deposit their energy very locally with a very high density of ionisations. More radiobiological studies are needed to determine which volume site is more relevant and therefore which kind of study better reflects the reality. / Le cadre général de ce travail est l'utilisation de sources scellées de photons de basse énergie en curiethérapie. Nous avons travaillé en particulier avec deux isotopes : I-125 (énergie moyenne de 27 keV) et le Pd-103 (énergie moyenne de 21 keV). Les sources scellées se présentent sous la forme de grains cylindriques de 4.5 mm de long et 0.8 mm de rayon. Bien que les dimensions externes soient standards, la configuration interne peut être extrêmement différente d'un modèle à l'autre vu que les fabricants essaient d'améliorer les caractéristiques dosimétriques de leurs sources. Ces isotopes sont utilisés comme implants permanents pour les tumeurs de la prostate mais peuvent aussi être utilisés dans le traitement des tumeurs ophtalmiques. Comparés aux sources de plus haute énergie, ils offrent l'avantage physique d'une manipulation plus facile du point de vue de la radioprotection et d'une réduction de la dose aux tissus sains avoisinants. Lors d'un traitement clinique, il est important de pouvoir rapporter précisément certains paramètres qui peuvent avoir un impact sur les résultats du traitement pour le patient. Ces paramètres sont par exemple la dose prescrite, les doses à différents organes, le degré d'uniformité atteint sur le volume cible ou certaines information dose-volume. Les systèmes de calcul de dose pour la curiethérapie permettent également de mieux en mieux de conformer le traitement à l'anatomie du patient, comme en radiothérapie externe. Dans le cas des tumeurs de la prostate, il est possible depuis quelques années de vérifier la position des sources et de calculer les distributions de dose en temps réel lors de l'implantation grâce à l'utilisation d'une sonde ultrason. Il est clair que pour arriver à une bonne précision lors du traitement lui-même, il est indispensable de connaître les caractéristiques dosimétriques des sources avec un maximum de précision. Vu que, à ces basses énergies, les caractéristiques dosimétriques sont très dépendantes de la structure interne de sources, cela implique une étude minutieuse et complète de chaque nouveau modèle de grain mis sur le marché. Ceci forme le cadre de la première partie de ce travail qui est la détermination des caractéristiques dosimétriques de deux nouveaux modèles de grains produits par la compagnie IBt. Cette détermination a été réalisée grâce à l'utilisation de détecteurs thermoluminescents et de calculs de Monte Carlo avec deux codes : MCNP4C et EGSnrc. La comparaison de ces deux codes avec les mesures a permit de mettre en évidence la nécessité d'utiliser des bibliothèques de sections efficaces récentes lors des calculs et de prendre en compte la présence des détecteurs lors des mesures. Cependant, la dosimétrie est seulement une partie du problème de l'étude des sources de basse énergie utilisée en radiothérapie. Les milieux irradiés sont des tissus vivants complexes, composés de plusieurs sous-structures sur lesquelles les radiations n'ont pas toujours le même effet. Il est maintenant largement accepté que la partie la plus radiosensible d'une cellule vivante est la molécule d'ADN. Lorsque des photons interagissent avec la matière, ils produisent des électrons qui eux-mêmes déposent leur énergie sous forme d'ionisation et d'excitations des atomes. Ces ionisations, si elles se passent à proximité ou dans la molécule d'ADN, peuvent endommager celle-ci. Cette structure peut supporter une certaine quantité de dommages tout en restant fonctionnelle grâce à des mécanismes de réparation, mais ces mécanismes connaissent des limitations dans la manière dont ils peuvent gérer les cassures multiples des brins d'ADN. Si ces cassures sont trop proches dans le temps et l'espace, la réparation n'est plus possible et la cellule meurt. La densité d'ionisation augmente lorsque l'énergie du photon décroît. Il est donc attendu que les photons de basse énergie auront un efficacité accrue pour tuer les cellules et donc une plus haute Efficacité Biologique Relative que les photons de haute énergie. Pour étudier cet effet, il faut réduire le volume de matière considéré lors des évènements responsables du dépôt d'énergie afin d'atteindre des dimensions d'un ordre représentatif pour les cellules, microns ou même nanomètres. C'est pourquoi la partie de la physique des radiations dédiée à ce problème est appelée microdosimétrie. La seconde partie de ce travail est dédiée à une étude microdosimétrique théorique des ces sources de photons de basse énergie utilisant le concept de fonctions microdosimétriques pour des volumes de 1 µm et une analyse de cluster pour des volumes de l'ordre du nanomètre. Les deux types d'étude montrent que les photons de 20-30 keV sont effectivement plus efficaces radiobiologiquement que les photons de hautes énergie. Les résultats microdosimétriques donnent un rapport de yD relatif au Co-60 de 2.6 pour le Pd-103 et de 2.2 pour l'I-125, ce qui est concordant avec les valeurs expérimentales publiées par Wuu et al. (1996). L'analyse de cluster montre aussi que les électrons générés par les photons de 20-30 keV produisent plus de clusters d'ordre élevé que les électrons de 300 keV. L'ordre de cluster moyen pour des clusters de 10 nm est de 3.0 pour le 103Pd et de 3.3 pour l'I-125 comparé à 2.1 pour des électrons de 300keV. Dans ce cas, l'I-125 montre un ordre de cluster moyen plus élevé, ce qui est lié à une efficacité biologique potentiellement plus élevée. Ceci est expliqué par le fait que les photons de l'I-125 interagissent par effet Compton avec une probabilité de 51 %, cette probabilité tombe à 27 % dans le cas du Pd-103. Les interactions Compton génèrent un grand nombre d'électrons de très basse énergie qui déposent leur énergie très localement avec une très haute densité d'ionisations. D'autres études biologiques sont nécessaires pour déterminer quel volume est plus adapté et donc quel type d'étude reflète le mieux la réalité.

Solutions were developed to solve the special photon transport problems. A respective Monte Carlo code were implemented. The photon transport calculations were made for simulation of light distribution in tissues, for 3D estimation of dosis brachytherapy source, for determination of scattering coincidences in PET and for solving of one special problem for ambient dosimetry. The developed calculation methods are based on a purely statistical approach and can therefore universal applied. The efficiency in respect of precision, numerical effectiveness, as well as memory requirement were optimised and verified by the calculations done. The modification and enhancement can be easy realised thanks to modular, object orientated implementation of the program and enables development of new application fields in physics and medicine. / Zur Lösung speziellen Probleme des Photonentransports wurden grundlegende Ansätze gewonnen und in einem Transportprogramm umgesetzt. Die Photonen-Transportrechnungen zur Simulation des Lichttransports in trüben Medien, für die dreidimensionale Dosisberechnung für interstitielle Brachytherapiequellen, für die Ermittlung der Streustrahlungsverteilungen in PET-Scannern und für die Lösung eines speziellen Problems der Umgebungsdosimetrie wurden durchgeführt. Die entwickelten Berechnungsmethoden basieren auf einer rein probabilistischen Herangehensweise und lassen sich deshalb universell verwenden. Ihre Leistungsfähigkeit hinsichtlich der Genauigkeit, der numerischen Effektivität sowie des Bedarfs an Rechenresourcen wurde optimiert und konnte durch die durchgeführten Berechnungen bestätigt werden. Durch den modularen, objektorientierten Programmaufbau sind Modifikationen und Erweiterungen relativ einfach durchzuführen. Das ermöglicht eine Erschließung von neuen Anwendungsgebieten in der Physik und Medizin.

En dosimétrie, la distribution énergétique des photons émis par une source constitue un paramètre incontournable. Dans la gamme des basses et moyennes énergies (i.e. E<150 keV ici), le LNHB possède 5 tubes à rayons X ainsi que des grains de curiethérapie à l’iode 125, présentant des hauts débits de fluence. La détermination du spectre émis par calcul (déterministe ou Monte-Carlo) est limitée, dans la gamme d’énergie considérée, par les incertitudes élevées sur les bases de données ainsi que par les approximations du modèle. La mesure directe avec un détecteur au germanium ultra-pur a donc été retenue, bien que nécessitant de lourds moyens. De plus, le spectre mesuré est le produit de convolution du spectre émis recherché par la réponse du système. Une fois la réponse du détecteur modélisée, il est possible de « déconvoluer» la mesure, i.e. de remonter au spectre réellement émis en corrigeant (par stripping, model-fitting, inférence bayésienne…) les déformations spectrales induites par le processus de détection. Pour la curiethérapie, le modèle de grain-source a ainsi pu être ajusté. Pour les tubes à rayons X, les résultats obtenus avec différents codes Monte-Carlo et 4 logiciels déterministes ont été comparés à un spectre dit de référence obtenu par mesure et déconvolué. Ainsi l’impact sur certaines grandeurs dosimétriques de la méthode utilisée a pu être quantifié. / In the field of dosimetry, the knowledge of the whole photon fluence spectrum is an essential parameter. In the low-to-medium energy range (i.e. E<150 keV), the LNHB possess 5 X-ray tubes and iodine-125 brachytherapy seeds, both emitting high fluence rates. The performance of calculation (either Monte Carlo codes or deterministic software) is flawed by increasing uncertainties on fundamental parameters at low energies, and modelling issues. Therefore, direct measurement using a high purity germanium is preferred, even though it requires a time-consuming set-up and mathematical methods to infer impinging spectrum from measured ones (such as stripping, model-fitting or Bayesian inference…). Concerning brachytherapy, the knowledge of the seed’s parameters has been improved. Moreover, various calculated X-ray tube fluence spectra have been compared to measured ones, after unfolding. The results of all these methods have then be assessed, as well as their impact on dosimetric parameters.

科学研究費補助金 研究種目:一般研究(C) 課題番号:04680228 研究代表者:青山 隆彦 研究期間:1992-1993年度