Monte Carlo'da Mutlak Doz: Gy/history'den Gy/mAs'e

Kısa cevap: Bu dönüşüm mutlak bir simülasyonla değil, kalibre bir iyon odasıyla yapılan tek bir ölçüme çapraz-kalibrasyonla elde edilir. Aşağıda bunun neden böyle olduğunu, formülü, düzeneği, oda seçimini ve başvuru belgelerini adım adım ele alıyoruz.

MC çıktısı neden "history başına"dır?

MC kodu, tanımladığınız kaynaktan bir parçacık başlatır (bir "history"), onu ve ürettiği tüm ikincil parçacıkları madde içinde izler, biriken enerjiyi skorlama hacmindeki kütleye böler ve tüm history'ler üzerinden ortalar. Sonuç, birim başlangıç parçacığı başına soğurulan dozdur.⁴ Kod, sizin bir mAs'te tüpten kaç foton çıktığını bilmez; çünkü bu, tüpün mutlak verimine, anot açısına, filtrasyona ve geometriye bağlı fiziksel bir büyüklüktür ve simülasyona girdi olarak verilmez.

Bir incelik: DOSXYZnrc bir faz-uzayı (phase-space) kaynağı kullandığınızda dozu "incident particle başına" verir; ama asıl fiziksel bağ, o faz-uzayını üreten hedefe düşen başlangıç elektron sayısıdır (BEAMnrc'de N_incident). Parçacık geri-dönüşümü (recycling/splitting) kullanıldığında "history" ile "gerçek başlangıç parçacığı" karışabilir; bu, yanlış normalizasyonun en yaygın kaynağıdır.⁶ Bu yüzden normalizasyonu her zaman fiziksel olarak ölçülebilir bir niceliğe (hava kermasına) bağlamak en güvenli yoldur.

İki yaklaşım

Mutlak doza geçmenin ilkesel olarak iki yolu vardır:

(a) Mutlak (ab initio) simülasyon: Tüpün bir mAs'te ürettiği foton sayısını ilk prensiplerden hesaplayıp simülasyona koymak. Bu, mutlak X-ışını verimini bilmeyi gerektirir; pratikte çok zor, hataya açık ve nadiren yapılır.
(b) Ölçüme çapraz-kalibrasyon: Kalibre bir odayla, bilinen bir geometride hava kermasını (mAs başına) ölçmek; aynı geometriyi MC'de simüle etmek; ikisinin oranından bir normalizasyon faktörü türetmek. Görüntüleme MC dozimetrisinde standart yaklaşım budur ve AAPM TG-195'in de benimsediği yoldur.¹

Mutlak doza geçiş: kalibre odayla ölç → aynı geometriyi simüle et → oranı (CF) bul → aynı demetin tüm MC skorlarına uygula.¹

Çapraz-kalibrasyon adımları

1. Referans ölçümü. Kalibre bir iyon odasını tanımlı bir referans noktaya (tipik olarak izosantr ya da belirli bir kaynak-oda uzaklığında, serbest-havada) yerleştirin ve hava kermasını ölçün. Aynı ışınlamanın mAs'ini kaydedin. Ham okumaya standart kV dozimetri düzeltmelerini uygulayın:²³

K_air = M · N_K · k_TP · k_pol · k_s · k_Q

Burada M ham okuma, N_K odanın (ilgili demet kalitesi için) hava-kerma kalibrasyon katsayısı, k_TP sıcaklık-basınç, k_pol polarite, k_s yeniden birleşme (rekombinasyon) ve k_Q demet kalitesi düzeltmesidir. Sonucu mAs'e bölün: K_air/mAs [mGy/mAs].

2. Aynı geometrinin simülasyonu. MC'de ölçümün birebir aynısını kurun: aynı kVp, aynı toplam filtrasyon, aynı kaynak-nokta geometrisi. Aynı referans hacimde hava kermasını (veya küçük bir hava/su vokselinde dozu) skorlayın → K_MC [Gy/history]. kV enerjilerde havada yüklü parçacık dengesi (CPE) geçerli olduğundan çarpışma kerması ≈ doz kabul edilebilir; akıdan kermaya geçişte kütle enerji-soğurma katsayıları (μ_en/ρ) kullanılır.¹⁰

3. Normalizasyon faktörü. İki niceliği oranlayın:

CF = (K_air / mAs) ⁄ (K_MC / history) → [history / mAs]

Bu faktör fiziksel olarak "bir mAs kaç başlangıç history'sine karşılık geliyor" sorusunun cevabıdır.

4. Uygulama. Artık aynı demetle yapılan her MC skorunu mutlak doza çevirebilirsiniz:

D_mutlak [Gy/mAs] = D_MC [Gy/history] × CF [history/mAs]

Böylece "Gy/mAs" sabiti, tek bir sihirli sayı değil, sizin tesisinize ve demetinize özgü, ölçümle türetilmiş bir normalizasyon faktörüdür.

İyon odası ve demet kalitesi

Oda seçimi. kV/tanısal enerjilerde, ilgili demet kalitesi için hava-kerma kalibrasyon katsayısı (N_K) taşıyan, birincil standarda izlenebilir bir oda kullanın:²

Nokta ölçümü için küçük hacimli, havalandırmalı bir oda: genel amaçlı thimble/Farmer tipi (≈0,6 cc) ya da 0,125 cc; çok düşük enerjilerde (mamografi bandı) ince pencereli "soft-X" odası.
CTDI benzeri entegre ölçüm için 100 mm CT kalem odası.

Kritik nokta: odanın kalibrasyonu, ölçtüğünüz demetin HVL'sine uyan demet kalitesinde olmalıdır. Aksi halde N_K yanlış enerjide uygulanır ve mutlak doz sistematik olarak kayar.

Demet kalitesi (RQR). kV dozimetride standart demet kaliteleri IEC 61267'nin RQR serisiyle tanımlanır ve HVL ile karakterize edilir.⁷ Ölçümünüzün ve simülasyonunuzun aynı HVL'yi vermesi, spektrum eşleşmesinin en pratik kanıtıdır (aşağıya bakınız).

MC tarafında dikkat edilecekler

Normalizasyon faktörü ancak simülasyon fiziği doğruysa anlamlıdır:

Spektrumu doğrula. Tüp spektrumunu ya tam bir BEAMnrc tüp modeliyle ya da doğrulanmış bir spektrum üreticisiyle (ör. spekpy/SpekCalc, TASMICS) oluşturun; ardından simüle HVL'yi ölçülen HVL ile karşılaştırarak doğrulayın. Spektrum yanlışsa, μ_en/ρ enerjiye bağlı olduğu için hem ölçüm-simülasyon eşleşmesi hem de doz bozulur.¹
Kerma mı, doz mu? Havada CPE koşulunda çarpışma kerması dozla örtüşür; küçük hava hacminde track-length (iz-uzunluğu) kerma tahmincisi düşük varyansla sonuç verir. Odayı açıkça modelliyorsanız duvar/kovuk perturbasyonlarını dahil edin; modellemiyorsanız serbest-havada hava kermasını skorlayıp ölçümle eşleştirmek en temizidir.
Faz-uzayı sayımı. BEAMnrc faz-uzayı başlığındaki başlangıç parçacık sayısını (N_incident) doğru okuyun; geri-dönüşüm/bölme kullanıldıysa "history" ile "gerçek başlangıç parçacığı" farkını normalizasyona doğru yansıtın.⁶
İstatistik. Referans vokseldeki dozun istatistiksel belirsizliğini (tipik hedef <%1) raporlayın; CF'nin belirsizliği doğrudan buraya bağlıdır.

CBCT'ye özgü notlar

kV-CBCT (Elekta XVI, Varian OBI benzeri) için metodoloji birebir aynıdır; yalnızca birkaç ek özen gerekir:⁸⁹

Bow-tie filtre ve kolimasyon demeti güçlü biçimde şekillendirir; simülasyonda bunlar eksiksiz modellenmeli, aksi halde eksenel doz profili ve dolayısıyla normalizasyon noktasındaki değer kayar.
Dönme (gantry) ile ışınlama: Tek projeksiyonun mAs'i ile tam dönüşün toplam mAs'i karıştırılmamalı; CF'yi tanımlı bir referans ışınlamaya (ör. tek statik projeksiyon ya da tam tarama) sabitleyin ve tutarlı kullanın.
Downes ve ark. bir kV-CBCT ünitesi için MC modelini tam da bu mantıkla — ölçülen doza normalize ederek — kurar ve hasta dozimetrisine uygular; izlenecek somut bir örnektir.⁸

Klinik/araştırma notu

Amaç mutlak doz değil de yalnızca bağıl dağılım (ör. protokoller arası karşılaştırma) ise normalizasyona hiç gerek olmayabilir. Mutlak Gy/mAs gerektiğinde ise tek doğru yol, yukarıdaki ölçüm-temelli çapraz-kalibrasyondur.

Belirsizlik bütçesi

Mutlak dozun belirsizliği birkaç bileşenin birleşimidir: odanın N_K kalibrasyon belirsizliği (standart laboratuvardan, tipik %1,5–3, k=2), düzeltme faktörleri (k_TP, k_s, k_pol), mAs tekrarlanabilirliği, spektrum/HVL eşleşmesi ve MC istatistiği. Bunları ayrı ayrı raporlamak, sonucun bilimsel olarak savunulabilir olması için gereklidir; TRS-457 bu bileşenleri ayrıntılı ele alır.²

Kontrol listesi

Adım adım

1) Kalibre oda + doğru demet kalitesinde N_K · 2) Serbest-havada hava kermasını ölç, mAs'i kaydet, düzeltmeleri uygula · 3) Aynı kVp/filtrasyon/geometriyi simüle et, aynı noktada K_MC/history skorla · 4) HVL ile spektrumu doğrula · 5) CF = (K/mAs)/(K_MC/history) · 6) D_mutlak = D_MC × CF · 7) Belirsizlik bütçesini raporla. Çerçeve: TG-195 + TRS-457 (veya TG-61).

Özet

Monte Carlo history başına doz verir; mutlak Gy/mAs, kalibre bir iyon odasıyla ölçülen hava kermasına çapraz-kalibrasyonla elde edilen tesise özgü bir faktörle bulunur. Anahtar belgeler: AAPM TG-195, IAEA TRS-457, AAPM TG-61; kod tarafında EGSnrc/DOSXYZnrc/BEAMnrc kılavuzları.

Kaynaklar

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ı N_K, 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 (N_incident) 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

Not: Bu içerik eğitim amaçlıdır; klinik karar veya mevzuat uyumu için yetkili medikal fizik uzmanına ve güncel düzenlemelere başvurun.

← Tüm makaleler