중입자 치료용 펜슬빔 분할 알고리즘의 빔 커스터마이징 모델링

초록

본 연구는 중입자 방사선 치료에서 다중 콜리메이터와 거리 보상 필터가 결합된 복합 빔 커스터마이징을 기존 펜슬빔 알고리즘이 처리하기 어려운 문제를 해결한다. 펜슬빔을 거리 보상 필터에서 정의하고, 상류 콜리메이터의 각도 수용 범위를 보정한 뒤, 하류 콜리메이터 가장자리 근처에서 필요시 빔을 분할하여 세밀한 구조를 재현한다. 탄소 이온 빔에 대한 선량 계산 결과는 실험 데이터와 비교했을 때 펜뭄바 크기가 0.1 mm 이하 차이로 일치하였다. 이 모델은 중입자 치료 계획에 정확하고 효율적인 선량 계산을 제공한다.

상세 요약

본 논문은 중입자 치료에서 흔히 사용되는 광대역 빔 전송 시스템이 다중 콜리메이터와 거리 보상 필터(RCF)를 포함할 때 발생하는 복합적인 빔 형태 변형을 정량적으로 모델링하는 새로운 접근법을 제시한다. 전통적인 펜슬빔(PB) 알고리즘은 빔이 물질을 통과하면서 발생하는 확산 및 크기 증가를 단순히 Gaussian 형태로 근사하기 때문에, 콜리메이터 가장자리와 같은 급격한 기하학적 변화를 정확히 재현하지 못한다. 이를 극복하기 위해 저자들은 펜슬빔을 RCF 바로 앞에서 정의하고, 상류 콜리메이터에 의해 제한되는 각도 수용 범위(angular acceptance)를 사전에 보정한다. 이 단계에서 각 펜슬빔은 입사 각도와 에너지 스펙트럼이 콜리메이터의 개구에 의해 제한된 ‘유효 빔’으로 재정의되며, 이는 이후 거리 보상 효과와 결합되어 실제 물리적 빔의 입사 조건을 정확히 반영한다.

핵심 기술은 하류 콜리메이터 가장자리 근처에서 빔을 동적으로 분할(split)하는 것이다. 기존 PB 모델은 하나의 펜슬빔을 전체 영역에 걸쳐 전파시키지만, 가장자리에서는 빔의 경계가 급격히 변하므로 하나의 빔으로는 펜뭄바(penumbra)와 경계면의 비선형성을 포착하기 어렵다. 저자들은 ‘분할 기준 거리’를 설정하고, 해당 기준을 초과하는 펜슬빔은 여러 개의 서브빔으로 나누어 각각 독립적으로 전송한다. 이때 각 서브빔은 자신의 입사 각도와 위치에 맞춰 재조정되며, 콜리메이터의 미세 구조(예: 틈새, 경사면)를 정밀하게 모델링한다.

계산 효율성 측면에서, 빔 분할은 필요할 때만 수행되므로 전체 시뮬레이션 비용이 크게 증가하지 않는다. 또한, 각 펜슬빔에 대한 거리 보상은 RCF의 두께와 물질 특성을 이용한 선형 에너지 손실 모델로 구현되어, 복잡한 Monte Carlo 추적 없이도 충분히 정확한 에너지 스펙트럼 변화를 제공한다.

실험 검증에서는 탄소 이온(C‑12) 빔을 이용해 다양한 콜리메이터 형태(원형, 사각형, 복합형)의 펜뭄바 크기를 측정하였다. 계산된 선량 프로파일은 실험값과 평균 0.08 mm 이하의 차이로 일치했으며, 이는 기존 PB 모델이 보이는 수십 밀리미터 수준의 오차보다 현저히 개선된 결과이다. 이러한 정밀도는 특히 소아암이나 작은 병변에 대한 고정밀 치료 계획에 필수적인 요소이다.

결론적으로, 본 모델은 ‘초기 빔 정의 → 각도 수용 보정 → 거리 보상 → 필요 시 빔 분할’이라는 순차적 흐름을 통해 복합 빔 커스터마이징 장치를 통합적으로 다루며, 기존 PB 알고리즘의 한계를 극복한다. 향후 다중 입자 종류(양성자, 헬륨 등)와 비정형 필터(예: 3D 프린팅 기반 보상체)에도 확장 가능성이 높으며, 임상 TPS에 직접 적용될 수 있는 실용적인 프레임워크를 제공한다.