Geant4-DNA models for the track-structure simulation of lithium and carbon ions in water: a new approach for ions heavier than helium

Wanneer artsen kanker behandelen met bundels deeltjes, gebruiken ze vaak protonen of zwaardere ionen zoals koolstof. Deze deeltjes zijn krachtige instrumenten omdat ze hun energie zeer precies afgeven, waardoor tumoren worden vernietigd terwijl het omliggende gezonde weefsel wordt gespaard. Om deze behandelingen veilig te plannen, vertrouwen wetenschappers op computersimulaties die volgen hoe deze minuscule deeltjes door het lichaam bewegen, dat grotendeels uit water bestaat. Decennialang hebben deze simulaties goed gewerkt voor protonen en heliumionen, maar ze hadden moeite met zwaardere ionen zoals lithium en koolstof. Het probleem is dat deze zwaardere ionen, naarmate ze afremmen, niet alleen botsen met atomen; ze wisselen voortdurend elektronen uit met de watermoleculen waar ze langs passeren, wat hun elektrische lading verandert in een complexe dans. Zonder een gedetailleerde kaart van deze interacties kunnen de simulaties niet nauwkeurig voorspellen waar de energie terechtkomt of hoeveel schade het aanricht aan levende cellen.

Een team van onderzoekers heeft dit gat nu opgevuld door een nieuwe set regels te creëren voor de Geant4-DNA software, een krachtig open-source programma dat wordt gebruikt om straling op de kleinste schaal te simuleren. Ze ontwikkelden specifieke modellen voor lithium- en koolstofionen die bijhouden elke keer dat een ion een elektron verliest of wint terwijl het door water reist. Door te berekenen hoe deze ladingveranderingen het pad en het energieverlies van het deeltje beïnvloeden, creëerde het team een veel realistischer beeld van wat er gebeurt wanneer deze ionen door biologisch weefsel bewegen. Hun werk stelt wetenschappers in staat om de reis van deze deeltjes te simuleren met een detailniveau dat voorheen alleen voor lichtere ionen was voorbehouden, wat de deur opent naar betere kankerbehandelingen en nauwkeurigere veiligheidsbeoordelingen voor astronauten die in de ruimte reizen.

De kern van dit nieuwe werk ligt in het begrijpen van hoe ionen interageren met water bij zeer lage snelheden. Wanneer een snel bewegend ion het lichaam binnendringt, is het meestal gestript van zijn elektronen en draagt het een sterke positieve lading. Terwijl het vertraagt, begint het elektronen te kapen van de watermoleculen om zich heen, waardoor het minder geladen wordt, en soms zelfs weer elektronen verliest. Deze voortdurende verschuiving van de ladingsstatus is cruciaal omdat het bepaalt hoeveel energie het ion overdraagt aan het weefsel. De onderzoekers realiseerden zich dat bestaande modellen te simpel waren; ze behandelden alle zware ionen alsof het één enkel, gemiddeld type deeltje was, waarbij de specifieke stappen van dit proces van ladingswisseling werden genegeerd. Om dit op te lossen, bouwde het team een nieuwe aanpak die de waarschijnlijkheid berekent dat een ion een elektron opvangt of verliest op elk moment van zijn reis. Ze gebruikten een methode die de paden van de deeltjes en elektronen simuleert alsof het kleine biljartballen zijn, waarbij de bewegingsvergelijkingen worden opgelost om precies te zien hoe vaak deze interacties plaatsvinden.

Voor de processen waarbij het ion simpelweg elektronen uit wateratomen slaat of ze exciteert zonder ze te verwijderen, gebruikte het team een schalingsmethode. Ze namen het bekende gedrag van protonen en pasten dit aan op basis van de snelheid en lading van het lithium- of koolstofion. Dit stelde hen in staat om te schatten hoeveel energie de zwaardere ionen zouden verliezen zonder dat er voor elk scenario miljoenen nieuwe, tijdrovende berekeningen uitgevoerd hoefden te worden. Ze namen ook de effecten op van het stuiteren van het ion tegen watermoleculen zonder energie te verliezen, een proces dat elastische verstrooiing wordt genoemd, wat de richting van het deeltje verandert. Door deze verschillende fysieke processen te combineren, creëerden ze een complete toolkit die een lithium- of koolstofion kan volgen van de snelle binnenkomst tot aan een volledige stilstand.

