Assessment of the power deposition on the MEGAPIE spallation target using the GEANT4 toolkit

Assessment of the power deposition on the MEGAPIE spallation target using the GEANT4 toolkit 1 2 corresp first-author Corresponding author, abdesslam.lamrabet@usmba.ac.ma Corresponding author, abdesslam.lamrabet@usmba.ac.ma abdesslam.lamrabet@usmba.ac.ma 1 1 1 University Sidi Mohamed Ben Abdellah Faculty of Sciences Dhar El mahraz Laboratory of Integration of Systems and Advanced Technologies, Fes-30000 Morocco. University Sidi Mohamed Ben Abdellah Faculty of Sciences Dhar El mahraz Laboratory of Integration of Systems and Advanced Technologies, Fes-30000 Morocco. Ain Taoujdate B.P 151, El Hajeb Morocco. Ain Taoujdate B.P 151, El Hajeb Morocco. This work aims at evaluating the reliability of the GEANT4 (GEometry ANd Tracking 4) Monte Carlo (MC) toolkit in calculating the power deposition on the Megawatt Pilot Experiment (MEGAPIE), the first liquid-metal spallation target worldwide. A new choice of codes to study and optimize this target is provided. The evaluation of the GEANT4 toolkit is carried out in comparison with the MCNPX and FLUKA MC codes. The MEGAPIE is an international project led by the Paul Scherrer Institute (PSI) in Switzerland. It aims to demonstrate the safe operation of an intense neutron source to power the next generation of nuclear reactors, accelerator driven systems (ADSs). In this study, we used the GEANT4 MC toolkit to calculate the power deposited by fast protons on the MEGAPIE target. The calculation focuses on several structures and regions. The predictions of our calculations were compared and discussed with that of the MCNPX and FLUKA codes, adopted by the MEGAPIE project. The comparison shows that there is a very good agreement between our results and those of the reference codes. GEANT4 Simulation MEGAPIE target Spallation Power deposition 1. Introduction 1 The Megawatt Pilot Experiment (MEGAPIE) is an international project led by the Paul Scherrer Institute (PSI) in Switzerland and supported by many international research institutions, such as the Commissariat à l'énergie atomique in France (CEA), Forschungszentrum Karlsruhe in Germany (FZK), Centre national de la recherche scientifique in France (CNRS), Agenzia Nazionale per le Nuove tecnologie in Italy (ENEA), StudieCentrum voor Kernenergie in Belgium (SCK-CEN), Department of Energy in USA (DOE), Japan Atomic Energy Research Institute (JAERI), and Korea Atomic Energy Research Institute (KAERI), and the European commission [ 1 1 ] . The MEGAPIE initiative was launched in 1999 to design, build, and operate a high-power 1 MW liquid lead–bismuth (PbBi) spallation target [ 2 2 ] . This target is powered by a high-energy proton beam (575 MeV) emitted by an accelerator that consists of three units: the pre-accelerator (Cockroft–Walton; 860 KeV), cyclotron injector (272 MeV), and main ring cyclotron [ 3 3 ] . The power deposited by protons on the target is one of the crucial parameters for the design of this target. In fact, this parameter affects the choice of materials, dimensions, and heat removal. Moreover, the production of neutrons in the target is also related to the rate of power deposition, particularly in the spallation region. This work explores the power deposition in various areas and structures of the MEGAPIE target (Figure 1) using the GEANT4 [ 4 4 - 6 6 4-6 ] MC code. The calculation focuses on structures such as the liquid PbBi metal, main flow guide tube, lower target container, two walls of the lower target enclosure, helium gas of the insulation gap, safety hull D 2 O, and moderator D 2 O. In addition to these structures, two very critical zones are analyzed, that is, the target window and axial area of the spallation zone. Although GEANT4 has many advantages, it is very difficult to handle. Advanced knowledge of C ++ language is necessary for optimal utilization of GEANT4. In addition to the complexity of the code, the geometry and distribution of the primary particles are also complex. In fact, the primary particles constitute a proton beam, which has a nonregular, two-dimensional distribution. In general, the simulation of the MEGAPIE target by the GEANT4 code itself represents a challenge. We used the two MC codes MCNPX [ 7 7 ] and FLUKA [ 8 8 , 9 9 ] as references to evaluate the predictions of our calculations. This choice is based on the conclusions of a benchmark exercise led by the X9 research group of the MEGAPIE project [ 10 10 ] . Based on these conclusions, MCNPX and FLUKA were adopted as reference codes for the project. In addition, most neutronic studies regarding the MEGAPIE target have been carried out with these two codes. 2. Materials and Methods 1 2.1 GEANT4 toolkit 2 The GEANT4 toolkit is a C++ library designed to simulate the passage of particles using the MC method. It is a universal code that is widely used worldwide. It covers several physical areas: nuclear physics, particle physics, astrophysics, and medical physics. The GEANT4 is not a simple code, where the user only manipulates predefined components, but a very flexible toolkit full of features. It provides the users the possibility to fully customize their application. It is always allowed to modify available implementations and/or add others according to the problem and equipment needs. On the other hand, the GEANT4 toolkit is used as basis for several specific codes such as the Gate code [ 11 11 , 12 12 ] and MCADS model [ 13 13 - 16 16 13-16 ] . The Gate code is specifically used for medical applications, whereas the MCADS model, which was developed by the Frankfurt Institute for Advanced Studies in Germany (FIAS), is intended to simulate the spallation targets. In the present simulation, we used GEANT4 version 10.2 [ 17 17 ] . Several improvements were made in this version, in particular, with respect to the compatibility with parallel calculations. The GEANT4 code considers two types of parallelism: clustering [ 18 18 ] and multithreading [ 19 19 ] . Clustering refers to the parallelism with separate memory that requires multiprocessor equipment, while multithreading represents parallelism with shared memory, which requires multicore processors. In this work, we adopted the multithreading technique, which is based on the optimization approach of memory management. Based on this technique, it was possible to reduce the computation time by a factor that is almost equal to the number of cores. Our simulation was performed using an application with many input files. The physics, geometry, materials, source (proton beam), and methods of results extraction in these files are all described using C++ language or command lines. These elements are mentioned in our previous work [ 20 20 ] and explained in more detail in the following sections. To ensure the optimal functioning of the GEANT4 toolkit in this simulation, we built GEANT4 version 10.2 in an environment consisting of the scientific Linux version 7.4 [ 21 21 ] , class library for high energy physics (CLHEP) version 2.3.1.1 [ 22 22 ] , data analysis framework ROOT version 6.12 [ 23 23 ] , and Qt development frameworks version 5.6.2 [ 24 24 ] . 2.2 Modeling physics 2 As a toolkit, GEANT4 does not provide physics models to use by default but provides a wide variety of several models. It is up to the user to choose which model is appropriate for the problem. The choice is made according to the energy field and type of the particle followed. The energy field of the MEGAPIE target is lower than 1 GeV. Therefore, we chose the predefined physics list FTFP_INCLXX_HP to treat the physics processes that may be produced in the MEGAPIE target, which is one of the reference physics lists provided by GEANT4. It was already validated as the most appropriate physics list for this type of problem [ 25 25 ] . This list includes the following models: Fritiof model (FTF) [ 26 26 , 27 27 ] coupled with the precompound one, Liège intranuclear cascade model, C++ version (INCL++ or INCLXX) version 5.2.9.5 [ 28 28 , 29 29 ] , and high-precision model (HP). The FTF is used in GEANT4 to simulate the interactions between hadrons, nucleus, antibaryons, and antinucleus–nucleus. It is suitable for energies higher than 3 GeV. The precompound model is intended to manage the pre-equilibrium emission of protons, neutrons, and light ions [ 30 30 ] . The INCLXX model is used the most for spallation reactions. Its original version INCL [ 31 31 ] was written in FORTRAN77 and translated to C++ such that it can be implemented in GEANT4. The translation was started by Kaitaniemi [ 32 32 ] and finalized by Mancusi [ 29 29 ] . The Liège intranuclear cascade model is designed to govern medium-energy reactions such as spallation reactions. It is permanently maintained and extended to cover more energies. The INCL++ version 5.2.9.5 implemented in GEANT4 10.2.03 has 20 GeV as an upper limit. As an intranuclear cascade model, INCL needs to be coupled with a de-excitation model to describe the two stages of the spallation reactions. Such reactions are assumed to start and finish with an avalanche of independent binary collisions and de-excitation phase, respectively. The nucleus remaining at the end of the first stage relaxes by emitting low-energy particles or by fissioning [ 33 33 ] . To control the processes corresponding to this stage, we coupled INCL++ with the ABLA version 3 statistical de-excitation model [ 34 34 - 36 36 34-36 ] . This model is one of the best de-excitation models based on the International Atomic Energy Agency (IAEA) benchmark of spallation models [ 37 37 , 38 38 ] . As mentioned above, the FTFP_INCLXX_HP physics list also includes the high-precision model HP, which is appropriate for neutron tracking below 20 MeV. Finally, note that it is necessary to add the thermal neutron scattering model [ 39 39 , 40 40 ] to previously mentioned models if one is interested in thermal neutrons. The threshold energy of secondary production (or Cut) is one of the most important physics parameters. It has a great impact on the calculation of the power deposition, especially for narrow volumes. The smaller this energy is, the more accurate is the power deposition. However, this precision comes at a cost. The calculation time increases. Therefore, a judicious choice of this parameter is essential for an optimal calculation within a reasonable time. 2.3 Geometry and materials 2 The MEGAPIE target is a complex device with respect to the geometry and materials. It consists of two parts that include many components ( Figure 1) Figure 1) . The lower part in which the spallation reaction occurs is constituted by a guide tube for the proton beam (AlMg 3 ), double-walled enclosure (AlMg 3 ), safety hull (D 2 O) between the two walls, lower liquid-metal container (martensitic steel T91), guide tube for the main flow of the liquid metal, central rod with microdetectors, bypass duct, and filling and draining tubes. The last five elements are all made of austenitic stainless steel SS316L. The upper part of the target that is dedicated to heat removal mainly consists of two electromagnetic pumps (copper, ferromagnetic), upper enclosure, upper liquid-metal container, heat exchanger, two Diphyl-THT (DTHT) oil boxes, heavy water box, and expansion tank. All these parts are made of stainless steel SS316L. In addition to all previously mentioned elements, there are shielding blocks (tungsten) and the target head with its equipment. Figures of the geometry obtained with GEANT4 were published in previous work [ 19 19 ] . The dimensions of the above-mentioned elements and other details can be found elsewhere [ 41 41 , 42 42 ] . 2.4 Proton Beam 2 In the MEGAPIE target, a monoenergetic proton beam with an energy of 575 MeV and intensity of 1.74 mA induces the spallation reaction. As already mentioned, the target is supplied with this beam by an accelerator with three units. Each unit raises the energy of the protons to a given value until it reaches 575 MeV. The spatial distribution of the proton beam ( Figure 2 Figure 2 ) is a parameter that strongly impacts the functioning of the target. It affects the production of neutrons, power deposition, and life of the materials (mechanical stress). Consequently, the safety of the device is affected by this parameter. On the other hand, the results of the calculations are very sensitive to any change in the distribution of the proton beam. Therefore, a precise description of the proton beam profile is essential to obtain correct results. The definition of the source (proton beam) in the present simulation was a difficult task. The difficulty is due to the complexity of the beam distribution and the lack of a function describing it correctly. To overcome this difficulty, we proceeded with a superposition of many sources, each with its own distribution and probability/intensity. The description was made in a macro file and the implementation was accomplished using a dedicated class, that is, G4GeneralParticleSource. Note that previous trials were made to determine a function that can describe the proton beam profile. Equation (1) [ 10 10 ] represents one of the functions adopted by the scientific community in charge of the MEGAPIE project. However, based on γ mapping and the results obtained from the MCNPX simulation, Equation (1) does not accurately describe the proton beam profile [ 43 43 ] . <math display="block" id="M1"><mrow><msub><mi>I</mi><mrow><mtext>x</mtext><mo>,</mo><mtext>y</mtext></mrow></msub><mrow><mo>[</mo><mrow><mfrac><mrow><mtext>μA</mtext></mrow><mrow><msup><mrow><mtext>cm</mtext></mrow><mn>2</mn></msup></mrow></mfrac></mrow><mo>]</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mn>0</mn></msub></mrow><mrow><mn>2</mn><msub><mrow><mtext>πσ</mtext></mrow><mtext>x</mtext></msub><msub><mtext>σ</mtext><mtext>y</mtext></msub><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mfrac><mrow><mo>−</mo><msup><mi>c</mi><mn>2</mn></msup></mrow><mn>2</mn></mfrac></mrow></msup></mrow><mo>)</mo></mrow></mrow></mfrac><msup><mi>e</mi><mrow><mfrac><mrow><mo>−</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>x</mi><mrow><msub><mtext>σ</mtext><mtext>x</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>y</mi><mrow><msub><mtext>σ</mtext><mtext>y</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow></mrow></msup></mrow></math> 3. Results and Discussion 1 3.1 Power deposited on the target 2 The MEGAPIE target is designed in a way that the maximum power of protons (1 MW) can be dissipated in its lower part, particularly in the spallation region ( Figure 1) Figure 1) . The simulation shows that approximately 85.1% of this power is dissipated in this part of the target. Approximately 71% of the power is deposited in the liquid PbBi metal (actual target), whereas 14% of the power is divided between the confining structures (SS316L, T91, and AlMg 3 ), coolant (D 2 O), and moderator (D 2 O). The power deposited in each component depends on the position, dimensions, and constituent material. Table 1 Table 1 shows the power dissipated in the various elements of the MEGAPIE target. The power values were previously calculated with the FLUKA and MCNPX codes [ 44 44 ] ; the GEANT4 toolkit was used in this work. Table 1 Table 1 indicates that the GEANT4 results are all framed by those of MCNPX and FLUKA, except for the power deposited on the insulation gap "gas_Gap." This difference is quite normal because it is due to the type of gas used in this area of the target. The MCNPX and FLUKA calculations were based on argon (Ar). Conversely, helium (He) was used as gas for the GEANT4 calculation [ 40 40 ] . The fact that all our results are in good agreement with that of the references indicates that the GEANT4 code provides accurate results. This conclusion can mainly be justified by the fact that the version of GEANT4 used is very recent compared with those of FLUKA and MCNPX. Therefore, the GEANT4 code has the potential to provide good results. However, the results need to be compared with experimental data to support this conclusion. 1 1 Material 3 1 Power deposition (kW) 3 1 Relative difference (%) 1 1 1 1 GEANT4 1 1 FLUKA 1 1 MCNPX 1 1 GEANT4/FLUKA 1 1 GEANT4/MCNPX 1 1 MCNPX/FLUKA 1 1 PbBi 1 1 709.7 1 1 710.8 1 1 708.2 1 1 0.15 1 1 0.21 1 1 0.37 1 1 T91 1 1 8.72 1 1 8.60 1 1 8.94 1 1 1.40 1 1 2.46 1 1 3.95 1 1 AlMg 3 _in 1 1 4.56 1 1 4.48 1 1 4.60 1 1 1.79 1 1 0.87 1 1 2.68 1 1 AlMg 3 _ext 1 1 4.19 1 1 4.15 1 1 4.21 1 1 0.96 1 1 0.48 1 1 1.45 1 1 Gas_Gap 1 1 0.00076 1 1 0.012 1 1 0.016 1 1 – 1 1 – 1 1 33.33 1 1 D 2 O 1 1 3.92 1 1 3.83 1 1 3.96 1 1 2.35 1 1 1.01 1 1 3.39 1 1 Moderator 1 1 114.9 1 1 114.4 1 1 115.4 1 1 0.44 1 1 0.43 1 1 0.87 1 1 Guid_Tub 1 1 5.33 1 1 5.55 1 1 5.23 1 1 3.96 1 1 1.91 1 1 5.77 1 1 Total 1 1 851,4 1 1 851.8 1 1 850.6 1 1 0.05 1 1 0.09 1 1 0.14 3.2 Power deposited in the target window 2 The window represents the most critical zone of the MEGAPIE target. It constitutes the interface between the proton beam and target. Thus, it is the most irradiated part of the target. Consequently, the window is the part that is most vulnerable to temperature increase and structural stress. Moreover, it is subjected to mechanical constraints due to the weight of the liquid PbBi metal and the difference in the pressure between the PbBi and helium gas sides. Therefore, the life of the target and its operation safety are directly related to the quantity of the power deposited in its window. Figure 3 Figure 3 shows the distribution of the power deposited in this part calculated by GEANT4 and FLUKA as a function of x at y = 0. Gaussian curves were obtained, with a maximum at the center ( x = 0, y = 0). The agreement between the results of GEANT4 and that of the reference code FLUKA [ 44 44 ] is very clear based on these curves. Note that the difference does not exceed 4% (maximum) and on average is 0.3% ( Figure 4) Figure 4) . Table 2 Table 2 shows a comparison between the maximal power values deposited in the structures of the target window. These data were calculated with GEANT4, FLUKA, and MCNPX. Once again, the values computed by GEANT4 are framed by the reference values. This confirms the ability of GEANT4 to provide results with satisfactory precision. 1 1 Material 3 1 Maximal deposited power (w/cm 3 ) 1 1 1 1 GEANT4 1 1 FLUKA 1 1 MCNPX 1 1 T91 1 1 902,7 1 1 879,5 1 1 924,8 1 1 AlMg 3 _in 1 1 332,6 1 1 321,9 1 1 330,4 1 1 AlMg 3 _ext 1 1 313,6 1 1 309,4 1 1 321,1 3.3 Power deposited along the z-axis of the target 2 The assessment of the power along the z-axis of the target is both of neutronic and thermohydraulic interest. It allows the location of the Bragg peak and therefore the identification of the spallation region, which will be useful for the dimensioning of any reactor powered by a similar target. The determination of the power deposited in the axial zone of the target also allows to consider which heat removal system should be used. Figure 5 Figure 5 represents the map of the volume density of the power deposited in the central region of the target. The radius of this region is r < 0.3 cm; the calculations were made using GEANT4, FLUKA [ 44 44 ] , and MCNPX [ 44 44 ] . The power deposition reaches its maximum of ~1 kW/cm 3 approximately 2 cm from the window. It then decreases with the altitude z up to 26 cm, where a peak can be observed, that is, the Bragg peak. Above 26 cm, the deposited power significantly decreases until it is practically null at approximately 27.5 cm. The curves show that the results are in very good agreement. The curves in Figure 6 Figure 6 representing the local ratios between GEANT4/FLUKA, GEANT4/MCNPX, and MCNPX/FLUKA confirm this agreement. The only observed disparity is related to the position of the Bragg peak. A lag of ~0.3 cm was observed between the FLUKA and GEANT4 and MCNPX results. This difference is mainly due to the physical choices adopted by each code. 4. Conclusion 1 In this paper, we evaluated the power deposition on the MEGAPIE spallation target using a MC method based on the GEANT4 toolkit. We faced many challenges during the accomplishment of this work: 1) the definition of a very complicated geometry; 2) the definition of a two-dimensional and nonregular proton beam; 3) the choice of the suitable physics for the problem; and 4) the development of a program providing correct results. The implemented geometry and proton beam are very similar to that of the references, as indicated in Figures 1 Figures 1 and 2 2 . The physics are governed by the combination of physics models used in this work: FTFP_INCLXX_HP coupled with the ABLA model. The results are in very good agreement with those of the reference codes FLUKA and MCNPX. Furthermore, the difference between our results and that of FLUKA and MCNPX is smaller than that between the FLUKA and MCNPX data. It is remarkable that the GEANT4 predictions are closer to that of the MCNPX code than to that of FLUKA. In general, the differences are mostly due to geometry and proton beam approximations, especially due to physics models and methods of sampling adopted by each code. However, our results are framed by the results of the references, which indicates that the GEANT4 simulation might be optimal. Finally, the GEANT4 toolkit reproduced the geometry and power deposition of MEGAPIE experiment well. Therefore, GEANT4 can be considered to be one of the most reliable and competitive codes for the simulation of great physics experiments. References C. Latge , F. Groeschel , P. Agostini et al ., Megapie spallation target: design, implementation and preliminary tests of the first prototypical spallation target for future ADS. HAL: in2p3-00290482v1 Megapie spallation target: design, implementation and preliminary tests of the first prototypical spallation target for future ADS. HAL: in2p3-00290482v1 G.S. Bauer , M. Salvatores , G. Heusener , MEGAPIE a 1 MW pilot experiment for a liquid metal spallation target . J. Nucl. Mat . 296 , 17 - 35 ( 2001 ). doi: 10.1016/S0022-3115(01)00561-X http://doi.org/10.1016/S0022-3115(01)00561-X MEGAPIE a 1 MW pilot experiment for a liquid metal spallation target C. Fazio , F. Gröschel , W. Wagner et al ., The MEGAPIE-TEST project: Supporting research and lessons learned in first-of-a-kind spallation target technology . Nucl. Eng. Des . 238 ( 6 ) 1471 - 1495 ( 2008 ). doi: 10.1016/j.nucengdes.2007.11.006 http://doi.org/10.1016/j.nucengdes.2007.11.006 . The MEGAPIE-TEST project: Supporting research and lessons learned in first-of-a-kind spallation target technology S. Agostinelli , J. Allison , K. Amako , et al ., GEANT4 - a simulation toolkit , Nucl. Instrum. Meth. A 506 , 250 - 303 ( 2003 ). 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 Trans. Nucl. Sci . 53 , 270 - 278 ( 2006 ). doi: 10.1109/TNS.2006.869826 http://doi.org/10.1109/TNS.2006.869826 GEANT4 developments and applications J. Apostolakis , M. Asai , A.G. Bogdanov et al ., Geometry and physics of the GEANT4 toolkit for high and medium energy applications , Rad. Phy. Chem . 78 , 859 - 873 ( 2009 ). doi: 10.1016/j.radphyschem.2009.04.026 http://doi.org/10.1016/j.radphyschem.2009.04.026 . Geometry and physics of the GEANT4 toolkit for high and medium energy applications L.S. Waters , ed. , MCNPX Users's Manual Version 2.4.0, LA-CP-02-408 , Los Alamos National Laboratory , 2002 . MCNPX Users's Manual Version 2.4.0, LA-CP-02-408 A. Fassò et al ., FLUKA: Status and prospective for hadronic applications , Proc. Monte Carlo 2000 Conf. (2000) , A. Kling , F. Barao , M. Nakagawa et al ., eds, Springer , 955 - 960 ( 2001 ). FLUKA: Status and prospective for hadronic applications A. Fassò , A. Ferrari , P. R. Sala , Electron-photon transport in FLUKA: Status . In: A. Kling , F. Barao , M. Nakagawa et al , eds, Advanced Monte Carlo for Radiation Physics, Particle Transport Simulation and Applications . Springer , Berlin, Heidelberg , 159 - 164 ( 2001 ). doi: 10.1007/978-3-642-18211-2_27 http://doi.org/10.1007/978-3-642-18211-2_27 Electron-photon transport in FLUKA: Status L. Zanini et al ., Summary Report for MEGAPIE R&D Task Group X9(PSI Bericht, 05-12, ISSN 1019-06432005) . Summary Report for MEGAPIE R&D Task Group X9(PSI Bericht, 05-12, ISSN 1019-06432005) G. Santin , D. Strul , D. Lazaro , et al ., GATE: a GEANT4-based simulation platform for PET and SPECT integrating movement and time management , IEEE Trans. Nucl. Sci . 50 ( 5 ), 1516 - 1521 ( 2003 ). doi: 10.1109/TNS.2003.817974 http://doi.org/10.1109/TNS.2003.817974 GATE: a GEANT4-based simulation platform for PET and SPECT integrating movement and time management D. Strulab , G. Santin , D. Lazaro , et al ., GATE: a PET/SPECT general-purpose simulation platform , Nucl. Phys. B . 125 , 75 - 79 ( 2003 ). doi: 10.1016/S0920-5632(03)90969-8 http://doi.org/10.1016/S0920-5632(03)90969-8 GATE: a PET/SPECT general-purpose simulation platform Y. Malyshkin , I. Pshenichnov , I. Mishustin , et al ., Neutron production and energy deposition in fissile spallation targets studied with GEANT4 toolkit . Nucl. Instrum. Meth. B 289 . 79 - 90 ( 2012 ) doi: 10.1016/j.nimb.2012.07.023 http://doi.org/10.1016/j.nimb.2012.07.023 Neutron production and energy deposition in fissile spallation targets studied with GEANT4 toolkit Y. Malyshkin , I. Pshenichnov , I. Mishustin , et al . Modeling spallation reactions in tungsten and uranium targets with the GEANT4 toolkit . EPJ. Web. Conf . 21 , 10006 ( 2012 ). doi: 10.1051/epjconf/20122110006 http://doi.org/10.1051/epjconf/20122110006 . Modeling spallation reactions in tungsten and uranium targets with the GEANT4 toolkit Y. Malyshkin , I. Pshenichnov , I. Mishustin , et al ., Monte Carlo modeling of spallation targets containing uranium and americium , Nucl. Instrum. Meth. B 334 , 8 - 17 ( 2014 ). doi: 10.1016/j.nimb.2014.04.027 http://doi.org/10.1016/j.nimb.2014.04.027 Monte Carlo modeling of spallation targets containing uranium and americium I. Mishustin , Y. Malyshkin , I. Pshenichnov , et al ., Possible Production of Neutron-Rich Heavy Nuclei in Fissile Spallation Targets . In: W. Greiner (eds) Nucl. Phys: Present and Future. FIAS Int. Sci. Series. Springer, Cham , 151 - 161 ( 2015 ) doi: 10.1007/978-3-319-10199-6_15 http://doi.org/10.1007/978-3-319-10199-6_15 Possible Production of Neutron-Rich Heavy Nuclei in Fissile Spallation Targets GEANT4 User's Guide for Application Developers , http://GEANT4.web.cern.ch/GEANT4/UserDocumentation/UsersGuides/ForApplicationDeveloper/BackupVersions/V10.2/html/index.html http://GEANT4.web.cern.ch/GEANT4/UserDocumentation/UsersGuides/ForApplicationDeveloper/BackupVersions/V10.2/html/index.html . GEANT4 User's Guide for Application Developers J. EL Bakkali , T. El Bardouni , Validation of MC GEANT4 code for a 6 MV Varian linac . J. King. Saud. Univ. – Sc . 29 ( 1 ), 106 - 113 ( 2017 ). doi: 10.1016/j.jksus.2016.03.003 http://doi.org/10.1016/j.jksus.2016.03.003 Validation of MC GEANT4 code for a 6 MV Varian linac B. Schmidt , J. González-Dominguez , C. Hundt , et al ., C++11 Multithreading, Parallel Programming, Chap 4, 77–133 (2018) . doi: 10.1016/B978-0-12-849890-3.00004-6 http://doi.org/10.1016/B978-0-12-849890-3.00004-6 C++11 Multithreading, Parallel Programming, Chap 4, 77–133 (2018) A. Lamrabet , A. Maghnouj and J. Tajmouati , GEANT4 Modeling of the International MEGAPIE Experiment . Adv. Stud. Theo. Phys . 11 ( 12 ), 567 - 575 ( 2017 ), doi: 10.12988/astp.2017.773 http://doi.org/10.12988/astp.2017.773 GEANT4 Modeling of the International MEGAPIE Experiment https://www.scientificlinux.org/ https://www.scientificlinux.org/ CLHEP – A Class Library for High Energy Physics , http://proj-clhep.web.cern.ch/proj-clhep/ http://proj-clhep.web.cern.ch/proj-clhep/ . CLHEP – A Class Library for High Energy Physics ROOT a Data analysis Framework - https://root.cern.ch/ https://root.cern.ch/ ROOT a Data analysis Framework Qt | Cross-platform software developement for embedded & desktop - https://www.qt.io/ https://www.qt.io/ Qt | Cross-platform software developement for embedded & desktop S. L. Meo , M.A. Cortés-Giraldo , C. Massimi , et al ., GEANT4 simulations of the n_TOF spallation source and their benchmarking, Eur . Phys. J. A 51 , 160 ( 2015 ). doi: 10.1140/epja/i2015-15160-6 http://doi.org/10.1140/epja/i2015-15160-6 GEANT4 simulations of the n_TOF spallation source and their benchmarking, Eur B. Andersson , G. Gustafson , B. Nilsson-Almqvist , et al ., A model for low-pT hadronic reactions with generalizations to hadron-nucleus and nucleus-nucleus collisions . Nucl. Phys. B 281 , 289 - 309 ( 1987 ). doi: 10.1016/0550-3213(87)90257-4 http://doi.org/10.1016/0550-3213(87)90257-4 A model for low-pT hadronic reactions with generalizations to hadron-nucleus and nucleus-nucleus collisions H. Pi , An event generator for interactions between hadrons and nuclei — FRITIOF version 7.0 . Comp. Phys. Comm . 71 , 173 - 192 ( 1992 ). doi: 10.1016/0010-4655(92)90082-A http://doi.org/10.1016/0010-4655(92)90082-A An event generator for interactions between hadrons and nuclei — FRITIOF version 7.0 A. Boudard , J. Cugnon , J.-C. David , et al ., New potentialities of the Liège intranuclear cascade model for reactions induced by nucleons and light charged particles . Phys. Rev. C 87 , 014606 ( 2013 ), doi: 10.1103/PhysRevC.87.014606 http://doi.org/10.1103/PhysRevC.87.014606 New potentialities of the Liège intranuclear cascade model for reactions induced by nucleons and light charged particles D. Mancusi , A. Boudard , J. Cugnon , et al ., Extension of the Liège intranuclear-cascade model to reactions induced by light nuclei . Phys. Rev. C 90 , 054602 ( 2014 ), doi: 10.1103/PhysRevC.90.054602 http://doi.org/10.1103/PhysRevC.90.054602 Extension of the Liège intranuclear-cascade model to reactions induced by light nuclei J.M. Quesada , V. Ivanchenko , A. Ivanchenko , et al ., Recent developments in pre-equilibrium and de-excitation models in Geant4 , Prog. Nucl. Sci. Technol . 2 , 936 - 941 ( 2011 ). doi: 10.15669/pnst.2.936 http://doi.org/10.15669/pnst.2.936 Recent developments in pre-equilibrium and de-excitation models in Geant4 A. Boudard , J. Cugnon , S. Leray , et al ., Intranuclear cascade model for a comprehensive description of spallation reaction data . Phys. Rev. C 66 , 044615 ( 2002 ). doi: 10.1103/PhysRevC.66.044615 http://doi.org/10.1103/PhysRevC.66.044615 Intranuclear cascade model for a comprehensive description of spallation reaction data P. Kaitaniemi , A. Boudard , S. Leray , et al ., INCL intra-nuclear cascade and ABLA de-excitation models in Geant4 , Prog. Nucl. Sci. Technol . 2 , 788 - 793 ( 2011 ) doi: 10.15669/pnst.2.788 http://doi.org/10.15669/pnst.2.788 INCL intra-nuclear cascade and ABLA de-excitation models in Geant4 D. Filges and F. Goldenbaum , Handbook of Spallation Research. Theory, Experiments and Applications . Wiley-VCH , ( 2009 ). doi: 10.1002/9783527628865 http://doi.org/10.1002/9783527628865 Handbook of Spallation Research. Theory, Experiments and Applications J.J. Gaimard and K.H. Schmidt , A reexamination of the abrasion-ablation model for the description of the nuclear fragmentation reaction . Nucl. Phys. A 531 ( 3-4 ), 709 - 745 ( 1991 ), doi: 10.1016/0375-9474(91)90748-U http://doi.org/10.1016/0375-9474(91)90748-U A reexamination of the abrasion-ablation model for the description of the nuclear fragmentation reaction A.R. Junghans , M. de Jong , H.-G. Clerc , et al ., Projectile-fragment yields as a probe for the collective enhancement in the nuclear level density . Nucl. Phys. A 629 ( 3-4 ) 635 - 655 ( 1998 ). doi: 10.1016/S0375-9474(98)00658-7 http://doi.org/10.1016/S0375-9474(98)00658-7 Projectile-fragment yields as a probe for the collective enhancement in the nuclear level density J. Benlliure , A. Grewe , M. de Jong , et al ., Calculated nuclide production yields in relativistic collisions of fissile nuclei . Nucl. Phys. A 628 ( 3 ), 458 - 478 ( 1998 ). doi: 10.1016/S0375-9474(97)00607-6 http://doi.org/10.1016/S0375-9474(97)00607-6 Calculated nuclide production yields in relativistic collisions of fissile nuclei http://www-nds.iaea.org/spallations http://www-nds.iaea.org/spallations . J.C. David , Spallation reactions: A successful interplay between modeling and applications . Eur. Phys. J. A ( 2015 ) 51 : 68 , doi: 10.1140/epja/i2015-15068-1 http://doi.org/10.1140/epja/i2015-15068-1 Spallation reactions: A successful interplay between modeling and applications A.R. García , E. Mendoza , and D. Cano-Ott , Validation of the Thermal Neutron Physics in geant4 . Department of Energy , Spain , ( 2013 ). Validation of the Thermal Neutron Physics in geant4 https://www-nds.iaea.org/geant4/ https://www-nds.iaea.org/geant4/ A. Cadiou , A. Guertin , T. Kirchner , et al , Final Summary Report on Target Design, 13-47 ( 2004 ), HAL:in2p3-00025250v1 . Final Summary Report on Target Design, 13-47 T. Kirchner et al ., MEGAPIE target design and dimensioning. proceedings of the 4 th MEGAPIE Technical Review Meeting, March 18-19, 2003, Paris, France, report FZKA 6876. HAL:hal-01082506v1 . MEGAPIE target design and dimensioning. proceedings of the 4 th MEGAPIE Technical Review Meeting, March 18-19, 2003, Paris, France, report FZKA 6876. HAL:hal-01082506v1 L. Zanini , H.U. Aebersold , K. Berg , et al ., Neutronic and Nuclear Post-Test Analysis of MEGAPIE, Part I, Chap 2 , PSI Bericht Nr . 08 - 04 , ( 2008 ), ISSN 1019-0643 Neutronic and Nuclear Post-Test Analysis of MEGAPIE, Part I, Chap 2 Y. Foucher , Nuclear Assessment of the MEGAPIE Target , in Proc. 4th MEGAPIE Tech. Rev. Meet., FZK. Sci. Rep. FZKA 6876 , 149 - 161 ( 2003 ). Nuclear Assessment of the MEGAPIE Target The calculations for this simulation were carried out using the national grid of calculation "MaGrid" managed by the National Center of Scientific and Technological Research (CNRST) in Morocco. The authors are grateful to the staff of "MaGrid," in particular Ms. Bouchra RAHIM, a computer engineer, for her availability and assistance with computer work. 10.1007/s41365-019-0590-6 2018-05-02 2018-06-13 2018-06-23 2019-03-12 2019-04-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2019 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-019-0590-6 http://dx.doi.org/10.1007/s41365-019-0590-6 . 2019 Nuclear Science and Techniques 04 30