프로톤 빔렛 모델을 활용한 정밀 선량 계산 기반

초록

본 논문은 GEANT4 기반 몬테카를로 시뮬레이션과 베테-볼크 방정식 적분을 결합한 분석적 프로톤 빔렛 모델을 제시한다. 1차·2차 프로톤, 반동핵 및 중성자 재활용을 모두 고려하고, 빔렛의 축방향 및 횡방향 선량 분포를 정확히 재현한다. 빌드업 현상은 2차 반응 프로톤(≈65 %)과 랜다우 꼬리·가우시안 범위 스트래깅으로 설명한다. 모델 파라미터는 시뮬레이션으로 사전 결정하거나, 실제 측정 데이터에 최소한의 자유 파라미터(주로 초기 에너지 스프레드)만을 이용해 보정한다. 중간 에너지에 대한 파라미터 예측이 가능함을 확인하였다.

상세 분석

이 연구는 임상에서 사용되는 프로톤 치료계의 선량 계산 정확도를 향상시키기 위해, 전통적인 경험적 모델을 넘어선 물리 기반의 분석 모델을 구축하였다. 핵심은 베테-볼크 방정식(BBE)의 정밀 적분을 통해 1차 프로톤의 에너지 손실을 연속적으로 추적하고, GEANT4 시뮬레이션으로부터 얻은 2차 입자(반응 프로톤, 재결합 프로톤, 무거운 반동핵)의 생성·전달 특성을 파라미터화한 점이다. 특히, 빌드업 현상이 고에너지 프로톤 빔에서 관찰되는 이유를 두 가지 메커니즘으로 분리하였다. 첫 번째는 핵반응에 의해 생성되는 2차 프로톤으로, 전체 빌드업의 약 65 %를 차지한다는 정량적 근거를 제시한다. 두 번째는 입자들의 에너지 손실 과정에서 발생하는 랜다우 꼬리와 가우시안 형태의 범위 스트래깅이다. 이 두 요소를 각각 수식화함으로써, 깊이선량곡선(D(z))의 초기 상승 구간을 물리적으로 설명한다. 모델 파라미터는 크게 두 그룹으로 나뉜다. 하나는 시뮬레이션으로부터 직접 추출된 고정 파라미터(예: 평균 전자밀도, 평균 핵반응 단면적 등)이며, 다른 하나는 빔라인 고유 특성(초기 에너지 스프레드, 입사각, 스팟 크기 등)을 반영하기 위해 실험 데이터에 맞춰 조정되는 자유 파라미터이다. 특히, 초기 에너지 스프레드가 빔렛의 범위 불확실성을 주도한다는 점을 강조하고, 몇 개의 대표 깊이선량곡선만으로도 중간 에너지에 대한 파라미터를 선형 보간하거나 경험적 관계식으로 예측할 수 있음을 보였다. 횡방향 확산은 다중 스캐터링 이론을 적용해 Gaussian 핵심과 장거리 꼬리(비탄성 산란)를 결합한 복합 모델로 구현했으며, 이는 실제 측정된 횡단면 프로파일과 높은 일치를 보였다. 전체적으로, 이 모델은 기존의 TPS(치료 계획 시스템)에서 요구하는 계산 속도와 정확도 사이의 균형을 만족시키면서, 물리적 근거에 기반한 파라미터 조정이 가능하도록 설계되었다.