Reaction initial state simulation.

Allowed projectiles and bombarding energy range for interaction with nucleon and nuclear targets

The Geant4 parton string models are capable to predict final states (produced hadrons which belong to the scalar and vector meson nonets and the baryon (antibaryon) octet and decuplet) of reactions on nucleon and nuclear targets with nucleon, pion and kaon projectiles. The allowed bombarding energy s>5  GeV is recommended. Two approaches, based on diffractive excitation or soft scattering with diffractive admixture according to cross-section, are considered. Hadron-nucleus collisions in the both approaches (diffractive and parton exchange) are considered as a set of the independent hadron-nucleon collisions. However, the string excitation procedures in these approaches are rather different.

MC initialization procedure for nucleus

The initialization of each nucleus, consisting from A nucleons and Z protons with coordinates ri and momenta pi, where i=1,2,...,A is performed. We use the standard initialization Monte Carlo procedure, which is realized in the most of the high energy nuclear interaction models:

  • Nucleon radii ri are selected randomly in the rest of nucleus according to proton or neutron density ρ(ri). For heavy nuclei with A>16 [GLMP91] nucleon density is

    ρ(ri)=ρ01+exp[(riR)/a]

    where

    ρ034πR3(1+a2π2R2)1.

    Here R=r0A1/3  fm and r0=1.16(11.16A2/3)  fm and a0.545 fm. For light nuclei with A<17 nucleon density is given by a harmonic oscillator shell model [B61], e. g.

    ρ(ri)=(πR2)3/2exp(r2i/R2),

    where R2=2/3r2=0.8133A2/3  fm2. To take into account nucleon repulsive core it is assumed that internucleon distance d>0.8  fm;

  • The initial momenta of the nucleons are randomly choosen between 0 and pmaxF, where the maximal momenta of nucleons (in the local Thomas-Fermi approximation [DA74]) depends from the proton or neutron density ρ according to

    pmaxF=c(3π2ρ)1/3

    with c=0.197327 GeV fm;

  • To obtain coordinate and momentum components, it is assumed that nucleons are distributed isotropicaly in configuration and momentum spaces;

  • Then perform shifts of nucleon coordinates rj=rj1/Airi and momenta pj=pj1/Aipi of nucleon momenta. The nucleus must be centered in configuration space around 0, i. e. iri=0 and the nucleus must be at rest, i. e. ipi=0;

  • We compute energy per nucleon e=E/A=mN+B(A,Z)/A, where mN is nucleon mass and the nucleus binding energy B(A,Z) is given by the Bethe-Weizsäcker formula [BA69]:

    B(A,Z)=0.01587A+0.01834A2/3+0.09286(ZA2)2+0.00071Z2/A1/3,

    and find the effective mass of each nucleon meffi=(E/A)2p2i.

Random choice of the impact parameter

The impact parameter 0bRt is randomly selected according to the probability:

P(b)db=bdb,

where Rt is the target radius, respectively. In the case of nuclear projectile or target the nuclear radius is determined from condition:

ρ(R)ρ(0)=0.01.

Bibliography

B61

Elton L. R. B. Nuclear Sizes. Oxford University Press, Oxford, 1961.

BA69

Mottelson B. R. Bohr A. Nuclear Structure. W. A. Benjamin, New York, vol. 1 edition, 1969.

DA74

Feshbach H. DeShalit A. Theoretical Nuclear Physics. Wyley, vol. 1: nuclear structure edition, 1974.

GLMP91

M E Grypeos, G A Lalazissis, S E Massen, and C P Panos. The 'cosh' or symmetrized woods-saxon nuclear potential. Journal of Physics G: Nuclear and Particle Physics, 17(7):1093, 1991. URL: http://stacks.iop.org/0954-3899/17/i=7/a=008.