Muon Photonuclear Interaction

The inelastic interaction of muons with nuclei is important at high muon energies ($E \geq 10$  GeV), and at relatively high energy transfers $\nu$ ( $\nu / E \geq 10^{-2}$). It is especially important for light materials and for the study of detector response to high energy muons, muon propagation and muon-induced hadronic background. The average energy loss for this process increases almost lineary with energy, and at TeV muon energies constitutes about 10% of the energy loss rate.

The main contribution to the cross section $\sigma ( E, \nu$) and energy loss comes from the low $Q^{2}$-region ( $Q^{2} \ll 1~{\rm GeV}^{2}$). In this domain, many simplifications can be made in the theoretical consideration of the process in order to obtain convenient and simple formulae for the cross section. Most widely used are the expressions given by Borog and Petrukhin [1], and Bezrukov and Bugaev [2]. Results from these authors agree within 10% for the differential cross section and within about 5% for the average energy loss, provided the same photonuclear cross section, $\sigma_{\gamma N}$, is used in the calculations.

Differential Cross Section

The Borog and Petrukhin formula for the cross section is based on:

• Hand's formalism [3] for inelastic muon scattering,
• a semi-phenomenological inelastic form factor, which is a Vector Dominance Model with parameters estimated from experimental data, and
• nuclear shadowing effects with a reasonable theoretical parameterization [4].

For $E \geq 10$ GeV, the Borog and Petrukhin cross section (${\rm cm}^{2}$/g GeV), differential in transferred energy, is

$$\sigma(E,\nu ) = \Psi(\nu ) \Phi( E,v ) ,$$ (9.20)

$$\Psi(\nu) = \frac{\alpha}{\pi} \frac{A_{\rm eff} N_{AV}}{A} \sigma_{\gamma N}(\nu) \frac{1}{\nu} ,$$ (9.21)

$$\Phi(E,v) = v-1 + \left[1-v+\frac{v^{2}}{2} \left(1+\frac... ...\Lambda}\left(1+\frac{\Lambda}{2M}+\frac{Ev}{\Lambda}\right)},$$ (9.22)

where $\nu$ is the energy lost by the muon, $v = \nu /$E, and $\mu$ and $M$ are the muon and nucleon (proton) masses, respectively. $\Lambda$ is a Vector Dominance Model parameter in the inelastic form factor which is estimated to be $\Lambda^{2}=0.4 \;{\rm GeV}^{2}$.

For $A_{\rm eff}$, which includes the effect of nuclear shadowing, the parameterization [4]

$$A_{\rm eff} = 0.22 A + 0.78 A^{0.89}$$ (9.23)

is chosen.

A reasonable choice for the photonuclear cross section, $\sigma_{\gamma N}$, is the parameterization obtained by Caldwell et al. [5] based on the experimental data on photoproduction by real photons:

$$\sigma_{\gamma N} = ( 49.2 + 11.1 \ln K + 151.8/ \sqrt{K} ) \cdot 10^{-30} {\rm cm}^{2} \quad K~~\mbox{in GeV} .$$ (9.24)

The upper limit of the transferred energy is taken to be $\nu_{\rm max}$ = $E - M/2$. The choice of the lower limit $\nu_{\rm min}$ is less certain since the formula 9.20, 9.21, 9.22 is not valid in this domain. Fortunately, $\nu_{\rm min}$ influences the total cross section only logarithmically and has no practical effect on the average energy loss for high energy muons. Hence, a reasonable choice for $\nu_{\rm min}$ is 0.2 GeV.

In Eq. 9.21, $A_{\rm eff}$ and $\sigma_{\gamma N}$ appear as factors. A more rigorous theoretical approach may lead to some dependence of the shadowing effect on $\nu$ and $E$; therefore in the differential cross section and in the sampling procedure, this possibility is forseen and the atomic weight $A$ of the element is kept as an explicit parameter.

The total cross section is obtained by integration of Eq. 9.20 between $\nu_{\min}$ and $\nu_{\max}$; to facilitate the computation, a $\ln (\nu )$-substitution is used.

Sampling

Sampling the Transferred Energy

The muon photonuclear interaction is always treated as a discrete process with its mean free path determined by the total cross section. The total cross section is obtained by the numerical integration of Eq. 9.20 within the limits $\nu_{\min}$ and $\nu_{\max}$. The process is considered for muon energies $1 \rm {GeV} \leq T \leq 1000 \rm {PeV}$, though it should be noted that above 100 TeV the extrapolation (Eq. 9.24) of $\sigma_{\gamma N}$ may be too crude.

The random transferred energy, $\nu_{p}$, is found from the numerical solution of the equation :

$$P = \int_{\nu_{p}}^{\nu_{\rm max}} \sigma(E,\nu ) d \nu \le... ...t_{\nu_{\rm min}}^{\nu_{\rm max}} \sigma(E,\nu) d\nu \right. .$$ (9.25)

