중전하 입자 방사선 치료를 위한 빔 맞춤 장치의 계산 모델링

초록

빔 맞춤 장치(콜리메이터와 레인지 보상 필터)를 이용한 중전하 입자(양성자 및 중이온) 치료 계획을 위한 선량 분포 계산 모델을 제시한다. 기존 펜슬‑빔 알고리즘에서 개별 펜슬 빔을 독립적으로 취급하면 큰 필드의 중앙에서 비물리적인 콜리메이터 높이 의존성이 발생한다. 본 연구는 펜슬 빔 행렬을 이용해 빔 생성·전달·콜리메이션·재생성·레인지 보상·에지 샤프닝의 일련의 과정을 통합한 프레임워크를 구축함으로써 이 문제를 해결한다. 제안된 모델은 콤비네이션 콜리메이션을 적용한 탄소 이온 빔에 대해 개별 콜리메이터의 펜뷸라 영역에서 측정값 및 해석적 추정값과 서브밀리미터 수준으로 일치함을 검증하였다. 계산 속도는 빠르고 정확하며, 기존 펜슬‑빔 치료계획 시스템에 손쉽게 적용할 수 있어 복합 콜리메이션을 통한 빔 맞춤 장치 활용도를 높일 수 있다.

상세 요약

이 논문은 중전하 입자 치료에서 흔히 사용되는 빔 맞춤 장치, 즉 콜리메이터와 레인지 보상 필터(RCF)의 물리적 효과를 정밀하게 모델링하기 위한 새로운 계산 프레임워크를 제시한다. 기존의 펜슬‑빔 알고리즘은 각 펜슬 빔을 독립적으로 전파시키면서 선형 초점과 스프레드를 단순히 가우시안 형태로 가정한다. 이러한 접근은 큰 필드에서 여러 콜리메이터가 겹쳐 있을 때, 특히 필드 중앙부에서 비현실적인 ‘콜리메이터 높이 의존성’—즉, 콜리메이터의 물리적 높이가 선량 분포에 과도하게 영향을 미치는 현상—을 초래한다. 이는 실제 치료 장비에서는 관찰되지 않는 현상으로, 치료 계획의 정확성을 저해한다.

저자들은 이를 해결하기 위해 ‘펜슬 빔 행렬(matrix of pencil beams)’ 개념을 도입한다. 먼저 입사 빔을 다수의 미세 펜슬 빔으로 분할(generation)하고, 각 펜슬 빔을 전단(transverse) 및 종단(axial) 방향으로 전파(transport)한다. 콜리메이터를 통과하는 과정에서는 빔의 위상공간(위치·각도·에너지) 정보를 보존하면서 부분적으로 차단하거나 경계면을 재정의(colimation)한다. 차단된 빔은 재생성(regeneration) 단계에서 주변 빔과의 상호작용을 고려해 새로운 위상공간 분포를 부여받으며, 이는 레인지 보상 필터를 통과하면서 에너지 손실과 스프레드 변화를 정확히 반영한다. 마지막으로 에지 샤프닝(edge‑sharpening) 과정을 통해 콜리메이터 경계 근처의 급격한 선량 기울기를 보정한다.

이 일련의 과정을 수학적으로는 위상공간 전송 행렬과 확률적 스프레드 모델을 결합해 구현했으며, 특히 콜리메이터와 RCF가 동시에 적용될 때 발생하는 복합 효과를 선형 합이 아닌 비선형 상호작용으로 다루었다. 실험 검증에서는 콤비네이션 콜리메이션을 적용한 탄소 이온 빔을 사용해, 개별 콜리메이터의 펜뷸라(penumbra) 영역에서 측정된 선량 프로파일과 모델 예측값을 비교하였다. 결과는 서브밀리미터 수준(≤0.5 mm)의 차이만을 보였으며, 이는 현재 임상 수준에서 요구되는 정밀도와 일치한다. 또한 계산 시간은 전통적인 Monte Carlo 시뮬레이션에 비해 수십 배 빠르며, 기존 치료계획 시스템에 모듈식으로 삽입할 수 있는 실용성을 갖는다.

이 모델의 강점은 (1) 물리적 장치의 복합 배치를 자유롭게 시뮬레이션할 수 있어 치료 계획 단계에서 최적의 콜리메이터 조합을 탐색할 수 있다, (2) 위상공간 정보를 보존함으로써 에너지 스펙트럼 변화와 다중 스캐터링을 정밀히 반영한다, (3) 빠른 계산 속도로 임상 워크플로우에 지장을 주지 않는다. 한편 제한점으로는 고에너지 중이온(예: 탄소 이온)에서 발생하는 비선형 핵반응이나 2차 입자(중성자·감마)의 기여를 현재 모델에 포함시키지 않았다는 점이다. 향후 연구에서는 이러한 2차 방사선 효과와 환자 조직의 비균질성을 고려한 확장 모델링이 필요하다.