No abstract

The new version (INCL4.6) of the Liège intranuclear cascade (INC) model for the description of spallation reactions is presented in detail. Compared to the standard version (INCL4.2), it incorporates several new features, the most important of which are: (i) the inclusion of cluster production through a dynamical phase space coalescence model, (ii) the Coulomb deflection for entering and outgoing charged particles, (iii) the improvement of the treatment of Pauli blocking and of soft collisions, (iv) the introduction of experimental threshold values for the emission of particles, (v) the improvement of pion dynamics, (vi) a detailed procedure for the treatment of light-cluster induced reactions taking care of the effects of binding energy of the nucleons inside the incident cluster and of the possible fusion reaction at low energy. Performances of the new model concerning nucleoninduced reactions are illustrated by a comparison with experimental data covering total reaction cross-sections, neutron, proton, pion and composite double differential cross-sections, neutron multiplicities, residue mass and charge distributions, and residue recoil velocity distributions. Whenever necessary, the INCL4.6 model is coupled to the ABLA07 de-excitation model and the respective merits of the two models are then tentatively disentangled. Good agreement is generally obtained in the 200 MeV-2 GeV range. Below 200 MeV and down to a few tens of MeV, the total reaction cross section is well reproduced and differential cross sections are reasonably well described. The model is also tested for light-ion induced reactions at low energy, below 100 MeV incident energy per nucleon. Beyond presenting the update of the INCL4.2 model, attention has been paid to applications of the new model to three topics for which some particular aspects are discussed for the first time. The first topic is the production of clusters heavier than alpha particle. It is shown that the energy spectra of these produced clusters are consistent with coalescence. The second topic regards the longitudinal residue recoil velocity and its fluctuations. Excellent results are obtained for these quantities. It addition, it is shown that the distributions of these quantities display typical random-walk characterics, at least for not too large mass losses. They are interpreted as a direct consequence of the independence of successive binary collisions occurring during the cascade process. The last topic concerns the total reaction cross section and the residue production cross sections for low energy incident light ions. It is shown that our new model can give a rather satisfactory account of these cross sections, offering so an alternative to fusion models and the advantage of a single model for the progressive change from fusion to pre-equilibrium mechanisms.

The purpose of this paper is twofold. First, we present the extension of the Liège intranuclear-cascade model to reactions induced by light ions. We describe here the ideas upon which we built our treatment of nucleus-nucleus reactions and we compare the model predictions against a vast set of heterogeneous experimental data. In spite of the discussed limitations of the intranuclear-cascade scheme, we find that our model yields valid predictions for a number of observables and positions itself as one of the most attractive alternatives available to GEANT4 users for the simulation of light-ion-induced reactions. Second, we describe the C++ version of the code, which is physicswise equivalent to the legacy version, is available in GEANT4, and will serve as the basis for all future development of the model.

We present a statistical-model description of fission, in the framework of compound-nucleus decay, which is found to simultaneously reproduce data from both heavy-ion-induced fusion reactions and proton-induced spallation reactions at around 1 GeV. For the spallation reactions, the initial compound-nucleus population is predicted by the Liège intranuclear cascade model. We are able to reproduce experimental fission probabilities and fission-fragment mass distributions in both reactions types with the same parameter sets. However, no unique parameter set was obtained for the fission probability. The introduction of fission transients can be offset by an increase of the ratio of level-density parameters for the saddle-point and ground-state configurations. Changes to the finite-range fission barriers could be offset by a scaling of the Bohr-Wheeler decay width as predicted by Kramers. The parameter sets presented allow accurate prediction of fission probabilities for excitation energies up to 300 MeV and spins up to 60h.

It is a well-established fact that intranuclear-cascade models generally overestimate the cross sections for one-proton removal from heavy, stable nuclei by a high-energy proton beam, but they yield reasonable predictions for one-neutron removal from the same nuclei and for one-nucleon removal from light targets. We use simple shell-model calculations to investigate the reasons for this deficiency. We find that a refined description of the neutron skin and of the energy density in the nuclear surface is crucial for the aforementioned observables, and that neither ingredient is sufficient if taken separately. As a by-product, the predictions for removal of several nucleons are also improved by the refined treatment.

International audience$^{12}$C ion therapy has had growing interest in recent years for its excellent dose conformity. However at therapeutic energies, which can be as high as 400 MeV/u, carbon ions produce secondary fragments. For an incident 400 MeV/u 12C ion beam, $\sim 70 \%$ of the beam will undergo fragmentation before the Bragg Peak. The dosimetric and radiobiological impact of these fragments must be accurately characterised, as it can result in increasing the risk of secondary cancer for the patient as well as altering the relative biological effectiveness. This work investigates the accuracy of three different nuclear fragmentation models available in the Monte Carlo Toolkit Geant4, the Binary Intranuclear Cascade (BIC), the Quantum Molecular Dynamics (QMD) and the Liege Intranuclear Cascade (INCL++). The models were benchmarked against experimental data for a pristine 400 MeV/u $^{12}$C beam incident upon a water phantom, including fragment yield, angular and energy distribution. For fragment yields the three alternative models agreed between $\sim 5$ and $\sim 35 \%$ with experimental measurements, the QMD using the “Frag” option gave the best agreement for lighter fragments but had reduced agreement for larger fragments. For angular distributions INCL++ was seen to provide the best agreement among the models for all elements with the exception of Hydrogen, while BIC and QMD was seen to produce broader distributions compared to experiment. BIC and QMD performed similar to one another for kinetic energy distributions while INCL++ suffered from producing lower energy distributions compared to the other models and experiment

Neutron production and transport in the spallation target of the n TOF facility at CERN has been simulated with GEANT4. The results obtained with different models of high-energy nucleon-nucleus interaction have been compared with the measured characteristics of the neutron beam, in particular the flux and its dependence on neutron energy, measured in the first experimental area. The best agreement at present, within 20% for the absolute value of the flux, and within few percent for the energy dependence in the whole energy range from thermal to 1 GeV, is obtained with the INCL++ model coupled with the GEANT4 native de-excitation model. All other available models overestimate by a larger factor, of up to 70%, the n TOF neutron flux. The simulations are also able to accurately reproduce the neutron beam energy resolution function, which is essentially determined by the moderation time inside the target/moderator assembly. The results here reported provide confidence on the use of GEANT4 for simulations of spallation neutron sources.

Total fission cross sections of 181 Ta induced by protons at different relativistic energies have been measured at GSI Darmstadt using the inverse kinematic technique. These data contribute to solve inconsistencies in previously reported measurements, but also help to benchmark state-of-the-art reaction codes. The energy range covered with these measurements allowed us to investigate the onset and temperature dependence of dissipative and transient effects at small deformation.