프로톤 깊이선량 시뮬레이션의 물리학적 불확실성 분석

초록

본 연구는 Geant4 툴킷에 구현된 다양한 물리 모델과 파라미터를 65 MeV 이하의 프로톤 에너지 범위에서 평가한다. 모델 선택·파라미터 설정에 내재된 인식론적(지식 부족) 불확실성이 선량 분포, 특히 브래그 피크 위치와 폭에 미치는 영향을 정량화하고, 방사선 치료 시뮬레이션에 미치는 잠재적 체계오차를 제시한다.

상세 분석

본 논문은 저에너지(≤ 65 MeV) 프로톤 치료에 필수적인 깊이선량 시뮬레이션을 대상으로, Geant4에서 제공되는 핵심 물리 모델들의 구현 방식과 입력 파라미터가 갖는 인식론적 불확실성을 체계적으로 조사한다. 먼저 전자기 상호작용 모델(표준, 저에너지 전자기, Penelope)과 핵반응 모델(Bertini‑Cascade, Binary‑Cascade, INCL‑XX) 사이의 차이를 비교한다. 각 모델은 입자 감쇠, 다중산란, 비탄성 핵반응 등에 대한 교차섹션과 평균 여기 에너지(I) 값을 다르게 적용한다는 점에서 근본적인 불확실성을 내포한다.

특히 평균 여기 에너지(I)는 물질별로 실험 데이터가 부족하거나 서로 다른 문헌값이 존재해, 스톱핑 파워 계산에 직접적인 편향을 초래한다. 저에너지 영역에서는 전자기 모델이 제공하는 에너지 손실률(dE/dx)이 실험값과 3 % 이상 차이날 수 있으며, 이는 브래그 피크의 깊이 오차로 이어진다. 핵반응 모델의 경우, 핵산란 단면적의 불확실성이 5 %~10 % 수준으로 보고되었으며, 이는 프래그멘테이션에 의한 부수선량 생성에 큰 영향을 미친다.

다중산란 모델(Gaussian, Urban, WentzelVI) 역시 경로 길이와 각도 분포를 다르게 예측한다. 논문은 동일한 입사 조건에서 Urban 모델이 평균 경로 길이를 2 % 과소평가하고, WentzelVI가 1 % 과대평가함을 보여준다. 이러한 차이는 특히 얇은 조직 경계에서 선량 경계면이 흐려지는 현상을 야기한다.

시뮬레이션 결과는 브래그 피크 위치가 모델에 따라 최대 1.2 mm, 피크 폭이 0.8 mm까지 변동함을 보여준다. 임상적으로는 1 mm 수준의 위치 오차가 치료 계획의 정확성을 저해할 수 있다. 또한, 부수선량(플루오르, 알파 등)의 발생률이 모델에 따라 15 % 차이를 보이며, 이는 장기 부작용 위험 평가에 중요한 변수이다.

저자들은 이러한 인식론적 불확실성을 정량화하기 위해 민감도 분석과 베이지안 모델 평균화 기법을 적용하였다. 결과적으로, 평균 여기 에너지와 핵산란 단면적이 전체 불확실성 기여도의 60 % 이상을 차지함을 확인했다. 따라서 실험적 교정 데이터 확보가 가장 시급한 과제로 제시된다.

마지막으로, 방사선 치료 시뮬레이션에서 모델 선택 가이드라인을 제시한다. 고정밀도가 요구되는 경우, 저에너지 전자기 모델과 Binary‑Cascade를 결합하고, 평균 여기 에너지 값을 최신 ICRU 보고서 값으로 교정할 것을 권고한다. 반면, 빠른 검증이 목적이라면 표준 전자기 모델과 Bertini‑Cascade를 사용할 수 있다.