李長楷 馬英杰 湯曉斌 謝 芹 耿長冉 陳 達

1（成都理工大學(xué)核技術(shù)與自動化工程學(xué)院 成都 610059）

2（南京航空航天大學(xué)核科學(xué)與工程系 南京 210016）

基于加速器7Li(p,n)反應(yīng)的硼中子俘獲治療中子源的優(yōu)化設(shè)計

李長楷1,2馬英杰1湯曉斌2謝 芹2耿長冉2陳 達2

1（成都理工大學(xué)核技術(shù)與自動化工程學(xué)院 成都 610059）

2（南京航空航天大學(xué)核科學(xué)與工程系 南京 210016）

7Li(p,n)反應(yīng)以中子產(chǎn)額大、反應(yīng)閾能低等優(yōu)點成為硼中子俘獲治療加速器驅(qū)動中子源所用中子反應(yīng)的候選類型之一。本文重點研究了該中子產(chǎn)生反應(yīng)作為加速器驅(qū)動中子源的中子產(chǎn)額及其能譜特性，并對產(chǎn)生的高能中子束流進行慢化，使其滿足BNCT治療要求。首先采用蒙特卡羅程序MCNPX2.5.0模擬加速器7Li(p,n)反應(yīng)過程，得到1.9?3.0 MeV能量入射質(zhì)子的中子產(chǎn)額及其能譜，并詳細研究了質(zhì)子入射能量為2.5 MeV的最佳條件下產(chǎn)生的中子束流特性；進而提出中子束流的慢化設(shè)計方案，并對慢化所得超熱中子束品質(zhì)進行分析研究。模擬計算結(jié)果表明，10 mA流量的2.5 MeV能量入射質(zhì)子所產(chǎn)生的中子束經(jīng)過慢化處理后，可以很好地滿足硼中子俘獲治療的中子束流要求。

硼中子俘獲治療，蒙特卡羅方法，超熱中子，中子慢化

1936年，Locher提出硼中子俘獲法(Boron Neutron Capture Therapy，BNCT)治療惡性腫瘤，其基本原理是利用熱中子和10B發(fā)生熱中子俘獲反應(yīng)：

其反應(yīng)產(chǎn)物7Li和4He均屬于重粒子，具有射程短、高LET(傳能線密度)、高REB(相對生物效應(yīng))等優(yōu)點，配合親腫瘤硼化合物使10B集中于腫瘤組織，可以局部性、選擇性殺死腫瘤細胞，而對正常細胞造成較小傷害[1]。

熱中子由于散射、俘獲等原因無法有效治療位于體內(nèi)的惡性腫瘤。超熱中子(0.4 eV?10 keV)經(jīng)表層組織慢化，到達腫瘤區(qū)域時已成熱中子，且在表層組織能量沉積較少，因此BNCT的理想治療中子束流一般選超熱中子束[2]。

目前，國內(nèi)對BNCT研究所用的中子源為反應(yīng)堆中子源，但由于反應(yīng)堆安裝復(fù)雜、造價和維護費用高及核安全方面的管制，因此，反應(yīng)堆中子源難以在醫(yī)院普及使用。相對于目前廣泛使用的加速器驅(qū)動中子源以安全性高、公眾接受度高、成本相對較低、慢化處理簡單等特點，成為BNCT最具未來應(yīng)用潛力的中子源[3,4]。目前，國內(nèi)關(guān)于BNCT加速器中子源的研究還很少，國際上對BNCT加速器中子源所使用的中子產(chǎn)生反應(yīng)主要包括7Li(p,n)、9Be(p,n)、9Be(d,n)、13C(d,n)以及D-T和D-D裂變反應(yīng)等，其中7Li(p,n)反應(yīng)具有中子產(chǎn)額大、反應(yīng)閾能低、產(chǎn)出中子平均能量低等優(yōu)點，成為首選方案；但由于該反應(yīng)產(chǎn)生的中子束流平均能量較BNCT所需超熱中子能量依然偏高，且存在γ光子污染，因此對病人照射前仍需對中子束流進行慢化處理[5?8]。

