Laser pulse-length effects in trident pair production

TL;DR

This paper investigates how the length and shape of laser pulses influence the trident pair production process, highlighting the role of pulse envelopes and bandwidth effects within a quantum electrodynamics framework.

Contribution

It isolates the impact of laser pulse envelopes on trident pair production using a Furry picture approach, providing a clear analysis of pulse length effects and bandwidth phenomena.

Findings

Pulse envelope Fourier transform encodes pulse length information.

Monochromatic limit can be derived transparently.

Bandwidth effects for short pulses can be separated and studied.

Abstract

Laser pulses facilitate multiphoton contributions to the trident pair production $e_L^- \to e_L^- + e_L^- + e_L^+$ , where the label $L$ indicates a laser field dressed electron ($e^-$) or positron ($e^+$). We isolate the impact of the pulse envelope in the trident S matrix element, formulated within the Furry picture, in leading order of a series expansion in the classical non-linearity parameter $a_0$. Generally, the Fourier transform of the envelope carries the information on the pulse length, which becomes an easily tractable function in the case of a $\cos^2$ pulse envelope. The transition to a monochromatic laser wave can be handled in a transparent manner, as also the onset of bandwidth effects for short pulses can be factorized out and studied separately.

Figures

Equations

$$\displaystyle\begin{minipage}{95.39693pt} \includegraphics[width=433.62pt]{pertC1} \end{minipage}\ +\begin{minipage}{95.39693pt} \includegraphics[width=433.62pt]{pertC2} \end{minipage}$$

$$+$$

$$S_{fi}=e^2\int d^{4}x\int d^{4}y$$

$$\displaystyle-(p_{1}\leftrightarrow p_{2})\bigg{]},$$

$$S_{fi}$$

$$\displaystyle\times\delta^{(4)}(p_{1}+p_{2}+p_{3}-p-(r+s)k)-(p_{1}\leftrightarrow p_{2})\bigg{]}$$

$${\Gamma }_0^{\mu }(r,q_1,q_2)$$

$${\Gamma }_1^{\mu }(r,q_1,q_2)$$

$${\Gamma }_2^{\mu }(r,q_1,q_2)$$

$${\Delta }^{\mu }(r,\ast )=\sum _{l=0}^2{B}_l(r,\ast )J_l^{\mu }(\ast ),$$

$$J_0^{\mu }(q_1,q_2)$$

$$J_1^{\mu }(q_1,q_2)$$

$$J_2^{\mu }(q_1,q_2)$$

$$B_l(s,\ast )={\int }_{-\infty }^{\infty }d\varphi f^l(\varphi )\exp \left\{i\sum _{n=0}^2{\alpha }_n(s,\ast ){\int }_0^{\varphi }d{\varphi }^{\prime }{f}^n({\varphi }^{\prime })\right\},$$

$${\alpha }_0(s;q_1,q_2)$$

$${\alpha }_1(s;q_1,q_2)$$

$${\alpha }_2(s;q_1,q_2)$$

$$B_0(s,\ast )$$

$$\mathcal{G}$$

$$g_{\mu \nu }{\Delta }^{\mu }(r;C){\Delta }^{\nu }(s;BW)$$

$$+\pi \mathcal{G}\delta (r)M_{11}(s)+\pi \mathcal{G}\delta (s)M_{12}(r)$$

$$+M_2(r,s),$$

$$M_0$$

$$M_{11}(s)$$

$$M_{12}(r)$$

$$M_2(r,s)$$

$${\tilde{B}}_0(s)$$

$${\tilde{B}}_j(s)$$

$$\mathcal{G}$$

$$F(s)={\int }_{-\infty }^{\infty }d\varphi f(\varphi )\exp \{is\varphi \}$$

$$M_0$$

$$M_{11}(s)$$

$$M_{12}(r)$$

$$M_2(r,s)$$

$$S_s$$

$$S_r$$

$$S_s$$

$$S_r$$

$$r=s=\frac{p^{+}-p_1^{+}-p_2^{+}-p_3^{+}}{k^{+}}\ne 0,$$

$$S_{fi}$$

$$M(s;BW)$$

Full text

1 Introduction

