SERIA: MEDYCYNA NUKLEARNA

Time-of-Flight PET: Matematyka, fizyka i zastosowania

Δx = c·Δt/2 - równanie, które zmniejszyło szum PET 2-krotnie i zredukowało dawkę promieniowania o połowę. Deep dive w fizyk TOF.

📅 1983-2026
⚡ TOF: 6ns → 214ps
📐 Matematyka + Fizyka
📖 ~7500 znaków

1. Zasada PET - koincydencja bez TOF

Zacznijmy od fundamentów. Pozytron e⁺ emitowany z jądra anihiluje z elektronem e⁻ tkanki:

e⁺ + e⁻ → 2γ (511 keV każdy, kąt ~180°) Zachowanie energii: 2mec² = 2 × 511 keV Zachowanie pędu: p⃗γ1 + p⃗γ2 = 0 → przeciwne kierunki (180°)

Klasyczny PET bez TOF: Dwa detektory rejestrują fotony w koincydencji (okno czasowe ~2-12 ns). Wiemy, że anihilacja nastąpiła GDZIEŚ na linii łączącej dwa detektory (Line of Response, LOR), ale NIE WIEMY GDZIE.

Problem lokalizacji w PET bez TOF

Pacjent średnicy D = 40 cm. Anihilacja może być w dowolnym punkcie na LOR. Każdy voxel 4×4×4 mm na LOR ma tę samą prawdopodobność a priori.

Liczba voxeli na LOR (L = 40 cm, voxel 4 mm): N = 400 mm / 4 mm = 100 voxeli Prawdopodobieństwo dla każdego voxela: p = 1/100 = 1% Szum obrazu (variance) ∝ √N = √100 = 10

Rekonstrukcja: miliony LOR-ów przecina się → FBRT (Filtered Back-Projection) lub OSEM (Ordered Subset Expectation Maximization). Ale bez dodatkowej informacji o pozycji, obraz ma fundamentalny szum statystyczny.

2. Time-of-Flight - genialny pomysł

Pomysł: Zmierzyć RÓŻNICĘ czasu detekcji fotonów!

Jeśli foton γ₁ dotrze do detektora A w czasie t₁, a foton γ₂ do detektora B w czasie t₂, to różnica Δt = t₁ - t₂ informuje nas o POŁOŻENIU anihilacji na LOR!

Matematyka TOF - wyprowadzenie

Układ współrzędnych: detektor A w x=0, detektor B w x=D Anihilacja w x=x₀ (0 ≤ x₀ ≤ D) Odległości: d₁ = x₀ (anihilacja → detektor A) d₂ = D - x₀ (anihilacja → detektor B) Czasy lotu fotonów (c = 30 cm/ns): t₁ = d₁/c = x₀/c t₂ = d₂/c = (D-x₀)/c Różnica czasów: Δt = t₁ - t₂ = (x₀/c) - ((D-x₀)/c) Δt = (x₀ - D + x₀) / c Δt = (2x₀ - D) / c Rozwiązujemy dla x₀: x₀ = (D + c·Δt) / 2

Kluczowa formuła TOF:

