Kurze Antwort: Diese Umrechnung erhält man nicht durch eine absolute Simulation, sondern durch Kreuzkalibrierung gegen eine einzelne Messung mit einer kalibrierten Ionisationskammer. Im Folgenden arbeiten wir Schritt für Schritt durch, warum das so ist, die Formel, den Aufbau, die Kammerwahl und die Referenzdokumente.
Warum ist die MC-Ausgabe „pro History“?
Der MC-Code startet ein Teilchen aus Ihrer definierten Quelle (eine „History“), verfolgt es und alle Sekundärteilchen durch die Materie, teilt die deponierte Energie durch die Masse im Scoring-Volumen und mittelt über alle Histories. Das Ergebnis ist die Energiedosis pro Startteilchen.4 Der Code weiß nicht, wie viele Photonen Ihre Röhre pro mAs aussendet, denn das ist eine physikalische Größe, die von der absoluten Röhrenausbeute, dem Anodenwinkel, der Filterung und der Geometrie abhängt — und sie wird nicht als Simulationseingabe bereitgestellt.
Eine Feinheit: Bei Verwendung einer Phasenraumquelle gibt DOSXYZnrc die Dosis „pro einfallendem Teilchen“ an, doch die eigentliche physikalische Verbindung ist die Zahl der auf das Target treffenden Startelektronen, die diesen Phasenraum erzeugt haben (Nincident in BEAMnrc). Bei Teilchen-Recycling/Splitting können „History“ und „echtes Startteilchen“ verwechselt werden — die mit Abstand häufigste Quelle einer falschen Normierung.6 Deshalb ist es am sichersten, die Normierung an eine physikalisch messbare Größe (die Luftkerma) zu binden.
Zwei Ansätze
Grundsätzlich gibt es zwei Wege zur absoluten Dosis:
- (a) Absolute (ab-initio-)Simulation: aus ersten Prinzipien die Zahl der Photonen berechnen, die die Röhre pro mAs erzeugt, und diese in die Simulation einspeisen. Das erfordert die Kenntnis der absoluten Röntgenausbeute; in der Praxis ist es sehr schwierig, fehleranfällig und wird selten gemacht.
- (b) Kreuzkalibrierung gegen die Messung: die Luftkerma (pro mAs) mit einer kalibrierten Kammer in einer bekannten Geometrie messen; dieselbe Geometrie in MC simulieren; aus dem Verhältnis einen Normierungsfaktor ableiten. Dies ist der Standardansatz in der bildgebenden MC-Dosimetrie und der von AAPM TG-195 übernommene Weg.1
Schritte der Kreuzkalibrierung
1. Referenzmessung. Platzieren Sie eine kalibrierte Ionisationskammer an einem definierten Referenzpunkt (typischerweise im Isozentrum oder in einem festgelegten Fokus-Kammer-Abstand, frei in Luft) und messen Sie die Luftkerma. Notieren Sie das mAs derselben Belichtung. Wenden Sie die üblichen kV-Dosimetrie-Korrekturen auf den Rohmesswert an:23
Dabei ist M der Rohmesswert, NK der Luftkerma-Kalibrierfaktor der Kammer (für die relevante Strahlenqualität), kTP die Temperatur-Druck-, kpol die Polaritäts-, ks die Rekombinations- und kQ die Strahlenqualitätskorrektur. Durch mAs teilen: Kair/mAs [mGy/mAs].
2. Dieselbe Geometrie simulieren. Bilden Sie die Messung in MC exakt nach: gleiche kVp, gleiche Gesamtfilterung, gleiche Fokusgeometrie. Scoren Sie die Luftkerma (oder die Dosis in einem kleinen Luft-/Wasservoxel) im gleichen Referenzvolumen → KMC [Gy/history]. Bei kV-Energien gilt in Luft das Sekundärelektronengleichgewicht (CPE), sodass Stoßkerma ≈ Dosis ist; die Umrechnung von Fluenz in Kerma nutzt Massen-Energieabsorptionskoeffizienten (μen/ρ).10
3. Normierungsfaktor. Bilden Sie das Verhältnis der beiden Größen:
Physikalisch beantwortet dieser Faktor die Frage „wie viele Start-Histories einem mAs entsprechen“.
4. Anwenden. Nun können Sie jeden mit dem gleichen Strahl erzeugten MC-Score in absolute Dosis umrechnen:
Die „Gy/mAs“-Konstante ist somit keine einzelne magische Zahl, sondern ein anlagen- und strahlenspezifischer, aus der Messung abgeleiteter Normierungsfaktor.
Kammer und Strahlenqualität
Kammerwahl. Verwenden Sie bei kV-/diagnostischen Energien eine Kammer mit einem Luftkerma-Kalibrierfaktor (NK) für die relevante Strahlenqualität, rückführbar auf ein Primärnormal:2
- Für eine Punktmessung eine kleine, belüftete Kammer: Allzweck-Fingerhut-/Farmer-Typ (≈0,6 cc) oder 0,125 cc; bei sehr niedrigen Energien (Mammographieband) eine dünnfenstrige „Soft-X“-Kammer.
- Für eine CTDI-ähnliche integrierte Messung die 100-mm-CT-Stiftkammer.
Kritischer Punkt: Die Kalibrierung der Kammer muss bei einer Strahlenqualität erfolgen, die zur HVL des gemessenen Strahls passt. Andernfalls wird NK bei der falschen Energie angewandt und die absolute Dosis ist systematisch verschoben.
Strahlenqualität (RQR). In der kV-Dosimetrie werden die Standard-Strahlenqualitäten durch die RQR-Reihe der IEC 61267 definiert und über die HVL charakterisiert.7 Dass Ihre Messung und Ihre Simulation dieselbe HVL liefern, ist der praktischste Nachweis der spektralen Übereinstimmung (siehe unten).
Die Monte-Carlo-Seite
Ein Normierungsfaktor ist nur dann aussagekräftig, wenn die Simulationsphysik stimmt:
- Validieren Sie das Spektrum. Erzeugen Sie das Röhrenspektrum entweder mit einem vollständigen BEAMnrc-Röhrenmodell oder einem validierten Spektrumgenerator (z. B. spekpy/SpekCalc, TASMICS) und validieren Sie es durch Vergleich der simulierten HVL mit der gemessenen HVL. Ist das Spektrum falsch, werden — weil μen/ρ energieabhängig ist — sowohl die Übereinstimmung von Messung und Simulation als auch die Dosis verfälscht.1
- Kerma oder Dosis? Unter CPE in Luft fällt die Stoßkerma mit der Dosis zusammen; ein Weglängen-Kerma-Schätzer in einem kleinen Luftvolumen liefert ein varianzarmes Ergebnis. Modellieren Sie die Kammer explizit, so beziehen Sie Wand-/Kavitätsperturbationen ein; falls nicht, ist das Scoren der Luftkerma frei in Luft und der Abgleich mit der Messung am saubersten.
- Phasenraum-Zählung. Lesen Sie die Zahl der Startteilchen (Nincident) im BEAMnrc-Phasenraumkopf korrekt aus; wurde Recycling/Splitting verwendet, spiegeln Sie den Unterschied zwischen „History“ und „echtem Startteilchen“ korrekt in der Normierung wider.6
- Statistik. Geben Sie die statistische Unsicherheit der Dosis im Referenzvoxel an (typisches Ziel <1 %); die Unsicherheit von CF hängt direkt davon ab.
CBCT-spezifische Hinweise
Für kV-CBCT (Elekta XVI, Varian OBI–ähnlich) ist die Methodik identisch; nur einige zusätzliche Sorgfaltspunkte kommen hinzu:89
- Der Bow-Tie-Filter und die Kollimation formen den Strahl stark; sie müssen in der Simulation vollständig modelliert werden, sonst verschiebt sich das axiale Dosisprofil — und damit der Wert am Normierungspunkt.
- Rotationsbelichtung (Gantry): Verwechseln Sie nicht das mAs einer einzelnen Projektion mit dem Gesamt-mAs einer vollen Rotation; fixieren Sie CF auf eine definierte Referenzbelichtung (z. B. eine einzelne statische Projektion oder einen vollständigen Scan) und verwenden Sie es konsistent.
- Downes et al. bauen ein MC-Modell einer kV-CBCT-Einheit genau mit dieser Logik — Normierung auf die gemessene Dosis — und wenden es auf die Patientendosimetrie an; ein konkretes Beispiel zum Nachmachen.8
Unsicherheitsbudget
Die Unsicherheit der absoluten Dosis ist eine Kombination mehrerer Komponenten: die NK-Kalibrierunsicherheit der Kammer (vom Normallabor, typisch 1,5–3 %, k=2), die Korrekturfaktoren (kTP, ks, kpol), die mAs-Reproduzierbarkeit, die Spektrum-/HVL-Übereinstimmung und die MC-Statistik. Diese getrennt zu berichten ist erforderlich, damit das Ergebnis wissenschaftlich verteidigbar ist; TRS-457 behandelt diese Komponenten im Detail.2
Checkliste
Quellen
- Sechopoulos I, Ali ESM, Badal A, Badano A, Boone JM, Kyprianou IS, Mainegra-Hing E, McMillan KL, McNitt-Gray MF, Rogers DWO, Samei E, Turner AC. Monte Carlo reference data sets for imaging research: Executive summary of the report of AAPM Research Committee Task Group 195. Medical Physics 42(10):5679–5691, 2015. Görüntülemede Monte Carlo çıktısının history başına verilmesi ve mutlak doza normalizasyonu için başvuru kaynağı.
- International Atomic Energy Agency. Dosimetry in Diagnostic Radiology: An International Code of Practice. Technical Reports Series No. 457 (TRS-457). IAEA, Viyana, 2007. kV hava-kerma dozimetrisi; iyon odası tipleri, kalibrasyon katsayısı NK, RQR demet kaliteleri ve düzeltme faktörleri.
- Ma C-M, Coffey CW, DeWerd LA, Liu C, Nath R, Seltzer SM, Seuntjens JP. AAPM protocol for 40–300 kV x-ray beam dosimetry in radiotherapy and radiobiology (TG-61). Medical Physics 28(6):868–893, 2001. Havada ve fantomda kV dozimetri formalizmi; hava-kerma yaklaşımı ve geri-saçılma/kütle enerji-soğurma oranları.
- Kawrakow I, Mainegra-Hing E, Rogers DWO, Tessier F, Walters BRB. The EGSnrc Code System: Monte Carlo Simulation of Electron and Photon Transport. NRC Report PIRS-701, National Research Council Canada, Ottawa. Taşıma fiziği ve skorlama konvansiyonları.
- Walters B, Kawrakow I, Rogers DWO. DOSXYZnrc Users Manual. NRC Report PIRS-794revB, National Research Council Canada. Doz çıktısının 'incident particle başına' (Gy/parçacık) verilmesi ve normalizasyonu.
- Rogers DWO, Walters B, Kawrakow I. BEAMnrc Users Manual. NRC Report PIRS-509(A)revL, National Research Council Canada. Faz-uzayı normalizasyonu ve hedefe düşen başlangıç elektron sayısının (Nincident) izlenmesi.
- International Electrotechnical Commission. Medical diagnostic X-ray equipment – Radiation conditions for use in the determination of characteristics. IEC 61267. Standart RQR demet kaliteleri (HVL ile tanımlı).
- Downes P, Jarvis R, Radu E, Kawrakow I, Spezi E. Monte Carlo simulation and patient dosimetry for a kilovoltage cone-beam CT unit. Medical Physics 36(9):4156–4167, 2009. Bir kV-CBCT ünitesi için MC modelinin ölçülen doza normalizasyonu — çalışılmış örnek.
- Ding GX, Duggan DM, Coffey CW. Characteristics of kilovoltage x-ray beams used for cone-beam computed tomography in radiation therapy. Physics in Medicine and Biology 52(6):1595–1615, 2007. CBCT kV demetinin karakterizasyonu ve kalibrasyonu.
- ICRU Report 90. Key Data for Ionizing-Radiation Dosimetry: Measurement Standards and Applications. ICRU/Oxford University Press, 2016. Kütle enerji-soğurma katsayıları (μen/ρ) ve etkileşim verileri için güncel referans (akıdan kermaya geçişte kullanılır). Alternatif: NIST XCOM/tables.
- İlişkili DoseSave yazıları: CTDIvol, DLP ve SSDE · Yarı Değer Katmanı (HVL) · Radyasyon Birimleri: Gy, Sv, kerma · Kalite Kontrolün Temelleri