High-intensity laser beams in the optical regime are customarily generated by the chirped pulse amplification (cf. [1]). Intensities up to $10^{22}\,\nicefrac{\mathrm{W}}{\mathrm{cm^{2}}}$ are achievable nowadays in several laboratories [2], yielding a classical non-linearity parameter of $a_{0}=\mathcal{O}(10-50)$ in the focal spot111The relation of $a_{0}$ vs. the laser peak intensity $I_{L}$ and frequency $\omega$ reads $a_{0}\simeq 7.5\frac{\text{eV}}{\omega}\sqrt{\frac{I_{L}}{10^{20}\text{W}/\text{cm}^{2}}}$, cf. [3].. Ongoing projects [4, 5, 6] of $10\,\mathrm{PW}$ class lasers envisage even larger values of $a_{0}$. Due to higher frequencies in XFEL beams, $\omega=\mathcal{O}(10\,\mathrm{keV})$, the parameter $a_{0}$ stays significantly below unity, despite similar intensities of $\mathcal{O}(10^{22}\,\nicefrac{\mathrm{W}}{\mathrm{cm^{2}}})$ when tight focusing is attained [7]. Given such a variety of laser facilities, the experimental exploration of nonlinear QED effects became feasible and is currently further promoted. Elementary processes are under consideration with the goal of testing QED in the strong-field regime. Most notable is the nonlinear Compton process $e_{L}^{-}\to e_{L}^{-}+\gamma$, also w.r.t. the subsequent use of the high energy photons ($\gamma$), up to prospects of industrial applications. While in the pioneering theoretical studies [8, 9] the higher harmonics, related to multi-photon effects, i.e. the simultaneous interaction of the electron with a multitude of photons, in monochromatic laser beams have been considered, the study of laser pulses revealed a multitude of novel structures in the $\gamma$ spectrum [3, 10, 11, 12, 13, 14, 15, 16, 17, 18].

The nonlinear Breit-Wheeler process, $\gamma\to e_{L}^{-}+e_{L}^{+}$ [19, 20, 3, 21, 13, 22, 23, 24, 25, 26, 27, 28, 29], as cross channel of the nonlinear Compton process, is in contrast a threshold process – sometimes termed a genuine quantum process – since the probe photon $\gamma$ energy in combination with the laser must supply sufficient energy to produce a $e^{+}e^{-}$ pair. When considering the seminal SLAC experiment E-144 [30, 31] as a two-step process (first step: generation of a high-energy photon $\gamma$ by Compton backscattering [32], second step: Breit-Wheeler process $\gamma+L\to e^{+}e^{-}$), also the nonlinear Breit-Wheeler process has been identified with the simultaneous interaction of up to five photons in the elementary subprocess.

Strictly speaking, the mentioned two-step process is only a part of trident pair production $e_{L}^{-}\to e_{L}^{-}+e_{L}^{+}+e_{L}^{-}$, as stressed in [33, 34]. Since the trident process is the starting point of seeded QED avalanches, expected to set in at high-intensities, it is currently a subject of throughout analyses [35, 36, 37], also for benchmarking PIC codes [38].

Given the high repetition rate of the European XFEL [39] a potentially interesting option is to combine it with a synchronized electron beam of about $50\,\mathrm{MeV}$ (to operate slightly above the threshold) in order to facilitate a high-statistics search for the dark photon. Such a dark photon (also dubbed $U$ boson or hidden photon) is a candidate for Dark Matter beyond the standard model of particle physics; it is a possible extension which enjoys intense theoretical [40, 41, 42] and experimental [43, 44, 45, 46] considerations. A corresponding theoretical analysis of the trident process can be found in [47]. In fact, the trident process – in a perturbative QED language – includes sub-diagrams of the type $\gamma^{*}\to e^{+}e^{-}$, i.e. an intermediate (virtual) photon which decays into a $e^{+}e^{-}$ pair. Via kinetic mixing, that virtual photon may “temporarily” couple to a dark photon $A^{\prime}$, e.g. $\gamma^{*}\to A^{\prime}\to\gamma^{*}$, thus signalizing its presence as a peak of the invariant mass distribution of $e^{+}e^{-}$. The peak would be at the mass of the dark photon and its width is related to the kinetic mixing strength.

We briefly mention the trident option of the LUXE project [48, 49] at DESY/Hamburg, which however is primarily dedicated to explore the “boiling of the vacuum” by means of the nonlinear Breit-Wheeler process in the Ritus corner, i.e. a kinematical region with a nonperturbative field strength dependence and coupling constant $\left|e\right|$ dependence analog to the Schwinger pair creation probability.

