\documentstyle[times,pramana,epsf,floats,amsmath,amssymb,amsthm,slashed]{ias}

\input{mydefs}

\begin{document}

\title{The stimulated Breit-Wheeler process as a source of background $e^{+}e^{-}$ pairs at the international
linear collider}

\author{A. Hartin}
\address{John Adams Institute, Oxford University Physics, Oxford UK, OX1 3RH}

\abstract{Passage of beamstrahlung photons through the bunch fields at the interaction point of the ILC
determines background pair production. The number of background pairs per bunch crossing due to the Breit-Wheeler, 
Bethe-Heitler and Landau-Lifshitz processes is well known. However the Breit-Wheeler process also 
takes place in and is modified by the bunch fields. A full QED calculation of this Stimulated Breit-Wheeler process reveals
cross section resonances due to the virtual particle reaching the mass shell. The one loop electron self energy 
in the bunch field is also calculated and included as a radiative correction. The bunch
field is considered to be a contant crossed electromagnetic field with associated bunch field photons.
Resonance is found to occur whenever the energy of contributed bunch field photons is equal to the
beamstrahlung photon energy. The stimulated Breit-Wheeler cross section exceeds the ordinary
Breit-Wheeler cross section by several orders of magnitude and a significantly different pair background may
result.}

\maketitle
\section{INTRODUCTION}

The design of the International linear Collider includes intense energetic electron and
positron bunch collisions in order to maximise centre of mass energy and
collider luminosity. Associated with these intense bunches are strong electromagnetic fields
which affect the physics processes at and near the interaction point. Of
particular concern is unwanted pair production which causes background
radiation in detectors and other collider components. The production of background pairs from
two real beamstrahlung photons is included in background calculations. However at present the effect of the
strong bunch fields on this process is ignored. This paper seeks to address that issue. 

In order to perform a full QED calculation of the Stimulated Breit-Wheeler cross section, 
a semi-classical approach is adopted in which the Dirac equation is solved exactly in a classical 
external plane wave field. Both the fermion field operators and the electron propagator are modified. The
Stimulated Breit-Wheeler process is a second order intense field QED process, and it is known that this class of
processes contains resonant cross-sections \cite{oleinik,bos}. In order to calculate the contribution from the
resonances radiative corrections have to be included.
The persistence of the bunch fields is very short being of order $10^{-14} \; s$. However the resonances 
are introduced in the virtual particle exchange which is of order $10^{-21} \; s$. There is more than enough
time for the resonant background pair process to occur and the calculation presented here is pertinent to
all intense charged bunch collisions.

\section{DIRAC EQUATION SOLUTIONS IN THE BUNCH FIELDS}

The electron and positron bunches at the ILC are ultra-relativistic
and the electromagnetic fields produced by ultra-relativistic beams approximate well 
to a constant, crossed plane wave propagating, say, along the $x$ axis with 4-momentum $k$ and described by 4-potential
$a_\mu(k.x)$. For the default Scheme 1 set of ILC bunch parameters including dimensions $\sigma_x \sigma_y
\sigma_z$, the number of charges $N$, relativistic parameters $\gamma ,\beta$) \cite{tor} and with $B_c$ being the Schwinger
critical field, the numerical values of the bunch field intensity $\nu^2$ and bunch field photon energy
$\omega$ can be established

\begin{equation}
\begin{array}{l}
  \dfrac{\omega}{m} = \dfrac{\hbar}{mc^2} \dfrac{e}{m} \dfrac{5}{6} B_c
  \dfrac{Nr^2_e \gamma}{\alpha \sigma_z (\sigma_x + \sigma_y)} = 0.0603 \\[10pt]
  \nu^2 = \dfrac{e^2}{\omega^2
  m^2 \beta} \dfrac{5}{6} B_c \dfrac{Nr^2_e \gamma}{\alpha \sigma_z
  (\sigma_x + \sigma_y)} = 1.000 
\end{array}
\end{equation}

Dirac equation solutions in plane electromagnetic waves are well known~\cite{volk}. The solutions are a product of
a bispinor $u_p$ and an exponential function of the bunch field 4-potential. The fourier transform of this
solution yields a set of $r$ Volkov field operators of the fermion 4-momentum $p$ each corresponding to a 
contribution of $r$ bunch field photons. The 4-potential of a contant crossed electromagnetic field introduces 
Airy functions into the Volkov solution as long as there is azimuthal symmetry (which is assumed here) 
about the axis of propagation direction of the bunch fields.