x₀ = D/2 + c·Δt/2 Interpretacja: - Jeśli Δt = 0 (jednocześnie) → x₀ = D/2 (środek LOR) - Jeśli Δt > 0 (t₁>t₂, A później) → x₀ > D/2 (bliżej A) - Jeśli Δt < 0 (t₁

Rozdzielczość przestrzenna TOF

W praktyce Δt ma niepewność σΔt (timing resolution). To przekłada się na niepewność pozycji:

σx = c · σΔt / 2 gdzie: c = 299,792,458 m/s ≈ 30 cm/ns σΔt = timing resolution FWHM (typ. 200-600 ps) Przykłady: TOF 600 ps (wczesne LSO systemy, 2006): σx = 30 cm/ns × 0.6 ns / 2 = 9 cm FWHM TOF 400 ps (LYSO+PMT, 2010): σx = 30 cm/ns × 0.4 ns / 2 = 6 cm FWHM TOF 214 ps (Siemens Biograph Vision, LYSO+SiPM, 2018): σx = 30 cm/ns × 0.214 ns / 2 = 3.2 cm FWHM TOF 150 ps (LaBr₃+SiPM, prototypy 2023): σx = 30 cm/ns × 0.15 ns / 2 = 2.25 cm FWHM

Interpretacja: Im lepsza rozdzielczość czasowa, tym precyzyjniej wiemy GDZIE na LOR nastąpiła anihilacja.

3. SNR gain - redukcja szumu

To kluczowa korzyść TOF: zmniejszenie szumu statystycznego obrazu. Wyprowadzenie:

Matematyka SNR gain

PET bez TOF: Każdy voxel na LOR ma równe prawdopodobieństwo p = 1/N Liczba voxeli: N = D / Δxvoxel Variance ∝ N PET z TOF: Tylko voxele w obrębie ±σx wokół x₀ mają znaczące p Effective number of voxels: Neff = 2σx / Δxvoxel Variance ∝ Neff SNR gain: Gain = √(N / Neff) = √(D / 2σx)
TOF SNR Gain = √(D / 2σx) = √(D / (c·σΔt)) Przykłady dla D = 40 cm: TOF 600 ps (σx = 9 cm): Gain = √(40 / 18) = √2.22 ≈ 1.49× (49% better SNR) TOF 400 ps (σx = 6 cm): Gain = √(40 / 12) = √3.33 ≈ 1.83× (83% better) TOF 214 ps (σx = 3.2 cm): Gain = √(40 / 6.4) = √6.25 = 2.5× (150% better!) TOF 150 ps (σx = 2.25 cm): Gain = √(40 / 4.5) = √8.89 ≈ 2.98× (198% better!)

Kluczowy wniosek: TOF 214 ps daje 2.5× lepszy SNR niż PET bez TOF. Można to wykorzystać na 3 sposoby:

  • Skrócić czas badania 2.5× (z 20 min → 8 min)
  • Zmniejszyć dawkę FDG 2.5× (z 370 MBq → 148 MBq)
  • Poprawić jakość obrazu 2.5× (lepsze CNR, wykrywalność małych zmian)

4. Historia TOF - dlaczego tak późno?

1983

Ter-Pogossian, Mullani (Washington University) budują pierwszy TOF PET z CsF scyntylatorem. τ=5 ns, TOF 1.2 ns → σx=18 cm. Gain tylko 1.2×. Rozczarowanie - zbyt mała poprawa.

1986-2000

Era BGO - TOF niemożliwy. BGO ma τ=300 ns. Timing resolution 3-6 ns. σx=45-90 cm - GORZEJ niż średnica pacjenta! TOF nie ma sensu. Wszyscy rezygnują z TOF.

1992

Melcher syntetyzuje LSO (τ=40 ns). Timing nagle skacze do 500-800 ps. Ale potrzeba czasu na rozwój technologii...

2006

Philips Gemini TF - pierwszy komercyjny TOF PET! LSO+PMT, TOF 585 ps (σx=8.8 cm), gain 1.5×. Klinika potwierdza: lepsze obrazy, krótszy czas, niższa dawka. Rewolucja.

2010-2015

Wszyscy producenci wprowadzają TOF: Siemens mCT (530 ps), GE Discovery (550 ps). TOF staje się standardem premium PET.

2018

Siemens Biograph Vision - LYSO+SiPM, TOF 214 ps! Gain 2.5×. To zmienia klinikę: 10 min whole-body, 185 MBq FDG.

2022

United Imaging uEXPLORER - total-body PET (194 cm FOV) z TOF 430 ps. Czułość 40× wyższa niż konwencjonalny PET. 2-min badanie!

❓ Dlaczego BGO uniemożliwiało TOF?

Problem 1: Długi czas zaniku (τ=300 ns)
Impuls światła z BGO rozciągnięty na 300-600 ns. Jitter czasowy ~3-6 ns. Dla D=40 cm:

σx = 30 cm/ns × 5 ns / 2 = 75 cm FWHM Gain = √(40 cm / 150 cm) = √0.27 ≈ 0.52× - POGORSZENIE!

Problem 2: Niska wydajność światła (8,200 ph/MeV)
Dla 511 keV: ~4,200 fotonów. Z QE PMT 25% → ~1,000 fotoelektronów. Poisson noise √1000 ≈ 32. Timing jitter rośnie jako 1/√Nph.

Wniosek: BGO fizycznie nie nadaje się do TOF. Potrzeba było LSO (τ=40 ns, 27k ph/MeV).

✓ Dlaczego LSO/LYSO umożliwiło TOF?

1. Szybki czas zaniku (τ=40 ns vs 300 ns):
Impuls światła krótki → timing precision 400-600 ps (vs 3-6 ns w BGO). To 7.5× lepsza rozdzielczość czasowa!

2. Wysoka wydajność światła (27k ph/MeV vs 8k):
Dla 511 keV: ~13,800 fotonów vs 4,200. Więcej fotonów → niższy Poisson noise → lepszy timing.

3. SiPM rewolucja (2010+):
QE 60% vs PMT 25% → 2× więcej fotoelektronów
Timing jitter <100 ps vs PMT 1-3 ns → 3-10× lepsze
Rezultat: LYSO+SiPM osiąga 214 ps!

5. Współczesne systemy TOF (2018-2026)

System Rok Scyntylator Detektor TOF σx SNR gain
Philips Gemini TF 2006 LSO PMT 585 ps 8.8 cm 1.51×
Siemens mCT Flow 2011 LSO PMT 530 ps 8.0 cm 1.58×
GE Discovery MI 2014 LYSO SiPM 375 ps 5.6 cm 1.88×
Philips Vereos 2016 LYSO SiPM 310 ps 4.7 cm 2.07×
Siemens Vision 2018 LYSO SiPM 214 ps 3.2 cm 2.50×
United uEXPLORER 2022 LYSO SiPM 430 ps 6.5 cm 1.76×

Uwaga: uEXPLORER ma "gorszy" TOF (430 ps) niż Vision (214 ps), ale jego 40× wyższa czułość (dzięki 194 cm FOV) całkowicie dominuje. Czułość > TOF dla total-body.

6. Wpływ kliniczny TOF

Jak TOF zmienia praktykę kliniczną?

Czas badania
20 → 10 min

2× szybciej dzięki gain 2×

Dawka FDG
370 → 185 MBq

50% redukcja dozowania

Wykrywalność
10 → 5 mm

Detekcja przerzutów 2× mniejszych

BMI limit
30 → 40

Otyłych pacjentów można badać

Przykłady kliniczne:

  • Onkologia: Wykrycie mikroprzerzutów <5 mm w wątrobie, płucach - niemożliwe w PET bez TOF. Zmiana staging 15-20% pacjentów.
  • Pediatria: Dzieci 10-20 kg otrzymują 50-100 MBq (vs 370 MBq dorośli). TOF umożliwia diagnostyczną jakość przy tak niskich dawkach.
  • Chorzy otyłość (BMI>35): Tłuszcz osłabia fotony. Non-TOF PET daje szumiące obrazy. TOF kompensuje osłabienie.
  • Monitorowanie terapii: Powtarzane badania co 3 miesiące - kumulatywna dawka! TOF zmniejsza lifetime exposure.
  • Badania przesiewowe (screening): Jeśli kiedyś PET stanie się narzędziem screeningu (kontrowersja), TOF jest kluczowy dla niskich dawek.

7. Granica TOF - jak szybko może być?

Fizyczne ograniczenia timing resolution

1. Czas zaniku scyntylatora (τ): Fundamentalny limit. Nawet dla pojedynczego fotonu, timestamp ma jitter ~τ.

Timing jitter ≥ τ / √Nph LSO: τ=40 ns, Nph=13,800 Jitter ≥ 40 ns / √13,800 ≈ 340 ps LaBr₃: τ=16 ns, Nph=32,000 Jitter ≥ 16 ns / √32,000 ≈ 89 ps Teoretyczne minimum dla VISIBLE light: τ ≥ 10 ns → jitter ≥ 50 ps

2. Photon transport time: Foton światła musi przejść przez kryształ do detektora. Dla LYSO 2 cm grubości, n=1.81:

ttransport = n × L / c = 1.81 × 20 mm / 300 mm/ns ≈ 120 ps

3. Detektor jitter: SiPM ma Single Photon Time Resolution (SPTR) 50-100 ps. Nie można zejść poniżej bez kriogeniki.

Wniosek: Praktyczna granica TOF dla room temperature PET to 100-150 ps FWHM. Systemy prototypowe osiągnęły 150 ps (LaBr₃+SiPM, 2023). Zejście poniżej 100 ps wymaga albo kriogeniki, albo Cherenkov light (inny mechanizm).

Przyszłość: TOF <100 ps?

  • Prompt photons: Pierwsze ~100 fotonów z LSO emitowane w <1 ns. Sampling tylko początek impulsu → 100-150 ps możliwe.
  • Cherenkov light: Elektrony >263 keV w krysztale emitują światło Cherenkov (instantaneous!). PbF₂ + LaBr₃ hybrid → TOF <100 ps demonstrowany.
  • Machine learning: CNN do timestamp estimation z raw waveform SiPM. Gain 20-30% w timing resolution (180 ps → 140 ps).
  • Cryo-PET: Chłodzenie detektorów do -20°C zmniejsza dark noise, pozwala na niższy threshold → lepszy timing. TOF 120 ps możliwe.

2030 prediction: Komercyjne systemy osiągną 150 ps, prototypy 100 ps. Gain 3× vs non-TOF. To równoważnik 9× redukcji szumu - jeden badanie FDG z 40 MBq!

Podsumowanie: TOF - game changer PET

Era Lata Scyntylator TOF Status TOF
Wczesny TOF 1983-1990 CsF, BaF₂ 1-2 ns Prototyp, marginal gain
Era BGO 1990-2005 BGO 3-6 ns TOF niemożliwy
Renesans TOF 2006-2015 LSO+PMT 400-600 ps Komercyjny, gain 1.5-2×
Era SiPM 2016-2026 LYSO+SiPM 210-400 ps Standard, gain 2-2.5×
Ultra-fast (prototypy) 2023+ LaBr₃+SiPM 100-150 ps R&D, gain 3×

Esencja: Time-of-Flight to najpotężniejsza innowacja w PET od lat 90. Prosta formuła Δx=c·Δt/2 zredukowała szum obrazu 2-2.5×, umożliwiając 50% redukcję dawki lub czasu badania. Z niemożliwego w erze BGO (300 ns), przez możliwy w LSO (40 ns), do standardu z SiPM (214 ps) - TOF jest przykładem jak fizyka materiałów (scyntylatory) + technologia (SiPM) zmieniają medycynę. Matematyka wygrała.

Bibliografia

  1. Ter-Pogossian MM, Mullani NA, Ficke DC, et al. Photon time-of-flight-assisted positron emission tomography. J Comput Assist Tomogr. 1981;5(2):227-239.
  2. Budinger TF. Time-of-flight positron emission tomography: status relative to conventional PET. J Nucl Med. 1983;24(1):73-76.
  3. Moses WW, Derenzo SE. Prospects for time-of-flight PET using LSO scintillator. IEEE Trans Nucl Sci. 1999;46(3):474-478.
  4. Surti S, Kuhn A, Werner ME, et al. Performance of Philips Gemini TF PET/CT scanner with special consideration for its time-of-flight imaging capabilities. J Nucl Med. 2007;48(3):471-480.
  5. Conti M. Focus on time-of-flight PET: the benefits of improved time resolution. Eur J Nucl Med Mol Imaging. 2011;38(6):1147-1157.
  6. Karp JS, Surti S, Daube-Witherspoon ME, Muehllehner G. Benefit of time-of-flight in PET: experimental and clinical results. J Nucl Med. 2008;49(3):462-470.
  7. Lecoq P. Pushing the limits in time-of-flight PET imaging. IEEE Trans Radiat Plasma Med Sci. 2017;1(6):473-485.
  8. Gundacker S, Auffray E, Jarron P, Meyer T, Lecoq P. On the comparison of analog and digital SiPM readout in terms of expected timing performance. Nucl Instrum Methods Phys Res A. 2015;787:6-11.
  9. van Dam HT, Borghi G, Seifert S, Schaart DR. Sub-200 ps CRT in monolithic scintillator PET detectors using digital SiPM arrays and maximum likelihood interaction time estimation. Phys Med Biol. 2013;58(10):3243-3257.
  10. Gundacker S, Turtos RM, Kratochwil N, et al. Experimental time resolution limits of modern SiPMs and TOF-PET detectors exploring different scintillators and Cherenkov emission. Phys Med Biol. 2020;65(2):025001.
  11. Cherry SR, Jones T, Karp JS, et al. Total-body PET: maximizing sensitivity to create new opportunities for clinical research and patient care. J Nucl Med. 2018;59(1):3-12.
  12. Schaart DR, Charbon E, Frach T, Schulz V. Advances in digital SiPMs and their application in biomedical imaging. Nucl Instrum Methods Phys Res A. 2016;809:31-52.
  13. Conti M, Bendriem B. The new opportunities for high time resolution clinical TOF PET. Clin Transl Imaging. 2019;7:139-147.
  14. Lecoq P, Morel C, Prior JO, et al. Roadmap toward the 10 ps time-of-flight PET challenge. Phys Med Biol. 2020;65(21):21RM01.
  15. Nemallapudi MV, Gundacker S, Lecoq P, Auffray E, Ferri A, Gola A, Piemonte C. Sub-100 ps coincidence time resolution for positron emission tomography with LSO:Ce codoped with Ca. Phys Med Biol. 2015;60(12):4635-4649.

📚 Seria: Medycyna nuklearna

Artykuł #5 • TOF: impossible (BGO) → possible (LSO) → standard (LYSO+SiPM) • Δx=c·Δt/2