本文主要采用蒙特卡羅程序MCNPX2.5.0對BNCT加速器中子源7Li(p,n)反應(yīng)過程進行模擬計算，獲得超熱中子束流的中子產(chǎn)額、能譜等信息；提出中子束流的慢化優(yōu)化設(shè)計方案，并對慢化所得超熱中子束品質(zhì)進行研究，通過跟蹤2×109個粒子的輸運，整個過程計算誤差在5%以內(nèi)，從而證明加速器中子源應(yīng)用于BNCT的可操作性。

1 研究方法

1.1計算模型及計數(shù)方法

考慮到7Li(p,n)反應(yīng)閾能為1.88 MeV，因此本文重點研究1.9?3.0 MeV的入射質(zhì)子轟擊靶核所得的中子能譜。計算模型為：質(zhì)子束為半徑5 cm的柱狀束；靶為圓盤狀厚靶，半徑8 cm，厚0.3 cm，經(jīng)計算驗證該厚度和半徑的厚靶可獲得較大中子產(chǎn)額。MCNPX的計數(shù)方法：在靶外部設(shè)一計數(shù)球面，用Ft卡配合余弦卡將整個計數(shù)球面從質(zhì)子入射方向按角度分為18個子區(qū)間，每個子區(qū)間對應(yīng)的角度為10°，用F1卡對每個子區(qū)間計數(shù)計算出微分中子產(chǎn)額角分布；設(shè)18個楔體將整個計數(shù)球面按面積均分為18個子區(qū)間，對子區(qū)間內(nèi)出射中子通量計數(shù)得出按立體角分布的微分中子產(chǎn)額分布；同時，設(shè)不同能量箱計算出微分中子產(chǎn)額按能量分布，為保證精確度對能量計數(shù)設(shè)3500個子區(qū)間。蒙特卡羅計算所用模型如圖1所示。

圖1 質(zhì)子束轟擊鋰靶計算模型Fig.1 Model of proton beam bombarding lithium target.

為研究7Li(p,n)反應(yīng)產(chǎn)生的高能中子束流經(jīng)過慢化得到的超熱中子束品質(zhì)，本文通過檢測、對比不同快中子慢化材料、慢中子和光子吸收材料以及反射材料后設(shè)計出一種新型慢化裝置方案：快中子慢化和慢中子吸收材料選用fluental材料[9]，此材料是由芬蘭開發(fā)的一種專利材料，設(shè)計初衷是作為反應(yīng)堆BNCT慢化劑，其成分是質(zhì)量分數(shù)為69%的AlF3、質(zhì)量分數(shù)為30%的Al和質(zhì)量分數(shù)為1%的LiF；光子吸收材料選用PbF2。裝置采用2 cm鉛層作為鋰襯底，裝置末端采用9 cm鉛板作為γ屏蔽材料；慢化物質(zhì)存放在2 mm厚的鋁制圓筒里，圓筒內(nèi)半徑為12 cm，外圍由18 cm厚的PbF2反射層包圍；整個裝置外圍由溶鋰重水屏蔽，其中鋰濃度為0.1 g·cm?3。慢化裝置結(jié)構(gòu)如圖2所示。

1.2慢化反射裝置設(shè)計標準[10]

根據(jù)BNCT系統(tǒng)要求，超熱中子場應(yīng)符合以下指標：超熱中子注量率φepi至少要達到109cm?2·s?1以上；超熱中子注量率與熱中子注量率的比值φepi/φth要達到100以上；快中子劑量率與超熱中子注量率的比值Df/φepi不能大于2×10?13Gy·cm2；γ劑量率與超熱中子注量率的比值Dγ/φepi不能大于2×10?13Gy·cm2。

圖2 中子束流慢化/屏蔽裝置結(jié)構(gòu)示意圖Fig.2 Moderating and shielding device.

2 計算結(jié)果與分析

2.1加速器7Li(p,n)反應(yīng)中子束流特性

2.1.1 模擬計算方法的可行性驗證

