그리드‑도스 스프레딩 알고리즘: 중입자 치료용 고속 선량 계산

초록

본 논문은 기존 펜슬‑빔(PB) 방식의 복잡성을 줄이고, 광빔(BB) 방식보다 정확한 빔 블러링을 제공하는 그리드‑도스 스프레딩(GDS) 알고리즘을 제안한다. 물리 계산을 격자‑대‑격자 에너지 전송으로 분리함으로써 연산 속도를 크게 향상시켰으며, 탄소이온 전치료에서 실제 임상 사례(전립선)와 각도‑입사 양성자 펜슬 빔을 통해 정확도와 적용 가능성을 검증하였다.

상세 분석

GDS 알고리즘은 PB 알고리즘의 핵심인 전단면 가우시안 커널을 3차원 타원형 커널로 근사한다. 이를 위해 먼저 모든 격자점에서 테르마(T)와 전단 스프레드(σ_t)를 사전 계산하고, 이 값을 3개의 1차원 가우시안 함수(σ_x, σ_y, σ_z)로 분해한다. 이렇게 하면 복잡한 레이 트레이싱 적분을 한 번만 수행하고, 이후의 컨볼루션은 단순히 격자‑대‑격자 가중합으로 처리된다. 커널 변형으로 인해 수직·수평(공평) 빔에서는 PB와 수학적으로 동일하지만, 빔이 격자축에 비스듬히 입사할 경우 타원형 변형이 발생한다. 논문에서는 30° 입사 양성자 펜슬 빔에 대해 변형에 따른 선량 오차가 최대 3 mm 이하임을 실험적으로 확인하였다. 이는 임상에서 허용되는 위치 오차 범위 내에 해당한다.
연산 속도 측면에서는 GDS가 BB 대비 1.4배 느리게 동작했으며, 전통적인 PB 알고리즘이 수십 배 정도 느린 것에 비해 현저히 개선된 것이다. 이는 테르마와 스프레드를 격자에 미리 저장하고, “interaction point of view” 방식으로 에너지 전송을 수행함으로써 가능해졌다. 또한, GDS는 기존 BB 코드와 거의 동일한 구조에 삽입할 수 있어, 기존 임상 시스템에 대한 적용 장벽이 낮다.
정확도 검증에서는 전립선 치료 사례에서 BB가 나타낸 비현실적인 급격한 선량 급변을 GDS가 자연스럽게 완화시켰으며, 전형적인 4 mm 전단 스프레드 조건에서 50 % 등선 위치 차이가 격자 해상도 수준(≈2 mm) 이하로 유지되었다. 각도 입사 펜슬 빔에서는 저선량 구역에서 약간의 오차가 발생했지만, 브래그 피크 근처에서도 3 mm 이내의 차이로 임상 허용 기준을 만족한다.
요약하면, GDS는 물리 모델의 정확성을 크게 손상시키지 않으면서 연산 효율성을 크게 높인 실용적인 대안이며, 특히 실시간 적응 방사선 치료나 스캔 빔 최적화와 같은 고속 계산이 요구되는 상황에 적합하다.