MCHIT 조직 모사 기반 양성자·중이온 치료용 Monte Carlo 모델

본 연구에서는 GEANT4 툴킷(버전 8.2)을 기반으로 한 Monte Carlo Model for Heavy‑ion Therapy(MCHIT)를 이용해 조직과 유사한 매질에서 핵자와 원자핵의 전파를 조사하였다. 모델은 투사체 핵의 파편화와 생성된 핵 파편의 2차 상호작용을 모두 고려한다. 물 및 PMMA 팬텀에 단일 에너지 탄소 이온 빔을 조사한 기존 실

MCHIT 조직 모사 기반 양성자·중이온 치료용 Monte Carlo 모델

초록

본 연구에서는 GEANT4 툴킷(버전 8.2)을 기반으로 한 Monte Carlo Model for Heavy‑ion Therapy(MCHIT)를 이용해 조직과 유사한 매질에서 핵자와 원자핵의 전파를 조사하였다. 모델은 투사체 핵의 파편화와 생성된 핵 파편의 2차 상호작용을 모두 고려한다. 물 및 PMMA 팬텀에 단일 에너지 탄소 이온 빔을 조사한 기존 실험 데이터와 비교하여 모델 예측을 검증하였다. MCHIT는 (1) 물과 PMMA에서의 깊이‑선량 분포, (2) 특정 전하를 가진 파편에 대한 선량, (3) 탄소 이온 빔에 의해 생성된 양전자 방출 핵 파편의 분포, (4) 빔 방향에 대해 다양한 각도에서 측정된 2차 중성자 에너지 스펙트럼을 잘 재현한다. 또한 1차 핵과 다양한 투사체 파편에 대한 방사선 횡단면 선량 프로파일을 계산·논의하여 생물학적 선량 분포 평가에 활용 가능성을 제시한다. 브래그 피크 근처 횡단면 주변에서는 2차 핵 파편이 제공하는 선량이 1차 핵이 제공하는 선량과 비교해도 손색이 없음을 보였다.

상세 요약

MCHIT는 GEANT4 8.2 버전을 기반으로 조직‑유사 매질에서 중이온 치료 빔의 물리적 거동을 정밀하게 시뮬레이션하도록 설계된 Monte Carlo 프레임워크이다. 핵심적인 특징은 투사체(예: 12C) 핵이 물질을 통과하면서 발생하는 파편화 과정을 상세히 모델링하고, 생성된 파편(헬륨, 리튬, 베릴륨 등)의 2차 핵반응까지 추적한다는 점이다. 이를 위해 물리 리스트는 전자·광자 상호작용뿐 아니라 저에너지 핵반응 모델(Bertini, Binary Cascade 등)과 고에너지 파편화 모델(Quantum Molecular Dynamics) 등을 복합적으로 적용하였다.

검증 단계에서는 물과 PMMA(폴리메틸메타크릴레이트) 두 종류의 표준 팬텀에 대해 단일 에너지(예: 290 MeV/u) 탄소 이온 빔을 조사한 실험 데이터를 활용하였다. 깊이‑선량 곡선에서는 브래그 피크 위치와 피크 전후의 선량 구배가 실험과 거의 일치했으며, 이는 입자 손실 메커니즘과 에너지 전이 모델링이 정확함을 의미한다. 또한 전하별 파편 선량을 측정한 결과, 1가~6가 파편에 대한 선량 비율이 실험값과 통계적 오차 범위 내에서 재현되었다. 이는 파편화 단면과 파편의 에너지 스펙트럼이 적절히 구현되었음을 보여준다.

양전자 방출 핵(예: 11C, 10C)의 공간 분포는 PET 기반 치료 검증에 핵심적인데, MCHIT는 이러한 핵의 생성 및 감쇠를 정확히 예측하여 빔 경로와 브래그 피크 위치를 비침습적으로 확인할 수 있게 한다. 또한 다양한 방위각(0°, 30°, 60° 등)에서 측정된 2차 중성자 에너지 스펙트럼을 재현함으로써 방사선 방호 설계와 환자 주변 장비 보호에 필요한 2차 방사선 환경을 정량화한다.

특히 횡단면 선량 프로파일을 상세히 분석한 결과, 브래그 피크 근처 주변부(반경 5–10 mm)에서는 2차 파편이 제공하는 선량이 1차 핵이 제공하는 선량과 비슷한 수준에 도달한다는 사실이 밝혀졌다. 이는 전통적인 치료 계획에서 주로 1차 입자만을 고려하는 경우, 주변 조직에 대한 과소평가 위험이 있음을 시사한다. 따라서 생물학적 효과를 정밀히 예측하려면 파편에 의한 선량 분포와 각 파편의 상대적 생물학적 효율(RBE)을 함께 고려해야 한다.

한계점으로는 GEANT4 물리 리스트의 파편화 단면 데이터가 아직 완전하지 않아 고에너지 영역에서 일부 파편 비율이 실험보다 약간 과소/과대 평가될 수 있다. 또한 모델은 조직의 이질성(예: 뼈와 연부 조직의 경계)을 단순 균일 매질로 가정했기 때문에 실제 환자 해부학에 바로 적용하려면 추가적인 변형이 필요하다. 향후 연구에서는 최신 GEANT4 버전의 개선된 핵반응 모델을 통합하고, 환자‑특이 CT 데이터와 연계한 voxel‑level 시뮬레이션을 수행함으로써 임상 적용성을 높일 수 있을 것이다.

요약하면, MCHIT는 중이온 치료 빔의 물리적 현상을 다중 스케일에서 정밀히 재현하며, 특히 파편에 의한 주변 선량 기여를 정량화함으로써 치료 계획의 정확도와 방사선 방호 설계에 중요한 도구가 될 잠재력을 입증하였다.


📜 논문 원문 (영문)

🚀 1TB 저장소에서 고화질 레이아웃을 불러오는 중입니다...