經(jīng)模擬計算得出2.5 MeV能量入射質(zhì)子產(chǎn)生的中子總產(chǎn)額為0.90×1012n·mC?1；中子最大能量0.782 MeV；中子平均能量0.332 MeV；中子角度區(qū)間為0°?180°；中子平均角度為62.7°。以上各參數(shù)已有的研究數(shù)據(jù)[11]如下：中子總產(chǎn)額為0.89× 1012n·mC?1；中子最大能量0.786 MeV；中子平均能量0.326 MeV；中子角度區(qū)間為0°?180°；中子平均角度為62.9°。對比本文數(shù)據(jù)與已有實驗數(shù)據(jù)看出，本文所用方法與實驗數(shù)據(jù)有較高的吻合度，從而驗證了該方法的可行性。

2.1.2 1.9?3.0 MeV能量入射質(zhì)子產(chǎn)生中子特征分析

為全面獲取7Li(p,n)反應(yīng)產(chǎn)出中子能譜特性，對單位立體角中子產(chǎn)額、單位弧度角中子產(chǎn)額和單位能量中子產(chǎn)額等微分中子產(chǎn)額進行模擬計算。圖3(a)、(b)分別給出單位立體角中子產(chǎn)額和單位弧度角中子產(chǎn)額隨角度的分布(0°≤θ≤180°)，其中，質(zhì)子入射方向為0°。

由圖3(a)，單位立體角中子產(chǎn)額在10°時達到最大，從90°開始曲線明顯變緩，可見在質(zhì)子入射方向上的中子產(chǎn)額面積密度較大。由圖3(b)，單位角度中子產(chǎn)額在40°達到最大，從90°開始變緩。綜合兩圖可見，對于能量大于1.9 MeV的入射質(zhì)子，其產(chǎn)出中子在90°以上的分布還是很可觀的。因此，為充分利用這些“向后”方向的中子，在設(shè)計慢化裝置時可將靶置于慢化材料之中，并在后面再加上反射層。

圖3 不同能量入射質(zhì)子單位立體角中子產(chǎn)額dY/d? (a)和單位弧度角中子產(chǎn)額dY/dθ (b)隨角度分布Fig.3 Different neutron yield distribution dY/d? (a) and dY/dθ (b) as a function of angle for incidence proton.

圖4 為單位能量區(qū)間的中子產(chǎn)額隨中子能量的分布情況。由圖4，隨著入射質(zhì)子能量的升高，高能量的出射中子逐漸出現(xiàn)且額度逐漸增多；當入射質(zhì)子能量大于2.3 MeV時，出射中子能譜有兩個峰，第二個峰形成的原因是7Li(p,n)在2.25 MeV處有個共振峰(圖5)，同時，中子總產(chǎn)額隨入射質(zhì)子能量分布(圖6)中，入射質(zhì)子能量從2.1 MeV提高到2.3MeV，中子產(chǎn)額有個較陡增幅也歸因于此?？紤]利用此共振峰使中子產(chǎn)額較大，以縮短治療時間，且不致使產(chǎn)出中子能量過高，本文最終選擇2.5MeV入射質(zhì)子所產(chǎn)生的中子束流作為實驗對象。

圖4 不同能量入射質(zhì)子微分中子產(chǎn)額dY/dE隨能量分布Fig.4 Different neutron yield distribution dY/dE as a function of energy for incidence proton.

圖5 7Li(p,n)反應(yīng)截面Fig.5 7Li(p,n) reaction cross section.

圖6 不同入射能量質(zhì)子對應(yīng)的中子產(chǎn)額Fig.6 Neutron yield for incidence proton with different energy.

2.1.3 2.5 MeV能量入射質(zhì)子產(chǎn)出中子特性

圖7給出微分中子產(chǎn)額隨能量和角度的分布，由19×3500個數(shù)據(jù)經(jīng)matlab繪出，可以精細反應(yīng)出微分中子產(chǎn)額隨能量和出射角度的分布以及中子出射能量和中子出射角度二者之間的分布關(guān)系；同時，由圖7可見，高能中子的分布角度較小。

圖7 2.5 MeV入射質(zhì)子能量多對應(yīng)的微分中子產(chǎn)額Fig.7 Different neutron yield for incidence proton with energy 2.5 MeV.

2.2中子束慢化模擬計算