While most of the above quoted papers focus on nonlinear effects in strong laser pulses, that is the impact of multiphoton contributions, we aim here at the study of apparent multiphoton effects due to bandwidth effects in weak and moderately strong laser pulses with $a_{0}<1$. The analysis of the Breit-Wheeler pair production in [22, 50] revealed that in such a regime interesting features appear for short and ultra-short pulses. For instance, despite of $a_{0}<1$ a significant subthreshold pair production is enabled. Roughly speaking, in short pulses the frequency spectrum contains high Fourier components, thus enabling the subthreshold pair creation. In that respect, we are going to study the relevance of the pulse duration for the trident process. In contrast to the elementary one-vertex processes, the trident process as a two-vertex process obeys a higher complexity, similar to the two-photon Compton scattering [51, 52, 53].

Our paper is organised as follows. In section 2, the matrix element is evaluated with emphasis on a certain regularisation required to uncover the perturbative limit. The weak-field expansion is presented in section 3, where the Fourier transform of the background field amplitude is highlighted as a central quantity. The case of a $\cos^{2}$ envelope is elaborated in dome detail in section 4, where also numerical examples are exhibited. The conclusion can be found in section 5.

2 Matrix element in the Furry picture

The leading-order tree level Feynman diagram of the trident process is exhibited in figure 1(a). The corresponding $S$ matrix reads

[TABLE]

where $\psi(x;p)$ stands for the Volkov solution of Dirac’s equation with a classical external electromagnetic (laser) field $A_{\mu}(\phi)=a_{0}\frac{m}{\left|e\right|}\,\varepsilon_{\mu}f(\phi)$, $\overline{\psi}$ its adjoint, and $D_{\mu\nu}$ is the photon propagator. The $p_{1}\leftrightarrow p_{2}$ term ensures the antisymmetrization of two identical fermions (mass $m$, charge $\left|e\right|$) in the final state. The laser field $A_{\mu}$ and its polarisation four-vector $\varepsilon_{\mu}$ and phase $\phi=k\cdot x$ is specialized further on below. The momenta $p,p_{1,2,3}$ and $k$ are four-vectors as well, and $\gamma^{\mu}$ stands for Dirac’s gamma matrices. Transforming the photon propagator into momentum space, $D_{\mu\nu}(x,y)=\int d^{4}k^{\prime}/\left(2\pi\right)^{4}\,e^{-ik\cdot(x-y)}D_{\mu\nu}(k^{\prime})$ and employing the Feynman gauge, $D_{\mu\nu}(k^{\prime})=\frac{g_{\mu\nu}}{k^{{}^{\prime}2}+i\epsilon}$, and a suitable splitting of phase factors of the Volkov solution, e.g. $\psi(x;p)=\left(1+e\frac{\not{k}\not{A}}{2p\cdot k}\right)u_{p}\,e^{-i\hat{S}_{p}(k\cdot x)}\,e^{-ip\cdot x}$, one can cast the above matrix element in the form

[TABLE]

upon Fourier transform of the vertex function $\left(1+\overline{\Omega}_{q_{2}}(k\!\cdot\!x)\right)\gamma^{\mu}\left(1+\Omega_{q_{1}}(k\!\cdot\!x)\right)\\ \times e^{-i(\hat{S}_{q_{1}}(k\cdot x)-\hat{S}_{q_{2}}(k\cdot x))}=\int\frac{dr}{2\pi}\Gamma^{\mu}(r;q_{1},q_{2})e^{-ir\,k\cdot x}$. A key for that is $\Omega_{p}(\phi=k\!\cdot\!x)=e\frac{\not{k}\not{A}}{2p\cdot k}$ as well as $\hat{S}(\phi\!=\!k\!\cdot\!x;p)=\frac{1}{2k\cdot p}\int_{0}^{k\cdot x}d\phi^{\prime}\left(2e\,p\!\cdot\!A(\phi^{\prime})-e^{2}A^{2}(\phi^{\prime})\right)$ as the nonlinear Volkov phase part. The quantities $u(p)$ and $v(p)$ are a free-field Dirac bispinors, with spin indices suppressed for brevity. In intermediate steps, one meets the local energy-momentum balance $p_{1}-p+k^{\prime}-sk=0$ and $p_{2}+p_{3}-k^{\prime}-rk=0$, which combine to the overall conservation in $\delta^{(4)}$. The corresponding representation of the matrix element in momentum space is exhibited in figure 1(b)-left, where we exposed the interaction with $s$ and $r$ laser photons marked by the crosses. In figure 1(b)-right, we suppress these explicit representations of the laser background field. Note that the momentum space diagrammatics differs from the notation in [54, 55]. We introduce the short-hand notations $C$ and $BW$ to mean momentum dependences on $(p_{1},-p)$ and $(p_{2},p_{3})$, respectively. For the currents $\Delta^{\mu}(r,C)=\overline{u}(p_{1})\Gamma^{\mu}u(p)$ and $\Delta^{\mu}(s,BW)=\overline{u}(p_{2})\Gamma^{\mu}v(p_{3})$ in (2), we note the decomposition, emerging from inserting the Volkov solution, $\Gamma^{\mu}(r,*)=\Gamma_{0}^{\mu}(r,*)+\Gamma_{1}^{\mu}(r,*)+\Gamma_{2}^{\mu}(r,*)$ with $*$ meaning $C$ or $BW$ and, by using a generic momentum pair $\left(q_{1},q_{2}\right)$ for them,

