A description of the IntraNuclear Cascade (INC), preequilibrium, evaporation, fission, coalescence, and Fermi breakup models used by the latest versions of our CEM03.03 and LAQGSM03.03 event generators is presented, with a focus on our most recent developments of these models. The recently developed "S" and "G" versions of our codes, that consider multifragmentation of nuclei formed after the preequilibrium stage of reactions when their excitation energy is above 2A MeV using the Statistical Multifragmentation Model (SMM) code by Botvina et al. ("S" stands for SMM) and the fission-like binary-decay model GEMINI by Charity ("G" stands for GEMINI), respectively, are briefly described as well. Examples of benchmarking our models against a large variety of experimental data on particle-particle, particle-nucleus, and nucleus-nucleus reactions are presented. Open questions on reaction mechanisms and future necessary work are outlined.

This paper presents part of an internal LANL Progress Report on LAQGSM03.03, an upgrade of the Los Alamos version of the Quark-Gluon String Model event generator for MCNPX/6 and MARS15 transport codes and on its validation and testing against a large variety of recent measurements. We present here an analysis with LAQGSM03.03 of the recent PHENIX mid-rapidity spectra of pi+, pi-, K+, K-, p, and p-bar produced in ultra-relativistic p + p interactions at sqrt(s) = 200 GeV; GSI cross sections for the fragmentation of Pb208 at 1 GeV/nucleon on Be9; fragmentation cross sections of Si28 on H, C, Al, Cu, Sn, and Pb at energies from 290 to 1200 MeV/nucleon measred recently at HIMAC and BNL; recent HIMAC data on B, Be, Li, and He production cross sections from fragmentation of C12 on H, C, Al, Cu, Sn, and Pb at 290 and 400 MeV/nucleon; BNL data on fragmentation cross sections of Fe56 on H, C, Al, Cu, and Pb targets at 1.05 GeV/nucleon; recent pi+ and pi- spectra from 6.4, 12.3, and 17.5 GeV/c p + Be9 from the E910 BNL measurements; and fragmentation cross sections of Ca40, Ca48, Ni58, and Ni64 on Be9 and Ta181 at 140 MeV/nucleon, and of Kr86 at 64 MeV/nucleon on the same targets measured recently at NSCL-MSU and RARF-RIKEN, respectively.

An improved version of the Cascade-Exciton Model (CEM) of nuclear reactions contained in the code CEM2k and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) are merged with the well-known sequential-binary-decay model GEMINI by Charity. We present some results on proton-induced fragmentation, fission-product yields and on particle spectra predicted by these extended versions of CEM2k and LAQGSM. We show that merging CEM2k and LAQGSM with GEMINI allows us to describe many fission and fragmentation reactions in addition to the spallation and evaporation reactions which are already described well by these codes. Nevertheless, the current version of GEMINI we use does not provide a completely satisfactory description of some complex-particle spectra, fragment emission, and spallation yields for some reactions, and is not yet a universal tool for applications. Our results show that GEMINI contains a powerful model to describe evaporation/fission/fragmentation reactions and often provides better results when compared to other models, especially for emission of heavy fragments from reactions on medium-heavy nuclei (where most other models simply fail), but it must be further extended and improved in order to properly describe arbitrary reactions.

The improved Cascade-Exciton Model (CEM) code CEM2k+GEM2 and the Los Alamos version of the Quark-Gluon String Model code LAQGSM are extended to describe photonuclear reactions. First, we incorporate into CEM2k+GEM2 new evaluations of elementary cross sections based on the latest experimental data and also make several improvements in the description of the de-excitation of nuclei remaining after the cascade stage of reactions induced by arbitrary projectiles. Next, for photonuclear reactions we include in CEM2k+GEM2 a normalization to evaluated experimental absorption cross sections based on the recent systematics by Kossov. Then, we extend our high-energy code LAQGSM by adding the photonuclear mode which was ignored in all its previous versions, and add to it the photonuclear part from our improved CEM2k+GEM2. In this work we present a short description of the photonuclear mode as incorporated into our codes, show several illustrative results, and point out some unresolved problems.

An improved version of the Cascade-Exciton Model (CEM) of nuclear reactions implemented in the codes CEM2k and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) has been developed recently at LANL to describe reactions induced by particles and nuclei at energies up to hundreds of GeV/nucleon for a number of applications. We present several improvements to the intranuclear cascade models used in CEM2k and LAQGSM developed recently to better describe the physics of nuclear reactions. First, we incorporate the photonuclear mode from CEM2k into LAQGSM to allow it to describe photonuclear reactions, not previously modeled there. Then, we develop new approximations to describe more accurately experimental elementary energy and angular distributions of secondary particles from hadron-hadron and photon-hadron interactions using available data and approximations published by other authors. Finally, to consider reactions involving very highly excited nuclei (E* > 2-3 MeV/A), we have incorporated into CEM2k and LAQGSM the Statistical Multifragmentation Model (SMM), as a possible reaction mechanism occurring after the preequilibrium stage. A number of other refinements to our codes developed recently are also listed.

Double-differential cross sections of neutron production at angles from 0 to 110 degrees from many reactions induced by light and medium nuclei on targets from 12-C to 208-Pb, at several incident energies from 95 to 600 MeV/A have been measured recently at the Institute of Physical and Chemical Research (RIKEN) Ring Cyclotron in Japan and at the Heavy-Ion Medical Accelerator of the National Institute of Radiological Science in Chiba, Japan using the time-of-flight technique. We have analyzed all these new measurements using the Quantum Molecular Dynamics (QMD) model, the Oak Ridge intranuclear cascade model HIC, the ISABEL intranuclear cascade model from LAHET, and the Los Alamos version of the Quark-Gluon String Model code LAQGSM03. On the whole, all four models used here describe reasonably well most of the measured neutron spectra, although different models agree differently with data from specific reactions and some serious discrepances are observed for some reactions. We present here some illustrative results from our study, discuss possible reasons for some of the observed discrepancies and try to outline ways to further improve the tested codes in order to address these problems.