由于加速器7Li(p,n)反應(yīng)產(chǎn)生的中子束流還需進行慢化，因此本文提出了一種新型慢化裝置設(shè)計方案，并對慢化裝置出口2 cm、4 cm、6 cm、8 cm處的超熱中子注量率φepi、熱中子注量率φth、快中子劑量率Df以及γ劑量率Dγ進行計數(shù)，慢化材料厚度選擇15?23 cm。

慢化裝置的快中子慢化材料和熱中子吸收材料為fluental材料，γ屏蔽材料為鉛板。分別計算了超熱中子注量率φepi、超熱中子注量率與熱中子注量率比值φepi/φth、快中子劑量與超熱中子注量率比值Df/φepi、γ劑量與超熱中子注量率比值Dγ/φepi隨慢化材料厚度的變化，如圖8所示。

圖8 φepi (a)、φepi/φth (b)、Df/φepi (c)和Dγ/φepi (d)隨慢化材料厚度的變化Fig.8 φepi (a), φepi/φth (b), Df/φepi (c) and Dγ/φepi (d) as a function of moderating material thickness.

由圖8(a)，慢化材料厚度小于21 cm時距出口2?4 cm處的超熱中子輻射場可滿足要求，6 cm、8cm處滿足要求的慢化材料最大厚度分別為17.4 cm和19.4 cm。由圖8(b)，距出口8 cm 處15?23 cm厚度慢化材料均可滿足要求，2 cm、4 cm、6 cm處為滿足要求對應(yīng)的慢化材料厚度區(qū)間分別為16?21 cm、17?20.4 cm、18.2?19.7 cm。由圖8(c)、(d)可以看出，距出口2?8 cm各點及15?23 cm各厚度的慢化材料均可滿足要求。綜上分析，最后選定慢化劑厚度為17 cm，滿足要求的最大超熱中子束通量治療點為距出口4 cm處。

2.3慢化后超熱中子束流技術(shù)指標參數(shù)

經(jīng)計算得出慢化材料厚度17 cm時距出口4 cm處的各項技術(shù)指標參數(shù)為：φepi為1.98×109n·cm?2·s?1；φepi/φth為100；Df/φepi為7× 10?14Gy·cm2；Dγ/φepi為1.4×10?14Gy·cm2。可以很好符合BNCT所需超熱中子束標準。

圖9為慢化后出射中子隨角度的分布。可以看出，中子出射角集中在50°左右。

圖9 慢化后中子出射角分布Fig.9 Neutron emergence angle as a function of different moderator thickness.■ 15 cm, ● 17 cm, ▲ 19 cm, ▼ 21 cm, ? 23 cm

3 結(jié)語

利用蒙特卡羅程序MCNPX的模擬計算，研究了7Li(p,n)反應(yīng)作為加速器驅(qū)動中子源的中子產(chǎn)額及其能譜特性；并提出慢化處理方案，確保產(chǎn)生的超熱中子束滿足BNCT治療要求；從理論上驗證了基于加速器10 mA流量2.5 MeV入射質(zhì)子轟擊鋰靶所得中子束流經(jīng)過慢化可作BNCT中子源使用。7Li(p,n)反應(yīng)閾能較低，通過提高入射質(zhì)子流量可以為加速器小型化提供可能。研究結(jié)果將有助于推動我國基于加速器驅(qū)動硼中子俘獲治療技術(shù)的發(fā)展，可作為未來研制7Li(p,n)反應(yīng)BNCT加速器中子源的重要理論依據(jù)。

1 Blue T E, Yanch J C. Accelerator-based epithermal neutron sources for boron neutron capture therapy of brain tumors[J]. Journal of Neuro-oncology, 2003, 62(1): 19?31

2 Hiroyuki Miyamaru, Isao Murata. Neutron and gamma-ray dose evaluation on accelerator neutrons source using p-Li reaction for BNCT[J]. Progress in Nuclear Science and Technology, 2011, 1: 533?536

3 劉云鵬, 陳達, 湯曉斌, 等. 質(zhì)子放療過程中射束與人眼體作用的蒙特卡羅模擬計算[J]. 原子能科學(xué)技術(shù), 2010, 44: 1?3