[TABLE]

which can be combined to arrive at the notation such as those of [34, 36]:

[TABLE]

with the currents

[TABLE]

where $\Psi(q)=u(q)$ for $C$ and $\Psi(q)=v(q)$ for $BW$, respectively. The phase integrals in equation (4) read

[TABLE]

with

[TABLE]

Note that in (6), $l$ is an index (label) on the l.h.s, while it is a power on the r.h.s., as $n$ too. An important step is the isolation of the divergent part in $\Gamma_{0}^{\mu}$ or $B_{0}$. We note $\lim_{A\to 0}\Gamma^{\mu}(s,*)=2\pi\delta(s)\gamma^{\mu}$ and regularise $\Gamma^{\mu}_{0}$ by inserting a damping factor $e^{-\epsilon\left|\phi\right|}$ and performing the limit $\epsilon\to 0$, similar to the method in [56], which results in $\Gamma^{\mu}_{0}(s,q_{1},q_{2})=\mathcal{G}\pi\,\delta(s)\gamma^{\mu}+\overset{12}{\Gamma^{\mu}_{0}},\text{ with }\overset{12}{\Gamma^{\mu}_{0}}=-\mathcal{P}\frac{1}{s}\int_{-\infty}^{\infty}d\phi\,\gamma^{\mu}e^{is\phi}e^{-i(\hat{S}_{q_{1}}(\phi)-\hat{S}_{q_{2}}(\phi))}\frac{\partial}{\partial\phi}(\hat{S}_{q_{1}}(\phi)-\hat{S}_{q_{2}}(\phi))$ or

[TABLE]

where $\mathcal{P}$ means the principal value in the variable $s$. One can exploit in reading (4) the crossing symmetry $\Delta^{\mu}(r,BW=(p_{2},p_{3}))=\Delta^{\mu}(r\to s,BW\to C=(p_{1},-p))$. Employing (8) with the short-hand notation $\tilde{B}_{0}=-\mathcal{P}\left[\frac{1}{s}\sum_{j=1}^{2}\alpha_{j}B_{j}(s,*)\right]$ and $\tilde{B}_{1,2}=B_{1,2}$ we arrive at

[TABLE]

where

[TABLE]

3 Weak-field expansion

With the argumentation given in the introduction we now attempt an expansion in powers of $a_{0}$. We note $J^{\mu}_{l}\propto a_{0}^{l}$ and $\alpha_{l}\propto a_{0}^{l}$; again, $l$ is a label (power) on the l.h.s. (r.h.s.). The leading-order terms of the phase integrals (6) thus become

[TABLE]

Denoting the Fourier transform of $f(\phi)$ by

[TABLE]

we recognize that

[TABLE]

The two delta distributions in the $M_{0}$ term in (10) enforce for the overall momentum conservation in (2) a factor $\delta^{(4)}(p_{1}+p_{2}+p_{3}-p)$ implying a zero contribution of the $M_{0}$ term (14a). In the spirit of the $a_{0}$ series expansion we neglect (14d) at all. (This term would give rise to on/off-shell contributions which require some care.) The remaining leading order terms in (14b) and (14c) generate the contributions

[TABLE]

