NIM B Beam Interactions with Materials & Atoms
Nuclear Instruments and Methods in Physics Research B 261 (2007) 321–325 www.elsevier.com/locate/nimb
Identification of materials hidden inside a container by using the 14 MeV tagged neutron beam Davorin Sudac
b
a,*
, Silvia Pesente b, Giancarlo Nebbia b, Giuseppe Viesti b, Vladivoj Valkovic a
a Rudjer Boskovic Institute, P.O. Box 180, 10002 Zagreb, Croatia Dipartimento di Fisica dell’Universita di Padova, Vai Marzolo 8, I-35131 Padova, Italy
Available online 5 April 2007
Abstract The results of the experiments aiming to confirm the presence of explosive inside the container by using the 14 MeV tagged neutron beam are presented. Measurements were performed with paper, sugar, flour, fertilizer, tobacco and explosive (Semtex1a) as target material placed in the center of an empty container. Additional measurements were done with paper and explosive placed in the center of the container filled with the iron matrix of 0.2 gcm3 density and with the paper target shielded by the 5.1 cm thick iron shield. The results of time of flight measurements and gamma ray spectra obtained by 14 MeV tagged neutron beam have showed that investigated materials could be well distinguished in the triangle plot with coordinates being the number of counts in the carbon peak, the number of counts in oxygen peak and the number of counts in transmitted neutron peak. By using such presentation we have been able to separate paper from Semtex1a, both hidden inside the 0.2 gcm3 iron matrix. We have also been able to confirm the presence of 64.4 kg of paper behind the 5.1 cm thick iron shield corresponding to the range of 300 keV X-rays. 2007 Elsevier B.V. All rights reserved. PACS: 25.40.Fq; 29.40. M; 82.80.Ej Keywords: Explosive detection; Fast neutrons; Container inspection
1. Introduction Method of inspection of cargo containers by using 14 MeV tagged neutron (TN) beam by detection of associated alpha particles has been recently discussed by several groups [1,2]. The aim of the proposed technique is to determine the chemical composition of the suspicious volume within the inspected container, which was previously identified by the X-ray imaging. Such a two sensors technique has been improved with the results promising it to become the powerful tool for container inspection ([3–5]). In the experiment presented, the 14 MeV tagged neutron tech-
*
Corresponding author. Tel.: +385 1 4561 161/1531; fax: +385 1 4680 239. E-mail address:
[email protected] (D. Sudac). 0168-583X/$ - see front matter 2007 Elsevier B.V. All rights reserved. doi:10.1016/j.nimb.2007.03.088
nique has been used to confirm the presence of the explosive Semtex1a within the inspected container. 2. Methods and results Fig. 1 shows the experimental set-up with the investigated target placed in the middle of the container. Target was the iron box dimensions 40 · 40 · 66 cm3 and mass 9.2 kg, filled with the different materials in different measurements: 86.2 kg of graphite, 64.4 kg of paper, 78 kg of sugar, 75 kg of flour, 75 kg of fertilizer, 23.2 kg of tobacco, 100 kg of sand and 100 kg of Semtex1a. Four 300 · 300 NaI(Tl) were put at the transmission position in the cone of the tagged neutrons beam. Each detector was shielded with the 5 cm thick lead. The alpha detector (YAP:Ce) was 8 cm away from the tritium target. A quadratic collimator, 1.8 cm long, was placed in front of YAP:Ce
322
D. Sudac et al. / Nucl. Instr. and Meth. in Phys. Res. B 261 (2007) 321–325 1000
2340
0.1 ns/ch 800
3'' x 3'' NaI (4 peac.)
400
1, 2
400
axis
Target
180
neutron source
1180
3, 4
beam
600 1130
Number of counts/ch
2350
1270
Tagged neutrons
600
400 Gammas from carbon 200
435
floor
0 20
40
60
80
Fig. 1. Experimental set-up inside the container (all units are in mm).
120
140
160
180
300 0.032 MeV/ch 250 1. escape
Number of counts/ch
scintilator. For data processing, a standard NIM electronics was used. The four fast outputs from the 4-segmented PMT connected to the YAP were fed through the octal constant fraction discriminator and fan-in fan-out to the STOP of the time to amplitude converters (TAC). Outputs from the NaI(Tl) detectors were fed through the constant fraction discriminators and ‘‘or’’ logic units to the START of the TACs. Slow signals from the gamma detectors were fed through the amplifiers, quad linear gates and delay amplifiers to the acquisition system. Neutron beam intensity was kept low, 2.5 · 106 n/s in 4p, in order to avoid the problems experienced in previous experiments [6]. Iron matrix was made from iron boxes filled with the iron wire. The axis of the tagged neutron cone was directed through the container approximately parallel with the container bottom (see Fig. 1). Vertical and horizontal profiles of the tagged neutron beam were measured with the 7.62 cm times 7.62 cm NE-213 neutron detector in coincidence with the alpha counter. The experimental diameter of the neutron spot in the middle of the container was (37 ± 2) cm and (35 ± 2) cm for the horizontal and the vertical profile, respectively. As a QC/QA measures the carbon target was used. The 86.2 kg of graphite with the volume of 30 · 30 · 40 cm3 was put in the middle of the container. Fig. 2 shows the corresponding time spectrum (a) and the c-ray spectrum (b) obtained with the appropriate time window and subtracted random background. The main peaks, corresponding to 4.4 MeV excited state as well as the first and the second escape peaks are clearly seen. Fig. 3 shows the time spectra from various materials taken with the same number of alphas. Integral number of the counts under the tagged neutron peaks gives the information on the densities of the various objects (and the surroundings matrices) while the time window set on the gamma peaks give the chemical compositions of the corresponding objects. All measurements with sand, sugar, tobacco, flour, paper, paper behind the 5 cm iron shield and fertilizer were done under the same conditions described in this article. The measurements with the explosive Semtex1a and paper in 0.2 g/cm3 iron matrix were done under the similar conditions, described in [3], and
100
Time (ns)
2. escape
200
4.4 MeV 150
100
50
0 0
1
2
3
4
5
6
7
8
9
10
11
12
Energy (MeV)
Fig. 2. The time spectrum (a) and the gamma ray spectrum (b) of the 86.2 kg of carbon. The gamma ray spectrum was obtained with the appropriate time window on gamma peak. (Target volume = 30 · 40 · 40 cm3, 2 · 108 tagged neutrons, Elapsed time = 17997 s, Average intensity = 0.3 · 107 n/s).
appropriately scaled (scaling factor was found with the help of paper which was measured in both cases). Fig. 4 shows the dependence of integral numbers of TN on density for various materials. Density of the materials was found by dividing the mass with the volume of the iron box. Experimentally obtained data for empty iron box, tobacco, paper, fertilizer, flour, sugar and sand were fitted according to the following formula: N ¼ AeBq
ð1Þ
where A and B are fitting parameters, N is integrated number of the tagged neutrons and q is density. Appropriately scaled Semtex1a nicely follows the exponential fit. Semtex1a in iron matrix, paper in iron matrix and paper behind the iron shield does not follow the exponential fit as was expected. Fig. 5(a) shows the gamma ray spectra obtained for paper and paper behind the iron shield, respectively. Fig. 5(b) shows the gamma ray spectra obtained for sugar and sand, respectively. The gamma ray spectra obtained for flour and tobacco are shown in Fig. 6(a) and gamma
D. Sudac et al. / Nucl. Instr. and Meth. in Phys. Res. B 261 (2007) 321–325 1000
323
300 0.033 MeV/ch
100
250
Carbon
3.3 3.9
Number of counts/ch
Tagged neutrons
Number of counts/ch
Oxygen
paper paper behind the iron shield flour
Gammas from target
Gammas from iron shield
200 4.4 Oxygen 150 5.6 5.1
100
6.1
50
10
0 200
400
600
800
1000
1200
1400
1600
1800
0
1
2
3
5
6
7
8
250
10000
0.033 MeV/ch
0.1 ns/ch
Oxygen
Tagged neutrons
3.8
3.9
Carbon
200
1000
Empty iron box tobacco sugar fertilzer sand
Gammas
100
Numbers of counts/ch
Number of counts/ch
4
Energy (MeV)
Channel
3.3
4.4 Oxygen
150 2.7
5.6 100
5.1
6.1
10
50
1
0 200
400
600
800
1000
1200
1400
1600
1800
0
1
Channel
2
3
4
5
6
7
8
Energy (MeV)
Fig. 3. The time spectra from various materials taken with the same number of alphas.
Fig. 5. The gamma ray spectra from paper (bold black, elapsed time = 39267 s, 4 · 108 tagged neutrons) and paper behind the iron shield (black, elapsed time = 50130 s, 5 · 108 tagged neutrons) (a) and sugar (bold black, elapsed time = 39456 s, 4 · 108 tagged neutrons) and sand (black, elapsed time = 20109 s, 2 · 108 tagged neutrons) (b), respectively.
Integral number of tagged neutrons
1e+6
ray spectrum obtained for fertilizer is shown in Fig. 6(b), respectively. Oxygen and carbon peaks are seen in each spectrum except in the case of fertilizer. In the case of fertilizer the nitrogen, sulphur and phosphorus peaks are visible. In the case of tobacco carbon peaks are visible together with the broad features around 6 MeV which is characteristic for oxygen.
1e+5
1e+4
3. Triangle diagram 1e+3 -0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
Density (g/cm3) Exponential fit empty iron box, tobacco, paper, fertilizer, flour, sugar, sand semtex1a paper behind 5 cm iron shield semtex1a in 0.2 g/cm 3 iron matrix paper in 0.2 g/cm3 iron matrix
Fig. 4. The integral number of the tagged neutrons versus density.
In the c-ray spectra of Semtex1a [3], it can be seen that the nitrogen peak at 5.1 MeV is overlapped with the second escape peak from the 6.1 MeV oxygen line while the peak at the 2.3 MeV is not recognizable at all. On the other side, some home-made explosives do not contain nitrogen, e.g. acetone peroxide (TATP). For that reason we used density as a parameter in triangle diagram instead the nitrogen.
324
D. Sudac et al. / Nucl. Instr. and Meth. in Phys. Res. B 261 (2007) 321–325 250 0.033 MeV/ch
Number of couts/ch
Carbon
3.9
200
4.4
C*
150
O*
Oxygen
5.1
100
5.6 6.1
50
Ro*
0 0
1
2
3
4
5
6
7
8
sand sugar tobacco flour paper fertilizer Semtex1a
Energy (MeV) 300 0.033 MeV/ch 2.2 ;2.3 Nitrogen, Sulfur, Phosphorus
Number of counts/ch
250
Semtex1a in iron matrix paper in iron matrix paper behind the iron shield
200 4.7; 4.6 ... Sulfur, Phosporus 150 1.escape
Fig. 7. Triangle diagram.
5.8; 6.0 Sulfur
100
50
0 0
1
2
3
4
5
6
7
8
Energy (MeV)
Fig. 6. The gamma ray spectra from flour (bold black, elapsed time = 41558 s, 4 · 108 tagged neutrons) and tobacco (black, elapsed time = 42230 s, 4 · 108 tagged neutrons) (a) and fertilizer (elapsed time = 41345 s, 4 · 108 tagged neutrons) (b), respectively.
Fig. 7 shows the triangle diagram. Axes are defined as: 2
Ro ¼ 2 N N
C ¼ 2 N N
2 2 þ CC þ OO 2 C C
þ
N N
2 C C
þ
2 O O
O ¼ 2
N N
þ
2 C C
þ
2
ð2Þ
O O
2 O O
where N is the average integral number of tagged neutrons obtained for sand, sugar, tobacco, flour, paper, Semtex1a and fertilizer, C is the sum of the 4.4 MeV main peak and his first escape. O is the sum of the 6.1 MeV peak, his first and the second escape. C and O are average values. The vertical axes in triangle diagram is Ro* and the horizontal axes is 0.57735 (1 + C* O*). Each peak used in
the calculations is found to be over the critical limit defined as [7]: rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi m ð3Þ Lc ¼ 1:645 B 1 þ 2 where B is the background under the peak and m is the number of channels covered by the peak. Most of the peaks used in the calculations are over the detection limit defined as: rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi m LD ¼ 2:71 þ 3:29 B 1 þ ð4Þ 2 4. Conclusions Semtex1a and Semtex1a in iron matrix are well separated in the triangle diagram from the other materials, even from the paper which has similar carbon to oxygen ratio. Semtex1a and Semtex1a in iron matrix are placed on the same positions in triangle diagram, as well as the paper (paper in iron matrix is a little bit off probably because of the scaling). This means that at the least for the presented cases the surrounding matrices do not influence the position of the point in the triangle diagram. Better statistic will improve the results and lower the error bar, which is important in order to identify various materials. Even with the limited statistics, the paper behind the 5 cm iron shield was successfully identified, which can not be done with the commercial 300 keV X-ray system. And finally, it was showed that with the presented set-up and the neutron beam intensity of 108 n/s, 100 kg of Semtex1a behind the 5 cm iron shield can be found in 20 min.
D. Sudac et al. / Nucl. Instr. and Meth. in Phys. Res. B 261 (2007) 321–325
Acknowledgement This work was supported by NATO Science Programme-Science for Peace Project Sfp – 980526, ‘‘Control of Illicit Trafficking in Threat Materials’’. References [1] S. Pesente et al., Nucl. Instr. and Meth. A 531 (2004) 657. [2] S. Blagus, D. Sudac, V. Valkovic, Nucl. Instr. and Meth.B 213 (2004) 434.
[3] [4] [5] [6]
325
D. Sudac, V. Valkovic, in: Proceedings of SPIE, Vol. 6204-02, 2006. B. Perot et al., in: Proceedings of SPIE, Vol. 6213-06, 2006. A.V. Kuznetsov et al., in: Proceedings of SPIE, Vol. 6213-21, 2006. D. Sudac, S. Blagus, V. Valkovic, Proceedings of the Sixth Topical Meeting on Industrial Radiation and Radioisotope Measurement Applications IRRMA6. [7] G. Gilmore, J. Hemingway, Practical Gamma Ray Spectrometry, John Wiley & Sons.