LIU Yunpeng, CHEN Da, TANG Xiaobin, etal. Monte-Carlo simulation of proton radiotherapy for human eye[J]. Atomic Energy Science and Technology, 2010, 44: 1?3

4 湯曉斌, 謝芹, 耿長冉, 等. 基于MCNP的超臨界水堆堆芯建模及中子通量計算[J]. 科技導(dǎo)報, 2012, 30(20): 39?43

TANG Xiaobin, XIE Qin, GENG Changran, etal. Core modeling and neutron flux calculation for supercritical water reactor using MCNP[J]. Science & Technology Review, 2012, 30(20): 39?43

5 Tang X B, Geng C R, Xie Q, etal. The simulation of effective dose of human body from external exposure by Monte Carlo methods[J]. First International Workshop on Complexity and Data Mining, 2011: 152?155

6 Verbeke J M, Vujic J L, Leung K N. Neutron beam optimization for boron neutron capture therapy using the DD and DT high-energy neutron sources[J]. Nuclear Technology, 2000, 129(2): 257?278

7 Friedman M, Cohen D, Paul M, etal. Simulation of the neutron spectrum from the7Li(p,n) reaction with a liquid-lithium target at soreq applied research accelerator facility[J]. Nuclear Instruments and Methods in Physics Research, 2013, A698: 117?126

8 Miyamaru H, Murata I. Neutron and gamma-ray dose evaluation on accelerator neutron source using p-Li reaction for BNCT[J]. Journal of Nuclear Science and Technology, 2011, 1: 533?536

9 Allen D A, Beynon T D. A design study for an accelerator-based epithermal neutron beam for BNCT[J]. Physics in Medicine and Biology, 1995, 40(5): 807?809

10 江新標, 陳達, 朱廣寧, 等. 西安脈沖堆超熱中子束的理論設(shè)計[J]. 核動力工程, 2001, 22(3): 199?203

JIANG Xinbiao, CHEN Da, ZHU Guangning, etal. Design of an epithermal neutron beam for Xi'an pulsed reactor[J]. Nuclear Power Engineering, 2001, 22(3): 199?203

11 Lee C L, Zhou X L. Thick target neutron yields for the7Li(p,n)7Be reaction near threshold[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 1999, 152(1): 1?11

CLCTL929

Research of accelerator-based neutron source for boron neutron capture therapy

LI Changkai1,2MA Yingjie1TANG Xiaobin2XIE Qin2GENG Changran2CHEN Da2

1(College of Nuclear Technology and Automation Engineering, Chengdu University of Technology, Chengdu 610059, China)2(Department of Nuclear Science and Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

Background:7Li(p,n) reaction of high neutron yield and low threshold energy has become one of the most important neutron generating reactions for Accelerator-based Boron Neutron Capture Therapy(BNCT). Purpose: Focuses on neutron yield and spectrum characteristics of this kind of neutron generating reaction which serves as an accelerator-based neutron source and moderates the high energy neutron beams to meet BNCT requirements. Methods: The yield and energy spectrum of neutrons generated by accelerator-based7Li(p,n) reaction with incident proton energy from 1.9 MeV to 3.0 MeV are researched using the Monte Carlo code-MCNPX2.5.0. And the energy and angular distribution of differential neutron yield by 2.5-MeV incident proton are also given in this part. In the following part, the character of epithermal neutron beam generated by 2.5-MeV incident protons is moderated by a new-designed moderator. Results: Energy spectra of neutrons generated by accelerator-based7Li(p,n) reaction with incident proton energy from 1.9 MeV to 3.0 MeV are got through the simulation and calculation. The best moderator thickness is got through comparison. Conclusions: Neutron beam produced by accelerator-based7Li(p,n) reaction, with the bombarding beam of 10 mA and the energy of 2.5 MeV, can meet the requirement of BNCT well after being moderated.

Boron neutron capture therapy, Monte Carlo method, Epithermal neutron, Neutron moderator

TL929

10.11889/j.0253-3219.2013.hjs.36.090203

江蘇省科技支撐計劃(1006-KFA13001)資助

李長楷，男，1987年出生，成都理工大學(xué)在讀碩士研究生，從事加速器硼中子俘獲治療研究

2013-07-15，

2013-08-24