plus the corresponding exchange terms upon $p_{1}\leftrightarrow p_{2}$. Introducing light-front coordinates222We use the definition $q^{\pm}=\frac{1}{2}\left(q^{0}\pm q^{3}\right)$ and $q^{\perp}=(q^{1},q^{2})$ for the light-front coordinates of a four-vector $q^{\mu}=(q^{0},q^{1},q^{2},q^{3})$. and the light-front delta distribution $\delta^{\mathrm{lf}}(p)=\delta(p^{-})\delta^{(2)}(p^{\perp})$ as well as choosing $k^{+}$ as the non-zero component of the laser four-momentum $k^{\mu}$ yields

[TABLE]

where $M_{11}$ refers to (11b) and $M_{12}$ to (11c) and one has to use in both cases

[TABLE]

and the principal value can be dropped due to the last inequality. For given entry channel parameters, equation (17) implies the dependence $s(E_{2,3},\cos\theta_{2,3},\varphi_{2,3})$ with spherical coordinates $E_{i},\cos\theta_{i},\varphi_{i}$ for the particles $1,2$ and $3$ (cf. figure 1(a)) in the exit channel. Fixing $E_{3},\cos\theta_{2,3}$ and $\varphi_{3}$ yields the contour plot $s(E_{2},\varphi_{2})$. An example is exhibited in figure 3. The locus of the apparent one- (two-) photon contribution with $s=1$ ($=2$) is highlighted by fat curves.

The final result is the leading-order matrix element

[TABLE]

with

[TABLE]

where the tildes indicate here that the factor $a_{0}$ is scaled out. These structures are suggestive: $M(BW)$ may be read as the coupling of the free Compton current $J_{0}^{\mu}(C)$ to the modified Breit-Wheeler current (in parenthesis of (19a)) and a free Breit-Wheeler current $J_{0}^{\mu}(BW)$ to a modified Compton current (in parenthesis of (19b)), both depicted in the middle panels of figure 1(c). The interaction with the external field is encoded in the modified currents. For practical purposes we replace the modified currents by the proper Volkov currents $\Delta^{\mu}(C)$ and $\Delta^{\mu}(BW)$ when evaluating numerically the matrix elements for $a_{0}\ll 1$. The differential probability is

[TABLE]

with three-body phase space element

[TABLE]

the Heaviside step-function $\Theta$ and the differential cross section is

[TABLE]

with the light-front volume $V_{\mathrm{lf}}$ and the normalisation factor $N=\frac{a_{0}^{2}m^{2}}{2e^{2}}\int_{-\infty}^{\infty}g^{2}(\phi)\,d\phi$ (cf. [16]). The matrix element $\mathcal{M}$ is given by (18) but without the pre-factor:

[TABLE]

We emphasize that the weak-field limit (18) in (20) corresponds to the standard perturbative tree level QED diagrams depicted in the figure 2, supposed

[TABLE]

$s=1$ selects then the admissible kinematics, and the phase space in (20) becomes five-dimensional. The decomposition (8) is essential for catching the proper weak-field perturbative limit. This is obvious when considering Møller or Bhabha scattering in an ambient background field as the cross channels of the trident process: For $A\to 0$ the standard perturbative QED must be recovered.

4 The case of a $\mathbf{\cos^{2}}$ envelope

We now specialise the linearly polarised external e.m. field $A^{\mu}$ with $\varepsilon^{\mu}=\left(0,1,0,0\right)$ to a $\cos^{2}$ envelope times the oscillating part yielding

[TABLE]

where $\bigsqcap\left(\frac{\phi}{2\Delta\phi}\right)$ is a box profile of width $2\Delta\phi$ centred at $\phi=0$. We leave a discussion of the carrier envelope phase $\phi_{\mathrm{cep}}$ for separate work, i.e. put $\phi_{\mathrm{cep}}=0$. Then (13) can be integrated analytically with the result

[TABLE]

where $\operatorname{sinc}(x)=\frac{\sin(x)}{x}$ for $x\neq 0$ and $\operatorname{sinc}(0)=1$ is the cardinal sine function. With a proper normalisation, (20) is proportional to $\frac{\left|F(s)\right|^{2}}{\Delta\phi}$. $F(s)$ is a real function due to the even symmetry of the special field (25), but in general it aquires also an imaginary part. We therefore keep the notation $\left|F(s)\right|^{2}$.

Given the properties of the sinc functions entering (26) we find

[TABLE]

