放射線治療利用醫用直線加速器產生高能光子射束殺死腫瘤細胞，同時高能光子也與機頭高原子序物質經由巨偶極共振(giant dipole resonance)產生光中子反應(□,n) (photoneutron reaction)，所產生的中子能量分佈從熱中子至快中子能量均有；高能光子射束伴隨中子汙染，使病人接受不必要的劑量。
利用TLD600與TLD700之雙TLD方法搭配鎘差法可以在光子、中子混合場中，分辨光子劑量與熱中子劑量，因此一般常將雙TLD放置於Rando假體以評估醫用直線加速器產生之熱中子劑量。實際上醫用直線加速器所產生之中子為熱中子至快中子皆有分佈，在假體內因空間有限無法放置鎘片，使雙TLD方法測得中子劑量有一部份來源為快中子貢獻而造成誤差。
本研究利用多球體搭配金箔測量醫用直線加速器照射不同假體深度下之中子能譜，利用此中子能譜以蒙地卡羅軟體計算出快中子劑量與熱中子劑量佔所有中子劑量之比例，將此修正因子用於實際量測值可以發現，在假體表面、4公分、10公分、15公分、20公分之誤差分別為99%、45%、21%、12%、12%，從結果可以發現在假體表面之中子能譜主要以快中子能譜為主，故此時以雙TLD方法(不搭配鎘差法)評估之熱中子劑量誤差最大，而在假體較深處，隨著中子進入假體產生作用使中子緩速，中子能量越來越接近熱中子，誤差也越來越小。

In radiation therapy, the high energy photon beams generated by the medical linear accelerators (Linac) were usually used to kill the tumors. Thorough the giant dipole resonances, (□,n) photoneutron reactions occurred by the high energy photons interacted with materials of high atomic number.. The energy distribution of induced photoneutrons is from thermal to fast energy. The neutrons induced by high energy photons increase the doses of patients.
Dual-TLD method coupling with Cadmium-difference method is usually used to distinguish the doses of thermal neutrons form photons in a mixed radiation field. TLD-600 and TLD-700 were easily used to put inside a Rando phantom to investigate the dose contributed by the thermal neutrons induced by a high energy photon beam. Due to the hole space inside the Rando phantom is small to fill a set of dual-TLD with Cd cover, in this case, the estimated neutrons would include neutrons with energies of thermal and above thermal, so the estimated thermal neutron doses would include a contribution for some fast neutrons.
In this study, a system of Bonner sphere spectrometer with gold foils was used to measure the neutron spectra (at different depth of phantom) induced by a Linac (15 MV photon beam). By means of the measured neutron spectra and the Monte Carlo method, the dose ratios contributed by thermal and fast neutrons were assessed and relative correction factors were calculated. In the results, the dose errors of the thermal neutrons due to not subtracting the doses contributed by the higher energy neutrons were 99 %, 45 %, 21 %, 12 % and 12 % for 0, 4, 10, 15 and 20 cm depth in phantom, respectively. On the surface of the phantom, the fast neutrons dominated in the neutron spectrum, therefore, dose of thermal neutrons estimated by dual-TLD method (without using cadmium-difference method) had a maximum dose error. In deeper depths, due to large amount of the fast neutrons were slowing down, estimated dose errors of thermal neutrons decreased.

1. Naseri, A. and A. Mesbahi, “A review on photoneutrons characteristics in radiation therapy with high-energy photon beams.”, Reports of Practical Oncology & Radiotherapy, 2010. 15(5): p. 138-144
2. International atomic energy agency(IAEA), “Handbook on photonuclear data for applications Cross-sections and spectra”, (2000) IAEA-TECDOC-1178
3. Vega-Carrillo, H.R., et al., “Neutron spectrum and doses in a 18 MV LINAC”,. Journal of Radioanalytical and Nuclear Chemistry, 2010. 283(1): p. 261-265.
4. Hsu, F.Y., et al., “Estimation of photon and neutron dose distributions in the THOR BNCT treatment room using dual TLD method.”, Radiation Measurements, 2008. 43(2-6): p. 1089-1094.
5. Hsu, F.-Y., et al., “Dose estimation of the neutrons induced by the high energy medical linear accelerator using dual-TLD chips.”, Radiation Measurements, 2010. 45(3-6): p. 739-741.
6. Esposito, A., et al., “Determination of the neutron spectra around an 18MV medical LINAC with a passive Bonner sphere spectrometer based on gold foils and TLD pairs.“, Radiation Measurements, 2008. 43(2-6): p. 1038-1043.
7. Pena, J., et al., “Monte Carlo study of Siemens PRIMUS photoneutron production.”, Physics in Medicine and Biology, 2005. 50(24): p. 5921-5933.
8. 陳俊丞，“醫用加速器中子汙染特性研究”，國立清華大學，博士論文，民國94年。
9. T. Bonner, M., et al.,”A new type of neutron spectrometer”, Nuclear Instruments and Methods,1960, 9(1):p.1-12
10. Alevra A. V. , et al., “Neutron spectrometry in mixed fields: multisphere spectrometers.”, Radiation Protection Dosimetry ,2003.107(1-3):p.37-72
11. Bedogni, R., et al., “Determination and validation of a response matrix for a passive Bonner sphere spectrometer based on gold foils.”, Radiation Measurements, 2008. 43(2-6): p. 1104-1107.

12. Rogus Ronald D., et al., ”Mixed field dosimetry of epithermal neutron beams for boron neutron capture therapy at the MITR-II research reactor.”, Medical Physics, 1994 21(10):p.1611-1625
13. Marcel Reginatto, Paul Goldhagen, “MAXED, A computer core for maximum entropy deconvolution of multisphere neutron spectrometer
14. Roberts, N.J., “Investigation of combined unfolding of neutron spectra using the UMG unfolding codes”. Radiation Protection Dosimetry, 2007. 126(1-4): p. 398-403.
15. Mares Vladimir, et al., “Evaluation of response matrix of a Bonner sphere spectrometer with LiI detector from thermal energy to 100 MeV.”, Nuclear Instruments and Methods in Physics Research A,1994 (337)p.461-473
16. 張瑜珊，”清華大學開放式水池反應器隻超熱中子束應用於硼中子捕獲治療之混核輻射場劑量研究”，國立清華大學，碩士論文，民國95年
17. Králík, M., et al., “Measurement with Bonner spheres spectrometer in pulsed neutron fields.”, Radiation Measurements, 2010. 45(10): p. 1245-1249.
18. Oongaro Carla , et al., “Analysis of photoneutron spectra produced in medical accelerators”, Physics in Medicine and Biology, 2000. (45): L55-61.
19. 朱葦翰，“多球體偵檢器度量醫用直線加速器環境中子能譜與劑量研究”，國立清華大學，碩士論文，民國99年。