We have analyzed the recent JINR measurements on nuclide production cross sections from interaction of 660 MeV proton beams with radioactive targets of enriched 129I (85% 129I and 15% 127I), 237Np, and 241Am with eleven different models, realized in eight transport codes and event-generators: LAHET (Bertini, ISABEL, INCL+ABLA, and INCL+RAL options), CASCADE, CEM95, CEM2k, LAQGSM+GEM2, CEM2k+GEM2, LAQGSM+GEMINI, and CEM2k+GEMINI. We found out that all these models have problems in a correct description of many of these cross sections, though some of these models describe very well most of the recent measurements done at GSI using inverse kinematics, as well as many other reactions. None of the tested here models is able to reproduce well all the JINR data and all of them should be further improved. Development of a better universal evaporation/fission model should be of a highest priority. We conclude that it is impossible to make a correct choice between fission and fragmentation reaction mechanisms analyzing only measurements on product cross sections; addressing this question would require analysis of two- or multi-particle correlation measurements.

The CEM2k and LAQGSM codes have been recently developed at Los Alamos National Laboratory to simulate nuclear reactions induced by particles and nuclei for a number of applications. We have benchmarked our codes against most available measured data at projectile energies from 10 MeV/A to 800 GeV/A and have compared our results with predictions of other current models used by the nuclear community. Here, we present a brief description of our codes and show illustrative results obtained with CEM2k and LAQGSM for A+p and A+A spallation, fission, and fragmentation reactions measured recently at GSI compared with predictions by other models. Further necessary work is outlined.

Two NaI (85% 129-I and 15% 127-I) targets were exposed to a beam of 660-MeV protons. Cross sections for formation of 76 residual nuclei were obtained by the induced activity method. The results are compared with other experimental data on 127-I and theoretical calculations by eleven models contained in the codes LAHET3 (using the Bertini+Dresner, ISABEL+Dresner, INCL+Dresner, and INCL+ABLA options), CASCADE, CEM95, CEM2K, LAQGSM+GEM2, CEM2k+GEM2, LAQGSM+GEMINI, and CEM2k+GEMINI. Most of the models describe spallation products with masses close to the target reasonably well while the reliability of the codes differs greatly in the deep spallation and fission/fragmentation regions. The difficulties in describing products with A=40-80 by all of the codes tested here except for CEM2k+GEMINI and LAQGSM+GEMINI is related to the neglect of fission (and fragmentation) processes for targets as light as 129-I.

The ITEP proton synchrotron U10 was used to irradiate isotopically-enriched 208Pb and natU thin targets with 1.0 GeV protons and 197Au thin targets with 0.8 GeV protons. More than 400 cross sections of the nuclides produced were measured using the direct spectrometry method with a high-resolution Ge detector. The measured cross sections are compared with similar data obtained at GSI for kinematically inverse reactions of 1 GeV/nucleon 208Pb, 1 GeV/nucleon 238U, and 0.8 GeV/nucleon 197Au interacting with a hydrogen target and with the ZSR data on natPb and 197Au irradiated with 1 and 0.8 GeV protons, respectively. Our results are on average ~ 10 to 20% higher than the GSI data. The measured data are analyzed with the LANL codes CEM2k+GEM2 and LAQGSM+GEM2 and with the INCL intranuclear cascade code from Liege merged with the GSI evaporation/fission code ABLA.

The LAQGSM code has been recently developed at Los Alamos National Laboratory to simulate nuclear reactions for proton radiography applications. We have benchmarked our code against most available measured data both for proton-nucleus and nucleus-nucleus interactions at incident energies from 10 MeV to 800 GeV and have compared our results with predictions of other current models used by the nuclear community. Here, we present a brief description of our code and show illustrative results obtained with LAQGSM for neutron spectra measured recently by Nakamura's groups for reactions induced by light and medium nuclei on targets from C to Pb at several incident energies from 95 to 600 MeV/nucleon and with the recent GSI measurements of spallation, fission, and fragmentation yields from A+p and A+A reactions at incident energies near and below 1 GeV/nucleon. Further necessary work is outlined.

We have benchmarked the QGSM code and event generators of the MARS and LAHET3 codes as potential candidates for high-energy programs to be used in simulations for the Proton Radiography (PRad) Project. We have compiled from the literature experimental data on spectra of particles emitted from proton-induced reactions at incident energies from 30 GeV to 70 GeV on different nuclei and have performed calculations for all reactions for which we found data with these three codes without any modifications and using only default parameters and standard inputs. Our results (514 plots) show that all three codes describe reasonably most of the studied reactions, though all of them should be further improved before becoming reliable tools for PRad. We present here our conclusions concerning the relative roles of different reaction mechanisms in the production of specific secondary particles. We comment on the strengths and weaknesses of QGSM, MARS, and LAHET3 and suggest further improvements to these codes and to other models.

The CEM2k and LAQGSM codes have been recently developed at Los Alamos National Laboratory to simulate nuclear reactions for a number of applications. We have benchmarked our codes against most available measured data at incident particle energies from 10 MeV to 800 GeV and have compared our results with predictions of other current models used by the nuclear community. Here, we present a brief description of our codes and show illustrative results to show that CEM2k and LAQGSM can be used as reliable event generators for space-radiation-shielding, cosmic-ray-propagation, and other astrophysical applications. Finally, we show the use of our calculated cross sections together with experimental data from our LANL T-16 compilation to produce evaluated files which we use in the GALPROP model of galactic particle propagation to better constrain the size of the CR halo.