A brief description of the IntraNuclear cascade, preequilibrium, evaporation, fission, coalescence, and Fermi breakup models used by the last versions of our CEM and LAQGSM event generators is presented, with a focus on our latest development of all 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 overviewed as well. Examples of benchmarking our models against a large variety of experimental data on particle-particle, particle nucleus, and nucleus-nucleus reactions, involving some very recent measurements of interest to accelerator radiation induced activation, are discussed.

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.

Cross-sections for radioactive nuclide production in 56Fe(p,x) reactions at 300, 500, 750, 1000, 1500, and 2600 MeV were measured using the ITEP U-10 proton accelerator. In total, 221 independent and cumulative yields of products of half-lives from 6.6 min to 312 days have been obtained via the direct-spectrometry method. The measured data have been compared with the experimental data obtained elsewhere by the direct and inverse kinematics methods and with calculations by 15 codes, namely: MCNPX (INCL, CEM2k, BERTINI, ISABEL), LAHET (BERTINI, ISABEL), CEM03 (.01, .G1, .S1), LAQGSM03 (.01, .G1, >.S1), CASCADE-2004, LAHETO, and BRIEFF. Most of our data are in a good agreement with the inverse kinematics results and disprove the results of some earlier activation measurements that were quite different from the inverse kinematics measurements. The most significant calculation-to-experiment differences are observed in the yields of the A<30 light nuclei, indicating that further improvements in nuclear reaction models are needed, and pointing out as well to a necessity of more complete measurements of such reactions.

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.

The work is aimed at experimental determination of the independent and cumulative yields of radioactive residual nuclei produced in intermediate-energy proton-irradiated thin targets made of highly isotopic enriched and natural lead (206,207,208,natPb) and 209Bi. 5972 radioactive product nuclide yields have been measured in 55 thin targets induced by 0.04, 0.07, 0.10, 0.15, 0.25, 0.6, 0.8, 1.2, 1.4, 1.6, and 2.6 GeV protons extracted from the ITEP U-10 proton synchrotron. The measured data have been compared with data obtained at other laboratories as well as with theoretical simulations by seven codes. We found that the predictive power of the tested codes is different but is satisfactory for most of the nuclides in the spallation region, though none of the codes agree well with the data in the whole mass region of product nuclides and all should be improved further.

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.

The nuclei fraction in cosmic rays (CR) far exceeds the fraction of other CR species, such as antiprotons, electrons, and positrons. Thus the majority of information obtained from CR studies is based on interpretation of isotopic abundances using CR propagation models where the nuclear data and isotopic production cross sections in p- and alpha-induced reactions are the key elements. This paper presents an introduction to the astrophysics of CR and diffuse gamma rays and discusses some of the puzzles that have emerged recently due to more precise data and improved propagation models. Merging with cosmology and particle physics, astrophysics of CR has become a very dynamic field with a large potential of breakthrough and discoveries in the near future. Exploiting the data collected by the CR experiments to the fullest requires accurate nuclear cross sections.

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.

Recent measurements of the CR antiproton flux have been shown to pose a problem for conventional propagation models. In particular, models consistent with secondary/primary nuclei ratio in CR produce too few antiprotons, while matching the ratio and the antiproton flux requires ad hoc assumptions. This may indicate an additional local CR component or new phenomena in CR propagation in the Galaxy. We discuss several possibilities which may cause this problem.

Accurate evaluation of production cross sections of light elements is important for models of CR propagation, galactic chemical evolution, and cosmological studies. However, experimental spallation cross section data are scarce and often unavailable to CR community while semi-empirical systematics are frequently wrong by a significant factor. We use all available data from LANL nuclear database together with modern nuclear codes to produce evaluated production cross sections of isotopes of Li, Be, B suitable for astrophysical applications.

The origin and evolution of isotopes of the lightest elements H2, He3, Li, Be, B in the universe is a key problem in such fields as astrophysics of CR, Galactic evolution, non-thermal nucleosynthesis, and cosmological studies. One of the major sources of these species is spallation by CR nuclei in the interstellar medium. On the other hand, it is the B/C ratio in CR and Be10 abundance which are used to fix the propagation parameters and thus the spallation rate. We study the production and Galactic propagation of isotopes of elements Z<6 using the numerical propagation code GALPROP and updated production cross sections.

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.

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux AND B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z<=28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three "radioactive clock" isotopes in CR, 10Be, 26Al, 36Cl, are all consistent and indicate a halo size z_h~4 kpc based on the most accurate data taken by the ACE spacecraft.

In a recent paper (Moskalenko et al., 2002), it has been shown that the flux of secondary cosmic ray (CR) antiprotons appears to be contradictory to measurements of secondary to primary nuclei ratios in cosmic rays when calculated in the same Galactic propagation model. The contradiction appears as a value of the diffusion coefficient necessary to match the secondary ratios antiproton/proton and B/C. In particlular, it was shown that the reacceleration models designed to match secondary to primary nuclei ratios produce too few antiprotons. It is, however, clear that some reacceleration is unavoidable in the turbulent interstellar medium. Here we discuss an idea of how to improve reacceleration model by allowing for the damping of interstellar turbulences on the small scale by cosmic rays, mostly protons. This would lead to increase in the mean free path lenghts at low energies, the well-known phenomena empirically discovered in the Leaky Box models, thus producing less secondary nuclei. Antiprotons will remain almost non-affected due to their high energy threshold of production cross section.

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., boron/carbon) produce too few antiprotons. Matching both the secondary to primary nuclei ratio and the antiproton flux requires artificial breaks in the diffusion coefficient and the primary injection spectrum suggesting the need for other approaches. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux AND B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z<=28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three "radioactive clock" isotopes in CR, 10Be, 26Al, 36Cl, are all consistent and indicate a halo size z_h~4 kpc based on the most accurate data taken by the ACE spacecraft.

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.

We use a new version of our numerical model for particle propagation in the Galaxy to study radioactive secondaries. For evaluation of the production cross sections we use the Los Alamos compilation of all available experimental cross sections together with calculations using the improved Cascade-Exciton Model code CEM2k. Using the radioactive secondary ratios 26Al/27Al, 36Cl/Cl, 54Mn/Mn, we show how the improved cross-section calculations together with the new propagation code allow us to better constrain the size of the CR halo.

A brief survey of nuclide abundances in the solar-system and in cosmic rays and of the believed mechanisms of their synthesis is given. The role of spallation processes in nucleosynthesis is discussed. A short review of recent measurements, compilations, calculations, and evaluations of spallation cross sections relevant to nuclear astrophysics is given as well. It is shown that in some past astrophysical simulations, old experimental nuclear data and theoretical cross sections that are in poor agreement with recent measurements and calculations were used. New astrophysical simulations using recently measured and reliably calculated nuclear cross sections, further researches in obtaining better cross sections, and production of evaluated spallation cross sections libraries for astrophysics are suggested.