Abstract. We study primary and secondary reactions induced by 600 MeV proton beams in monolithic cylindrical targets made of natural tungsten and uranium by using Monte Carlo simulations with the Geant4 toolkit [1][2][3]. Bertini intranuclear cascade model, Binary cascade model and IntraNuclear Cascade Liège (INCL) with ABLA model [4] were used as calculational options to describe nuclear reactions. Fission cross sections, neutron multiplicity and mass distributions of fragments for 238 U fission induced by 25.6 and 62.9 MeV protons are calculated and compared to recent experimental data [5]. Time distributions of neutron leakage from the targets and heat depositions are calculated.
MotivationOne of the most important applications of high-energy accelerators is to produce neutrons in protoninduced spallation reactions on extended targets made of heavy materials. In this way neutron spallation sources operate, see, e.g., [6] and [7], providing neutrons for neutron scattering experiments and material research. Neutrons produced by energetic protons can be also used to maintain a chain of nuclear fission reactions in a subcritical assembly of an accelerator-driven system (ADS) for nuclear waste incineration [8,9]. In all such cases the neutron flux around the spallation target has to be thoroughly quantified. In particular, the time dependence of the neutron flux on the surface of the spallation target is of a special concern.
SimulationsNeutron production by 600 MeV proton beam in monolithic non-fissile and fissile cylindrical targets are studied my Monte Carlo Simulations. A dedicated software called MCADS (Monte Carlo for Accelerated Driven Systems) was created in FIAS. It is based on the Geant4 toolkit, which is widely used nowadays in nuclear and particle physics.We consider spallation targets made of natural tungsten, uranium and also of 238 U. The latter material was used in validating Geant4 fission models. The radius of these cylindrical targets was taken as 10 cm and the length as (1) 1 cm for the fission studies making multiple interactions of beam protons improbable, and (2) 20 cm for the neutron yields and energy deposition calculations ensuring that all beam protons are absorbed inside the target.Visualization of a history of one proton and all secondary particles is shown on Fig. 1.