# Design and Simulation of a New Model for Treatment by NCT

Authors: Seyed Alireza Mousavi Shirazi; Dariush Sardari
Journal: Science and Technology of Nuclear Installations
 (2012)

---

## Abstract

In this investigation, neutron capture therapy (NCT) through high energy neutrons using Monte Carlo method has been studied. In this study a new method of NCT for a sample liver phantom has been defined, and interaction of 12 MeV neutrons with a multilayer spherical phantom is considered. In order to reach the desirable energy range of neutrons in accord with required energy in absence of eligible clinical neutron source for NCT, this model of phantom might be utilized. The neutron flux and the deposited dose in the all components and different layers of the mentioned phantom are computed by Monte Carlo simulation. The results of Monte Carlo method are compared with analytical method results so that by using a computer program in Turbo-Pascal programming, the deposited dose in the liver phantom has been computed.


---

## Body

## 1. Introduction

Neutron capture therapy (NCT) has been one of the most important methods for treatment of cancers in recent years. This method of radiation therapy is applicable in treatment of liver cancer. During clinical practice, it is always essential to stop absorption of additional dose by normal tissue. On the other hand, measurement and assessment of the absorbed dose and its calibration is an important matter [1, 2]. Thus computation and modeling of the deposited dose by Monte Carlo method before practical treatment is recommended. An appropriate software tool for this purpose is MCNP4C code. It is a particular engineering solution when BNCT facilities such as low energy neutron source are not available.In this paper for simulation by MCNP4C code, a phantom is considered so that it has been encased by polyethylene sphere with 20 cm radius. This sphere is covered with a layer of cadmium which has 100μm thickness [3, 4]. The cadmium layer has high absorption cross-section for thermal neutrons and helps to reentrance the scattered neutrons from surface of the sphere to the phantom [5]. The polyethylene sphere is surrounded by a graphite shell which has 25 cm radius and 5 mm thickness (according to moderation ratio: ξ(Σs/Σa)). This layer serves as a reflector to reduce escaping the fast neutrons [6].In the present work, neutrons are emitted from an external source, and after passing through polyethylene and slowing down, their deposited energy in the phantom’s materials is computed by the MCNP4C code. The F6 tally in the MCNP4C code is applied. The absorbed energy in the liver is computed through analytical computations as well. It includes generation of random numbers along with using the neutron diffusion equation [7]. The outcomes of two methods are finally compared.The results of computation provide assurance that whether the absorbed dose in cancerous and healthy tissues is in accord with requirements [8, 9]. Then neutron therapy would be implemented on the patient. In neutron capture therapy it is suggested to use monoenergetic neutrons. For this purpose, one choice is D-Be source which produces 14 MeV neutrons [10].As Am-Be source produces a wide spectrum of neutron energies, thus, is not suitable for this purpose [11]. Liver tissue includes substances such as water, glycogen, and heavy molecules like protein. Figure 1 shows interaction of incident neutron and depositing the energy during neutron irradiation.Figure 1
The emitted neutron from source and entering into sphere including tissue.


## 2. Monte Carlo Calculations

To simulate neutron therapy and related dosimetry, a phantom is considered as shown in Figure2.Figure 2
The sphere including the liver tissue for the neutron capture therapy.Figure3 depicts the geometry which has been defined by MCNP4C code, and Figure 4 is a schematic view of the physical shape of the phantom.Figure 3
The geometry of the phantom defined by MCNP4C code.Figure 4
The schematic view (3D) of the related geometry.Laboratory experiments have already rendered the precise molecular composition of liver tissue according to Table1 [12].Table 1
The components and the structural materials of a liver tissue (relating to male sex).
Mass percent
Material
69.69%
Water
0.35%
Glycogen (C24H42O21)
29.9%
Protein and Glucose (C44189H71252N12428O14007S321 and C6H12O6)The composition and geometrical data belonging to liver tissue have been given to MCNP4C as input data [13, 14]. In the simulation, neutron slowing down has been taken into account as well [15]. The input data to MCNP4C code are the radius of polyethylene sphere (20 cm), the thickness of cadmium layer (100 μm), radius and thickness of graphite reflector (25 cm and 5 mm, resp.), and the weight of liver tissue that is 190 g.


