Abstract
Currently, the dose distribution calculation used by commercial treatment planning systems (TPSs) for high-dose rate (HDR) brachytherapy is derived from point and line source approximation method recommended by AAPM Task Group 43 (TG-43). However, the study of Monte Carlo (MC) simulation is required in order to assess the accuracy of dose calculation around three-dimensional Ir-192 source. In this study, geometry factor was calculated using segmented sources integration method by dividing microSelectron HDR Ir-192 source into smaller parts. The Monte Carlo code (MCNPX 2.5.0) was used to calculate the dose rate $\dot{D}(r,\theta)$ at a point ($r,\theta$) away from a HDR Ir-192 source in spherical water phantom with 30 cm diameter. Finally, anisotropy function and radial dose function were calculated from obtained results. The obtained geometry factor was compared with that calculated from line source approximation. Similarly, obtained anisotropy function and radial dose function were compared with those derived from MCPT results by Williamson. The geometry factor calculated from segmented sources integration method and line source approximation was within 0.2% for $r{\geq}0.5$ cm and 1.33% for r=0.1 cm, respectively. The relative-root mean square error (R-RMSE) of anisotropy function obtained by this study and Williamson was 2.33% for r=0.25 cm and within 1% for r>0.5 cm, respectively. The R-RMSE of radial dose function was 0.46% at radial distance from 0.1 to 14.0 cm. The geometry factor acquired from segmented sources integration method and line source approximation was in good agreement for $r{\geq}0.1$ cm. However, application of segmented sources integration method seems to be valid, since this method using three-dimensional Ir-192 source provides more realistic geometry factor. The anisotropy function and radial dose function estimated from MCNPX in this study and MCPT by Williamson are in good agreement within uncertainty of Monte Carlo codes except at radial distance of r=0.25 cm. It is expected that Monte Carlo code used in this study could be applied to other sources utilized for brachytherapy.
고선량률 근접치료에 사용되는 상업용 선원과 치료계획 시스템들은 AAPM TG 43에서 권고하는 점 및 선 선원에 의해 선량분포를 계산한다. 하지만, 근접치료용 선원에 대한 인체 내의 정확한 선량계산을 위해서 3차원 부피의 선원을 고려하는 MC 기반의 선량계산 방법이 필요하다. 본 연구에서는 microSelectron HDR Ir-192 선원을 작은 부분으로 분할하여 계산하는 미소선원 적분법을 이용하여 기하학적 인수를 계산하였다. 또한, 범용 방사선 수송코드인 MCNPX를 사용하여 30 cm 직경의 구형 물 팬텀 내에서 선원의 선량률을 계산하여 비등방성함수와 반경선량함수를 구하였다. 그 결과를 MC 기반 광자 수송코드인 MCPT를 사용하여 계산한 Williamson의 결과와 비교 및 분석하였다. 미소선원 적분법과 선 선원 근사법에 따른 기하학적 인수는 $r{\geq}0.5cm$에서는 0.2% 이내에서 일치하였고 r=0.1 cm일 때 1.33%의 차이를 보였다. 본 연구에서 계산된 비등방성함수와 반경선량함수가 Williamson의 계산된 결과의 차이는 비등방성함수의 경우 r=0.25 cm에 서 2.33%의 가장 큰 R-RMSE를 보였고 $r{\geq}0.5cm$에서는 1% 미만의 R-RMSE를 보였다. 반경선량함수의 경우는 r=0.1~14.0 cm에서 0.46%의 R-RMSE를 보였다. 미소선원 적분법과 선 선원 근사법으로 계산한 기하학적 인수는 $r{\geq}0.1cm$에서 잘 일치하지만 3차원의 Ir-192 선원을 적용하여 계산한 미소선원 적분법이 실제 기하학적 인수를 잘 반영할 것으로 생각된다. r=0.25 cm에서 비등방성함수를 제외하고는 MCPT와 MCNPX의 몬테칼로 코드를 이용하여 얻어진 비등방성함수와 반경선량함수는 각각의 몬테칼로 코드에 대한 불확실성 이내에서 잘 일치함을 확인하였다. 따라서 MCNPX 전산모사 결과를 통해 TG-43의 선량 계산식에 사용된 인자를 Williamson 등의 결과와 비교 및 검증함으로써, 추후 다른 종류의 선원에 대해서도 Monte Carlo 기반의 연구가 가능할 것으로 기대된다.