\begin{equation}
\label{eq2}
\begin{array}{c}
  \Psi^{V}_{p,r}(x) = E_p(x)\; u_p\; e^{-i(p+Q)x-irkx} \\[8pt]
  E_p (x) = [ 1 + \frac{e}{2 (kp)} \slashed{k} \slashed{a} (kx)] Q^{-1/3} Ai(z)\\[8pt]
  \text{where} \quad Q = \frac{e^2 a^2}{2 (kp)} k \quad \text{and} \quad z=-Q^{-1/3}(r+Q) 
\end{array}
\end{equation} 

\section{ELECTRON SELF ENERGY IN THE BUNCH FIELDS}

The class of second order external field QED processes are known to have resonant cross sections due to the
virtual particle reaching the mass shell. Therefore radiative corrections to the Stimulated
Breit-Wheeler cross section calculation considered in this paper are probably necessary. 
The only radiative correction considered in this instance is the
Electron Self Energy which also occurs in the presence of the bunch field (Fig.~\ref{f1}). 


\begin{figure}[htbp]
\epsfxsize=15mm
\centerline{\epsfbox{eself.eps}}
\caption{Electron Self Energy Feynman diagrams.}
\label{f1}
\end{figure}

The Electron Self Energy in a constant crossed field $\Sigma(p)$ has been calculated before, and contains the 
usual UV and IR divergences. In addition, an extra pole is introduced in the limit of vanishing bunch 
field intensity~\cite{Ritus72}. The UV divergence appears in the photon propagator which is unmodified by the
bunch field. The regularisation and renormalisation procedures are therefore the same as the non external field
case and the corrected electron propagator to be used for the Stimulated Breit-Wheeler scattering matrix is

\begin{equation}
\label{eq3}
  [\slashed{p} -m-\Sigma(p)+\Sigma(p\rightarrow m)]^{-1}
\end{equation}

The Electron Self Energy can be calculated after inclusion of the electron propagator in the external field.
This is the non-external field electron propagator sandwiched between Volkov $E_p$ functions. The integrations
over internal 4-momenta $k'$ and $q'$ are transformed using Heaviside step functions and a shift to 4
dimensional spherical coordinates. The imaginary and real parts are then separated using dispersion 
relations \cite{kallen}. Renormalisation introduces an infinitesimal shift to the lower bounds of the remaining
integrations and the, thus avoided, IR divergences are removed by further transformations in 
integration variables. The imaginary part of the Self Energy finally is rendered as

\begin{figure}[htbp]
\epsfxsize=7cm
\centerline{\epsfbox{eesaz_vary_rho_nu.eps}}
\caption{imaginary part of the Electron Self Energy in the bunch field.}
\label{fig2}
\end{figure}

