Neutron Time-Of-Flight spectrometer based on HIRFL for studies of spallation reactions related to ADS project

Neutron Time-Of-Flight spectrometer based on HIRFL for studies of spallation reactions related to ADS project 1 2 first-author 1 corresp Corresponding author, zqchen@impcas.ac.cn Corresponding author, zqchen@impcas.ac.cn zqchen@impcas.ac.cn 1 1 1 2 1 2 1 1 1 2 1 1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China A Neutron Time-Of-Flight (NTOF) spectrometer, based at the Heavy Ion Research Facility in Lanzhou (HIRFL) was developed for studies of neutron production of proton induced spallation reactions related to the ADS project. After the presentation of comparisons between calculated spallation neutron production double-differential cross sections and the available experimental data, a detailed description of the NTOF spectrometer is given. Test beam results show that the spectrometer works well and data analysis procedures are established. The comparisons of the test beam neutron spectra with those of GEANT4 simulations are presented. Time-Of-Flight spectrometer Neutron production cross section Spallation reaction ADS project GEANT4 I. INTRODUCTION 1 The interest in spallation reactions [ 1 1 , 2 2 ] has been increasing with the development of Accelerator Driven Systems (ADS) and the applications of Spallation Neutron Sources (SNS) in many fields, such as the transmutation of nuclear waste, nuclear energy generation, neutron sources for material irradiation, and neutron scattering science. The spallation reaction is defined as interactions between a light projectile (e.g. proton) with GeV range energy and a heavy target nucleus, which is split to a large number of hadrons (mostly neutrons) or fragments. Many nuclear models, such as the intra-nuclear cascade-evaporation (INC/E) [ 3 3 ] model and the quantum molecular dynamics (QMD) model [ 4 4 ], have been developed to study spallation reactions. These nuclear models are coded for thin-target simulations. In order to perform calculations for practical applications, which may involve complex geometries and multiple composite materials, nuclear models need to be embedded into a sophisticated transport code. Nuclear models are combined with a Monte Carlo transportation code, like MCNP [ 5 5 ], GEANT4 [ 6 6 , 7 7 ], and FLUKA [ 8 8 , 9 9 ], nucleon meson transport codes (NMTC) [ 10 10 ] are widely utilized for designing the facilities of engineering applications of the spallation reaction. However, the nuclear models embedded into the transportation codes need to be validated by experimental measurements. Several laboratories [ 11 11 - 14 14 11-14 ] worldwide have constructed specific experimental setups for measuring the double-differential cross sections and the spectra of neutrons produced in proton induced spallation reactions. In their experiments, a time-of-flight technique with organic liquid scintillators (or plastic scintillators) or/and a recoil proton measurement combined with a magnetic spectrometer has been utilized for neutron measurements. Those experimental measurements have made great contributions to the improvement of nuclear models. A Neutron Time-Of-Flight (NTOF) spectrometer, which is based at the Heavy Ion Research Facility in Lanzhou (HIRFL), is designed for further investigation of reaction mechanisms and neutron productions of proton induced spallation reactions related to the ADS project. In this work, the prediction ability of three nuclear models embedded in the GEANT4 for double-differential cross sections of spallation neutron production is examined by comparing the calculated results and the available experimental data [ 15 15 ]. The configuration and test beam results of the NTOF spectrometer are also presented and discussed in detail. II. THEORETICAL CALCULATIONS 1 GEANT4 is a Monte Carlo transport code developed in CERN for simulating the passage of particles through matter. It is widely used in particle and nuclear physics, accelerator physics, medical science, astrophysics, and aerospace studies. In GEANT4, users have abundant choices of physics models to handle the interactions of particles with matter over an extended energy range. Spallation reactions are usually described in two stages: the intra-nuclear cascade and de-excitation. At the first stage, the incident particle transfers its kinetic energy to target nucleons by elastic collisions and a cascade of nucleon-nucleon collisions. Some of the particles that obtained enough energy to escape from the nucleus are emitted. The residual nucleus before the de-excitation stage can emit particles with the low energy, which are called pre-equilibrium particles, before the de-excitation stage. The energies of pre-equilibrium particles are generally greater than the energies of particles emitted during the later stage. In the de-excitation stage, the excited nuclei lose their energy through the evaporation of neutrons, protons, or light charged particles. If the nuclei do not have enough energy to evaporate the particles, it may emit gamma rays and decay to a stable state. In this work, the BIC [ 16 16 ], INCL [ 17 17 ], and QMD [ 18 18 ] models embedded in GEANT4 were used to calculate the first stage of the reactions. A statistical model is used for the second stage. The dynamical model is linked to the de-excitation handler provided by GEANT4. As a result, the double differential cross section of neutron production can be given at each stage. In Fig. 1 Fig. 1 , the experimentally measured neutron production double-differential cross sections from the reaction on the Tungsten target induced by 800 MeV protons at detection angles of 30 to 150 are compared with the simulated results of BIC, INCL, and QMD models embedded in the GEANT4. The experimental data were taken from the EXFOR database [ 15 15 ]. Similar comparisons are given in Fig. 2 Fig. 2 for 800 MeV to 1600 MeV protons at a 30 detection angle. Both of the spectra are multiplied by a factor of 10 -n ( n =0, 2, 4, 6) from top to bottom. Good agreements are achieved among the INCL calculations and experimental measurements in the entire angular and energy range examined. The BIC and QMD calculations slightly underestimate the experimental results in some parts of the neutron spectrum, especially in the high energy tail. These discrepancies are related to the neutron detection angle and the incident proton energy. III. NTOF SPECTROMETER 1 The NTOF spectrometer is capable of measuring neutrons in wide energy and angular ranges. It consists of a beam pick-up detector and 10 individual neutron detection modules. The schematic view of the NTOF spectrometer is shown in Fig. 3 Fig. 3 . The beam pick-up detector is composed of a BC404 plastic scintillator (5 cm×5 cm×0.2 cm detector with a dual-PMT readout at both ends and is located upstream from the target. It is used to obtain the time and position of the incident particle beam. Each of the neutron detection modules is composed of a thin plastic scintillator detector (veto detector) and organic liquid scintillator detector. Two different sizes of cylindrical BC501A organic liquid scintillator detectors are placed at different angles. The larger BC501A detectors have a 12.7 cm diameter and are 12.7 cm in length. They are used for detecting higher energy neutrons with longer neutron flight paths. The smaller BC501A detectors have a size of 5.08 cm in diameter and 5.08 cm in length and they are used to measure low energy neutrons. In front of the individual BC501A detectors, a BC404 plastic scintillator (15 cm×15 cm×0.3 cm) coupled to a 9813KB PMT are mounted as a veto detector to distinguish charged particles from non-charged particles (neutrons and gamma-rays). For the separation of gamma-rays from neutrons, the pulse shape discrimination (PSD) property of the organic liquid scintillator detector is utilized. Neutron kinetic energy is calculated from the time-of-flight (TOF) spectrum between the organic liquid scintillator detector and the beam pick-up detector. The neutron flight path from the target to neutron detector is about 1.5 m. The 400 ns TDC range is used in the experiment for measuring above 0.1 MeV neutrons. A VME-based data acquisition system (DAQ) for recording the experimental data on an event by event basis has been developed. The DAQ software is based on the CERN ROOT framework and runs on a Linux operating system. The main functions of the DAQ include data acquisition, online monitoring and offline data analysis with a graphical user interface (GUI), which provide convenient operation. IV. RESULTS AND DISCUSSION 1 A test experiment for the whole system of the NTOF spectrometer was performed by measuring neutrons, gamma rays, and charged particles in the energy spectrum, production yield, and angular distribution from a Tungsten target with an 16 O beam at a bombarding energy of 400 MeV/u. In the experiment, the combination of the light output spectra of the veto detector and the BC501A scintillator detector was utilized for the separation of charged particle events from non-charged particle (neutron and gamma-ray) ones, as shown in Fig. 4 Fig. 4 . The anode signal of the BC501A scintillator detector was divided into three pulses. Two of them were fed into charge-to-digital converters (QDCs), which had different gate widths (total and slow). QDCs with total and slow gates were used to eliminate gamma-ray events with PSD by the two gate integration method. A typical result of PSD is shown in Fig. 5 Fig. 5 . The TOF spectra were measured between the beam pick-up detector and the BC501A detector. In Fig. 6 Fig. 6 , a typical TOF spectrum, excluding the charged particle events, is given. The prompt gamma peak of the TOF spectrum was used as a time reference of the spectrum. The energy calibration of the organic liquid scintillator detectors was accurately determined by comparing the experimental light output of standard gamma sources with a GEANT4 simulated one [ 19 19 ]. In Fig. 7 Fig. 7 , the comparisons of the experimental neutron production yield at detection angles ranging from 10 to 60 using lager BC501A detectors with GEANT4 results are shown. The neutron production yield is converted from TOF spectra with normalizing by unit solid angle and incident ion numbers. The experimental neutron spectra shape are well reproduced by the simulations. The experimental data have been normalized to the simulated one at 10. The experimental results show that the whole system of the NTOF spectrometer works well, and the data analysis procedure is established. V. CONCLUSION 1 For nuclear engineering facilities, sophisticated simulation tools based on nuclear reaction models are required. In this work, the spallation neutron production cross sections of the BIC, INCL, and QMD reaction models embedded in GEANT4 are compared with the available experimental data. The simulated results with INCL models are in good agreement with the experimental data in the entire energy and angular ranges studied. The predictions of the BIC and QMD models show slight discrepancies with the experimental results. In the design of the ADS project, accurate nuclear data were required. For this reason, the NTOF spectrometer was developed at the Institute of Modern Physics, Chinese Academy of Sciences. A test run for the whole spectrometer system was performed by the experiment of 400 MeV/u 16 O bombarded on Tungsten target. The experimental results show that the experimental apparatus works well and the data analysis method is established. Further experiments will be performed in future for the improvement and validation of the nuclear models embedded in Monte Carlo transportation codes. These transportation codes will be the main simulation tools for the design of the spallation target of ADS project in China. References B Slowinski . Spallation reactions and accelerator-driven systems . Appl Energ , 2003 , 75 : 129 - 136 . DOI: 10.1016/S0306-2619(03)00025-4 http://doi.org/10.1016/S0306-2619(03)00025-4 Spallation reactions and accelerator-driven systems D Ding and S Fu . Modern Physics , 2001 , 13 : 20 - 25 . DOI: xx http://doi.org/xx Modern Physics J Cugnon , C Volent and S Vuillier . Improved intranuclear cascade model for nucleon-nucleus interactions . Nucl Phys A , 1997 , 620 : 475 - 509 . DOI: 10.1016/S0375-9474(97)00186-3 http://doi.org/10.1016/S0375-9474(97)00186-3 Improved intranuclear cascade model for nucleon-nucleus interactions J Aichelin , G Peilert , A Bohnet , et al . Quantum molecular dynamics approach to heavy ion collisions: Description of the model, comparison with fragmentation data, and the mechanism of fragment formation . Phys Rev C , 1988 , 37 : 2451 - 2468 . DOI: 10.1103/PhysRevC.37.2451 http://doi.org/10.1103/PhysRevC.37.2451 Quantum molecular dynamics approach to heavy ion collisions: Description of the model, comparison with fragmentation data, and the mechanism of fragment formation J F Briesmeister . MCNP a general monte carlo n-particle transport code. LA-12625 , 1993 . MCNP a general monte carlo n-particle transport code. LA-12625 S Agostinelli , J Allison , K Amako , et al . GEANT4-a simulation toolkit . Nucl Instrum Meth A , 2003 , 506 : 250 - 303 . DOI: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 GEANT4-a simulation toolkit J Allison , K Amako , J Apostolakis , et al . Geant4 developments and applications . IEEE T Nucl Sci , 2006 , 53 : 270 - 278 . DOI: 10.1109/TNS.2006.869826 http://doi.org/10.1109/TNS.2006.869826 Geant4 developments and applications G Battistoni , F Cerutti , A Fassò , et al . The FLUKA code: description and benchmarking . AIP Conf Proc , 2007 , 896 : 31 - 49 . DOI: 10.1063/1.2720455 http://doi.org/10.1063/1.2720455 The FLUKA code: description and benchmarking A Ferrari , P R Sala , A Fasso , et al . FLUKA: a multi-particle transport code. CERN-2005-010: INFN/TC_05/11; SLAC-R-773 , Geneva , 2005 . FLUKA: a multi-particle transport code. CERN-2005-010: INFN/TC_05/11; SLAC-R-773 W A Coleman and T W Wamstrong . The nucleon-meson transport code NMTC. Technical report, ORNL-4606 , 1970 . The nucleon-meson transport code NMTC. Technical report, ORNL-4606 M M Meier , D A Clark , C A Goulding , et al . Differential neutron production Cross Sections and Neutron Yields from stopping-length targets for 113-MeV protons . Nucl Sci Eng , 1989 , 102 : 310 - 321 . DOI: 10.2172/6455417 http://doi.org/10.2172/6455417 Differential neutron production Cross Sections and Neutron Yields from stopping-length targets for 113-MeV protons K Ishibashi , H Takada , T Nakamoto , et al . Measurement of neutron-production double-differential cross sections for nuclear spallation reaction induced by 0.8, 1.5 and 3.0 GeV protons . J Nucl Sci Tech , 1997 , 34 : 529 - 537 . DOI: 10.1080/18811248.1997.9733705 http://doi.org/10.1080/18811248.1997.9733705 Measurement of neutron-production double-differential cross sections for nuclear spallation reaction induced by 0.8, 1.5 and 3.0 GeV protons X Ledoux , F Borne , A Boudard , et al . Spallation neutron production by 0.8, 1.2, and 1.6 GeV protons on Pb targets . Phys Rev Lett , 1999 , 82 : 4412 - 4415 . DOI: 10.1103/PhysRevLett.82.4412 http://doi.org/10.1103/PhysRevLett.82.4412 Spallation neutron production by 0.8, 1.2, and 1.6 GeV protons on Pb targets Yu V Trebukhovsky , Yu E Titarenko , V F Batyaev , et al . Double-differential cross sections for the production of neutrons from Pb, W, Zr, Cu, and Al targets irradiated with 0.8, 1.0, and 1.6 GeV protons . Phys Atom Nucl+ , 2005 , 68 : 3 - 15 . DOI: 10.1134/1.1858552 http://doi.org/10.1134/1.1858552 Double-differential cross sections for the production of neutrons from Pb, W, Zr, Cu, and Al targets irradiated with 0.8, 1.0, and 1.6 GeV protons B Pritychenko . Exfor Experimental Nuclear Reaction Data Retrievals . http://www.nndc.bnl.gov/exfor/exfor00.htm http://www.nndc.bnl.gov/exfor/exfor00.htm Exfor Experimental Nuclear Reaction Data Retrievals G Folger , V N Ivanchenko and J P Wellisch . The binary cascade . Eur Phys J A , 2004 , 21 : 407 - 417 . DOI: 10.1140/epja/i2003-10219-7 http://doi.org/10.1140/epja/i2003-10219-7 The binary cascade P Kaitaniemi , A Boundard , S Leray , et al . INCL intra-nuclear cascade and ABLA De-Excitation models in Geant4 . Prog Nucl Sci Tech , 2011 , 2 : 788 - 793 . INCL intra-nuclear cascade and ABLA De-Excitation models in Geant4 K Niita , S Chiba , T Maruyama , et al . Analysis of the (N, xN’) reactions by quantum molecular dynamics plus statistical decay model . Phys Rev C , 1995 , 52 : 2620 - 2635 . DOI: 10.1103/PhysRevC.52.2620 http://doi.org/10.1103/PhysRevC.52.2620 Analysis of the (N, xN’) reactions by quantum molecular dynamics plus statistical decay model S Zhang , Z Chen , R Han , et al . Study on gamma response function of EJ301 organic liquid scintillator with GEANT4 and FLUKA . Chinese Phys C , 2013 , 37 : 126003 . DOI: 10.1088/1674-1137/37/12/126003 http://doi.org/10.1088/1674-1137/37/12/126003 Study on gamma response function of EJ301 organic liquid scintillator with GEANT4 and FLUKA 10.13538/j.1001-8042/nst.26.030502 Supported by National Natural Science Foundation of China (No. 11075189), 100 Persons Project (Nos. 0910020BR0 and Y010110BR0) and ADS Project 302 of the Chinese Academy of Sciences (No. XDA03030200) 2014-08-28 2014-10-13 2015-06-20 2015-06-20 2021-04-11 2015 Nuclear Science and Techniques 3 26