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.

The Los Alamos version of the Quark-Gluon String Model implemented in the code LAQGSM03 was extended to describe photonuclear reactions at energies up to tens of GeV. We have incorporated 56 channels to consider gamma+p elementary interactions during the cascade stage of reactions and 56 channels for gamma+n interactions, in addition to absorption of photons on two-nucleon pairs, into LAQGSM03 and have tested the extended code referred to as LAQGSM03.01 on a number of measured high-energy photonuclear reactions. Here, we present the photonuclear reaction model of LAQGSM03.01 and show examples of its results compared with available experimental data.

This paper presents part of an internal LANL Progress Report on completion of the "S" and "G" versions of the improved Cascade-Exciton Model (CEM03.01) and the Los Alamos Quark-Gluon String Model (LAQGSM.03.01) codes. The "S" versions consider fragmentation of compound nuclei produced after the preequilibrium stage of reactions for excitation energies above 2A MeV using the Statistical Multifragmentation Model (SMM) by Botvina et al. ("S" stands for SMM), while the "G" versions describe evaporation/fission stages of reactions using the fission-like binary-decay model GEMINI of Charity et al. ("G" stands for GEMINI) instead of using the the Generalized Evaporation Model GEM2 of Furihata incorporated into the standard versions of these codes. We present here an analysis of the recent 660 MeV p + 129I and 3.65 GeV p + 112Sn JINR measurements, of the new COSY data on 1.2 GeV p + (13 nuclei from Al to Th), of the 300 MeV and 1 GeV p + 56Fe data measured at GSI in inverse kinematics, and of the new GSI data on 1 GeV/nucleon 124Xe and 136Xe + Pb. To better understand the mechanisms of fragment production, we discuss several calculated but not-yet-measured kinematic characteristics of products of these reactions, which are predicted to be quite different by SMM, GEMINI, and GEM2. We find these kinematic quantities to be potentially useful in differentiating these reaction mechanisms if they can be measured in future experiments.

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.

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.

We measure cross sections for residual nuclide formation in the mass range 6 < A < 97 caused by bombardment with protons and deuterons of 3.65 GeV/nucleon energy of enriched tin isotopes (112-Sn, 118-Sn, 120-Sn, 124-Sn). The experimental data are compared with calculations by the codes FLUKA, LAHET, CEM03, and LAQGSM03. Scaling behavior is observed for the whole mass region of residual nuclei, showing a possible multifragmentation mechanism for the formation of light products (6 < A < 31). Our analysis of the isoscaling dependence also shows a possible contribution of multifragmentation to the production of heavier nuclides, in the mass region 39 < A < 81.

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.

In order to meet the challenges of new accelerator and space projects and further improve modelling of radiation effects in microscopic objects, heavy-ion interaction and transport physics have been recently incorporated into the MARS15 Monte Carlo code. A brief description of new modules is given in comparison with experimental data.

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.

The improved Cascade-Exciton Model code CEM2k and the Los Alamos version of the Quark-Gluon String Model code LAQGSM, previously merged with the Generalized Evaporation Model code of Furihata (GEM2) were further modified to provide reliable proton-induced fission cross sections for applications. By adjusting two parameters in GEM2 for each measured reaction, we were able to describe very well with CEM2k+GEM2 and LAQGSM+GEM2 all available experimental fission cross sections induced by protons with energies from 20 MeV to 10 GeV both for subactinide and actinide targets. We also successfully tested our approach on several reactions induced by neutrons, pions, and photons.

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.