thus making (24) explicit. Since $\frac{\left|F(s)\right|^{2}}{\Delta\phi}>0$ for $s\neq 1$, in particular for $s>1$, at finite values of $\Delta\phi$, we see that this signals bandwidth effects, despite $a_{0}\ll 1$. Such effects have been observed in [50] for the Breit-Wheeler pair production below the threshold and the Compton process as well [22, 23]. Since $s$ is a continuous variable one must not identify it with a “photon number”; instead, $s$ could be interpreted as fraction of energy or momentum in units of $\omega=\left|\vec{k}\right|$ participating in creating a final state different from the initial state (see [9] for discussions of that issue). Even more, $s>1$ does not mean proper multiphoton effects due to our restriction on leading order in $a_{0}$, rather one could speak an “apparent multiphoton effects” caused by finite bandwidth of the pulse.

The dependence of $\frac{\left|F(s)\right|^{2}}{\Delta\phi}$ on $s$ for various values of $\Delta\phi$ is displayed in figure 4. The function is symmetric, $\left|F(s)\right|^{2}=\left|F(-s)\right|^{2}$, with main maxima at $s=\pm 1$ and the envelopes are

[TABLE]

We use the absolute value of the analytic signal of $F(s)$, i.e. the Fourier transform of the amplitude function (25) constrained to the positive half-line [57]. The first side maxima are located between $1\pm 2\frac{\pi}{\Delta\phi}$ and $1\pm 3\frac{\pi}{\Delta\phi}$ and their heights are $7\times 10^{-4}$ of the respective main maximum, meaning that their contribution is not negligible, in particular for a kinematical situation where $s>1$ (cf. figure 3)

From the definition of $s$ in equation (17), we note (i) the dependence $s(E_{2},E_{3},\theta_{2,3},\varphi_{2,3})$ and (ii) the U shape of $s(E_{2})$ when keeping constant the other variables of the momenta in polar coordinates. Denoting the minimum of $s(E_{2})$ by $s_{\mathrm{min}}$, then for kinematical situations, where $s_{\mathrm{min}}>1$ the side peaks of $\frac{\left|F(s)\right|^{2}}{\Delta\phi}$ with spacing of about $\frac{\pi}{\Delta\phi}$ become relevant. In fact, $\frac{\left|F(s)\right|^{2}}{\Delta\phi}$ carries much of the energy dependence of the differential cross section. As an example, we exhibit in figure 5 (bottom) the differential cross section $d^{6}\sigma/\left(dE_{2}\,d\cos\theta_{2}\,d\varphi_{2}\,dE_{3}\,d\cos\theta_{3}\,d\varphi_{3}\right)$ (multiplied by $m^{4}$ to make it dimensionless) as a function of $E_{2}$ and compare it with $\frac{\left|F(s(E_{2}),\Delta\phi)\right|^{2}}{\Delta\phi}$ scaled by a factor $5.2\times 10^{-7}$. Note the near-perfect match.333Supposed the match of the cross section (22) and $\frac{\left|F(s,\Delta\phi)\right|^{2}}{\Delta\phi}$ scaled by one common factor continues over a wide region in phase space one could envisage a greatly simplified numerical procedure by evaluation (or even estimating) (22) only a few times and continue it with $\frac{\left|F(s,\Delta\phi)\right|^{2}}{\Delta\phi}$.

The selected kinematics implies $s_{\mathrm{min}}\approx 1.1$, i.e. in the limit $\Delta\phi\to\infty$, this setting would be kinematically forbidden, but bandwidth effects for finite values of $\Delta\phi$ enable the selected kinematics. In fact, increasing $\Delta\phi$ causes $(i)$ a rapid dropping of the differential cross section and $(ii)$ make the oscillatory pattern more dense. At the heart of the behavior is essentially the quantity $\frac{\left|F(s,\Delta\phi)\right|^{2}}{\Delta\phi}$ of figure 4; the additional $E_{2}$ dependence is fairly smooth. To translate $\frac{\left|F(s,\Delta\phi)\right|^{2}}{\Delta\phi}$ of figure 4 into the bottom panels of figure 5 one needs explicitly $s(E_{2})$ which can be inferred from figure 3.

