Leslie M. Kerby

MCNP6 fragmentation of light nuclei at intermediate energies

Fragmentation reactions induced on light target nuclei by protons and light nuclei of energies around 1 GeV/nucleon and below are studied with the latest Los Alamos Monte Carlo transport code MCNP6 and with its cascade-exciton model (CEM) and Los Alamos version of the quark-gluon string model (LAQGSM) event generators, version 03.03, used as stand-alone codes. Such reactions are involved in different applications, like cosmic-ray-induced single event upsets (SEU’s), radiation protection, and cancer therapy with proton and ion beams, among others; therefore, it is important that MCNP6 simulates them as well as possible. CEM and LAQGSM assume that intermediate-energy fragmentation reactions on light nuclei occur generally in two stages. The first stage is the intranuclear cascade (INC), followed by th e second, Fermi breakup disintegration of light excited residual nuclei produced after INC. Both CEM and LAQGSM account also for coalescence of light fragments (complex particles) up to 4 He from energetic nucleons emitted during INC. We investigate the validity and performance of MCNP6, CEM, and LAQGSM in simulating fragmentation reactions at intermediate energies and discuss possible ways of further improving these codes.

Production of energetic light fragments in extensions of the CEM and LAQGSM event generators of the Monte Carlo transport code MCNP6

We extend the cascade-exciton model (CEM), and the Los Alamos version of the quark-gluon string model (LAQGSM), event generators of the Monte-Carlo N-particle transport code version 6 (MCNP6), to describe production of energetic light fragments (LF) heavier than 4He from various nuclear reactions induced by particles and nuclei at energies up to about 1 TeV/nucleon. In these models, energetic LF can be produced via Fermi break-up, preequilibrium emission, and coalescence of cascade particles. Initially, we study several variations of the Fermi break-up model and choose the best option for these models. Then, we extend the modified exciton model (MEM) used by these codes to account for a possibility of multiple emission of up to 66 types of particles and LF (up to 28Mg) at the preequilibrium stage of reactions. Then, we expand the coalescence model to allow coalescence of LF from nucleons emitted at the intranuclear cascade stage of reactions and from lighter clusters, up to fragments with mass numbers A < 8, in the case of CEM, and A < 13, in the case of LAQGSM. Next, we modify MCNP6 to allow calculating and outputting spectra of LF and heavier products with arbitrary mass and charge numbers. The improved version of CEM is implemented into MCNP6. Finally, we test the improved versions of CEM, LAQGSM, and MCNP6 on a variety of measured nuclear reactions. The modified codes give an improved description of energetic LF from particle- and nucleus-induced reactions; showing a good agreement with a variety of available experimental data. They have an improved predictive power compared to the previous versions and can be used as reliable tools in simulating applications involving such types of reactions.

An energy-dependent numerical model for the condensation probability, γj

Abstract The “condensation” probability, γ j , is an important variable in the preequilibrium stage of nuclear spallation reactions. It represents the probability that p j excited nucleons (excitons) will “condense” to form complex particle type j in the excited residual nucleus. It has a significant impact on the emission width, or probability of emitting fragment type j from the residual nucleus. Previous formulations for γ j were energy-independent and valid for fragments up to 4 He only. This paper explores the formulation of a new model for γ j , one which is energy-dependent and valid for up to 28 Mg, and which provides improved fits compared to experimental fragment spectra.

MCNP6 simulation of light and medium nuclei fragmentation at intermediate energies

MCNP6, the latest and most advanced LANL Monte Carlo transport code, representing a merger of MCNP5 and MCNPX, is actually much more than the sum of those two computer codes; MCNP6 is available to the public via RSICC at Oak Ridge, TN, USA. In the present work, MCNP6 was validated and verified (V&V) against different experimental data on intermediate-energy fragmentation reactions, and results by several other codes, using mainly the latest modifications of the Cascade-Exciton Model (CEM) and of the Los Alamos version of the Quark-Gluon String Model (LAQGSM) event generators CEM03.03 and LAQGSM03.03. It was found that MCNP6 using CEM03.03 and LAQGSM03.03 describes well fragmentation reactions induced on light and medium target nuclei by protons and light nuclei of energies around 1 GeV/nucleon and below, and can serve as a reliable simulation tool for different applications, like cosmic-ray-induced single event upsets (SEU’s), radiation protection, and cancer therapy with proton and ion beams, to name just a few. Future improvements of the predicting capabilities of MCNP6 for such reactions are possible, and are discussed in this work.

Production of Energetic Light Fragments in Spallation Reactions