\begin{equation}
\label{eq4}
\begin{array}{l}
  \Im \Delta \varepsilon^R_p (\rho) = \frac{e^2 m^2}{16 \pi \varepsilon_p}
  \dint^{\infty}_{\nths 0+\epsilon} dr \dint^{u_r}_{\nths 0+\epsilon} \frac{3du}{2(\rho + u^3 r^{3/2})^2} \\[10pt]
   \times \left[ - 4 \rho r u^2 Ai(z)^2 - 2 \nu^2 ( 1-u ) \left( 2\rho + 
  \frac{r^3 u^6}{\rho + u^3 r^{3/2}} \right)\right. \\[10pt] \left.
  \times\left( ruAi(z)^2 + \frac{1}{1-u} Ai' (z)^2 \right) \right]
   \; \text{where} \;  \rho = \frac{\nu^2}{(kp)}
\end{array}
\end{equation}

The final integrations are straightforward and the imaginary part of the Self Energy bears an almost linear
relationship with the parameter $\rho$ which is a function of the scalar product of virtual particle and bunch
field photon 4-momenta and the bunch field intensity $\nu^2$ (see Figure \ref{fig2}).

\section{STIMULATED BREIT-WHEELER CROSS-SECTION}

The cross section of the Stimulated Breit-Wheeler process is calculated from the matrix element whose form is
written down with the aid of the Feynman diagrams (Figure \ref{fig3} and Equation \ref{eq5}).

\begin{figure}[htbp]
\epsfxsize=8cm
\centerline{\epsfbox{pprod.eps}}
\caption{Stimulated Breit-Wheeler Feynman diagrams.}
\label{fig3}
\end{figure}


\begin{equation}
\label{eq5}
\begin{array}{l}
  S_{f i}^e = - \dfrac{e^2}{2} \dint^{\infty}_{\nths -\infty} d^4 x_1 d^4 x_2 d^4 q 
  \; \dfrac{1}{p^2 - m^2}  \bar{u}(p_-) \\[10pt]
  \quad\quad\quad\quad \times\left(\mathcal{A}e^{-ik_1 x_2 - ik_2 x_1} 
   +\mathcal{B}e^{-ik_2 x_2- ik_1 x_1} \right) v(p_+) \\[10pt]
  \text{with} \mathcal{A}= \bar{E}_{p_{-}} (x_2) \gamma_{\mu} E_p (x_2) (
  \slashed{p} + m) \bar{E}_{p} (x_1) \gamma_{\nu} E_{-p_{+}} (x_1)  \\[10pt]
   \text{and} \; \mathcal{B}= \bar{E}_{p_{-}} (x_2) \gamma_{\nu} E_p (x_2) 
  ( \slashed{p}+ m) \bar{E}_{p} (x_1) \gamma_{\mu} E_{-p_{+}} (x_1)  
\end{array}
\end{equation}

Squaring the matrix element yields $\sim 100000$ scalar product terms and four integrations over contributions from 
bunch field photons. Simplifications of the problem were sought.
Beamstrahlung photons are usually emitted in forward directions, so all photons in the pair production process
were considered collinear. Furthermore, a centre of mass-like reference frame (Equation~\ref{eq6}) was chosen
which permitted scalar products of one of the beamstrahlung photons and bunch field photons to be zero. Two of
the integrations over bunch field photons were so eliminated and the square of the matrix element contained
products of two Airy functions.

\begin{equation}
\label{eq6}
\begin{array}{rl}
  \sql{k}_1 + \sql{k}_2 + l \sql{k} &= \squ{q}_- + \squ{q}_+ = 0 \\
  \omega_1 + \omega_2  + l \omega &= \varepsilon_{q_-} + \varepsilon_{q_+} 
\end{array}
\end{equation}

The full cross-section is too bulky to write down here. Nevertheless, numerical results were obtained using the
Feyncalc Mathematica package and custom written FORTRAN routines.

The appearance of cross-section resonances was established by examining the propagator
denominator. With $\theta_f$ being the angle which the produced electron makes with the propagation direction of
the bunch field photons, and $\omega,\omega_1$ being the beamstrahlung photon and bunch field photon energy
respectively, resonance was found to occur whenever the total energy of bunch field photons contributed to the
process was the same as the beamstrahlung photon energy

\begin{equation}
\label{eq7}
  (\omega_1 - r \omega)(\cos\theta_f + \omega_1 - r \omega)
\end{equation}

The numerical value of the resonant differential cross section was thereby established (Figure~\ref{fig4}) and integration over all
production angles gave the full cross section (Figure~\ref{fig5}). Comparison was made with the ordinary, non
external field Breit-Wheeler cross-sections The differential cross section for the
collinear case of the Stimulated Breit-Wheeler process is peaked in longitudinal directions, as is the ordinary Breit Wheeler
process. However the overall values are greatly enhanced. Indeed for the beamstrahlung photon energies
expected at ILC collisions the enhancement is between five and six orders of magnitude.

\begin{figure}[htbp]
\epsfxsize=9cm
\centerline{\epsfbox{diff_xsect_final.eps}}
\caption{The Stimulated Breit-Wheeler differential cross-section peak.}
\label{fig4}
\end{figure}

\begin{figure}[htbp]
\epsfxsize=9cm
\centerline{\epsfbox{full_xsect_final.eps}}
\caption{The full cross-section of the Stimulated Breit-Wheeler process.}
\label{fig5}
\end{figure}

\section{CONCLUSION}

The cross section of the Stimulated Breit-Wheeler process was calculated using solutions of the Dirac equation
in a constant crossed electromagnetic field and including as a radiative correction, the one loop 
Electron Self Energy in the bunch field. For default accelerator parameters at the ILC, the Breit-Wheeler cross
section used in background pair calculations is underestimated by several orders of magnitude. Background pair
generators used in simulation studies need to be reconfigured in order to included the amended cross section.
The extra background pairs that can be expected will be produced longitudinally and will be twice as energetic
due to the contribution from the bunch fields. The theoretical calculation presented here can be further refined
by allowing for non-azimuthal symmetry in which the pair production takes place at the periphery of the bunches
rather than on axis.

\begin{thebibliography}{9}   % Use for  1-9  references
%\begin{thebibliography}{99} % Use for 10-99 references

\bibitem{oleinik}
V.P. Oleinik, Soviet Phys JETP, 26 (1968) 1132

\bibitem{bos}
J. Bos et al, J Phys A, 12 (1979) 2573

\bibitem{tor}
T. Raubenheimer, "Suggested ILC Beam Parameter Range"

\bibitem{volk}
D.M. Volkov, Zeitschrift f{\"u}r Physik, 94 (1935) 250.

\bibitem{Ritus72}
V.I. Ritus, Ann Phys, 69 (1972) 555.

\bibitem{kallen} 
G. K{\"a}llen, Elementary Particle Physics, Addison-Wesley 1964



\end{thebibliography}

\end{document}