Displaying the differential cross section as contour plot over the $E_{2}-E_{3}$ plane (see figure 5-top) and keeping $\theta_{2,3}$ and $\varphi_{2,3}$ fixed as above, one observes a pronounced fringe pattern which become denser with increasing values of $\Delta\phi$. The fringe pattern occupies a finite region in the $E_{2}-E_{3}$ plane. At the origin of the fringe pattern is again the function $\frac{\left|F(s(E_{2}),\Delta\phi)\right|^{2}}{\Delta\phi}$. The fringe pattern make numerical integrations towards total cross section fairly challenging.

The relevance of the function $\frac{\left|F\right|^{2}}{\Delta\phi}$ continues of course for kinematical situations which are not forbidden in the limit $\Delta\phi\to\infty$. Examples are exhibited in figure 6. For the selected kinematic situation, $s_{min}=0.976$, meaning that with increasing values of $\Delta\phi$ the differential cross section does not drop but gets concentrated at two values of $E_{2}$ (at given $E_{3}$) which are allowed for $s=1$. In the limit $\Delta\phi\to\infty$, two delta peaks arise when $E_{3}$ is appropriately fixed. They are the result of cutting the strength distribution over the $E_{2}-E_{3}$ plane at $E_{3}=\text{const.}$, i.e. there is a sharp ring given by the solution of $E_{3}(E_{2})$. The tilt of the U shaped function $s(E_{2})$ at fixed other parameters to the right makes the region $s\approx 1$ larger at large but finite values of $\Delta\phi$. Correspondingly, the r.h.s peak structure is wider, as seen e.g. in the right bottom panel.

To compare with the standard perturbative QED (pQED) result, based on the evaluation of the Feynman diagrams depicted in figure 2, one has to perform the $\varphi_{2}$ integration. In fact, making $\Delta\phi$ larger and larger, the differential cross section $d^{5}\sigma/\left(dE_{2}\,d\cos\theta_{2}\,dE_{3}\,d^{2}\Omega_{3}\right)$ from equation (22) approaches the pQED result which is for a monochromatic photon beam, see figure 7. This figure illustrates how the pQED is approached for $a_{0}\ll 1$ and $\Delta\phi\to\infty$. It also demonstrates that for $\Delta\phi<100$ a significantly larger phase space beyond the perturbatively accessible (one-photon) region (indicated by $\longleftrightarrow$) is occupied due to bandwidth effects in laser pulses.

5 Summary

The length of laser pulses has a decisive impact on the pair production in the trident process. The rich phase space patterns, already found and analysed in some detail for non-linear one-vertex processes à la Compton and Breit-Wheeler, show up also in the two-vertex trident process. Even for weak laser fields, a region becomes accessible which would be kinematically forbidden in a strict perturbative, leading-order tree level QED approach. The key is the frequency distribution in a pulse which differs significantly from a monochromatic laser beam, which would mean an “infinitely long laser pulse”. By resorting to a special pulse model $\propto\cos^{2}\left(\pi\phi/2\Delta\phi\right)\cos(\phi)$, we quantify in a transparent manner the effect of the pulse duration $\Delta\phi$ and identify the transition to monochromatic laser fields, $\Delta\phi\to\infty$, for small values of the classical laser non-linearity parameter. To ensure the contact to a perturbative QED approach the proper decomposition of a phase factor is mandatory. To be specific, the first term in (8) is essential for the two-vertex process; for one-vertex processes it does not contribute. The effect of short laser pulses manifests itself in bandwidth effects mimicking apparent multiphoton contributions even for weak fields, where the intermediate photon is off-shell and and the process is of one-step nature. The rich pattern of the differential phase space distribution is traced back to the Fourier transform of the external e.m. field. In obvious further extensions of our approach, more general pulse envelope shapes should be studied, e.g. within the slowly varying envelope approximation. Going to larger values of the classical laser non-linearity parameter means checking whether an analog of the Fourier transform of the e.m. field can be isolated as crucial element of the phase space distribution of produced particles444In fact, for $a_{0}<0.01$ and large values of $\Delta\phi$ we find numerical agreement of (22) and the pQED result.. The final state phase distribution is important for planing corresponding experimental designs.

The authors gratefully acknowledge the collaboration with A. Otto and very useful discussions with F. Mackenroth. We appreciate conversations with our colleges B. King, A. M. Fedotov, G. Torgrimsson, A. Di Piazza, A. Ilderton, T. Heinzl, A. Hartin, A. I. Titov, D. Seipt, M. Bussmann and T. Nousch. The work is supported by R. Sauerbrey and T. E. Cowan w.r.t. the study of fundamental QED processes for HIBEF.