[논문 리뷰] Nuclear halo of a 177 MeV proton beam in water: theory, measurement and parameterization

B. Gottschalk, Ethan W. Cascio외 2명|arXiv (Cornell University)|2014. 09. 05.

Radiation Therapy and Dosimetry참고 문헌 37인용 수 9

한 줄 요약

이 연구는 수은에서 177 MeV 프로톤 펜치 빔의 핵 홀로를 측정하고 모델링하며, 모델에 의존하는 및 모델에 의존하지 않는 피팅을 통해 코어, 홀로, 어우라 선량 성분을 분리한다. 연구는 홀로 반경이 빔 범위의 약 1/3임을 발견하고, T(w)의 사용을 도전하며, 정확한 펜치 빔 모델링을 위해 브라그 피크 챔버가 불필요하다고 보여주며, 모델에 의존하지 않는 피팅이 더 뛰어난 이식성과 정확도(전역 RMS 오차 9%)를 제공함을 보여준다.

ABSTRACT

The dose distribution of a monoenergetic pencil beam in water consists of an electromagnetic "core", a "halo" from charged nuclear secondaries, and a much larger "aura" from neutral secondaries. These regions overlap, but each has distinct spatial characteristics. We have measured the core/halo using a 177MeV test beam offset in a water tank. The beam monitor was a fluence calibrated plane parallel ionization chamber (IC) and the field chamber, a dose calibrated Exradin T1, so the dose measurements are absolute (MeV/g/p). We performed depth-dose scans at ten displacements from the beam axis ranging from 0 to 10cm. The dose spans five orders of magnitude, and the transition from halo to aura is clearly visible. We have performed model-dependent (MD) and model-independent (MI) fits to the data. The MD fit separates the dose into core, elastic/inelastic nuclear, nonelastic nuclear and aura terms, and achieves a global rms measurement/fit ratio of 15%. The MI fit uses cubic splines and the same ratio is 9%. We review the literature, in particular the use of Pedroni's parametrization of the core/halo. Several papers improve on his Gaussian transverse distribution of the halo, but all retain his T(w), the radial integral of the depth-dose multiplying both the core and halo terms and motivating measurements with large "Bragg peak chambers" (BPCs). We argue that this use of T(w), which by its definition includes energy deposition by nuclear secondaries, is incorrect. T(w) should be replaced in the core term, and in at least part of the halo, by a purely electromagnetic mass stopping power. BPC measurements are unnecessary, and irrelevant to parameterizing the pencil beam.

연구 동기 및 목표

• 177 MeV 프로톤 펜치 빔의 공간 선량 분포를 수은에서 측정하고 모델링하여 코어, 홀로, 어우라 성분을 구분한다.
• 코어 및 홀로를 매개변수화하는 데 관용적으로 사용되는 T(w), 즉 깊이-선량의 반경 적분의 사용을 도전하며, 핵 2차 입자 기여를 잘못 포함한다고 주장한다.
• 25개 매개변수를 사용하는 모델에 의존하는 피팅보다 모델에 의존하지 않는 피팅(삼차 스퍼린)을 사용한 피팅이 더 이식성 있고 정확하다는 것을 보여준다.
• 브라그 피크 챔버가 정확한 펜치 빔 선량 측정에 개념적으로 잘못되었고, 특히 소형 필드에서 불필요하다는 것을 보여준다.
• 몬테카를로 시뮬레이션을 위한 기준점으로 제공하기 위해, 전자기적, 탄성, 비탄성 핵 반응 과정이 홀로 영역의 총 선량에 기여하는 정량적 기여도를 정확히 규명한다.

제안 방법

