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 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.

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.

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.

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 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.