De onderzoekers testten hun nieuwe modellen door de resultaten te vergelijken met bestaande gegevens uit experimenten en andere theoretische berekeningen. Ze keken naar hoeveel energie de ionen verloren per eenheid afstand, een waarde die bekend staat als remkracht (stopping power), en hoe ver de ionen reisden voordat ze tot stilstand kwamen, bekend als hun bereik (range). De simulaties toonden aan dat de nieuwe modellen zeer goed overeenstemmen met experimentele gegevens bij hoge snelheden en bij zeer lage snelheden. Echter, in een specifiek middelsnelheidsbereik, tussen 50 en 500 kiloelektronenvolt per atomaire massa-eenheid, voorspelden de nieuwe modellen iets minder energieverlies dan sommige oudere referentiedata. Dit middensegment is het moeilijkst te voorspellen omdat het precies de plek is waar de ionen snel van ladingsstatus veranderen. Het feit dat de nieuwe modellen, die deze ladingveranderingen individueel volgen, minder energieverlies voorspellen dan oudere modellen die een gemiddelde lading gebruiken, suggereert dat de oudere modellen de schade in dit specifieke gebied mogelijk overschatten.

Ondanks dit verschil in het middelsnelheidsbereik, is de algemene overeenstemming met de meest betrouwbare beschikbare experimentele gegevens sterk. Het team controleerde ook hoe de ionen hun energie verspreiden in een cirkel rond hun pad, een patroon dat radiale dosisverdeling wordt genoemd. Hier kwamen de nieuwe modellen bijna perfect overeen met de verwachte patronen, wat aantoont dat ze de energieafgifte in het weefsel rond het spoor van het deeltje correct beschrijven. Deze consistentie geeft wetenschappers het vertrouwen dat de modellen de essentiële fysica van de situatie correct vastleggen. De onderzoekers merkten op dat hoewel hun resultaten robuust zijn, het gebrek aan een grote hoeveelheid nieuwe experimentele gegevens voor deze specifieke ionen betekent dat er enige onzekerheid blijft bestaan. Zij benadrukten dat er meer metingen nodig zijn om de nauwkeurigheid van de modellen in dat lastige middelsnelheidsbereik volledig te bevestigen.

De impact van dit werk reikt verder dan alleen het verbeteren van een computerprogramma. Door deze gedetailleerde modellen beschikbaar te maken in Geant4-DNA, hebben de onderzoekers de simulatie van koolstoftherapie mogelijk gemaakt, een groeiende kankerbehandeling die koolstofionen gebruikt om tumoren gericht aan te pakken. Het ondersteunt ook onderzoek naar Boron Neutron Capture Therapy, een behandeling die berust op de interactie van neutronen met boronatomen binnen cellen, wat lithium en andere zware ionen als bijproducten kan produceren. Bovendien zijn de modellen essentieel voor ruimteverkenning, waarbij astronauten worden blootgesteld aan kosmische straling die zware ionen bevat die DNA kunnen beschadigen. Met deze nieuwe instrumenten kunnen wetenschappers nu de biologische effecten van deze zware ionen met veel grotere precisie simuleren, wat helpt bij het ontwerpen van betere stralingsschilden en behandelplannen. De software-infrastructuur die gebouwd is voor lithium en koolstof kan ook snel worden aangepast voor andere ionen, zoals zuurstof en neon, waardoor het vakgebied van de stralingsfysica kan blijven evolueren en de menselijke gezondheid kan beschermen in zowel ziekenhuizen als in de kosmos.