1 |
Extension de l'outil Monte Carlo généraliste Geant4 pour la simulation de la radiolyse de l'eau dans le cadre du projet Geant4-DNAKaramitros, Mathieu 23 November 2012 (has links)
Ce travail, réalisé dans le cadre du projet Geant4-DNA, consiste à concevoir un prototype pour la simulation des effets chimiques précoces des rayonnements ionisants. Le modèle de simulation étudié repose sur la représentation particule-continuum où toutes les molécules sont explicitement simulées et où le solvant est traité comme un continuum. La méthode proposée par cette thèse a pour but d'améliorer les performances de ce type de simulation. Elle se base sur (1) la combinaison d'une méthode de pas en temps dynamiques avec un processus de pont Brownien pour la prise en compte des réactions chimiques et afin d'éviter une simulation à pas en temps fixe, coûteuse en temps de calcul, et (2) sur la structure de données k-d tree pour la recherche du voisin le plus proche permettant, pour une molécule donnée, une localisation rapide du réactif le plus proche. La précision de l'algorithme est démontrée par la comparaison des rendements radiochimiques en fonction du temps et en fonction du transfert d'énergie linéaire avec des résultats d'autres codes Monte-Carlo et des données expérimentales. A partir de ce prototype, une tentative de prédiction du nombre et du type d'interactions radicaux-ADN a été entreprise basée sur d'une description simplifiée du noyau cellulaire. / The purpose of this work, performed under the Geant4-DNA project, is to design a prototype for simulating early chemical effects of ionizing radiation. The studied simulation model is based on the particle-continuum representation where all the molecules are explicitly simulated, and where the solvent is treated as a continuum. The method proposed by this thesis aims at improving the performance of this type of simulation. It is based on (1) a dynamical time steps method with a Brownian bridge process, to account for chemical reactions, which avoids the costly fixed time-step simulations, and (2) on the k-d tree data structure for quickly locating, for a given molecule, its closest reactants. The accuracy of the algorithm is demonstrated by comparing radiochemical yields over time and depending on the linear energy transfer with results obtained from other Monte Carlo codes and experimental data. Using this prototype, an attempt to predict the number and type of radical attacks on DNA has been performed using a simplified description of the cell nucleus.
|
Page generated in 0.1153 seconds