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.

An improved version of the Cascade-Exciton Model (CEM) of nuclear reactions realized in the code CEM2k and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) have been developed recently at LANL to describe reactions induced by particles and nuclei for a number of applications. Our CEM2k and LAQGSM merged with the GEM2 evaporation/fission code by Furihata have predictive powers comparable to other modern codes and describe many reactions better than other codes; therefore both our codes can be used as reliable event generators in transport codes for applications. During the last year, we have made a significant improvements to the intranuclear cascade parts of CEM2k and LAQGSM, and have extended LAQGSM to describe photonuclear reactions at energies to 10 GeV and higher. We have produced in this way improved versions of our codes, CEM03.01 and LAQGSM03.01. We present a brief description of our codes and show illustrative results obtained with CEM03.01 and LAQGSM03.01 for different reactions compared with predictions by other models, as well as examples of using our codes as modeling tools for nuclear applications.

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 MARS code is under continuous development and has recently undergone substantial improvements that further increase its reliability and predictive power in numerous shielding, accelerator, detector and space applications. The major developments and new features of the MARS15 (2004) version described in this paper concern an extended list of elementary particles and arbitrary heavy ions and their interaction cross-sections, inclusive and exclusive nuclear event generators, module for modelling particle electromagnetic interactions, enhanced geometry and histograming options, improved MAD-MARS Beam Line Builder, enhanced Graphical-User Interface, and an MPI-based parallelization of the code.

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 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 present the current status of the improved Cascade-Exciton Model (CEM) code CEM2k and of the Los Alamos version of the Quark-Gluon String Model code LAQGSM. To describe fission and light-fragment (heavier than He4) production, both CEM2k and LAQGSM have been merged with the GEM2 code of Furihata. We present some results on proton- and deuteron-induced spallation, fission, and fragmentation reactions predicted by these extended versions of CEM2k and LAQGSM. We show that merging CEM2k and LAQGSM with GEM2 allows us to describe many fission and fragmentation reactions in addition to the spallation reactions which are already relatively well described. Nevertheless, the standard version of GEM2 does not provide a completely satisfactory description of complex particle spectra, heavy-fragment emission, and spallation yields, and is not yet a reliable tool for applications. We conclude that we may choose to use a model similar to the GEM2 approach in our codes, but it must be significantly extended and further improved. We observe that it is not sufficient to analyze only A and Z distributions of the product yields when evaluating this type of model, as is often done in the literature; instead it is important to study all the separate isotopic yields as well as the spectra of light particles and fragments.

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.

A brief description of our improvements and refinements that led from the CEM95 version of the Cascade-Exciton Model (CEM) code to CEM97 and to CEM2k is given. The increased accuracy and predictive power of the code CEM2k are shown by several examples. To describe fission and light-fragment (heavier than 4He) production, the CEM2k code has been merged with the GEM2 code of Furihata. We present some results on proton-induced fragmentation and fission reactions predicted by this extended version of CEM2k. We show that merging CEM2k with GEM2 allows us to describe many fission and fragmentation reactions in addition to the spallation reactions which are already relatively well described. We have initiated another approach to describe fission, complex particles and fragment emission by developing further our CEM2k code addressing specifically these problems. In this effort, we have developed our own universal approximation for inverse cross sections, new routines to calculate Coulomb barriers and widths of emitted particles and to simulate their kinetic energy using arbitrary approximations for the inverse cross sections. To describe fission-fragment production, we have incorporated into CEM2k a thermodynamical model of fission by Stepanov. This extended version of CEM2k allows us to describe much better complex particle emission and many fission fragments, but it is still incomplete and needs further work.