The goal of radiation therapy is to eradicate tumour cells while minimising radiation dose to healthy tissues. Ions including protons and carbon ions have gained increasing interest for cancer treatment. Advantages of ion beam therapy are conformal dose distribution, and for ions heavier than protons increased biological effectiveness in cell killing, compared to conventional radiation therapy using photons. Despite these advantages, fundamental problems in ion beam therapy include accuracy of dose determination at the cellular level, and characterisation of the radiation quality at the microscopic scale. Due to the high density of interactions along ion tracks, inhomogeneity of dose and track parameters at the cellular level is one of the major concerns for ion beam therapy. The aim of the thesis is to develop computational tools for dosimetry of ion tracks at the molecular level. Event-by-event Monte Carlo track structure (MCTS) simulations were developed for full-slowing-down tracks of protons and carbon ions in water representing cellular environment. In Paper I, the extension of the MCTS code KURBUC_proton was carried out to energies up to 300 MeV, covering the entire proton energy range used in radiation therapy. Physical properties and microdosimetry of proton tracks were investigated and benchmarked with the experimental data. Papers II-V describe the development of the MCTS code for full-slowing-down tracks of carbon ions. In Papers II-IV, the classical trajectory Monte Carlo (CTMC) model was developed for the calculation of interaction cross sections for low and intermediate energy carbon projectiles of all charge states (C0 to C6+) in water. In Paper V, the calculated cross sections were implemented in a new MCTS code KURBUC_carbon simulating carbon ions of energies 1-104 keV/u in water. This development allows the investigation of track parameters in the Bragg peak region of carbon ion beams. Publication of the thesis and the published papers make contribution to the physics of ion interactions in matter, and provide a new and complete database of electronic interaction cross sections for low and intermediate energy carbon projectiles of all charge states in water. The MCTS codes for protons and carbon ions provide new tools for biophysical study, including microdosimetry, of ion tracks at cellular and subcellular levels, in particular in the Bragg peak region of these ions. / <p>At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 4: Submitted. Paper 5: Submitted.</p><p> </p>

Abstract: Monte Carlo simulations were used here to compare the radiation chemistry of pure water and aqueous bromide solutions after irradiation with two different types of radiation, namely, tritium β-electrons (~7.8 keV) and [superscript 60]Co γ-rays/fast electron (~1 MeV) or high energy protons. Bromide ions (Br-) are known to be selective scavengers of hydroxyl radicals •OH precursors of hydrogen peroxide H[subscript 2]O[subscript 2]. These simulations thus allowed us to determine the yields (or G-values) of H[subscript 2]O[subscript 2] in the radiolysis of dilute aqueous bromide solutions by the two types of radiations studied, the first with low linear energy transfer (LET) (~0.3 keV/μm) and the second with high LET (~6 keV/μm) at 25 °C. This study was carried out under a wide range of Br- concentrations both in the presence and the absence of oxygen. Simulations clearly showed that irradiation by tritium β-electrons favored a clear increase in G(H[subscript 2]O[subscript 2]) compared to [superscript 60]Co γ-rays. We found that these changes could be related to differences in the initial spatial distributions of radiolytic species (i.e., the structure of the electron tracks, the low-energy β-electrons of tritium depositing their energy as cylindrical “short tracks” and the energetic Compton electrons produced by γ-radiolysis forming mainly spherical “spurs”). Moreover, simulations also showed that the presence of oxygen, a very good scavenger of hydrated electrons (e-[subscript aq]) and H• atoms on the 10[superscript-7] s time scale (i.e., before the end of spur expansion), protected H[subscript 2]O[subscript 2] from further reactions with these species in the homogeneous stage of radiolysis. This protection against e-[subscript aq] and H• atoms therefore led to an increase in the H[subscript 2]O[subscript 2] yields at long times, as seen experimentally. Finally, for both deaerated and aerated solutions, the H[subscript 2]O[subscript 2] yield in tritium β-radiolysis was found to be more easily suppressed than in the case of cobalt-60 γ-radiolysis, and interpreted by the quantitatively different chemistry between short tracks and spurs. These differences in the scavengeability of H[subscript 2]O[subscript 2] precursors in passing from low-LET [superscript 60]Co γ-ray to high-LET tritium β-electron irradiation were in good agreement with experimental data, thereby lending strong support to the picture of tritium-β radiolysis in terms of short tracks of high local LET. / Résumé: Les simulations Monte Carlo constituent une approche théorique efficace pour étudier la chimie sous rayonnement de l'eau et des solutions aqueuses. Dans ce travail, nous avons utilisé ces simulations pour comparer l’action de deux types de rayonnement, à savoir, le rayonnement γ de [indice supérieur 60]Co (électrons de Compton ~1 Me V) et les électrons β du tritium (~ 7,8 keV), sur la radiolyse de l’eau et des solutions aqueuses diluées de bromure. Les ions Br- sont connus comme d’excellents capteurs des radicaux hydroxyles •OH, précurseurs du peroxyde d’hydrogène H[indice inférieur 2]O[indice inférieur 2]. Les simulations Monte Carlo nous ont donc permis de déterminer les rendements (ou valeurs G) de H[indice inférieur 2]O[indice inférieur 2] à 25 °C pour les deux types de rayonnements étudiés, le premier à faible transfert d'énergie linéaire (TEL) (~0,3 keV/μm) et le second à haut TEL (~6 keV/μm). L’étude a été menée pour différentes concentrations d’ions Br-, à la fois en présence et en absence d'oxygène. Les simulations ont montré que l’irradiation par les électrons β du tritium favorisait nettement la formation de H[indice inférieur 2]O[indice inférieur 2] comparativement aux rayons γ du cobalt. Ces changements ont pu être reliés aux différences qui existent dans les distributions spatiales initiales des espèces radiolytiques (i.e., la structure des trajectoires d'électrons, les électrons β du tritium déposant leur énergie sous forme de «trajectoires courtes» de nature cylindrique, et les électrons Compton produits par la radiolyse γ formant principalement des «grappes» de géométrie plus ou moins sphérique). Les simulations ont montré également que la présence d'oxygène, capteur d’électrons hydratés et d’atomes H• sur l'échelle de temps de ~10[indice supérieur -7] s (i.e., avant la fin des grappes), protégeait H[indice inférieur 2]O[indice inférieur 2] d’éventuelles réactions subséquentes avec ces espèces. Une telle «protection» conduit ainsi à une augmentation de G(H[indice inférieur 2]O[indice inférieur 2]) à temps longs. Enfin, en milieu tant désaéré qu’aéré, les rendements en H[indice inférieur 2]O[indice inférieur 2] obtenus lors de la radiolyse par les électrons β du tritium ont été trouvés plus facilement supprimés que lors de la radiolyse γ. Ces différences dans l’efficacité de capture des précurseurs de H[indice inférieur 2]O[indice inférieur 2] ont été interprétées par les différences quantitatives dans la chimie intervenant dans les trajectoires courtes et les grappes. Un excellent accord a été obtenu avec les données expérimentales existantes.