• 10개의 반경 오프셋(0–10 cm)에서 빔 축에서의 깊이-선량 스캔을 수행하였으며, 선량 측정을 위해 유도력 캘리브레이션된 이온실린 침대와 선량 캘리브레이션된 Exradin T1 침대를 사용하여 절대 선량 측정(MeV/g/p)을 수행하였다.
• 25개 매개변수를 사용하는 모델에 의존하는(MD) 피팅을 통해 선량을 전자기적(코어), 탄성 핵 반응, 비탄성 핵 반응, 어우라 성분으로 분해하였다.
• 깊이 및 반경 방향에 삼차 스퍼린을 사용하는 모델에 의존하지 않는(MI) 피팅을 적용하였으며, 이중선형 보간과 두 번째 도함수 밀도 기반 최적화된 초기 스퍼린 점을 사용하였다.
• 운동량 보존과 질량 중심 산란 각도를 기반으로 최대 홀로 반경을 추정하기 위해 운동역학 분석을 수행하였으며, 이는 빔 범위의 약 1/3과 일치함을 보였다.
• 코어 및 홀로 매개변수화에서 T(w)를 사용하는 대신, 순수하게 전자기적 질량 정지력인 S_em/ρ_water를 사용하여 핵 2차 입자 기여로 인한 오염을 제거하였다.
• 문헌 데이터, 몬테카를로 시뮬레이션(Gate/Geant4), 기존 매개변수화(Pedroni, Zhang, Zhu)와의 비교를 포함하여 결과를 검증하였으며, 브라그 피크 챔버 측정 결과와의 비교도 수행하였다.

실험 결과

연구 질문

• RQ1 177 MeV 프로톤 빔에서 수은 내 핵 홀로의 공간적 범위와 반경 의존성은 무엇이며, 빔 범위와 어떻게 관련되어 있는가?
• RQ2 깊이-선량의 반경 적분인 T(w)를 사용하는 것이 코어 및 홀로 모델링에 문제가 되는 이유는 무엇이며, 더 물리적으로 정확한 대안은 무엇인가?
• RQ3 모델에 의존하는 피팅과 모델에 의존하지 않는 피팅은 코어/홀로 선량 분포를 특성화하는 데 있어 정확도와 이식성 측면에서 어떻게 비교되는가?
• RQ4 브라그 피크 챔버는 펜치 빔 선량 측정에서 어느 정도 오차나 중복을 유발하며, 그 사용을 제거할 수 있는가?
• RQ5 전자기적, 탄성, 비탄성 핵 반응 과정이 홀로 영역의 총 선량에 기여하는 정량적 기여도는 무엇인가?

주요 결과

• 최대 홀로 반경은 탄성 및 준탄성 산란에 의한 운동역학적 제약과 일치하며, 산소 및 수소에서의 탄성 및 준탄성 산란 기반으로 약 빔 범위의 1/3이다.
• 모델에 의존하지 않는 피팅은 측정값과 피팅값 간 전역 RMS 비율이 9%로, 모델에 의존하는 피팅(15% RMS)을 뛰어넘었으며, 개념적으로 더 단순하고 더 이식성이 뛰어나다.
• 코어 및 홀로 매개변수화에 T(w)를 사용하는 것은 잘못된 것으로, 핵 2차 입자에 의해 에너지가 소모되는 기여를 포함하기 때문이다. 이를 S_em/ρ_water(순수 전자기적 정지력)로 대체함으로써 물리적 정확도가 향상된다.
• 브라그 피크 챔버는 펜치 빔 매개변수화에 불필요하며, 비코어 기여를 포함하기 때문에 측정치가 중복되고 때로는 오해의 소지가 있다.
• 홀로는 특히 큰 반경에서 총 선량에 크게 기여하며, 필드 크기 의존성 오차를 피하기 위해 소형 필드 선량 측정에서 고려되어야 한다.
• Pedroni의 원래 매개변수화는 빔 크기에 따라 조정된 후 코어 영역에서 개선된 일치를 보였지만, 홀로 영역에서는 거의 변화가 없었으며, 이는 T(w)-기반 모델의 한계를 드러낸다.