Here $P$ is the random uniform probability, with $\nu_{\rm max}= E-M/2$ and $\nu_{\rm min}=0.2$ GeV.

For fast sampling, the solution of Eq. 9.25 is tabulated at initialization time. During simulation, the sampling method returns a value of $\nu_{p}$ corresponding to the probability $P$, by interpolating the table. The tabulation routine uses Eq. 9.20 for the differential cross section. The table contains values of

$$x_p = \ln (\nu_p / \nu_{\rm max})/\ln (\nu_{\rm max}/\nu_{\rm min}),$$ (9.26)

calculated at each point on a three-dimensional grid with constant spacings in $\ln (T)$, $\ln(A)$ and $\ln(P)$ . The sampling uses linear interpolations in $\ln (T)$ and $\ln(A)$, and a cubic interpolation in $\ln(P)$. Then the transferred energy is calculated from the inverse transformation of Eq. 9.26, $\nu_{p}=\nu_{\rm max}(\nu_{\rm max}/\nu_{\rm min})^{x_{p}}$. Tabulated parameters reproduce the theoretical dependence to better than 2% for $T > 1$ GeV and better than 1% for $T > 10$ GeV.

Sampling the Muon Scattering Angle

According to Refs. [1,6], in the region where the four-momentum transfer is not very large ( $Q^{2} \leq 3 {\rm GeV}^{2}$), the $t$ - dependence of the cross section may be described as:

$$\frac{d \sigma }{dt} \sim \frac{(1- t / t_{\rm max}) } {t (1+t/ \nu^{2})(1+t/m^{2}_{0})} [(1-y)(1-t_{\rm min}/t)+y^{2}/2] ,$$ (9.27)

where $t$ is the square of the four-momentum transfer, $Q^{2} = 2 (EE' - PP'\cos\theta - \mu^{2})$. Also, $t_{\rm min} = (\mu y)^{2} / (1-y)$, $y=\nu /E$ and $t_{\rm max} = 2M\nu$. $\nu = E-E'$ is the energy lost by the muon and $E$ is the total initial muon energy. $M$ is the nucleon (proton) mass and $m_{0}^{2} \equiv \Lambda^{2} \simeq 0.4\;{\rm GeV}^{2}$ is a phenomenological parameter determing the behavior of the inelastic form factor. Factors which depend weakly, or not at all, on $t$ are omitted.

To simulate random $t$ and hence the random muon deflection angle, it is convenient to represent Eq. 9.27 in the form :

$$\sigma( t ) \sim f(t) g(t),$$ (9.28)

where

$$f(t) = \frac{1}{t(1+ t/t_{1} )} ,$$ (9.29)

$$g(t) = \frac{1-t/t_{\max}}{1+t/t_{2}} \cdot \frac{(1-y)(1-t_{\min}/t)+y^{2}/2} {(1-y)+y^{2}/2},$$

and

$$t_{1} = {\min} (\nu^{2}, m_{0}^{2}) \quad t_{2} ={\max} (\nu^{2}, m_{0}^{2}) .$$ (9.30)

$t_{P}$ is found analytically from Eq. 9.29 :

$$t_{P} = \frac{t_{\max}t_{1}} {\displaystyle(t_{\max}+t_{1})\l... ...min}+t_{1})} {t_{\min}(t_{\max}+t_{1})}\right]^{P}-t_{\max}} ,$$

where $P$ is a random uniform number between 0 and 1, which is accepted with probability $g(t)$. The conditions of Eq. 9.30 make use of the symmetry between $\nu^{2}$ and $m_{0}^{2}$ in Eq. 9.27 and allow increased selection efficiency, which is typically $\geq 0.7$. The polar muon deflection angle $\theta$ can easily be found from ^9.1

$$\sin^{2}(\theta /2) = \frac {t_{P} - t_{\rm min}} {4\,(EE'- \mu^{2}) - 2\,t_{\rm min}}.$$

The hadronic vertex is generated by the hadronic processes taking into account the four-momentum transfer.

Status of this document

12.10.98 created by R.Kokoulin, A.Rybin.
18.05.00 edited by S.Kelner, R.Kokoulin, and A.Rybin.
07.12.02 re-worded by D.H. Wright
30.08.04 correction of eq. 8.24 (to 1/sqrt) from H. Araujo

Bibliography

1. V.V.Borog and A.A.Petrukhin, Proc. 14th Int.Conf. on Cosmic Rays, Munich,1975, vol.6, p.1949.
2. L.B.Bezrukov and E.V.Bugaev, Sov. J. Nucl. Phys., 33, 1981, p.635.
3. L.N.Hand. Phys. Rev., 129, 1834 (1963).
4. S.J.Brodsky, F.E.Close and J.F.Gunion, Phys. Rev. D6, 177 (1972).
5. D.O. Caldwell et al., Phys. Rev. Lett., 42, 553 (1979).
6. V.V.Borog, V.G.Kirillov-Ugryumov, A.A.Petrukhin, Sov. J. Nucl. Phys., 25, 1977, p.46.