## 3. Analytical Calculations

It is essential to obtain data on neutron scattering cross-section and its angular distribution because of computing the neutron penetration in liver tissue (or in the corresponding equivalent material) and obtaining the absorbed energy in it. Since neutrons pass through the polyethylene sphere, the knowledge about neutron angular distribution after scattering is vital. Thus, using the random sampling techniques might help to compute the probability of those neutrons which might be scattered in a definite angle [16]. Since it might be comprehended, the atomic composition of the soft tissue is approximated as well with hydrocarbon materials.In the process of neutron transport in matter, it is subjected to three major kinds of interactions with carbon and hydrogen nuclei. These are elastic scattering, inelastic scattering, and radioactive capture. Each kind of interactions takes place with a specific probability with a magnitude proportional to related interaction cross-section according to following equations (1)–(4):
(1)P1=∑H(Elastic)∑tot,P2=∑C(Elastic)∑tot.
The probability of inelastic scattering in carbon (for the first excitation level in: 4.43 MeV):
(2)P3=∑C∑tot.
The probability of inelastic scattering in carbon (for the second excitation level in: 7.65 MeV):
(3)P4=∑C∑tot.
The probability of neutron absorption in hydrogen and carbon:
(4)P5=[∑C+∑H(absorption)]∑tot.
Collision of neutrons on carbon and hydrogen nuclei is caused to transfer some energy from neutrons to the target nucleus. The recoiled nucleus moves a short distance through the matter and deposits energy along its path. The problem is to compute the energy of recoiled nucleus. By using the Monte Carlo method, collision history of each neutron might be tracked down to the energy where it is either absorbed or escapes from the volume of material.The lost energy of neutron in each elastic scattering is described by (5):
(5)E1-E2(cosθ)P(cosθ)d(cosθ).E2 is also obtained by the following equation (6):
(6)E2=E1·A2+2Acosθ+1(A+1)2.
Thus the average energy of recoiled nucleus through assuming the isotropic scattering is
(7)ER¯=2AE1(A+1)2.
The inelastic scattering of neutron in polyethylene is important. There are two excited statuses in carbon nucleus, namely, 4.43 MeV and 7.65 MeV. Thus remarkable value of energy is absorbed in the recoiled nucleus, and ER is less than the energy of inelastic scattered neutrons [17].The transferred energy to a recoiled nucleus with mass numberA due to collision of a neutron with energy En is computed with (8) [18]:
(8)ER=2A(A+1)2Encos2ψ.
Since remarkable percentage of the soft tissue is hydrogen, the interaction of neutron with hydrogen must be studied as well. Approximately 85%–95% of neutron energy transferred to soft tissue is attributed to its interaction with hydrogen [19]. For instance, for En>8 MeV the (α,n) reaction makes the significant fraction of absorbed dose in the tissue. At 14 MeV the contribution of recoiled proton (hydrogen nucleus) is related to energy deposition phenomenon, that is, 66% of total, and the rest is due to α particle plus other heavy nuclei [20, 21].In fact, the penetration of neutron through polyethylene sphere and its interaction in liver tissue require precise knowledge about interaction cross-sections and angular distribution of the scattered neutron.For this purpose by using the Turbo-Pascal programming and considering several values forEn in the interval 0.001–12 MeV and providing the mass number (A) of the components of the phantom as the input to the program, the absorbed dose in the liver tissue (ER) and other components such as polyethylene sphere, air, and collimator are computed [Gy]. The number of neutron collisions is calculated according to (9):
(9)n=ln(En/ER)ξ.
In order to determine absorbed dose in the phantom or other components such as polyethylene, the numeric range of En (energy of incident neutrons) and ER (deposited energy) into the number of neutron collisions (n) is divided, and then its value is added to early energy (10):
(10)ER(new)=ER+En-ERn.
The Turbo-Pascal programming which has been developed in the present work computes the energy transferred from incident neutrons to the tissue based on given neutron energy, scattering angle (ψ), and mass number of the target nuclei.Figure5 shows a schematic view of the phantom which is comprised by four layers. A hypothetical track and slowing down of neutron in consecutive collisions are shown as well.Figure 5
The collision of emitted neutron from the external source to the polyethylene shell.Some useful equations in analyzing the energy transfer of neutron from high energy to the medium are according to (11)-(12) [22]:
(11)ER=Ene-nξ,(12)Lsl2=n(λtr)23,n=3Lsl2λtr2≈∑tr∑sl=∑a+∑s(1-(2/3A))∑sl.
Thus,
(13)n≈∑a+∑s(1-(2/3A))∑sl.
According to (11), (14) is obtained as follows:
(14)ER=Ene-nξ=Ene-((∑a+∑s(1-(2/3A)))/∑sl)×(2/(A+(2/3))).
The variables of A, En, and Σsl are as the inputs for Turbo-Pascal programming so that En is inputted: 12 MeV. Therefore according to (14) the values of absorbed dose in the phantom and other components using Turbo-Pascal programming are obtained.