Different reaction mechanisms contribute to the production of light fragments (LF) from nuclear reactions. Available models cannot accurately predict emission of LF from arbitrary reactions. However, the emission of LF is important for many applications, such as cosmic-ray-induced single event upsets, radiation protection, and cancer therapy with proton and heavy-ion beams, to name just a few. The cascade-exciton model (CEM) and the Los Alamos version of the quark-gluon string model (LAQGSM), as implemented in the CEM03.03 and LAQGSM03.03 event generators used in the Los Alamos Monte Carlo transport code MCNP6, describe quite well the spectra of fragments with sizes up to 4 He across a broad range of target masses and incident energies. However, they do not predict high-energy tails for LF heavier than 4 He. The standard versions of CEM and LAQGSM do not account for preequilibrium emission of LF larger than 4 He. The aim of our work is to extend the preequilibrium model to include such processes. We do this by including the emission of fragments heavier than 4 He at the preequilibrium stage, and using an improved version of the Fermi Break-up model, providing improved agreement with various experimental data.

Production of energetic light fragments in extensions of the CEM and LAQGSM event generators of the Monte Carlo transport code MCNP6 [Production of energetic light fragments in CEM, LAQGSM, and MCNP6]

We extend the cascade-exciton model (CEM), and the Los Alamos version of the quark-gluon string model (LAQGSM), event generators of the Monte-Carlo N-particle transport code version 6 (MCNP6), to describe production of energetic light fragments (LF) heavier than 4He from various nuclear reactions induced by particles and nuclei at energies up to about 1 TeV/nucleon. In these models, energetic LF can be produced via Fermi break-up, preequilibrium emission, and coalescence of cascade particles. Initially, we study several variations of the Fermi break-up model and choose the best option for these models. Then, we extend the modified exciton model (MEM) used by these codes to account for a possibility of multiple emission of up to 66 types of particles and LF (up to 28Mg) at the preequilibrium stage of reactions. Then, we expand the coalescence model to allow coalescence of LF from nucleons emitted at the intranuclear cascade stage of reactions and from lighter clusters, up to fragments with mass numbers A < 8, in the case of CEM, and A < 13, in the case of LAQGSM. Next, we modify MCNP6 to allow calculating and outputting spectra of LF and heavier products with arbitrary mass and charge numbers. The improved version of CEM is implemented into MCNP6. Finally, we test the improved versions of CEM, LAQGSM, and MCNP6 on a variety of measured nuclear reactions. The modified codes give an improved description of energetic LF from particle- and nucleus-induced reactions; showing a good agreement with a variety of available experimental data. They have an improved predictive power compared to the previous versions and can be used as reliable tools in simulating applications involving such types of reactions.

Production of Energetic Light Fragments in CEM, LAQGSM, and MCNP6

We extend the cascade-exciton model (CEM), and the Los Alamos version of the quark-gluon string model (LAQGSM), event generators of the Monte-Carlo N-particle transport code version 6 (MCNP6), to describe production of energetic light fragments (LF) heavier than 4He from various nuclear reactions induced by particles and nuclei at energies up to about 1 TeV/nucleon. In these models, energetic LF can be produced via Fermi break-up, preequilibrium emission, and coalescence of cascade particles. Initially, we study several variations of the Fermi break-up model and choose the best option for these models. Then, we extend the modified exciton model (MEM) used by these codes to account for a possibility of multiple emission of up to 66 types of particles and LF (up to 28Mg) at the preequilibrium stage of reactions. Then, we expand the coalescence model to allow coalescence of LF from nucleons emitted at the intranuclear cascade stage of reactions and from lighter clusters, up to fragments with mass numbers A < 8, in the case of CEM, and A < 13, in the case of LAQGSM. Next, we modify MCNP6 to allow calculating and outputting spectra of LF and heavier products with arbitrary mass and charge numbers. The improved version of CEM is implemented into MCNP6. Finally, we test the improved versions of CEM, LAQGSM, and MCNP6 on a variety of measured nuclear reactions. The modified codes give an improved description of energetic LF from particle- and nucleus-induced reactions; showing a good agreement with a variety of available experimental data. They have an improved predictive power compared to the previous versions and can be used as reliable tools in simulating applications involving such types of reactions.

Comparison of Expanded Preequilibrium CEM Model with CEM03.03 and Experimental Data, FY2013

Emission of light fragments (LF) from nuclear reactions is an open question. Different reaction mechanisms contribute to their production; the relative roles of each, and how they change with incident energy, mass number of the target, and the type and emission energy of the fragments is not completely understood. None of the available models are able to accurately predict emission of LF from arbitrary reactions. However, the ability to describe production of LF (especially at energies

Total reaction cross sections in CEM and MCNP6 at intermediate energies

\gtrsim 30

Preequilibrium Emission of Light Fragments in Spallation Reactions

MeV) from many reactions is important for different applications, such as cosmic-ray-induced Single Event Upsets (SEUs), radiation protection, and cancer therapy with proton and heavy-ion beams, to name just a few. The Cascade-Exciton Model (CEM) version 03.03 and the Los Alamos version of the Quark-Gluon String Model (LAQGSM) version 03.03 event generators in Monte Carlo N-Particle Transport Code version 6 (MCNP6) describe quite well the spectra of fragments with sizes up to