# Impact of the laser spatio-temporal shape on Breit-Wheeler pair production

A. Mercuri-Baron, M. Grech, F. Niel, A. Grassi, M. Lobet, A. Di Piazza, C. Riconda

► To cite this version:

A. Mercuri-Baron, M. Grech, F. Niel, A. Grassi, M. Lobet, et al. Impact of the laser spatio-temporal shape on Breit-Wheeler pair production. New Journal of Physics, 2021, 23 (8), pp.085006. 10.1088/1367-2630/ac1975. hal-03261268

HAL Id: hal-03261268

https://hal.science/hal-03261268

Submitted on 1 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SUBJECT • OPEN ACCESS

# Impact of the laser spatio-temporal shape on Breit–Wheeler pair production

To cite this article: A Mercuri-Baron et al 2021 *New J. Phys.* **23** 085006

View the [article online](#) for updates and enhancements.

## You may also like

- Laser pulse-length effects in trident pair production
U Hernandez Acosta and B Kämpfer

- Nonlinear Breit–Wheeler pair creation with bremsstrahlung rays
T G Blackburn and M Marklund

- PICSAR-QED: a Monte Carlo module to simulate strong-field quantum electrodynamics in particle-in-cell codes for exascale architectures
Luca Fedeli, Neül Zaïm, Antonin Sainte-Marie et al.
<figure><img src="image_2.png" /></figure>

CrossMark

OPEN ACCESS

RECEIVED

24 May 2021

REVISED

26 July 2021

ACCEPTED FOR PUBLICATION

30 July 2021

PUBLISHED

18 August 2021

Original content from
this work may be used
under the terms of the
Creative Commons
Attribution 4.0 licence.

Any further distribution
of this work must
maintain attribution to
the author(s) and the
title of the work, journal
citation and DOI.

<figure><img src="image_3.png" /></figure>

# Impact of the laser spatio-temporal shape on Breit-Wheeler pair production

A Mercuri-Baron<fnref n="1" />,*, M Grech<fnref n="2" />, F Niel<fnref n="1" />, A Grassi<fnref n="1" />, M Lobet<fnref n="3" />, A Di Piazza<fnref n="4" /> and C Riconda<fnref n="1" />

E-mail: anthony.mercuri@live.fr

**Keywords:** strong field quantum electrodynamics, pair creation, nonlinear Breit–Wheeler, OAM, Laguerre–Gauss laser pulse

## Abstract

The forthcoming generation of multi-petawatt lasers opens the way to abundant pair production by the nonlinear Breit–Wheeler process, i.e. the decay of a photon into an electron–positron pair inside an intense laser field. In this paper we explore the optimal conditions for Breit–Wheeler pair production in the head-on collision of a laser pulse with gamma photons. The role of the laser peak intensity versus the focal spot size and shape is examined keeping a constant laser energy to match experimental constraints. A simple model for the soft-shower case, where most pairs originate from the decay of the initial gamma photons, is derived. This approach provides us with a semi-analytical model for more complex situations involving either Gaussian or Laguerre–Gauss (LG) laser beams. We then explore the influence of the order of the LG beams on pair creation. Finally we obtain the result that, above a given threshold, a larger spot size (or a higher order in the case of LG laser beams) is more favorable than a higher peak intensity. Our results match very well with three-dimensional particle-in-cell simulations and can be used to guide upcoming experimental campaigns.

## 1. Introduction

Abundant electron–positron pair production is among the most exotic and striking processes of extremely high-intensity laser interaction with light and/or matter that will be made possible by the emerging class of multi-petawatt (multi-PW) laser systems [1]. Under conditions envisioned for forthcoming facilities such as Apollon [2], CoReLS [3], ELI [4], OMEGA-EP OPAL [5], XCELS [6] and zeus [7], the dominant process leading to pair production is nonlinear Breit–Wheeler, i.e. the decay of a high-energy photon, as it interacts with the strong laser field, in an electron–positron pair [8, 9].

In the last decade, various simulation campaigns have been conducted to help design the future experiments to efficiently drive this process, an effort that was largely made possible and supported by particle-in-cell (PIC) simulations accounting for quantum electrodynamics (QED) processes [10–15]. Several configurations have been considered. Electron-seeded electromagnetic/QED cascades in the collision of two counter-propagating laser pulses, a setup originally proposed by Bell and Kirk [16], have attracted particular attention [11, 14, 17–20] as such cascades were identified as a possible limitation on the attainable intensity of high power lasers [21]. This setup was further built upon considering either the use of plasma channels [22] or the collision of multiple laser pulses [23–25]. All these works considered multiple laser beams, but electromagnetic cascades were also predicted considering the direct irradiation of a solid target by a single, extremely intense laser pulse, a situation that leads to the production of dense pair plasmas [26] (see also [27] for an application to laboratory astrophysics). In these cases, the solid target, an overdense plasma acts as a mirror and pair production is efficiently achieved in the field of the incident and reflected laser light at the target front [28].

<footnote n="1">¹ LULI, CNRS, CEA, Sorbonne Université, École Polytechnique, Institut Polytechnique de Paris, F-75255 Paris, France</footnote>

<footnote n="2">² LULI, CNRS, CEA, Sorbonne Université, École Polytechnique, Institut Polytechnique de Paris, F-91128 Palaiseau, France</footnote>

<footnote n="3">³ Maison de la Simulation, CEA, CNRS, Université Paris-Sud, UVSQ, Université Paris-Saclay, F-91191 Gif-Sur-Yvette, France</footnote>

<footnote n="4">⁴ Max Planck Institute for Nuclear Physics, Saupfercheckweg 1, D-69117, Heidelberg, Germany</footnote>

<footnote>* Author to whom any correspondence should be addressed.</footnote>
All the previously discussed configurations considered cascades seeded by electrons initially at rest in the laboratory frame. The seed electrons, strongly accelerated in the electromagnetic fields of laser pulses with extreme intensity (typically ~ $10^{24}$ W cm$^{-2}$), gain energies at the GeV-level (as expected, at such intensities, from the ponderomotive scaling) and are the source (via nonlinear Compton scattering) of the high-energy photons triggering pair production. Another possible path toward pair production in the laboratory relies on an already existing high-energy electron or photon source. The two latter schemes are related, as photons with energy of the order of the electron one are produced via nonlinear Compton scattering of electrons in the intense laser fields [29–33]. Indeed, this configuration was exploited in the seminal E-144 experiment at SLAC [34], where the head-on collision of 46.6 GeV electron beams with TW laser pulses gave the first direct demonstration of nonlinear Breit–Wheeler pair production. However, the laser intensity in this experiment was only marginally relativistic, and very few positrons were produced (~ 100 positrons over nearly 22 000 electron-laser collisions).

Combining GeV-electron beams with PW and multi-PW lasers will allow to extend this study deep into the relativistic regime, as investigated in [35] that numerically demonstrated the possibility of producing extremely bright high-energy photon sources as well as positrons bunches for laser intensities $\gtrsim 10^{22}$ W cm$^{-2}$, in conditions relevant to the Apollon facility. This configuration was further studied theoretically and via PIC simulations in [36–38], while the possibility to produce pairs via nonlinear Breit–Wheeler with high-energy photons from a Bremsstrahlung source was presented in [39].

In this paper, we focus on the study of nonlinear Breit–Wheeler pair production following from the head-on collision of an extremely intense laser (intensities in the range $10^{21}$–$10^{25}$ W cm$^{-2}$) with a burst of gamma photons (with energies ranging from 100 MeV to few 10 s of GeV). In contrast to [39], where the authors focused on the optimization of the photon source, we consider here the high-energy photon burst as given (only its energy will vary), and do not discuss its origin (various sources of high-energy photons have been proposed [25, 39–46]). Rather, we aim to optimize the conditions of interaction with the colliding high-intensity laser. Motivated by experimental constraints, we investigate the optimal conditions for pair production varying the laser polarization, focusing, spatial or temporal profiles, always considering a fixed laser energy.

In particular, special attention will be paid to the use of Laguerre–Gauss (LG) beams, which have recently attracted the interest of the ultra-high intensity laser-plasma community [38, 47–50] for their ability to carry orbital angular momentum (OAM) [51]. These beams, that emerge as eigen-modes of the paraxial equation, have unique properties, such as a ring-shaped intensity distribution and the large OAM, that could have an impact on pair production. It was for instance demonstrated in [38] that the collision of an electron beam with a LG beam can lead to efficient production of gamma rays carrying large OAM, to enhanced secondary radiation emission and pair production. However, that study was performed at constant maximum intensity, so it is hard to distinguish the impact of the increased energy (up to $3 \times$ higher for the LG beam) from the role of the laser spatial profile itself. We will discuss this aspect in detail in this paper.

Because of the remarkable number of free parameters to study (polarization, focusing, spatial and temporal profile) and the complex spatio-temporal dependence of the electromagnetic fields of LG beams, it is useful to develop and validate a reduced analytical model for pair creation (see reference [52] for an analytical study of the Breit–Wheeler pair production process in a Gaussian beam). Starting from the probability for a given high-energy photon to decay in an electron–positron pair while it crosses an extremely intense laser pulse, we propose a simple, yet accurate, model to compute the number of pairs created during this interaction. The computation relies on three assumptions. (i) The electron–positron pair production rates are obtained in the locally constant cross-field approximation [53]. (ii) The newly created charged particles do not further contribute to pair production. (iii) The potentially complex electromagnetic fields explored by the high-energy photons as they cross the laser pulse are approximated, assuming that each photon sees a succession of single half-periods of monochromatic plane waves, the field amplitude of each half-period being computed so that both the temporal and spatial profiles of the laser pulse at focus are accounted for. As will be further discussed in this work, these assumptions are well satisfied for a wide range of parameters of interest for this study and currently accessible laser parameters. Our model goes further than previous studies [36, 53] as it accounts for arbitrary polarization and spatio-temporal profile of the laser beam, and it is not limited to small pair production probability and/or rate.

Three-dimensional PIC simulations performed with the code SMILEI [54] and embarking the relevant QED modules [15] are presented and validate our model. Assumptions (ii) and (iii) are shown to hold under conditions relevant to ultra-short laser pulses (up to few tens laser cycles) typical of the Ti:sapphire technology. Assumption (ii) holding is of particular importance as it demonstrates that, under conditions relevant to forthcoming ultra-short laser facilities, the high-energy-photon-laser interaction proceeds
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Dependence of $b_0(\chi_\gamma)$ on the photon quantum parameter $\chi_\gamma$ as defined by equation (3) (black line) and asymptotic behaviors given by equations (5) and (6) (grey dashed lines). (b) Dependence of $I_e(\chi_0)$ (equation (19)) on the maximum photon quantum parameter $\chi_0$ (as defined by equation (13)), for a linearly polarized (LP) background field ($\epsilon = 0$, black line) and a circularly polarized (CP) background field ($\epsilon = \pm 1$, blue line). Dashed lines show the asymptotic behavior given by equations (22) and (23) for the LP case. The inset highlights the values of $\chi_0$ for which LP ($\chi_0 < 2.52$) or CP ($\chi_0 > 2.52$) gives higher pair production probability and the relative difference of $I_e$ (in absolute value) between the two cases.</figcaption></figure>

essentially in a regime which we refer to as the soft-shower regime [55]. Furthermore, our PIC simulations allow to test our model against complex laser pulse geometries, e.g. considering LG beams, showing that assumption (iii)—which is at the core of the simplification of our model—provides accurate predictions even for non trivial electromagnetic field configurations.

The proposed model and performed study make it possible to clearly identify the laser parameters and high-energy photon energy that optimize pair production for a given laser energy. This work also highlights the conditions in which either stronger focusing or increased gamma photon energy does not improve the production of primary pairs.

The paper is structured as follows. In section 2, we introduce the key properties of the nonlinear Breit–Wheeler process, and we propose a reduced model for pair creation in the head-on collision between gamma photons and a plane wave with arbitrary polarization. The possibility to account for a time envelope of the laser field is also considered and the model is benchmarked against 1D PIC simulations. In section 3 we consider the effect of the spatial dependence of the laser field. Both LG and Gaussian beams are considered and we introduce the notion of a total cross-section for which we propose a reduced model, shown to be in very good agreement with 3D PIC simulations. The optimal conditions for pair creation with respect to the laser geometry and intensity are also discussed. In section 4 we then use our model to discuss pair production on upcoming Ti:sapphire laser facilities, before giving our conclusions in section 5.

## 2. Reduced model for pair production in a time-varying background field

This section provides the theoretical framework for describing nonlinear Breit–Wheeler pair production, and aims to establish a simple model for pair production in a time-varying background field. Throughout this section, we consider the head-on collision of high-energy gamma photons with a strong, optical, laser pulse (also referred to as the background field). In addition to exact expressions, we provide practical formulae relevant to forthcoming high-power laser facilities, focusing on optical laser systems with micrometric wavelength $\lambda$, with relativistic field strength $a_0 = eE_0/(m_e c\omega) \gg 1$ (with $\omega = 2\pi c/\lambda$ the background field angular frequency), corresponding to $I_0\lambda^2$ well beyond $10^{18} \text{ W cm}^{-2} \text{ µm}^{-2}$, with $I_0 = \epsilon_0 c^2/2$ the laser intensity. Note that SI units are used throughout this work: $e$, $m_e$ and $c$ denote the elementary charge, electron mass and speed of light in vacuum, respectively, $\epsilon_0$ the vacuum permittivity, $r_e = e^2/(4\pi\epsilon_0 m_e c^2)$ the classical radius of the electron, $\tau_e = r_e/c$ the time for light to cross $r_e$, $\hbar$ the reduced Planck constant and $\alpha = e^2/(4\pi\epsilon_0\hbar c)$ the fine-structure constant.

### 2.1. Pair production instantaneous rate and probability

Nonlinear Breit–Wheeler is the dominant QED process leading to electron–positron pair production from the collision of a high-energy photon with a strong (optical) laser pulse [56, 57]. As discussed in reference [53], for a relativistically intense background field, $a_0 \gg 1$, the local constant field approximation holds and the pair production probability can be computed by integrating the instantaneous pair production rate. For a high-energy photon propagating with velocity **c**, normalized energy $\gamma_\gamma = \hbar\omega_\gamma/(m_e c^2)$ and instantaneous
quantum parameter

$$ \chi_{\gamma} = \frac{\gamma_{\gamma}}{E_{S}} \sqrt{(\mathbf{E} + \mathbf{c} \times \mathbf{B})^{2} - (\mathbf{c} \cdot \mathbf{E})^{2}/c^{2}}, \qquad (1) $$

where $E_S = m_e^2 c^3 / (e\hbar)$ $\simeq 1.32 \times 10^{18}$ V m$^{-1}$ is the Schwinger field and $\mathbf{E}$ and $\mathbf{B}$ denote the electric and magnetic fields at the photon position. The instantaneous pair production rate is given by (see, e.g. references [58, 59]):

$$ W_{BW} = W_0 \frac{b_0(\chi_\gamma)}{\gamma_\gamma}, \qquad (2) $$

where $W_0 = 2\alpha^2/(3\tau_e)$. For a fixed value of the quantum parameter $\chi_\gamma$, the rate is inversely proportional to the gamma photon energy $\gamma_\gamma$, and depends non-trivially on $\chi_\gamma$ through the function

$$ b_0(\chi_\gamma) = \frac{\sqrt{3}}{2\pi} \int_0^1 \frac{d\xi}{\xi(1-\xi)} \left[ \frac{2}{3\chi_\gamma} \frac{1-2\xi}{(1-\xi)} K_{5/3}(\mu) + K_{2/3}(\mu) \right], \qquad (3) $$

with $\mu = 2/[3\chi_\gamma\xi(1-\xi)]$ and $K_n(x)$ the modified Bessel function of the second kind. This function, shown in figure 1(a), can be well approximated (within less than 1% from the numerically computed values for $\chi_\gamma \in [10^{-2}, 10^3]$) by<fnref n="5" />

$$ b_0(\chi_\gamma) \approx 0.242 \frac{K_{1/3}^2(4/(3\chi_\gamma))}{1 - 0.172/(1 + 0.295 \chi_\gamma^{2/3})}. \qquad (4) $$

The denominator in equation (4) provides a simple improvement of the well known Erber approximation [60], which also allows to recover the correct asymptotic behavior (dashed grey lines in figure 1(a)) of $b_0(\chi)$:

$$ b_0(\chi_\gamma) \rightarrow c_1 \chi_\gamma e^{-8/(3\chi_\gamma)} \quad \text{for } \chi_\gamma \ll 1, \qquad (5) $$

$$ b_0(\chi_\gamma) \rightarrow c_2 \chi_\gamma^{2/3} \quad \text{for } \chi_\gamma \gg 1, \qquad (6) $$

with $c_1 = (3/16)(3/2)^{3/2} \simeq 0.344$ and $c_2 = 3^{2/3}45\Gamma^4(2/3)/(56\pi^2) \simeq 0.569$ ($\Gamma(x)$ denoting the gamma function).

The quantity $W_{BW}$ defined by equation (2) is the instantaneous rate of pair production, so that considering, at a time $t_0$, a given number $N_0$ of identical high-energy photons entering a strong background field, the number of photons $N_\gamma$ remaining at later times $t > t_0$ satisfies the rate equation:

$$ \frac{d}{dt} N_{\gamma} = -W_{BW} N_{\gamma}, \qquad (7) $$

where both $N_\gamma$ and $W_{BW}$ depend on time (the latter when the high-energy photons explore electromagnetic fields with spatial and/or temporal variations). Note that equation (7) does not account for subsequent radiation from the produced pairs, which will be verified to be accurate for the investigated parameter range in section 3. The solution to this equation reads

$$ N_{\gamma}(t) = N_{0} \exp \left[ - \int_{t_0}^{t} W_{BW}(t') dt' \right], \qquad (8) $$

from which one can derive the probability $P(t_0, \Delta t) = 1 - N_\gamma(t_0 + \Delta t)/N_0$ for a given high-energy photon to decay into an electron–positron pair between time $t_0$ and $t_0 + \Delta t$:

$$ P(t_0, \Delta t) = 1 - \exp \left[ - \int_{t_0}^{t_0 + \Delta t} W_{BW}(t') dt' \right]. \qquad (9) $$

When the time-integrated rate (i.e. the integral term in equation (8)) is small, it can be assimilated to the probability $P(t_0, \Delta t)$, as $P(t_0, \Delta t) \rightarrow \int_{t_0}^{t_0 + \Delta t} W_{BW}(t') dt' \ll 1$. This leads some authors to refer to $W_{BW}$ as the pair production probability per unit time. While this is correct in the limit $P(t_0, \Delta t) \ll 1$, it leads to inconsistencies as the probability increases and the time-integrated rate assumes values greater than unity (see, also, references [53, 61–63]). Hence, throughout this work, we will distinguish pair production rate (as defined by equation (2)) and probability (as defined by equation (9)).

<footnote n="5">⁵ Note that $0.242 \sim 0.16 \times 3/2$, and $b_0(\chi) = 3\chi T(\chi)/2$, with $T(\chi)$ the function used by Erber [60] and Blackburn [36].</footnote>
## 2.2. Pair production probability in an arbitrarily polarized plane wave

Let us now focus on the case of a high-energy photon colliding head-on with an ultra-intense laser pulse described by a plane wave (propagating in the $\hat{z}$-direction), with electric field:

$$ \mathbf{E}(z, t) = \frac{E_0}{\sqrt{1 + \varepsilon^2}} \left[ \sin(\omega t - kz) \hat{\mathbf{x}} + \varepsilon \cos(\omega t - kz) \hat{\mathbf{y}} \right], \quad (10) $$

magnetic field $\mathbf{B} = \hat{\mathbf{z}} \times \mathbf{E}/c$, and polarization ellipticity $\varepsilon \in [-1, 1]$. It takes the particular values $\varepsilon = \pm 1$ for a CP wave ($\pm$ standing for opposite helicities), and $\varepsilon = 0$ for a LP wave ($\hat{\mathbf{x}}$ denoting by convention the direction of polarization). It is important to stress that our choice of parametrisation (equation (10)) ensures that the laser energy is the same for any value of $\varepsilon$. On the contrary, $E_0/\sqrt{1 + \varepsilon^2}$ measures the maximum laser electric field amplitude and varies with the ellipticity.

When a high-energy photon collides with such a laser field (we define $t_0 = 0$ as the time at which the photon enters the background field), the value of its quantum parameter

$$ \chi_{\gamma}(t) = \frac{\gamma_{\gamma}}{E_{S}} \left| \mathbf{E}(z_{\gamma}(t), t) + \mathbf{c} \times \mathbf{B}(z_{\gamma}(t), t) \right| \quad (11) $$

evolves in time as the photon explores regions with different electric and magnetic field amplitude (here $z_{\gamma}(t)$ denotes the photon time-dependent position). In the case of head-on collisions ($\mathbf{c} = -c\hat{\mathbf{z}}$), the time-dependent photon quantum parameter reduces to

$$ \chi_{\gamma}(t) = \chi_{0} \Psi_{\varepsilon}(2\omega t) \quad \text{with} \quad \Psi_{\varepsilon}(\varphi) = \sqrt{\frac{\sin^2\varphi + \varepsilon^2 \cos^2\varphi}{1 + \varepsilon^2}}, \quad (12) $$

and

$$ \chi_0 = 2\gamma_\gamma \frac{E_0}{E_S} \approx 0.801 \left( \frac{\hbar\omega_\gamma}{1 \text{ GeV}} \right) \sqrt{\frac{I_0}{10^{22} \text{ W cm}^{-2}}}. \quad (13) $$

Note that, with our convention on the background field parametrization (equation (10)), $\chi_0$ denotes the maximum quantum parameter for the LP case only. For arbitrary $\varepsilon$, the maximum achievable quantum parameter $\chi_m$ is a decreasing function of $|\varepsilon|$

$$ \chi_m = \frac{\chi_0}{\sqrt{1 + \varepsilon^2}}. \quad (14) $$

Note also that equation (12) takes a simple form when considering either a CP or a LP background field:

$$ \chi_{\gamma}(t) = \chi_{0}/\sqrt{2} \quad \text{for } \varepsilon = \pm 1, \quad (15) $$

$$ \chi_{\gamma}(t) = \chi_{0} |\sin(2\omega t)| \quad \text{for } \varepsilon = 0. \quad (16) $$

### 2.2.1. Decay probability after propagating through half a period of the background field

For arbitrary values of $\varepsilon$, $\chi_{\gamma}(t)$ is either a constant (for $\varepsilon = \pm 1$) or a periodic function of time with period $\tau/4$ (we recall $\tau = 2\pi/\omega$) corresponding to the time for the high-energy photon to cross half a wavelength of the background field. It is thus convenient to compute the probability $P_m$ of the high-energy photon to decay into a pair during this interval of time. Let us then consider a time $t_m$ at which the high-energy photon sees a maximum of the background electric field, and compute the probability $P_m$ for the high-energy photon to decay into an electron–positron pair in between the times $t_0 = t_m - \tau/8$ and $t_0 + \tau/4 = t_m + \tau/8$:

$$ P_m = P(t_0, \tau/4) = 1 - \exp\left(-R_m \frac{\tau}{4}\right) \quad \text{with} \quad R_m \frac{\tau}{4} = \frac{W_0}{2\omega\gamma\gamma} \int_0^\pi b_0(\chi_0 \Psi_\varepsilon(\varphi)) \, d\varphi, \quad (17) $$

where $R_m\tau/4$ is the time-integrate rate, while $R_m$ denotes the average rate.

Let us stress that the dimensionless quantity $W_0/(2\omega\gamma_\gamma)$ in equation (17) is of order 1 for optical background fields (with micrometric wavelength $\lambda$) and GeV-level high-energy photons:

$$ \frac{W_0}{2\omega\gamma_\gamma} = \frac{\alpha m_e c^2}{3 \hbar\omega} \frac{m_e c^2}{\hbar\omega_\gamma} \approx 0.512 \left( \frac{\lambda}{1\mu\text{m}} \right) \left( \frac{1 \text{ GeV}}{\hbar\omega_\gamma} \right). \quad (18) $$

The probability $P_m$, that measures the contribution of a single field maximum, depends on the integral quantity

$$ \mathcal{I}_{\varepsilon}(\chi_0) = \int_{0}^{\pi} b_0(\chi_0 \Psi_{\varepsilon}(\varphi)) \, d\varphi, \quad (19) $$
which, for a given polarization ellipticity $\epsilon$, is a function of $\chi_0$ only as shown in figure 1(b) for LP (black line) and CP (blue line).

This integral is the key object to compute the pair creation probability and can be calculated either numerically or by means of an approximated formula as described in the following. In appendix A, we generalize the approach proposed in [36] to arbitrary $\epsilon$ and $\chi$, and show that $\mathcal{I}_{\epsilon}(\chi_0)$ can be written in the approximate form:

$$ \mathcal{I}_{\epsilon}(\chi_0) \approx \pi b_0(\chi_m) \min\{F(s_{\epsilon}(\chi_m)), f(\epsilon)\} \quad \text{with} \quad \begin{cases} F(s) = \sqrt{2/\pi s} \operatorname{erf}(\pi\sqrt{2}/(4s)) \\ f(\epsilon) = \frac{1}{\pi} \int_0^{\pi} [\sin^2 \varphi + \epsilon^2 \cos^2 \varphi]^{1/3} d\varphi \end{cases}, \quad (20) $$

where $\chi_m$ is defined in equation (14) and $\operatorname{erf}(x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} dt$ is the error function. The function $F(s)$ emerges from the saddle point approximation used to compute the integral when the main contribution comes from the vicinity of the field maximum. It varies slowly with $s$, as $F(s) \approx \frac{2}{\sqrt{\pi}} s$ for $s \ll 1$ and $F(s) \to 1$ for $s \to +\infty$. In equation (20), $F(s)$ takes for argument $s_{\epsilon}(\chi_m)$, with

$$ s_{\epsilon}(\chi) = \sqrt{\frac{3}{2} \frac{c(\chi)}{\sqrt{1-\epsilon^2}}} \quad \text{and} \quad c(\chi) = \sqrt{\frac{2 b_0(\chi)}{3 \chi b'_0(\chi)}}, \quad (21) $$

where $b'_0(\chi)$ denotes the derivative of $b_0(\chi)$ and $c(\chi)$ is a slowly varying function of $\chi$, as $c(\chi) \approx \sqrt{\chi}/2$ for $\chi \ll 1$ and $c(\chi) \to 1$ for $\chi \to +\infty$.

Equation (20) is exact for CP background fields ($\epsilon = \pm 1$) for which $s_{\pm 1}(\chi_m) \to +\infty$, so that $\mathcal{I}_{\pm 1}(\chi_0) = \pi b_0(\chi_0/\sqrt{2})$. Moreover it allows to recover the asymptotic behavior of equation (19) in the limiting cases of small and large quantum parameter $\chi_0$:

$$ \mathcal{I}_{\epsilon}(\chi_0) \xrightarrow{\chi_0 \ll 1} \begin{cases} c_3(1-\epsilon^2)^{-1/2}(1+\epsilon^2)^{-1/4} \chi_0^{3/2} \exp\left(-\frac{8\sqrt{1+\epsilon^2}}{3\chi_0}\right) & \text{for } \chi_0 \ll 1-\epsilon^2, \\ c_4 \chi_0 \exp\left(-\frac{8\sqrt{2}}{3\chi_0}\right) & \text{otherwise,} \end{cases} \quad (22) $$

$$ \mathcal{I}_{\epsilon}(\chi_0) \xrightarrow{\chi_0 \gg 1} c_5 \frac{f(\epsilon)}{(1+\epsilon^2)^{1/3}} \chi_0^{2/3}, \quad (23) $$

where we used equations (5) and (6), and $c_3 = 27\sqrt{\pi}/64 \approx 0.529$, $c_4 = 9\sqrt{3\pi}/64 \approx 0.765$ and $c_5 = 3^{2/3}45\Gamma^4(2/3)/(56\pi) \approx 1.789$.

### 2.2.2. Total decay probability

As will be further discussed in section 2.3, considering parameters typical of the upcoming generation of laser facilities, the probability for a high-energy photon to decay into a pair after crossing a single half-wavelength of the background field will in general be small compared to 1. This will not however be the case when considering the cumulative effect of crossing several wavelengths.

Let us denote as $t_0 = 0$ the time at which the photon starts to interact with the background field$^6$. The probability for a high-energy photon to decay into a pair after $t = n\tau/4$ of interaction with the background field (or, equivalently, after crossing $n$ local maxima of the field amplitude) $P_{\text{tot}}(t = n\tau/4)$, can be derived from the number $N_{\gamma}$ of photons surviving after a time $t$:

$$ N_{\gamma}(t = n\tau/4) = N_0 \prod_{m=1}^{n} \bar{P}_m \quad \text{with} \quad \bar{P}_m = 1 - P_m = \exp(-R_m \tau/4). \quad (24) $$

Here we have used $m$ as a running index for compactness, and $P_m$ (given by equation (17)) and $\bar{P}_m$ denote the probabilities for the photon to decay and not to decay, respectively, during the $m$th interval. We can then write the total probability of producing a pair as

$$ P_{\text{tot}}(t = n\tau/4) = 1 - \prod_{m=1}^{n} \bar{P}_m = 1 - \exp(-Rt) \quad \text{with} \quad R = \frac{1}{n} \sum_{m=1}^{n} R_m, \quad (25) $$

$R$ denoting the average pair production rate.

<footnote>$^6$Throughout this paper, we consider background fields with a finite (e.g. $\sin^2$) temporal profile so that this time $t_0$ is uniquely defined.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Temporal evolution of the probability for a high-energy photon with normalized energy $\gamma_{\gamma} = \hbar\omega_{\gamma}/(m_e c^2) = 10^3$ to decay into an electron–positron pair while interacting with a plane wave background field in head-on collision. Considering a constant background field envelop with (a) $a_0 \simeq 82$ corresponding to $\chi_0 = 0.5$ and (b) $a_0 \simeq 660$ corresponding to $\chi_0 = 4$. Considering a background field with $\sin^2$ temporal profile with maximum field amplitude (c) $a_0 \simeq 82$ corresponding to $\chi_0 = 0.5$, (d) $a_0 \simeq 660$ corresponding to $\chi_0 = 4$. Solid lines correspond to the exact probability obtained by integrating numerically equation (9), dashed lines and dots to the theoretical prediction obtained from equation (25), for LP ($\epsilon = 0$, black) and CP ($\epsilon = \pm 1$, blue) cases. The red dash-dotted lines correspond to the theoretical prediction from equation (25) using the approximation of equation (20) for LP.</figcaption></figure>

Let us note that this formula is exact at $t = n\tau/4$ with $n \in \mathbb{N}$, and—as will be shown later—can also be used with good approximation at all times $t \gg \tau/4$. Moreover, in the case of a monochromatic plane wave (no temporal envelop), the average pair production rate simply reduces to $R = R_m$. When dealing with a pulse with a finite temporal envelop, however, equation (25) needs to be computed combining the contribution of successive $(\tau/4)$-long intervals, using the local maximum of the background field strength for each time interval. This approach gives very good results even considering ultra-short, few cycles laser pulses. To confirm this, figure 2 shows the temporal evolution of the decay probability for a high-energy photon with $\gamma_{\gamma} = 10^3$ colliding head-on with different background fields. In panels (a) and (b), the high-energy photon interacts with a background field with constant amplitude $a_0 \simeq 82$ ($\chi_0 = 0.5$, figure 2(a)) and $a_0 \simeq 660$ ($\chi_0 = 4$, figure 2(b)). The interaction lasts for a time $\tau$, i.e. the photon explores two wavelengths of the background field. In panels (c) and (d), the high-energy photon interacts with a background field with a $\sin^2$ temporal profile in intensity, a full-width-half-maximum (FWHM) of $5\tau$, and maximum field strength $a_0 \simeq 82$ and $a_0 \simeq 660$, respectively.

In all panels, the solid lines denote the *exact* probability computed by numerically integrating equation (9), considering either LP (black lines) or CP (blue lines) background field. We can see that, for $\chi_0 = 0.5$, CP produces less pairs than LP, while for $\chi_0 = 4$ CP is slightly more efficient (the difference between the two polarizations will be discussed in details in section 2.3.1). These probabilities are compared against the prediction of equation (25), computed using either the numerically evaluated $I_{\epsilon}(\chi_0)$ (dashed lines in panels (a) and (b), dots in panels (c) and (d)) or the approximation given by equation (20) to compute $R_m$ (red dot-dashed lines, only computed for LP as it is exact for CP). Note that here and in the following, whenever equation (20) is used, the improved Erber approximation equation (4) is exploited.

Panels (a) and (b) show an excellent agreement between the exact computations (solid lines) and equation (25) (dashed lines), both using the numerically integrated values of $I_{\epsilon}(\chi_0)$. A remarkably good agreement is also found when considering the fully analytical approximation of equation (20) for $I_{\epsilon}(\chi_0)$ (red dot-dashed lines). The probability is slightly overestimated in this case, as expected for $\chi_0 \sim 1$ for which equation (20) shows the greatest departure from the exact expression of $I_{\epsilon}(\chi_0)$.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Probability $P_m$ for a high-energy photon to decay into an electron–positron pair after crossing half-a-wavelength of a LP background field, as a function of the photon energy $\gamma_{\gamma} = \hbar\omega_{\gamma}/(m_e\epsilon^2)$ and background field amplitude $a_0$: (a) integrating numerically equation (17), (b) computed using the approximate expression equation (20), (c) extracted from one-dimensional PIC simulations. In all three panels, the solid white lines report the isocontours of the first panel. In panel (b), the isocontours obtained from equation (20) are shown in red dot-dashed lines. White dotted lines correspond to constant values of $\chi_0 = 0.1, 1, 16.5$ from bottom left to top right corners.</figcaption></figure>

Similarly, an excellent agreement between all three approaches is found in panels (c) and (d) considering a short ($5\tau$ FWHM in intensity) background field. In this case the value of $R_m$ in each $\tau/4$ interval is computed by taking the value of $I_\epsilon$ given by the local maximum of the field, i.e. the time variation of the envelope is considered slow during $\tau/4$. The validity of this assumption is confirmed by PIC simulations discussed in the following. In particular, figure 2 shows that the approach leading to the derivation of equation (25), which consists in treating the cumulative contributions of successive maxima of the background field, remains a very good approach even for ultra-short, few cycle, background fields.

Last, we note that the results of 1D PIC simulations performed with Smilei (see appendix B for details) are not distinguishable from the exact computations (solid lines) for all reported cases, and are therefore not shown in figure 2.

## 2.3. Guidelines for pair production optimizations

### 2.3.1. Influence of the background field polarization

Interesting insights into the role of the background field polarization on pair production can be gained from the asymptotic behavior provided by equations (22) and (23).

We recall that, motivated by experimental constraints, we compare LP and CP at fixed energy: this condition implies that the CP beam has a lower maximum amplitude. Because of this, for $\chi_0 \ll 1$, equation (22) shows that any departure from the LP case leads to a tremendous decrease of $I_\epsilon(\chi_0)$. Indeed the exponential cut-off of the pair production rate is strongly affected by the decrease of $\chi_m$ with $|\epsilon|$ as shown in equation (14). In contrast, for $\chi_0 \gg 1$, $f(\epsilon)/(1+\epsilon^2)^{1/3}$ is an increasing function of $|\epsilon|$: the reduction of the field peak amplitude is compensated by the fact that the absolute value of the amplitude is not time dependent (see equation (15)). Hence, in this range, increasing the background field ellipticity increases the pair production rate. However, the rate increase from the LP to the CP case is small (less than 12%), and it is thus of marginal use to consider CP for boosting pair production$^7$.

The differences between the LP and CP cases are clearly shown in the inset in figure 1(b), where $I_0(\chi_0)$ is found to be orders of magnitude larger than $I_{\pm 1}(\chi_0)$ for $\chi_0 \ll 1$ while $I_{\pm 1}(\chi_0)$ is only about 10% larger than $I_0(\chi_0)$ for $\chi_0 \gg 1$. Similar conclusions can be drawn from figure 2. For small values of the photon parameter, $\chi_0 = 0.5$ (left panels), the pair production probability is significantly larger considering a LP background field instead of a CP one. For the higher value of $\chi_0 = 4$ (right panels), CP increases, but only marginally, pair production. For this reason, we will focus on pair production in LP background fields, unless otherwise specified.

### 2.3.2. Influence of the background field amplitude and gamma-photon energy

The relative influence on pair production of the photon energy $\gamma_\gamma$ and the background field amplitude $a_0$ (or equivalently intensity $I_0$, as we consider here a given background field wavelength $\lambda = 0.8~\mu\text{m}$), is summarized in figure 3. It shows the probability $P_m$ for a high-energy photon to decay in a pair after crossing half a wavelength of a LP background field: (a) integrating numerically equation (17) which provides an exact value of $P_m$; (b) using the approximate but fully analytical expression given by

<footnote>$^7$ In particular as, in practice, implementing CP on multi-petawatt facilities would lead to a decrease of the delivered laser pulse energy.</footnote>
<figure><img src="image_7.png" /><figcaption>Figure 4. Position in the (γγ, a₀)-plane [equivalently in the (Eγ, I₀)-plane] at which RΔt = 1 [Ptot(Δt) ≈ 0.63] for the interaction with a single maximum of the incoming pulse (half-wavelength, i.e. pulse duration equal to τ/2) for which the interaction lasts for Δt = τ/4 (yellow line) and for the interaction with pulses with a step-like time profile (solid lines) of duration 5τ (green) and 100τ (blue), and pulses with a sin² time envelope with FWHM of 5τ (green dashed line) and 100τ (blue dashed line). The color-shaded areas above the plain lines correspond to RΔt > 1 (considering the step-like time profiles). The dotted black lines correspond to constant χ₀ = 0.1, 1, 4, 16.5. The symbols highlight the parameters of the simulations discussed in sections 3.4 and 3.5.</figcaption></figure>

equation (20); (c) extracted from 1D PIC simulations of the interaction of a flash of high-energy photons colliding head-on with a plane wave (details are given in appendix B).

Let us first discuss panel (a) and note that, in addition to the probability isocontours (solid white lines), the dotted white lines represent the contours of constant quantum parameter, χ₀ = 0.1, 1 and 16.5. As χ₀ ∝ γγa₀, these are straight lines with a −45° inclination.

In the limit χ₀ ≪ 1, i.e. the bottom-left of figure 3, Pm assumes very small values. In this limit, the time-integrated rate Rmτ/4 scales as γγ⁻¹χ₀³/² exp[−8/(3χ₀)] (equation (22)) and the exponential term, depending only on χ₀, gives the dominating contribution to Pm. As a result, the isocontours of Pm (solid white lines) behave as nearly straight lines, roughly parallel to the contours of constant χ₀. We wish to stress however that the dependence on γγ⁻¹ cannot be fully ignored: the isocontours of Pm are less steep than the isocontours of χ₀ which indicates that Pm increases faster with a₀ than with γγ. In this range of small χ₀, for which probability and time-integrated rate are equivalent, the importance of increasing a₀ to improve pair production was already pointed out in [36]. However, the probability variation with both a₀ and γγ needs to be examined in more details at higher χ₀.

Indeed, the dependence with γγ⁻¹ plays a more important role as χ₀ increases so that χ₀ cannot be considered as the only relevant parameter in order to optimize the pair production rate. This is clearly seen by considering the isocontours of Pm in (a₀, γγ)-plane, figure 3(a): the isocontour slope becomes shallow and eventually changes sign. Thus, for any value of Pm, a minimum field strength a₀ is needed in order to obtain the desired level of probability. A corollary is that, at constant a₀, increasing γγ increases the probability Pm up to a maximum value, beyond which a further increase of γγ would only decrease Pm. The minimum of each isocontour can thus be found by solving ∂Pm/∂γγ = 0, which one can recast in the form

$$ \frac{d}{\mathcal{I}_{\varepsilon}} d\chi_0 - \frac{\mathcal{I}_{\varepsilon}(\chi_0)}{\chi_0} = 0. \qquad (26) $$

Interestingly, this equation involves only χ₀ so that the minima of the probability isocontours lie on a straight line of constant χ₀. This a priori surprising result can be better understood noting that both a₀ and Pm are Lorentz invariants, the minimum value of a₀ to reach a given probability Pm can depend only on the Lorentz invariant involving the photon energy, i.e. χ₀. Solving numerically equation (26) for ε = 0 (LP case), we obtain<fnref n="8" /> χ₀ ≈ 16.5, which is reported by a dotted line in figure 3.

To conclude with figure 3, let us now turn to panels (b) and (c). Panel (b) reports the probability Pm computed from equation (20) and the corresponding isocontours (red dashed lines). Superimposing the isocontours (white lines) of the probability from panel (a), we can see an excellent agreement between the two approaches. This validates the approximate form, equation (20), which has the advantage to be

<footnote n="8">⁸ Solving equation (26) for arbitrary values of ε ∈ [−1, 1] leads to χ₀ in between 16.5 and 17.2, and corresponding values of Iε(χ₀) in between 5.3 and 6.0.</footnote>
completely analytical and can now be used to get quick yet precise estimates of the probability $P_m$. Last, panel (c) reports the probability extracted from 1D PIC simulations of the head-on collision of a flash of high-energy photons with a laser beam. The probability was extracted from the depletion of high-energy photons after crossing half a laser wavelength. Here again, a very good agreement is observed with the isocontours (white lines) obtained from the exact integration shown in panel (a). This last panel thus provides a cross-benchmark of our model and the Monte-Carlo module for nonlinear Breit-Wheeler pair production implemented in our PIC code SMILEI.

Our model can be used to predict the effectiveness of pair production for a given set of laser and high-energy photon parameters<fnref n="9" />. Most of the photons crossing the laser will be converted in pairs over a time $\Delta t$ if the quantity $R\Delta t$ in equation (25) becomes of order 1, for which<fnref n="10" /> $P_{\text{tot}}(\Delta t) \ge 0.63$. In figure 4, we highlight in yellow the region in the ($a_0, \gamma_\gamma$)-plane where the probability for a high-energy photon to decay into a pair after interacting with a single maximum (half-wavelength) of the laser pulse $P_{\text{tot}}(\tau/4)$ is equal or larger than 0.63 (the pulse duration is here $\tau/2$, but the time required for the photon to explore half the laser wavelength is $\Delta t = \tau/4$). The condition $P_{\text{tot}}(\Delta t) = 0.63$ is also shown for pulses with a step-like time profile (solid lines) and duration $5\tau$ (green, $\Delta t = 2.5\tau$) or $100\tau$ (blue, $\Delta t = 50\tau$), and pulses with a $\sin^2$ intensity time profile with FWHM $5\tau$ (green dashed line, $\Delta t = 5\tau$) and $100\tau$ (blue dashed line, $\Delta t = 100\tau$). Note that both step-like and $\sin^2$ time profiles correspond to laser pulses with the same energy and peak power, and lead to similar requirements to reach $P_{\text{tot}} \approx 0.63$. Comparing the blue and green lines, we see that at constant $\gamma_\gamma$ a longer pulse satisfies the condition of efficient conversion for a lower value of $a_0$. Hence, even though the condition $R\tau/4 = 1$ (yellow solid line) is quite stringent and achieving a high probability level over a time interval $\tau/4$ can be difficult, the condition $R\Delta t \gtrsim 1$ is significantly relaxed considering longer interaction time. Indeed, significant (order 1) pair production probability is expected on forthcoming multi-PW laser facilities (see section 4 for details).

Let us also note that, when $R\tau/4 \ge 1$ (yellow area), more than 63% of conversion is achieved after crossing a single wavelength of the background field, and more than 99% of the incident gamma photons are converted into pairs in less than 3 periods. This gives a very robust condition for abundant pair creation, and it determines a limit above which is not useful to further increase the laser amplitude to increase the number of primary pairs. As we will show in the following, this condition ($R\tau/4 \ge 1$) can also be invoked to find the optimal laser transverse shape and focal spot for a given laser energy. Indeed 3D PIC simulations, shown by the symbols in figure 4 and discussed in detail in the next section 3, show a different behaviour depending on whether their initial parameters lie in the yellow or green region of figure 4. In the case of the two simulations reported as squares in the figure, we will see that these differences are observed even though the simulation parameters correspond to the same value of $\chi_0$.

As a final remark, we point out that the probabilities considered so far are computed considering a given high-energy photon energy. In practice, we expect the high-energy photon beam to have a broad energy spectrum or at least a finite energy width (see e.g. references [25, 32, 39–46]). Accounting for such energy distributions is however straightforward as there is no coherence effects between the high-energy photons. The total probabilities can thus be computed independently for all relevant high-energy photon energies so that $P_{\text{tot}}(\gamma_\gamma)$ is actually a function of $\gamma_\gamma$. As an example, to compute the total number $N_{\text{pair}}$ of produced pairs, one only needs to integrate over the whole high-energy photon spectrum: $N_{\text{pair}} = \int d\gamma_\gamma P_{\text{tot}}(\gamma_\gamma) dN_\gamma/d\gamma_\gamma$, where $dN_\gamma/d\gamma_\gamma$ denotes the high-energy photon energy distribution.

### 3. Geometrical effects on pair production: the case of LG beams

This section aims to generalize the model developed in section 2.2.2, to overcome the plane wave approximation and account for the laser beam transverse profile. We start by reminding some properties of the LG beams, limiting ourselves to the critical aspects for the understanding of pair production. We then compare the predictions of our theoretical model (generalized to account for the laser transverse profile in section 3.3) with the results of 3D PIC simulations, whose initialization is discussed in section 3.2. A thorough comparison of the efficiency of pair production considering beams with the same energy and different field structures is presented in sections 3.4 and 3.5. We explore first the impact of the beam geometry and amplitude in a regime of efficient pair production (time-integrated rate $R_m\tau/4 > 1$) and then in a regime ($R_m\tau/4 < 1$) more relevant for current experimental facilities (as will be discussed in the following section 4). Finally in section 3.6, we discuss some properties of the produced pairs, such as the energy spectrum and the propagation direction.

<footnote n="9">⁹ In the case of large secondary pair production, our prediction will underestimate the number of pairs as our analysis is limited to the soft-shower regime.</footnote>

<footnote n="10">¹⁰ Note that in this limit, the probability cannot be assimilated to the time-integrated rate, as in the case of very weak pair production.</footnote>
<figure><img src="image_8.png" /><figcaption>Figure 5. Photon quantum parameter $\chi_{\gamma}$ in the head-on collision of photons of energy $\gamma_{\gamma} = 400$ with a LP LG beam at focus (i.e. $z = 0$ in equation (27)) with $a_0 = 2000$, $p = 0$ and (a) $\ell = 0$, (b) $\ell = 1$, (c) $\ell = 2$, (d) $\ell = 3$. (e) Maximum of the complex envelope $u_{p\ell}^{\text{max}}$ from equation (29) (black line) and for $p \in [0, 3]$ as a function of $|\ell|$. For $p \neq 0$, equation (28) is solved numerically.</figcaption></figure>

### 3.1. Properties of LG beams

The LG modes are solutions of the paraxial equation with cylindrical symmetry and a generalization of the Gaussian beam. Each mode is indexed by two integers: $p \ge 0$ for the radial order and $\ell$ for the azimuthal one. Hence, the notation $LG_{p\ell}$ will be used in the following to identify a specific mode.

Starting from equation (10), the electric field of an LG beam can be expressed using $E_0 = a_0 \Re(u_{p\ell})$, where the complex envelope $u_{p\ell}$ in cylindrical coordinates is [64]

$$u_{p\ell}(\rho, \phi, z) = C_{p\ell} \frac{w_0}{w(z)} \left(\frac{\sqrt{2}\rho}{w(z)}\right)^{|\ell|} L_p^{|\ell|} \left(\frac{2\rho^2}{w^2(z)}\right) \exp\left[-\frac{\rho^2}{w^2(z)}\right] \exp\left[-i\psi_{p\ell}(z) + i\ell\phi + i\frac{z\rho^2}{w^2(z)}\right], \quad (27)$$

where $C_{p\ell} = \sqrt{p!(p + |\ell|)!}$ is a normalization constant, $L_p^{|\ell|}(x)$ are the generalized Laguerre polynomials, $w(z) = w_0\sqrt{1 - z^2/z_R^2}$ is the beam waist with $w_0 = w(z = 0)$ the waist at focus and $z_R = \pi w_0^2/\lambda$ the Rayleigh length, and $\psi_{p\ell}(z) = (2p + |\ell| + 1) \arctan(z/z_R)$ is the generalized Gouy phase. Note that the LG modes described by equation (27) have all the same energy and that $LG_{00}$ corresponds to a Gaussian beam with maximum field amplitude $a_0$.

In order to infer how the structure of the LG modes might affect pair production we show in figures 5(a)–(d) the quantum photon parameter $\chi_{\gamma}$ for the head-on collision of photons of energy $\gamma_{\gamma} = 400$ with a LP LG beam at focus (i.e. $z = 0$) with $a_0 = 2000$, $p = 0$ and $\ell \in [0, 3]$. As can be seen in figures 5(a)–(d), for $\ell \neq 0$ the maximum of the field, and thus the peak values of the quantum parameter, are no longer along the beam axis (i.e. $x = y = 0$), in contrast with the Gaussian beam case. The radius at which the field has its maximum amplitude $\rho^{\text{max}}$ can be obtained by solving $\partial_{\rho}|u_{p\ell}|^2 = 0$, leading to:

$$\frac{2\rho^2}{w(z)^2} + (1 - \delta_{0p}) \frac{4\rho^2}{w(z)^2} \frac{L_{p-1}^{|\ell|+1} \left( \frac{2\rho^2}{w(z)^2} \right)}{L_p^{|\ell|} \left( \frac{2\rho^2}{w(z)^2} \right)} - |\ell| = 0, \quad (28)$$

with $\delta_{ij}$ the Kronecker delta. For $p = 0$, equation (28) can be solved analytically and gives $\rho^{\text{max}} = w(z)\sqrt{|\ell|/2}$. While the radial position of the field maximum for $p = 0$ increases with $|\ell|$, the amplitude decreases, as shown in figure 5(e), and tends asymptotically to zero for $|\ell| \gg 1$ as $|2\pi\ell|^{-1/4}$. The maximum amplitude of the envelope $u_{0\ell}$ can be computed inserting $\rho^{\text{max}}$ into equation (27), which gives

$$u_{0\ell}^{\text{max}} = |u_{0\ell}(\rho^{\text{max}}, \varphi, z)| = \frac{1}{\sqrt{1 + z^2/z_R^2}} \frac{|\ell|^{\ell/2} e^{-|\ell|/2}}{\sqrt{|\ell|!}}. \quad (29)$$

Figure 5(e) also shows $u_{p\ell}^{\text{max}}$ for $p \neq 0$, obtained by inserting the numerical solution of equation (28) into equation (27). The maximum field amplitude for the LG beams decreases with both $|\ell|$ and $p$, and as a consequence the same behaviour is observed for the maximum value of $\chi_{\gamma}$ (see figures 5(a)–(d)). Therefore, in the following we will consider only $|\ell| \le 5$ and $p = 0$, because of the higher field amplitude (hence higher values of $\chi_{\gamma}$) and because of the practical challenges in producing high order modes at relativistic intensity.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Typical set-up of the 3D PIC simulations, with $\lambda$ and $\tau$ the laser wavelength and the optical cycle, respectively.</figcaption></figure>

Note that in the considered cases with $\ell \neq 0$, $\chi_\gamma$ has $2|\ell|$ maxima in the transverse plane (see figures 5(b)–(d)) and, looking at the evolution in time, each maximum makes a full rotation around the laser propagation axis in $|\ell|\tau$. This means that, in a fixed transverse plane, a maximum rotates to the location of the next one in half a period.

If we consider the head-on collision of the LG beam with a gamma photon, the photon is crossing half a period of the field in $\tau/4$. Hence, the photon crosses a rotating local maximum in $\tau/4$. This observation is fundamental for generalization of the reduced model presented in the previous section 2.2.2, where we consider subsequent intervals of duration $\tau/4$ to infer the total probability.

As noted in figure 5, the LG beams not only have a non-trivial spatial distribution of the fields in the transverse plane, but also a maximum field amplitude decreasing with $|\ell|$ (following equation (29)). In order to distinguish the impact on pair production of these two aspects, we discuss the interaction with extended Gaussian (EG) beams. For each LG beam (i.e. each value of $|\ell|$), we define the associated EG beam as the Gaussian beam that has the same maximum field amplitude and carries the same total energy. As $u_{\ell\ell}^{\text{max}}$ decreases with $|\ell|$ (equation (29)) and the energy is kept constant, the waist of the EG beams increases with respect to the fundamental Gaussian (equivalent to $\text{LG}_{00}$) as $w_\ell = w_0 \sqrt{|\ell|!} e^{|\ell|^{1/2}} |\ell|^{-|\ell|/2} \ge w_0$.

## 3.2. Simulation set-up

In what follows, we present 3D PIC simulations performed with the open-source code SMILEI [54], considering the setup schematically represented in figure 6. Each simulation reproduces a volume of $40\lambda \times 40\lambda \times 14\lambda$ (in the x, y and z directions, respectively) with spatial resolution $\lambda/24$ and temporal resolution at 95% of the CFL condition [65].

An intense laser pulse (with wavelength $\lambda = 0.8~\mu\text{m}$) is injected from the left boundary $z = -5\lambda$ using the method presented in reference [66] that allows for the Maxwell-consistent injection of an electromagnetic wave with an arbitrary spatio-temporal profile. Here, this method is used to inject the laser pulse by prescribing its spatio-temporal profile in its focal plane [i.e. the (x, y)-plane at $z = 0$], as given by the real part of equation (27) at $z = 0$, multiplied by a temporal envelope so that the laser pulse has a $\sin^2$ temporal profile in intensity with FWHM $\tau_{\text{FWHM}} = 5\tau$, and total duration of $10\tau$. The laser beam collides head-on with a counter-propagating flash of high-energy photons with a finite longitudinal extension $L_\gamma = \lambda/6$, transverse extension equal to the simulation box, and a density $n_\gamma$ equal to the critical density $n_c = \epsilon_0 m_e \omega^2 / e^2$. The duration of the simulation is long enough for the two light beams to stop interacting, but short enough so that none of the photons or secondary particles escape from the simulation box. All results presented hereafter have been extracted at the end of the simulations.

In sections 3.4 and 3.5, we present two series of simulations: the first one with $a_0 = 2000$ and $\gamma_\gamma = 400$, the second with $a_0 = 400$ and $\gamma_\gamma = 2000$, both series corresponding to a reference quantum parameter $\chi_0 \simeq 6 \times 10^{-6} a_0 \gamma_\gamma \sim 4.85$ (for $\lambda = 0.8~\mu\text{m}$). For each series, we varied the transverse spatial profile of the high-intensity laser pulse in order to study the impact of the spatio-temporal field profile, comparing different LG beams (varying $\ell$, but keeping $p=0$) and the corresponding EG beams. Note however that within each series, $a_0$ is left unchanged and refers to the maximum field amplitude of the reference Gaussian beam (for which we consider $w_0 = 3\lambda$). As a result, within each set of simulations, the laser pulse energy is unchanged, equal to $\simeq 10$ kJ for $a_0 = 2000$ and $\simeq 410$ J for $a_0 = 400$. Conversely, the photon flash surface energy density $m_e c^2 \gamma_\gamma n_\gamma d$ is $\simeq 5$ mJ/($\lambda^2$ for the first series (with $\gamma_\gamma = 400$) and $\simeq 25$ mJ/($\lambda^2$ for the second one ($\gamma_\gamma = 2000$).
<figure><img src="image_10.png" /><figcaption>Figure 7. Probability map for a photon with energy $\gamma_{\gamma} = 400$ to produce a pair during the interaction with the LG$_{05}$ beam with FWHM of $5\tau$ and maximum amplitude of the reference Gaussian beam equal to $a_0 = 2000$, corresponding to $\chi_0 \simeq 4.85$. (a) Analytical prediction from equation (25). (b) Numerical results from the 3D PIC simulation.</figcaption></figure>

As highlighted in figure 4 (squares), the reference Gaussian case of the first series of simulations allows us to explore the physics in the regime of high pair production probability after a single $\tau/4$ interval (i.e. $R_m\tau/4 > 1$, yellow region). On the contrary, with the second series we investigate a regime where efficient pair production is achieved thanks to the cumulative effects of several wavelength (i.e. $R_m\tau/4 < 1$, but the reference Gaussian case is still above the green line in figure 4).

### 3.3. Reduced model for pair creation with arbitrary spatio-temporal laser beam profiles

In this section we generalize the model developed in section 2.2.2 to accurately describe pair production from complex, spatially structured laser beams.

To do so, we exploit the fact that we are here dealing with high-energy photons colliding head-on with the strong background field. Hence, the photon momenta are left unchanged as the photon crosses the background field, and their trajectories are straight lines so that the field seen by a given photon along its trajectory is known and depends only on the photon initial position $(x_0, y_0)$ in the plane transverse to the direction of propagation of both light beams$^{11}$. It follows that the quantum parameter of any photon is known throughout the photon interaction with the background field, and for a given background field, is determined uniquely by the photon initial position $(x_0, y_0)$. As a result, the total probability $P_{tot}$ for any photon to decay into a pair after interacting with the laser pulse is also determined by the photon initial position, and the total number of produced pairs at the end of the interaction between the flash of high-energy photons (with density $n_{\gamma}$ and longitudinal width $L_{\gamma}$) simply reads

$$N_{\text{pair}} = n_{\gamma}L_{\gamma}\sigma_{\text{tot}} \quad \text{with} \quad \sigma_{\text{tot}} = \int dx_0 dy_0 P_{\text{tot}}(x_0, y_0). \qquad (30)$$

A key quantity to model pair production is thus the total cross-section $\sigma_{tot}$, here defined as the integral of the total probability over the transverse plane.

Notice that, as discussed in the last paragraph of section 2, for a given gamma photon spectrum $\sigma_{tot}$ can be considered a function of $\gamma_{\gamma}$, and the total number of produced pairs can be computed integrating over the incident high-energy photon spectrum. Moreover, equation (30) can be easily generalised to account for a transverse profile of the incident high-energy photon beam. In particular, in the limit where this beam is much narrower than the laser beam, equation (30) reduces to $\sigma_{tot} \simeq \sigma_{\gamma}P_{tot}$, where $\sigma_{\gamma}$ is the transverse area of the high-energy photon beam, and $P_{tot}$ is computed at the (transverse) location of the high-energy photon beam.

To compute the total probability at any position $(x_0, y_0)$ we use equation (25) and proceed as described in section 2.2.2, reconstructing the total probability from successive time intervals of duration $\tau/4$ (as discussed in section 3.1, a photon at any position in the transverse plane explores in $\tau/4$ a local maximum of the field amplitude whatever value $\ell$ assumes). One then only needs to compute the maximum field amplitude satisfactorily. For simplicity, we do so measuring the field amplitude at focus, taking the absolute value of the field envelope from equation (27) at $z = 0$. This approach, which neglects the effects of diffraction, is justified whenever the duration of the interaction ($\sim \tau_{FWHM}$) does not exceed the time $z_R/c$

<footnote n="11">11 For analytical tractability, we consider all photons to be propagating exactly in the z-direction.</footnote>
<figure><img src="image_11.png" /><figcaption>**Figure 8.** (a) Number of produced pairs, and corresponding cross section $\sigma_{tot}/\lambda^2$ after the interaction of gamma photons with energy $\gamma_{\gamma} = 400$ and a counter-propagating LG with $p = 0$, $\ell \in [[-1,5]]$ and $a_0 = 2000$ [CP (blue), LP (black)], and the corresponding LP EG (green). Frozen simulations (dots), full simulations (squares) and theoretical predictions from equation (30) (dashed lines). The insets show the probability map for the $LG_{00}$, $LG_{05}$ and $EG_{05}$ cases. (b) Relative difference in the number of produced pairs between the frozen and full simulation results.</figcaption></figure>

for light to cross the laser pulse Rayleigh length. As shown in the following, this approximation proves satisfactory considering typical ultra-short (Ti:sapphire) light pulses.

The result of this procedure is shown in figure 7(a) for the case of a LG beam with $\ell = 5$, $p = 0$ and $a_0 = 2000$ colliding with a flash of photons with energy $\gamma_{\gamma} = 400$ (corresponding to $\chi_0 = 4.85$; all other parameters are specified in section 3.2). Panel (a) shows a map, in the transverse plane, of the probability $P_{tot}$ computed over the whole interaction time in the transverse plane. Integrating over the transverse plane provides the measure of the total cross-section $\sigma_{tot} \approx 91.4 \lambda^2$.

To test the impact of the approximations used in the model (phase and diffraction effects are neglected), we show in figure 7(b) the probability map extracted from a 3D PIC simulation, where all effects neglected in our model are taken into account. Here the probability is extracted at the end of the interaction as the ratio of the remaining photon linear density by the initial one. An excellent agreement is found with the theoretical prediction of panel (b), with a measured total cross-section $\sigma_{tot} \approx 91.0 \lambda^2$, which confirms the predictive capability of our simple model for the case of spatially structured beams.

### 3.4. Simulations in the regime of high pair production probability ($R_m\tau/4 > 1$)

In order to investigate the dependence on the laser structure of the number of produced pairs, we present here a study in the regime satisfying $R_m\tau/4 > 1$ for the reference Gaussian beam, in which we expect efficient pair production. We remind that we consider here photons with normalized energy $\gamma_{\gamma} = 400$ and a reference Gaussian beam with maximum amplitude $a_0 = 2000$, leading to $\chi_0 \approx 4.85$ (all other parameters have been specified in section 3.2).

To test our analytical predictions and the assumptions at the base of our model, we performed simulations (hereafter referred to as frozen cases) in which the produced pairs are prevented to further radiate. In this condition, no secondary pairs are produced and the simulations allow to directly test the proposed analytical model. We also performed simulations including the full system dynamics (hereafter referred to as full cases), meaning that the produced pairs interact with the laser fields and potentially radiate. This allowed us to verify that we are still in the soft-shower regime, and that only a small departure from predicted number of produced pairs is indeed observed.

Figure 8(a) shows the number of produced pairs at the end of the interaction for all the tested cases, LP LG (in black), CP LG (blue) and EG (green), dots refer to the frozen simulations and squares to the full ones. Considering a LP LG laser beam, pair production efficiency increases with $\ell$. This means that for the high value of $a_0$ considered here, it is preferable to increase $\ell$, even if it is decreasing the maximum amplitude (as shown in figure 5(e)), as it increases the characteristic transverse size of the beam, and so the
<figure><img src="image_12.png" /><figcaption>**Figure 9.** (a) Number of produced pairs and corresponding cross section $\sigma_{tot}/\lambda^2$, after the interaction of gamma photons of energy $\gamma_{\gamma} = 2000$ with a counter-propagating LG beam with $p = 0$, $\ell \in [-1, 5]$ and $a_0 = 400$ [CP (blue) and LP (black)], and the corresponding LP EG (green). Full simulations (squares) and theoretical predictions from equation (30) (dashed lines). The insets show the probability map for the $LG_{00}$, $LG_{05}$ and $EG_{05}$ cases. (b) Relative difference in the number of produced pairs between the frozen and full simulation results. The frozen cases are not shown in panel (a) as the difference with the full ones is always below 15%.</figcaption></figure>

interaction area where the probability is high (see inset in figure 8(a)). The EG beams (which have by construction the same maximum amplitude of the LG beams but Gaussian profile) are more efficient than the corresponding LG, up to $\ell = 4$, as they have a bigger region with substantial probability in their probability maps. However, a slightly greater number of pairs are produced in the simulation with $LG_{05}$ than $EG_{05}$ beam. This because the maximum field amplitude has dropped substantially and a large part of the interaction region in the EG case has low probability of pair production (see inset in figure 8(a)). Predictions form our model discussed in section 3.3 and based on equation (30), shown in dashed lines, are in very good agreement with both EG and LG simulations.

Based on the prediction from figure 1, for the considered $\chi_0 \simeq 4.85$ we would expect at first look a higher efficiency in the CP case than with LP, contrary to what is obtained from PIC simulations (blue symbols in figure 8(a)). However, given the spatial structure of the laser field in the transverse plane, in a large part of the interaction region the quantum number is $\chi < 2.5$, i.e. the value above which CP should be favoured based on the analysis shown in figure 1(b) obtained within the plane wave approximation. Hence, the slightly higher number of pairs produced in the CP case with respect to LP in the region where $\chi > 2.5$, cannot compensate the higher efficiency of LP in the region of $\chi < 2.5$. Note that the results using a CP beam are independent from the sign of $\ell$ (i.e. from the total angular momentum). This suggests that the total angular momentum of the laser has no effect on the number of pair produced. Note however that we do not account for the spin polarization effects. Since we are focusing on the soft shower regime where there are few secondary pairs, the dynamics of the pairs affected by spin effects is expected to play a minor role for pair production in our case. Moreover the dynamics of the produced pairs is not the main concern of this study. Nevertheless, in other configurations spin effects can impact the created pairs and gamma photons dynamics as shown in references [67, 68].

For completeness, we compare the results of the frozen cases with simulations reproducing the full dynamics of the system (squares in figure 8(a)). As highlighted in figure 8(b), the difference in the produced number of pairs is always below 25% of its value in the frozen case, and it decreases with $|\ell|$. This confirms that the majority of pairs comes from the conversion of primary gamma photons and not from subsequent radiation, as expected for a soft-shower. In this regime, our model, which accurately captures the physics of pair production for the frozen simulations (dashed lines), can be therefore used to predict pair production in realistic conditions.
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Positron distribution in Lorentz factor $\gamma$ (equiv. normalized energy) and angle $\theta = \arctan(p_x/p_z)$ for the first simulation series (left column) and for the second one (right column); LP (CP) in the top (bottom) row. The insets show the energy spectra for positrons propagating in the positive (black line) and in the negative (red line) directions.</figcaption></figure>

### 3.5. Simulations in the regime of moderate pair production probability ($R_m\tau/4 < 1$)

In this second series of simulations we explore the regime of $R_m\tau/4 < 1$, maintaining the same maximum $\chi_0$ as in the previous section by exchanging the values of field amplitude (now equal to $a_0 = 400$) and photons energy (now equal to $\gamma_\gamma = 2000$). Note that considering a laser duration ($\tau_{FWHM} = 5\tau$), we still expect efficient pair production (the reference plane wave simulation lies above the green dashed line in figure 4). This can be seen in the inset of figure 9(a) for the Gaussian case $LG_{00}$ where a significant region exhibits a probability of order 1.

As figure 9 shows, for the simulations with LG beams the number of produced pairs is optimized for $\ell = 3$. The increase for $|\ell| \le 3$ can be explained as in the previous section by geometrical considerations, i.e. the area with the field being large enough to reach $P_{\text{tot}}$ of order one is getting wider with $|\ell|$. In the contrary for $|\ell| > 3$, the laser field amplitude becomes too low and having a laser beam with a larger transverse extension does not improve pair production. This is due to the fact that the field is becoming too low, so that the effect of decreasing $P_{\text{tot}}$ on pair production is more important, and is not compensated anymore by the increase in transverse size.

Moreover, in this regime, all LG beams perform better than the corresponding EG, given that the region with high (even though smaller than in the previous section 3.4) probability is larger for the LG than EG, as shown in the inset of figure 9 for $\ell = 5$. Also, a similar behaviour as discussed in the previous section 3.4 is observed for the CP simulations and can be explained with the same reasoning. Last but not least, our model is reproducing with high accuracy the observed number of pair in this case too.

### 3.6. Energy and angular distribution of the produced pairs

The energy-angular distribution of the produced pairs in the $z-x$ plane (i.e. the plane formed by the laser propagation direction $\hat{z}$ and the polarization direction $\hat{x}$ for the LP case) recorded at the end of the interaction is shown in figure 10 considering the reference LP and CP Gaussian beams. The top row gives the positron distribution in $\theta = \arctan(p_x/p_z)$ and in energy for the LP beams and the bottom one for the CP beams. The insets show the energy spectra of particles moving in the positive (same as the laser, black line) and negative (same as the gamma flash, red line) z-directions.

At high $a_0$ (panels (a) and (c)) the produced pairs, initially generated with momentum in the negative z-direction, are predominantly moving in the same direction as the laser (corresponding to $\theta = 0^\circ$) within a
cone of aperture 40°, meaning that they have been slowed down and then pushed by the laser. They reach a maximum energy of ≈ 10⁴m_ec² for the LP case and ≈ 8 × 10³m_ec² for the CP beam, consistent with the decrease of the peak field amplitude in CP with respect to LP at constant energy.

The spatial distribution is drastically modified when considering $a_0 = 400$. Even if there still is a substantial amount of pairs propagating in the laser direction, the majority of them keep their original direction, along the initial gamma photons propagation axis and anti-parallel to the laser pulse one. For these pairs the maximum energy is independent from the polarization, as shown in figures 10(b)–(d), contrarily to the pairs that are slowed down and pushed by the laser (right side of the quadrants or black line in the insets).

This study, supported by complementary simulations (not shown), suggests that for $a_0 \gtrsim \gamma_\gamma$ most of the created particles will end up propagating along the laser direction, while for smaller values of the laser maximum amplitude a large fraction of pairs will keep their original direction. A complete characterisation of the pairs spectrum and directionality is beyond the scope of this work, and will be investigated in more details in the future, as it can give important information for the planning of upcoming experimental campaigns.

## 4. Discussion

The model developed in this work allows to predict the pair production capability of upcoming facilities. In figure 11, we discuss the pair production maximum probability and total cross-section as a function of the photon energy and laser amplitude. The considered parameter range is relevant to ultra-short Ti:sapphire ($\lambda = 0.8$ µm) lasers, where studies of high-field physics and QED are already planned [69–72]. Two types of facilities can be distinguished. First, the Apollon<fnref n="12" />, CoReLS<fnref n="13" /> and ZEUS<fnref n="14" /> facilities (shown in black in figure 11) are designed to deliver multiple light beams with duration in between 15 and 30 fs FWHM and peak power from 1 to 10 PW. Second, the FACET-II<fnref n="15" /> and LUXE<fnref n="16" /> experiments will couple conventional electron accelerators with high-intensity (100 TW-class) ultra-short laser pulses (in white in figure 11).

In figures 11(a) and (b) the accessible $a_0$ and typical $\gamma_\gamma$ envisioned with these facilities (following the reference articles) are reported. For Apollon only, we computed $a_0$ considering a 20 fs light pulse delivering 20, 60 and 150 J on target, corresponding to peak power of 1, 3 and 7.5 PW, respectively. The reported value of $a_0$ is then computed for a LP Gaussian beam with focusing aperture $f/3$ (where the beam waist for a focusing aperture $f/N$ is $w_0 \approx 0.90 N\lambda \approx 2.2$ µm), using

$$a_0 \approx \frac{151}{N} \sqrt{\frac{\mathcal{E}_{\text{laser}}}{10 \text{ J}}} \sqrt{\frac{25 \text{ fs}}{\tau_{\text{FWHM}}}} \approx \frac{239}{N} \sqrt{\frac{P_{\text{laser}}}{1 \text{ PW}}}. \quad (31)$$

Note in addition that, as most of these facilities will operate in the $\chi_0 \gtrsim 1$ regime, we report for the photon energy that of the electron beams either expected from laser-driven wakefield acceleration (for Apollon, CoReLS and ZEUS) or emerging from the linear accelerators (for FACET-II and LUXE). In the case of a broad (e.g. synchrotron-like) energy distribution for the high-energy photon source, the expected number of produced pairs can be computed integrating the reported values of $P_{tot}$ or $\sigma_{tot}$ over the full high-energy photon spectrum. As detailed below, this integration can be carried out in two rough but straightforward ways depending on the kind of facility (multi-PW or laser combined with a conventional accelerator) that is considered.

Let us now focus on panel (a), which shows the maximum total decay probability (from equation (25)) for a photon with energy $\gamma_\gamma$ colliding head-on with a LP Gaussian laser pulse with maximum field strength $a_0$ and duration $\tau_{FWHM} = 25$ fs (typical for the facilities listed above). This panel suggests the strategy to optimize pair production depending on the type of facility.

First, despite the extremely high electron beam energies, facilities such as LUXE and FACET-II (in white in figure 11) offer the possibility to probe the regime of moderate quantum parameter ($\chi_0 \lesssim 1$) for which pair production may be observed ($P_{tot} \approx 10^{-2}$) but will not be abundant. This limit might be overcome by increasing the laser power up to few 100's of TW (300 TW are envisioned at LUXE), which allows to enter the quantum regime ($\chi_0 > 1$) and increase $P_{tot}$ to values exceeding 0.1. On these facilities however, because

<footnote n="12">¹² The Apollon, located on the Plateau de Saclay, south of Paris, France, started operating at the PW-level [2, 73]; see also the Apollon website.</footnote>

<footnote n="13">¹³ CoReLS, in Gwanju, South Korea, is operating at up to 4PW [3]; see also the CoReLS website.</footnote>

<footnote n="14">¹⁴ A 3PW laser facility being built at the University of Michigan; see the ZEUS website.</footnote>

<footnote n="15">¹⁵ FACET-II is planned at SLAC National Accelerator Facility [74]; see also the FACET-II website.</footnote>

<footnote n="16">¹⁶ LUXE is an experiment proposed at the European XFEL [75]; see also the LUXE experiment website.</footnote>
<figure><img src="image_14.png" /></figure>

$P_{\text{tot}}$ will not assume large values, increasing $a_0$ (e.g. by focusing the laser pulse as much as technically possible) is one of the most promising path to achieve abundant pair production.

In contrast, provided that multi-GeV electron beams can be obtained from laser wakefield acceleration, abundant pair production ($P_{\text{tot}} \gtrsim 0.63$) is expected on all multi-PW laser facilities with the focalisation technique already in place (typically an aperture in between f/3 and f/4 was considered in figure 11).

Similar conclusions can be drawn from panel (b) where we examine in more details the influence of the laser pulse spatial profile, by looking at the total cross-section<fnref n="17" /> $\sigma_{\text{tot}}$ (as defined by equation (30)) normalized to $w_0^2$, considering a Gaussian laser beam with the same temporal properties as in panel (a).

For the FACET-II and LUXE facilities, the total cross-section assumes very small values but increases fast as $a_0$ is increased: e.g. for LUXE, $\sigma_{\text{tot}}$ increases by 3 orders of magnitude increasing $a_0$ by a factor 10. Clearly, operating at large field strength will be a bottleneck for achieving abundant pair production on these facilities.

In contrast, in the range of parameters covered by multi-PW facilities, the dependence of $\sigma_{\text{tot}}$ with $a_0$ is much weaker. As in addition, $\sigma_{\text{tot}} \simeq w_0^2$ (consistent with $P_{\text{tot}} \simeq 1$ at the center of the beam), it becomes interesting for these facilities to increase the laser transverse size rather than opt for tight focusing.

Let us then discuss in more details the impact of focusing for the three upgrades of the Apollon facility (at 1, 3 and 7.5 PW). In panel (c), we present the total cross section (in units of $\lambda^{-2}$) as a function of $\gamma_\gamma$, considering either the standard f/3 aperture (black lines) or the more challenging f/1.5 aperture (red lines). For each laser power, a tighter focusing reduces the beam waist but increases the maximum $a_0$.

In all cases, the cross section is rapidly increasing with the photon energy for $\gamma_\gamma \lesssim 2 \times 10^3$ ($\simeq 1$ GeV). It is almost flat as $\gamma_\gamma$ reaches $5 \times 10^3$ ($\simeq 2.5$ GeV), which suggests that it is not worth increasing the photon energy above this value in forthcoming experiments to optimize pair production in the soft-shower regime.

<footnote n="17">¹⁷ Since the Gaussian field profiles depend on the transverse coordinate through the ratio $\rho/w_0$ (equation (27)), $\sigma_{\text{tot}}/w_0^2$ is independent of $w_0$.</footnote>
As shown in panel (c), for a given laser power the cross section for $f/1.5$ is larger than for $f/3$ only for small values of $\gamma_\gamma$, while the opposite behaviour is observed at large $\gamma_\gamma$ (of the order of a few GeV). This means that such tight focus, which is technologically very challenging, is not necessary for pair creation in forthcoming experiments on Apollon.

To test this prediction, we have performed complementary 3D PIC simulations of the 3 PW case with $f/3$ and $f/1.5$ ($a_0 = 138.8$ and 276.7, respectively) interacting with a gamma flash with energy $\gamma_\gamma = 10^4$. These simulations give $\sigma_{tot} \simeq 26.8 \lambda^2$ for $f/3$, which is indeed larger than the value $\sigma_{tot} \simeq 19.5 \lambda^2$ obtained for the tight-focusing $f/1.5$ case. These simulation results are consistent with our reduced model predictions, but also noticeably above its predictions, which are $\sigma_{tot} \simeq 16.2 \lambda^2$ for $f/3$ and $\sigma_{tot} \simeq 6.5 \lambda^2$ for $f/1.5$ (as given by figure 11(c)).

One of the reasons of the discrepancy is that our calculations are performed at focus, neglecting diffraction. While this was a very good approximation for the parameters considered earlier (see section 3.3), it is not correct when considering tightly focused intense background fields and an interaction region larger than the Rayleigh length. The model can be easily generalized to include diffraction. Doing so improves our predictions to $\sigma_{tot} \simeq 22.2 \lambda^2$ for $f/3$ and $\sigma_{tot} \simeq 9.8 \lambda^2$ for $f/1.5$. These values are now significantly closer to our 3D PIC simulation results, and are also in excellent agreement (within 1%) with additional 3D PIC simulations in which secondary pair production is not accounted (frozen case). Secondary pairs, responsible for a further increase of about 20% for the $f/3$ case and by a factor $\times 2$ for the tightly-focused $f/1.5$ case, are not included in the model and are the source of the remaining discrepancy. This suggests that for the highest peak intensity considered here we are at the limit of the soft-shower regime and the model should be considered as providing a lower bound estimate for pair production.

Let us finally briefly comment on how the effect of a broad energy spectrum and transverse extension of the high-energy photon beam can be accounted for depending on the type of facility considered.

In the case of FACET-II and LUXE facilities, both $P_{tot}$ and $\sigma_{tot}$ are very steep functions of $\gamma_\gamma$, e.g. increasing by several (more than 4) orders of magnitude for the FACET-II configuration as $\gamma_\gamma$ is increased from $5 \times 10^3$ to $3 \times 10^4$. One can thus expect that, for these facilities, only the photons with the highest energy contribute to pair production. The total number of produced pairs can then be estimated considering only those photons with the highest energy (maximum $\gamma_\gamma$ as reported on panels (a) and (b)). If the high-energy photon beam is much narrower than the laser beam, the total number of produced pair can be computed from the number of photons and $P_{tot}$ evaluated at this maximum $\gamma_\gamma$. If instead the high-energy photon beam has a large transverse distribution, the final number of pairs can be calculated by taking the high-energy photon surface density and $\sigma_{tot}$ computed at the maximum $\gamma_\gamma$.

For the multi-PW laser facilities, and large background field strength $a_0$, we have seen that both $P_{tot}$ and $\sigma_{tot}$ (for the latest, see e.g. panel (c) for the Apollon facility) reach a plateau above a minimum value of the photon energy (e.g. $\gamma_\gamma \simeq 5 \times 10^3$ in the case of panel (c) discussed above). Hence, pair production will follow primarily from the decay of the high-energy photons in this plateau. A simple estimate of the number of produced pairs will then be obtained from either the number of all photons in the plateau region and a typical value of $P_{tot}$ (if the photon beam transverse size is small with respect to the laser beam waist) or from their surface density and $\sigma_{tot}$ obtained for this energy range (if the photon beam transverse size is large).

## 5. Conclusion

In summary, we have presented a systematic study of nonlinear Breit-Wheeler pair production in the head-on collision of high-energy gamma photons with an ultra-high intensity laser beam. Combining analytical modeling and PIC simulations embarking the relevant QED modules, we have evaluated the impact on the efficiency of pair production of the laser spatio-temporal profile (comparing e.g. Gaussian with different focal spots and LG pulses) and parameters such as polarization, intensity and duration. Motivated by experimental constraints, we have considered fixed laser energies while changing these parameters.

We have explored laser field strength and photon energy parameter ranges relevant for currently and/or upcoming high-power laser facilities, focusing in particular on ultra-short Ti:sapphire laser facilities. A reduced model was proposed that allows to describe pair production in the regime of soft-shower, where secondary pair production can be neglected. This model, that would in principle allow to predict a minimum number of produced pairs, is found to agree remarkably well with 3D PIC simulations over a broad range of parameters, highlighting the importance of the soft-shower regime for experiments on the forthcoming laser facilities.

The model also allows to distinguish two regimes of interaction depending on whether the probability for a photon to decay into a pair as it crosses a single half-wavelength of the laser pulse is large (of order 1)
or not. The two regimes are not however fully determined by the photon maximum quantum parameter $\chi_0$, so that the laser field strength parameter $a_0$ and photon energy $\gamma_\gamma$ do not play a symmetric role in the system dynamics. This was confirmed in 3D PIC simulations where the two regimes were investigated at $\chi_0 = 4.85$ considering laser beams with complex spatio-temporal profiles, LG beams in particular, and different polarizations. It was found that, for a fixed laser energy, using circular polarization or LG does not improve significantly pair production. As they are also quite difficult to achieve experimentally, their effect on pair production is too marginal to be interesting for applications. It was found however that, in the regime of high pair production probability, it can be preferable not to focus much the laser beam (keeping the same total energy with lower peak intensity) and maintain a Gaussian shape in order to maximise the area with high-enough fields for efficient pair production.

These findings help to draw guidelines for future experiments on ultra-high intensity facilities. We show in particular that the path to abundant pair production is different whether one considers 100 TW-class or multi-PW laser systems. In particular, for facilities such as FACET-II or LUXE, 10 to 100 TW-class lasers will be coupled to ultra-relativistic 10–20 GeV electron beams emerging from conventional electron accelerators. On these facilities, pair production will develop in the low probability regime and operating in a tightly-focused configuration to access the highest possible laser field strength will be a major experimental challenge. In contrast, using multiple laser beam multi-PW facilities will operate in the regime of high pair production probability. In this regime, the pair production efficiency (measured in our work in terms of a total cross section) depends more strongly on the laser pulse transverse size than on its maximal field strength so that operating with standard (not too tight) focusing aperture increases the number of produced pairs. Interestingly, producing high-energy photons (or electrons, e.g. through laser wakefield acceleration) of a few (~ 5) GeV shall be enough for these facilities as higher energies do no help increase significantly the pair production efficiency in this high probability regime.

To conclude, this work provides a deeper understanding of the optimal conditions for pair production in upcoming experimental campaigns exploring nonlinear Breit–Wheeler pair production. It focused on pair production seeded by high-energy photons, and it will be interesting for future works to investigate the differences with the case where ultra-relativistic electrons are used to seed the pair production process. Studying the transition from the soft-shower to the cascade/avalanche regime in which secondary pairs and radiations play a dominant role is another interesting perspective.

## Acknowledgments

The authors are grateful to F Amiranoff, P Audebert, S Corde, L Gremillet, L Lancia, R Nuter and C Thaury for stimulating discussions. This work used the open-source PIC code SMILEI, the authors are grateful to all SMILEI contributors and to the SMILEI-dev team for its support. Simulations were performed on the Irene-Joliot-Curie machine hosted at TGCC, France, using High Performance Computing resources from GENCI-TGCC (Grant No. 2020-x2016057678). This publication is also supported by the Collaborative Research Centre 1225 funded by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 273811115—SFB 1225.

## Data availability statement

All data that support the findings of this study are included within the article (and any supplementary files).

## Appendix A. Approximation of $\mathcal{I}_\epsilon(\chi_0)$

In this appendix, we detail how the integral $\mathcal{I}_\epsilon(\chi_0) = \int_0^\pi b_0(\chi_0 \Psi_\epsilon(\varphi)) d\varphi$ (equation (19)) can be efficiently approximated by equation (20). To do so, let us start by stressing that, in general, and in particular at small $\chi$, $b_0(\chi)$ is a steep function of $\chi$. One thus expects the principal contribution to the integral to come from the phase $\varphi$ around $\varphi_m = \pi/2$, for which $\chi_0 \Psi_\epsilon(\varphi) \approx \chi_m$ with $\chi_m = \chi_0/\sqrt{1 + \epsilon^2}$. To compute the contribution of this maximum analytically, we approximate the integrand in equation (19) around $\varphi_m$ as

$$b_0(\chi_0 \Psi_\epsilon(\varphi)) \approx b_0(\chi_m) \exp \left( -\frac{(\varphi - \varphi_m)^2}{2s_\epsilon^2(\chi_m)} \right), \quad (A.1)$$

where $s_\epsilon(\chi_m)$ is chosen so that the exact and approximated integrand have the same second derivative at $\varphi = \varphi_m$, which gives

$$s_{\epsilon}(\chi_{m}) = \sqrt{\frac{3}{2} \frac{c(\chi_{m})}{\sqrt{1 - \epsilon^{2}}}} \quad \text{with } c(\chi) = \sqrt{\frac{2b_{0}(\chi)}{3\chi b'_{0}(\chi)}}. \quad (A.2)$$
With the ansatz equation (A.1), the integral can be performed analytically, leading to

$$ \mathcal{I}_{\epsilon}(\chi_0) \approx \pi b_0(\chi_m) F(s_{\epsilon}(\chi_m)) \quad \text{with} \quad F(s) = \sqrt{2/\pi s} \operatorname{erf}(\pi \sqrt{2}/(4s)). \qquad (A.3) $$

The ansatz equation (A.1) gives a good approximation of the integrand for arbitrary values of $\epsilon$ as long as $\chi_m$ is small enough. It is also exact (for all $\chi_0$) in the cases $\epsilon = \pm 1$ for which $s_{\pm 1}(\chi_m) \to +\infty$. Hence, equation (A.3) provides a very good, fully analytical approximation of equation (19) for a wide range of $\epsilon$ and $\chi_0$.

The approximation however needs to be corrected for $\chi_0 \gg 1$ and $|\epsilon| < 1$. In this limit, the integrand in equation (19) can be approximated using equation (6) leading to

$$ \mathcal{I}_{\epsilon}(\chi_0) \approx \pi b_0(\chi_m) f(\epsilon) \quad \text{with} \quad f(\epsilon) = \frac{1}{\pi} \int_{0}^{\pi} [\sin^2 \varphi + \epsilon^2 \cos^2 \varphi]^{1/3} d\varphi. \qquad (A.4) $$

Equations (A.3) and (A.4) can be combined in equation (20). This form ensures the correct asymptotic behaviour of the integral equation (19) in both small and large $\chi_0$ limits. It departs from the exact expression (integrated numerically) by less than 20%, the error being maximum for intermediate values of $\chi_0 \sim 1$. This was tested over a broad range of $\chi_0 \in [0.1, 1000]$ for various polarization parameter $\epsilon \in \{0, 0.25, 0.5, 0.75, 1\}$.

## Appendix B. Details on 1D PIC simulations

A series of 1D3V (1 dimensional in space and 3 dimensional in velocity) PIC simulations, considering a plane wave colliding head-on with a flash of gamma photons, have been performed to produce figure 3(c). We simulated a box of length $2.5 \lambda$ with spatial resolution $\lambda/256$, for a simulation time of $1.75\tau$, with temporal resolution $\tau/512$. This is a long enough time to study the full interaction of the plane wave with a flash of gamma photons of extension $\lambda/2$ and extract the probability from the number of surviving photons. The probability map in figure 3(c) is obtained from the results of 4096 ($64 \times 64$) simulations performed over a the logarithmically spaced range of $a_0 \in [10, 4000]$ and $\gamma_\gamma \in [100, 4000]$.

## ORCID iDs

A Mercuri-Baron https://orcid.org/0000-0003-2011-6097

A Di Piazza https://orcid.org/0000-0003-1018-0458

## References

[1] Di Piazza A, Müller C, Hatsagortsyan K and Keitel C 2011 Extremely high-intensity laser interactions with fundamental quantum systems Rev. Mod. Phys. 84 1177–228

[2] Cros B et al 2014 Laser plasma acceleration of electrons with multi-PW laser beams in the frame of CILEX Nucl. Instrum. Methods Phys. Res. A 740 27–33

[3] Nam C H, Sung J H, Lee H W, Youn J W and Lee S K 2018 Performance of the 20 fs, 4 PW Ti: sapphire laser at CoReLs Conf. Lasers and Electro-Optics (Optical Society of America) p STu4O.3

[4] Extreme light infrastructure (ELI) https://eli-laser.eu

[5] Bromage J et al 2019 High Power Laser Sci. Eng. 7 e4

[6] Exawatt center for extreme light studies (XCELS) http://xcels.iapras.ru

[7] Zetawatt-equivalent ultrashort pulse laser system (ZEUS) https://zeus.engin.umich.edu

[8] Reiss H R 1962 J. Math. Phys. 3 59

[9] Nikishov A I and Ritus V I 1964 Quantum processes in the field of a plane electromagnetic wave and in a constant field Sov. Phys. J. Exp. Theor. Phys. 19 529

[10] Duclous R, Kirk J and Bell A 2010 Monte Carlo calculations of pair production in high-intensity laser-plasma interactions Plasma Phys. Control. Fusion 53 10

[11] Elkina N V, Fedotov A M, Kostyukov I Y, Legkov M V, Narozhny N B, Nerush E N and Ruhl H 2011 QED cascades induced by circularly polarized laser fields Phys. Rev. Spec. Top. Accel. Beams 14 054401

[12] Arber T D et al 2015 Contemporary particle-in-cell approach to laser-plasma modelling Plasma Phys. Control. Fusion 57 113001

[13] Gonoskov A et al 2015 Extended particle-in-cell schemes for physics in ultrastrong laser fields: review and developments Phys. Rev. E 92 023305

[14] Grismayer T, Vranic M, Martins J L, Fonseca R A and Silva L O 2016 Laser absorption via quantum electrodynamics cascades in counter propagating laser pulses Phys. Plasmas 23 056706

[15] Lobet M, d'Humières E, Grech M, Ruyer C, Davoine X and Gremillet L 2016 Modeling of radiative and quantum electrodynamics effects in PIC simulations of ultra-relativistic laser-plasma interaction J. Phys.: Conf. Ser. 688 012058

[16] Bell A R and Kirk J G 2008 Possibility of prolific pair production with high-power lasers Phys. Rev. Lett. 101 200403

[17] Nerush E N, Kostyukov I Y, Fedotov A M, Narozhny N B, Elkina N V and Ruhl H 2011 Laser field absorption in self-generated electron–positron pair plasma Phys. Rev. Lett. 106 035001
[18] Grismayer T, Vranic M, Martins J L, Fonseca R A and Silva L O 2017 Seeded QED cascades in counterpropagating laser pulses
*Phys. Rev.* E **95** 023210

[19] Tamburini M, Di Piazza A and Keitel C H 2017 Laser-pulse-shape control of seeded QED cascades *Sci. Rep.* **7** 5694

[20] Jirka M, Klimo O, Vranic M, Weber S and Korn G 2017 QED cascade with 10 PW-class lasers *Sci. Rep.* **7** 15302

[21] Fedotov A, Narozhny N and Korn G 2010 Limitations on the attainable intensity of high power lasers *Phys. Rev. Lett.* **105** 080402

[22] Zhu X-L, Chen M, Yu T-P, Weng S-M, He F and Sheng Z-M 2019 Collimated GeV attosecond electron-positron bunches from a
plasma channel driven by 10 PW lasers *Matter Radiat. Extremes* **4** 014401

[23] Gelfer E, Mironov A, Fedotov A, Bashmakov V, Nerush E, Kostyukov I and Narozhny N 2015 Optimized multibeam
configuration for observation of QED cascades *Phys. Rev. A* **92** 022113

[24] Vranic M, Grismayer T, Fonseca R and Silva L 2016 Electron-positron cascades in multiple-laser optical traps *Plasma Phys. Control. Fusion* **59** 014040

[25] Gonoskov A et al 2017 Ultrabright GeV photon source via controlled electromagnetic cascades in laser-dipole waves *Phys. Rev. X* **7** 041003

[26] Ridgers C P, Brady C S, Duclous R, Kirk J G, Bennett K, Arber T D, Robinson A P L and Bell A R 2012 Dense electron-positron
plasmas and ultraintense γ rays from laser-irradiated solids *Phys. Rev. Lett.* **108** 165006

[27] Lobet M, Ruyer C, Debayle A, d'Humières E, Grech M, Lemoine M and Gremillet L 2015 Ultrafast synchrotron-enhanced
thermalization of laser-driven colliding pair plasmas *Phys. Rev. Lett.* **115** 215003

[28] Kostyukov I Y and Nerush E N 2016 Production and dynamics of positrons in ultrahigh intensity laser-foil interactions *Phys. Plasmas* **23** 093119

[29] Mackenroth F and Di Piazza A 2011 Nonlinear Compton scattering in ultrashort laser pulses *Phys. Rev. A* **83** 032106

[30] Blackburn T G, Ridgers C P, Kirk J G and Bell A R 2014 Quantum radiation reaction in laser-electron beam collisions *Phys. Rev. Lett.* **112** 015001

[31] Niel F, Riconda C, Amiranoff F, Duclous R and Grech M 2018 From quantum to classical modeling of radiation reaction: a focus
on stochasticity effects *Phys. Rev. E* **97** 043209

[32] Niel F, Riconda C, Amiranoff F, Lobet M, Derouillat J, Pérez F, Vinci T and Grech M 2018 From quantum to classical modeling
of radiation reaction: a focus on the radiation spectrum *Plasma Phys. Control. Fusion* **60** 094002

[33] Artemenko I I, Krygin M S, Serebryakov D A, Nerush E N and Kostyukov I Y 2019 Global constant field approximation for
radiation reaction in collision of high-intensity laser pulse with electron beam *Plasma Phys. Control. Fusion* **61** 074003

[34] Burke D L et al 1997 Positron production in multiphoton light-by-light scattering *Phys. Rev. Lett.* **79** 1626-9

[35] Lobet M, Davoine X, d'Humières E and Gremillet L 2017 Generation of high-energy electron-positron pairs in the collision of a
laser-accelerated electron beam with a multipetawatt laser *Phys. Rev. Accel. Beams* **20** 043401

[36] Blackburn T G, Ilderton A, Murphy C D and Marklund M 2017 Scaling laws for positron production in laser-electron-beam
collisions *Phys. Rev. A* **96** 022128

[37] Vranic M, Klimo O, Korn G and Weber S 2018 Multi-Gev electron-positron beam generation from laser-electron scattering *Sci. Rep.* **8** 4702

[38] Chen Y, Li J, Hatsagortsyan K and Keitel C 2018 γ-ray beams with large orbital angular momentum via nonlinear Compton
scattering with radiation reaction *Phys. Rev. Lett.* **121** 074801

[39] Blackburn T G and Marklund M 2018 Nonlinear Breit-Wheeler pair creation with Bremsstrahlung γ rays *Plasma Phys. Control. Fusion* **60** 054009

[40] Ta Phuoc K, Corde S, Thaury C, Malka V, Tafzi A, Goddet J P, Shah R C, Sebban S and Rousse A 2012 All-optical Compton
gamma-ray source *Nat. Photon.* **6** 308-11

[41] Gong Z et al 2017 High-efficiency γ-ray flash generation via multiple-laser scattering in ponderomotive potential well *Phys. Rev.* E **95** 013210

[42] Capdessus R, King M, Del Sorbo D, Duff M, Ridgers C P and McKenna P 2018 Relativistic Doppler-boosted γ-rays in high fields
*Sci. Rep.* **8** 9155

[43] Vranic M, Grismayer T, Meuren S, Fonseca R A and Silva L O 2019 Are we ready to transfer optical light to gamma-rays? *Phys. Plasmas* **26** 053103

[44] Magnusson J et al 2019 Laser-particle collider for multi-Gev photon production *Phys. Rev. Lett.* **122** 253801

[45] Zhu X-L et al 2020 Extremely brilliant GeV γ-rays from a two-stage laser-plasma accelerator *Sci. Adv.* **6** eaaz7240

[46] Sampath A et al 2021 Extremely dense gamma-ray pulses in electron beam-multifoil collisions *Phys. Rev. Lett.* **126** 064801

[47] Mendonça J T and Vieira J 2014 Donut wakefields generated by intense laser pulses with orbital angular momentum *Phys. Plasmas* **21** 033107

[48] Nuter R, Korneev P, Dmitriev E, Thiele I and Tikhonchuk V T 2020 Gain of electron orbital angular momentum in a direct laser
acceleration process *Phys. Rev. E* **101** 053202

[49] Longman A, Salgado C, Zeraouli G, Apınaniz J I, Pérez-Hernández J A, Eltahlawy M K, Volpe L and Fedosejevs R 2020 Off-axis
spiral phase mirrors for generating high-intensity optical vortices *Opt. Lett.* **45** 2187-90

[50] Duff M et al 2020 High order mode structure of intense light fields generated via a laser-driven relativistic plasma aperture *Sci. Rep.* **10** 105

[51] Allen L 1992 Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes *Phys. Rev. A* **45** 8185-9

[52] Di Piazza A 2016 *Phys. Rev. Lett.* **117** 213201

[53] Meuren S, Hatsagortsyan K Z, Keitel C H and Di Piazza A 2015 Polarization-operator approach to pair creation in short laser
pulses *Phys. Rev. D* **91** 013009

[54] Derouillat J et al 2018 Smilei: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation *Comput. Phys. Commun.* **222** 351-73

[55] Fedotov A M 2016 Quantum regime of laser-matter interactions at extreme intensities (arXiv:1612.02038v1)

[56] Hu H, Müller C and Keitel C H 2010 Complete QED theory of multiphoton trident pair production in strong laser fields *Phys. Rev. Lett.* **105** 080401

[57] He Y, Blackburn T G, Toncian T and Arefiev A V 2020 Dominance of γ-γ electron-positron pair creation in a plasma driven by
high-intensity lasers *Communications Physics* **4** 139

[58] Ritus V I 1979 Quantum effects of the interaction of elementary particles with an intense electromagnetic field *PhD Thesis* P N Lebedev Physical Institute of the USSR Academy of Sciences
[59] Baier V N, Katkov V M and Strakhovenko V M 1998 *Electromagnetic Processes at High Energies in Oriented Single Crystals* (Singapore: World Scientific)

[60] Erber T 1966 High-energy electromagnetic conversion processes in intense magnetic fields *Rev. Mod. Phys.* **38** 626–59

[61] Di Piazza A, Hatsagortsyan K Z and Keitel C H 2010 *Phys. Rev. Lett.* **105** 220403

[62] Tamburini M and Meuren S 2020 arXiv:1912.07508

[63] Podszus T and Di Piazza A 2021 arXiv:2103.14637

[64] Siegman A E 1986 *Lasers* revised edn (Sausalito: University Science Books)

[65] Nuter R, Grech M, Gonzalez de Alaiza Martinez P, Bonnaud G and d'Humières E 2014 Maxwell solvers for the simulations of the laser-matter interaction *Eur. Phys. J. D* **68** 177

[66] Pérez F and Grech M 2019 Oblique-incidence, arbitrary-profile wave injection for electromagnetic simulations *Phys. Rev. E* **99** 033307

[67] Chen Y-Y, He P-L, Shaisultanov R, Hatsagortsyan K Z and Keitel C H 2019 Polarized positron beams via intense two-color laser pulses *Phys. Rev. Lett.* **123** 174801

[68] Seipt D and King B 2020 Spin- and polarization-dependent locally-constant-field-approximation rates for nonlinear Compton and Breit–Wheeler processes *Phys. Rev. A* **102** 052805

[69] Lee S K, Kim H T, Choi I W, Kim C M and Nam C H 2018 Exploration of strong field physics with multi-PW lasers *J. Korean Phys. Soc.* **73** 179–89

[70] Meuren S et al 2020 On seminal HEDP research opportunities enabled by colocating multi-petawatt laser with high-density electron beams (arXiv:2002.10051v1)

[71] Zhang P, Bulanov S S, Seipt D, Arefiev A V and Thomas A G R 2020 Relativistic plasma physics in supercritical fields *Phys. Plasmas* **27** 050601

[72] Grech M et al 2021 Investigating strong-field QED processes in laser-electron beam collisions at Apollon https://hal.archives-ouvertes.fr/hal-03229914

[73] Papadopoulos D N et al 2016 The Apollon 10 PW laser: experimental and theoretical investigation of the temporal characteristics *High Power Laser Sci. Eng.* **4** e34

[74] Yakimenko V et al 2019 Facet-II facility for advanced accelerator experimental tests *Phys. Rev. Accel. Beams* **22** 101301

[75] Abramowicz H et al 2019 Letter of intent for the LUXE experiment (arXiv:1909.00860)