## 4. Results

The derived graphs for absorbed dose in the liver tissue and other components of mentioned system (per emitted neutrons) by both Monte Carlo simulation with nps: 106 and analytical method (Turbo-Pascal programming) are as Figures 6, 7, 8, 9, 10, 11, and 12.Figure 6
The absorbed dose in the liver tissue.Figure 7
The absorbed dose in the polyethylene sphere.Figure 8
The absorbed dose in the cadmium layer.Figure 9
The absorbed dose in the air of sphere.Figure 10
The absorbed dose in the collimator.Figure 11
The absorbed doses in the liver tissue and other components obtained by MCNP4C.Figure 12
The absorbed doses in the liver tissue and other components obtained by Turbo-Pascal programming.In addition the neutron fluence [n/cm2] or on the other hand the number of neutrons which reach the liver tissue and other components is obtained as Figure 13. Meantime the X-axis means the neutron energies reached the phantom’s components that are low energies [eV].Figure 13
The neutron fluence in the phantom’s components obtained by MCNP4C.


## 5. Conclusion

The absorbed dose by liver tissue in the course of neutron capture therapy has been simulated as function of irradiation time and neutron physical data. This simulation might be carried out for various shapes and size of liver, and the results must be saved in a data bank. For each patient, in accord with given required dose for treatment by the clinic and also knowing the specifications of neutron source, the eligible irradiation time according to both the neutron fluence and energy of existing clinical neutron source through the present simulation is obtained. Meantime the neutron energy reached the eligible interval of energy on the tissue computed by (14).Therefore in case there is no eligible clinical neutron source for NCT, then this model of phantom might be utilized to reach the desirable energy range of neutrons according to the required energy in NCT.From Figures6–10 it is observed that within neutron energy range of 0.001 eV–8 MeV the computed dose by Turbo-Pascal programming is approximately similar to obtained results by MCNP code, and the derived graphs of both methods agree together as well for neutron energy below 8 MeV. Actually, through increasing incident neutron energy, the absorbed dose is increased. For neutron energy above 8 MeV, the results of these two approaches produce significant error. This is because of computation based on MCNP code that computation bases of this code are transport equation and neutron tracking. The equations in Turbo-Pascal programming include diffusion approximations. In addition for neutron energy greater than 8 MeV, (n,α) reactions become more probable to occur. In such nonelastic reactions, the angular distribution of incident neutron and the recoiled nucleus are complicated more than lower energies. This is observed in Figure 7 for polyethylene as well. The absorbed dose in polyethylene phantom is increased through neutron energy so that excess doses may be delivered to the liver tissue. Meantime the absorbed doses because of thermal and epithermal neutron irradiation are negligible.



---

*Source: 213640-2012-09-17.xml*