Note on the conjectured breakdown of QED perturbation theory in strong fields
NNote on the conjectured breakdown of QED perturbation theory in strong fields
A. Ilderton ∗ Centre for Mathematical Sciences, University of Plymouth, Plymouth, PL4 8AA, UK
Strong background fields require a non-perturbative treatment, which is afforded in QED by theFurry expansion of scattering amplitudes. It has been conjectured that this expansion breaks downfor sufficiently strong fields, based on the asymptotic growth of loop corrections with increasing“quantum nonlinearity”, essentially the product of field strength and particle energy. However, cal-culations to date have assumed that the background is constant. We show here, using general planewaves of finite duration, that observables at high quantum nonlinearity scale differently depend-ing on whether intensity or energy is large. We find that, at high energy, loop contributions toobservables tend to fall with increasing quantum nonlinearity, rather than grow.
I. INTRODUCTION
A strong electromagnetic field is characterised by a dimensionless coupling to matter which is larger than one. Hencethe interaction of the field with matter cannot be treated in perturbation theory. For QED processes in a strong fieldthe required semi-non-perturbative treatment is given by the Furry expansion of scattering amplitudes. Here the cou-pling between matter and generated/absorbed photons, still characterised by the fine-structure constant α , is treatedin perturbation theory as normal, while the coupling to the strong field is treated exactly. This amounts to replacingthe fermion propagator by a dressed propagator, while all other position space Feynman rules are unchanged [1].It has however been conjectured that for sufficiently strong fields, this semi-perturbative expansion also breaksdown [2–6], and the theory becomes “fully nonperturbative”. This conjecture is based on the identification of aneffective coupling parameter αχ / , where χ is essentially the product of field strength and probe particle energy (seebelow). As field strength increases, loop corrections appear to grow with higher powers of the effective coupling, andso when αχ / ∼ the Furry expansion breaks down. The physics of the regime αχ / ∼ is thus unknown, butfor experimental proposals for how to approach it see [7–9]. For a review of the conjecture see [10]. (For effectiveexpansion parameters in relation to muon g − see [11–13].)The calculations behind the conjecture have, however, been performed in constant crossed fields; these are the zerofrequency limit of plane waves, commonly used as a first model of laser fields at ultra-high intensity [14–19]. (We willtypically refer to field intensity, rather than strength.) Notably, the power of / which appears is tied to the Airyfunctions particular to the constant field case. This prompts the question of what happens in more general fields.Furthermore, the literature to date has tended to focus on (loop corrections to) the polarisation and mass operators,neither of which are observables. Indeed, the one-loop vertex correction has been argued to scale asymptotically aseither αχ / or just α , depending on the gauge used [6, 20].There are, then, several issues to address. First, we should only consider gauge invariant observables. Second, fieldswhich decay at infinity, rather than constant fields, should be the generic case; here we will work with plane wavepulses of finite duration and arbitrary intensity, which can be treated exactly in the Furry expansion [14, 15, 17, 18].Third, observables do not in general depend on χ alone, but on intensity and energy individually. A dependencesolely on χ is again particular to the constant field limit, but χ can also be made large by making the energy large.(Literature investigations of one-loop processes in pulsed fields already contain hints that the scalings attributed to theconstant field case may not be universal, see for example [21, 22].) Fourth, while loop corrections have been calculatedin the high intensity approximation, emissivity corrections and inclusive observables have not been discussed, and weexpect higher numbers of emissions to also become more important at high intensity. In this note we will presentan initial investigation of these issues using the few one-loop diagrams which have been calculated exactly in pulsedplane waves, rather than constant fields. We will contrast two high- χ limits, reached by high intensity or high energy.We introduce relevant parameters in Sect. II. In Sect. III we consider nonlinear Compton scattering at tree level,and forward scattering at one loop. In Sect. IV we consider photon helicity flip at one loop. In both cases we will showthat observables scale differently in the two high- χ limits, demonstrating that there is no universal high- χ behaviour.The high- χ (energy) limit, in particular, does not seem to suffer from a potential breakdown of perturbation theoryas the high- χ (intensity) limit does. A further consequence of our results is that, as we will show, currently used high-intensity approximations in numerical laser-plasma models are unable to properly reproduce high-energy, quantum,behaviour. We conclude in Sect. V. Closely related calculations have very recently appeared in [23], and we will findagreement with the conclusions presented there. ∗ [email protected] a r X i v : . [ h e p - ph ] M a r F = i F = N = FIG. 1. F is the Furry picture propagator in a background field, F its one-loop correction, and N tree level photon emission,“nonlinear Compton scattering”. II. PARAMETERS & INVARIANTS
Consider the interaction of some particle, momentum p µ , with a strong background field F µν ( x ) . A key parameteris the “quantum nonlinearity” χ of the particle, defined by [14] χ = em (cid:112) p µ F µν F νσ p σ , (1)in which e and m are, by convention, the electron charge and mass. We will only deal with electrons and photonshere. It is clear that χ looks, in general, like the product of particle energy and field strength, with the proportionalityfactor depending on field- and collision geometry. (For an electron, χ is equal to the ratio of the electric field in the restframe, to the Schwinger field E S = m /e .) Since χ is a composite parameter it can be made large by increasing theparticle energy, the field strength, or both. We emphasise that the physics differs depending on whether a particularvalue of χ is reached by going to high energy, or high intensity [24].Our focus here is on plane waves. For a wave characterised by the lightlike propagation direction n µ we have χ = n.pm | E | E S . (2)The plane wave will have some typical scale ω associated with it, be it a central frequency, inverse width, etc, and itis convenient to define dimensionless variables in terms of this scale . So, let k µ = ωn µ be a typical momentum vectorassociated to the wave, and φ = k.x the dimensionless “lightfront time” on which the wave depends. This allows usto define the dimensionless energy, b , and field intensity parameter, a , by b := k.pm , a := χb = | e E | mω . (3)(Note that χ = a b is independent of the chosen scale ω .) The intensity parameter a is the coupling between theplane wave and matter, and easily exceeds unity in modern laser experiments [25–27]. For, say, a head-on collisionbetween an electron and the wave we have b (cid:39) ω/m ) γ at high energy.Below we will compare the behaviour of processes at high χ reached via high intensity or high energy, meaning large a or large b respectively. There are many classic examples of processes which exhibit a very different dependence ona single (usually energy) invariant, compare for example the low [28] and high [29] energy behaviours of the photon-photon cross-section, see also [30] for a review. Here, though, we will examine how the two independent invariants a and b essentially compete to affect different behaviours of processes in strong fields. III. FORWARD SCATTERING AND NONLINEAR COMPTON SCATTERING
We begin with the diagrams in Fig. 1. These are the tree level fermion propagator, its one-loop correction, andthe tree level vertex, in the Furry picture. The loop diagram F (associated with the mass operator in the literatureon αχ / ) contributes to forward scattering of the electron, and to electron spin-flip, while the vertex yields e.g. treelevel photon emission from a field-accelerated electron, or “nonlinear Compton scattering" (NLC) [14, 31–35].Consider the electron forward scattering amplitude; this is degenerate with soft emission [36–38], so we shouldconsider not just exclusive but also inclusive processes. Let P ( b (cid:12)(cid:12) a ) be the probability of an exclusive scatteringprocess a → b as calculated in the Furry expansion, and let P ( (cid:94) b (cid:12)(cid:12) a ) be the inclusive probability. Then the probabilityof observing an electron scattered with no photon emission above the detector resolution ε min of the system is, writing e for electron and γ s for a soft photon, The physics is independent of this choice , but not of the physical scale ω . If, as for constant fields, no such scale is apparent, one cansimply use the electron rest mass. Any natural scale will emerge during the evaluation of the spacetime integrals to be introduced below. P ( (cid:94) e (cid:12)(cid:12) e ) = P ( e (cid:12)(cid:12) e ) + P ( e, γ s (cid:12)(cid:12) e ) + . . . = (cid:90) ε IR |F + i F + . . . | + ε min (cid:90) ε IR |N + . . . | + . . . = 1 − (cid:90) ε IR Im ( F (cid:63) F ) + ε min (cid:90) ε IR |N | + O ( α ) , (4)where ε IR is the infra-red (IR) cutoff. By the optical theorem the one-loop imaginary part is exactly equal to the treelevel probability of NLC. (For the explicit demonstration of this in plane waves see [39].) Hence we have P ( (cid:94) e (cid:12)(cid:12) e ) = 1 − (cid:90) ∞ ε IR |N | + ε min (cid:90) ε IR |N | + O ( α ) = 1 − ∞ (cid:90) ε min |N | + O ( α ) . (5)For plane waves we can take ε min → without introducing any IR divergence, assuming as usual that there is no DCmode [40]. The forward scattering probability is then equal to one minus the tree-level NLC probability, demonstratingthe standard IR result that the proper inclusion of higher loop and emissivity corrections can reduce probabilities [36].The diagrams necessary to study this in depth in the context of αχ / have not yet been calculated, however. (Wereturn to this in the conclusions.) Here, we observe that the tree-level NLC probability also determines the one-loopforward-scattering probability, and so we turn to the explicit calculation of NLC in two, distinct, high- χ limits.Consider an electron, momentum p µ , colliding with a plane wave depending on lightfront time k.x , so b = k.p/m .The wave may be described by a potential eA µ /m = a µ ( k.x ) which is transverse, k.a ( k.x ) = 0 , and where the profile a µ ( k.x ) is characterised by amplitude a , with a the intensity parameter [15, 18]. The electron emits a photon ofmomentum k (cid:48) µ and scatters. The total probability of this NLC scattering at tree level is given by an integral over thelightfront momentum fraction s := k.k (cid:48) /k.p of the emitted photon, and over two lightfront times θ and φ as [41] P ( e, γ (cid:12)(cid:12) e ) = − απb (cid:90) d s (cid:90) d φ ∞ (cid:90) d θ sin (cid:18) s θµ b (1 − s ) (cid:19)(cid:20) µ ∂µ∂θ + (cid:18)
12 + 14 s − s (cid:19) (cid:104) a (cid:48) (cid:105) θ (cid:21) , (6)where we define, here and throughout, a floating average (cid:104)·(cid:105) and Kibble’s (normalised) effective mass µ by [42, 43] (cid:104) f (cid:105) = 1 θ φ + θ/ (cid:90) φ − θ/ d( k.x ) f ( k.x ) , µ ( φ, θ ) = 1 − (cid:104) a (cid:105) + (cid:104) a (cid:105) . (7)Note that the probability (6) is a function of a and b individually. The s -integral can be performed analytically [41]. A. High intensity
The high-intensity limit is reached when a (cid:29) . An old result is that the “formation length” of quantum pro-cesses scales in this limit as /a , so that amplitudes can be calculated in the locally constant field approximation(“LCFA") [14]. (It is now known that this argument fails for the emission of photons with low energy [44] or lowlightfront energy [45–47].) The scaling of processes in constant fields, on which the αχ / conjecture is based, is sharedby the LCFA. We consider the slightly more general case of the LCFA because of its use in particle-in-cell codes, usedto plan and analyse intense laser experiments, for a review see [48]. Taking a (cid:29) , the LCFA for NLC is [14, 46, 47], P ( e, γ (cid:12)(cid:12) e ) (cid:39) − αb (cid:90) d φ (cid:90) d s Ai ( z ) + (cid:18) z + χ γ √ z (cid:19) Ai (cid:48) ( z ) , where z ( φ ) := (cid:18) χ e ( φ ) s − s (cid:19) . (8)Note that, aside from the /b prefactor, the LCFA (and constant field) expressions now depend only on χ . Theasymptotic scaling of the probability is obtained by expanding the Airy functions, which yields P ( e, γ (cid:12)(cid:12) e ) ∼ αb (cid:90) d φ χ / ( φ ) . (9)We are only interested in typical scalings, so “ ∼ ” indicates throughout that we neglect purely numerical factors,keeping track only of dependence on the important parameters b , a , χ and, below, pulse shape effects. For pulses, E / E S a / a - - - ϕ - - FIG. 2. The potential and electric field of the short pulse a ( φ ) = a φ exp( − φ ) , where φ = x + /τ for some width τ , as usedin the text. the integral in (9) generates a finite factor which, if we consider a short pulse characterised only by some width τ anduse the choice of scale described in Sect. II, is independent of parameters. For the constant field limit proper, χ isconstant and the φ -integral generates a (dimensionless) length factor. Either way we have, for some peak value of χ , P ( e, γ (cid:12)(cid:12) e ) ∼ αχ / b ∼ αχ / γ ( mτ ) . (10)The first expression gives the literature scaling. The pulse length is made explicit in the second expression. As χ increases with increasing intensity, i.e. at fixed b , the probability clearly grows. It easily exceeds unity, even for shortpulses, demonstrating the need for higher order corrections [41]. B. High energy
The high-energy, b → ∞ , limit of NLC may not be so well known, so we give a little more detail. Starting from (6),we rescale θ → bθ . Then we use the assumption that the plane wave has finite duration, for then µ ( φ, ∞ ) = 1 [43, 49],i.e. the effective mass at large distance equals the free mass. To leading order in b (cid:29) we can then replace µ → .It may be checked that the first term in square brackets of (6) only contributes sub-leading terms in the high-energylimit, so we turn to the second term. We Fourier transform each factor of a (cid:48) appearing. The overall d φ integral isthen trivial, setting the Fourier variables to be equal and opposite, such that we obtain ˜ a ( ν ) . ˜ a (cid:63) ( ν ) in the integrandfor ν the remaining Fourier variable. Writing t ≡ s/ (1 − s ) , the θ -integral reduces at this stage to ∞ (cid:90) d θθ sin νθ θt b = π Θ(2 b | ν | − t ) . (11)This puts a limit on the s –integral, which is most easily performed by changing variables to t . The final result for theleading order high energy limit is P ( e, γ (cid:12)(cid:12) e ) = − α b (cid:90) d ν π ˜ a ( ν ) . ˜ a (cid:63) ( ν ) (cid:18)
12 + log 2 b | ν | (cid:19) + . . . . (12)This describes the convolution of the plane wave intensity profile ˜ a ( ν ) . ˜ a (cid:63) ( ν ) with the term in large brackets, whichis just the high-energy limit of ordinary Compton scattering of an electron, momentum p µ , against a photon ofmomentum νk µ (see e.g. [51, §6.1]). The whole result is quadratic in a and is, interestingly, exactly the same as thehigh-energy limit of the wholly perturbative , small a , calculation. The remaining Fourier integral is finite for finitepulses, and it is clear that the asymptotic behaviour of the probability is P ( e, γ (cid:12)(cid:12) e ) ∼ αa b log b . (13) Interestingly, this method is similar to that needed to obtain the s → ( low lightfront energy) limit of the differential probability. Whilethis appears to be zero from (6), the correct nonzero value is only obtained after performing the lightfront time integrals [46, 50]. ϵ ′ µ ϵ µ FIG. 3. The one-loop diagram contributing to helicity flip and forward scattering at one loop (and, via the optical theorem,pair production at tree level).
This may be confirmed with an example. Consider the pulse in Fig. 2, with a ( φ ) = a φ exp( − φ ) and linear polari-sation. Then we can evaluate the Fourier integral in (12) exactly to find P ( e, γ (cid:12)(cid:12) e ) = αa b (cid:114) π (cid:0) log b + 12 log 2 + 32 − γ E (cid:1) + . . . . (14) C. Discussion
The expressions (10) and (13) give two different high- χ limits of both NLC at tree level, and electron forwardscattering at one loop. These two limits have different functional forms, and a comparison is facilitated by making itexplicit that the high intensity limit corresponds to increasing χ by increasing a , at fixed b , while the high energylimit corresponds to increasing χ by increasing b at fixed a . This allows us to write (9) and (13) as P ( e, γ (cid:12)(cid:12) e ) χ ∼ a →∞ ∼ αχ / b , P ( e, γ (cid:12)(cid:12) e ) χ ∼ b →∞ ∼ α a χ log χ . (15)Hence the nonlinear Compton probability falls with increasing χ at high energy, rather than increasing as it doesat high intensity. Furthermore, this implies that at high- χ reached via high energy, the loop correction to forward-scattering is smaller than the tree level contribution, again in contrast to the high intensity limit. This showsexplicitly that there is no universal high- χ behaviour of observables. Also, at least for the observables consideredhere, the increase in size of higher-loop corrections at high intensity does not appear at high energy.The LCFA is a commonly employed tool which allows for progress where analytic results are lacking, especially inthe consideration of inhomogeneous background fields with realistic spacetime structure. It is already known thatthe CCF and LCFA approximations, as currently employed, give incorrect results at low photon energy [44] and lowlightfront momentum [45–47], and we can now show that they also fail in the high-energy regime.The LCFA for nonlinear Compton is (aside from the /b prefactor) a function of χ alone, and the high- χ limit istherefore (9); it is impossible to obtain the correct high energy behaviour in (15) from the LCFA expression. Thescaling is different, and furthermore the logarithmic dependence in (15) is missed. The reason for this is that the LCFArequires, at the level of the integrated rates as considered here, not just a (cid:29) but a (cid:29) b [24, 52], so that makingthe assumption that the LCFA holds precludes the possibility of going to high-energy. Our example of NLC makesthis concrete; the high-intensity limit (10) comes from the first term in square brackets of (6), while the high-energylimit (13) comes from the second term. This shows straightforwardly that the two limits do not commute.It follows that approximations used in particle-in-cell codes break down at high energy. It is therefore important tounderstand, when modelling potential experiments, whether high χ is being reached via high energy or high intensity. IV. HELICITY FLIP
For our next example we consider photon helicity flip. In a plane wave a single photon can produce pairs, orotherwise only scatter forward (due to the many symmetries of plane waves), and may flip helicity due to loopinteractions with the wave. Let the photon have momentum k (cid:48) µ , so that b = k.k (cid:48) /m is the energy parameter, and let (cid:15) µ and (cid:15) (cid:48) µ be two orthogonal helicity polarisation vectors, see Fig. 3. Then the probability P ( γ (cid:48) | γ ) of helicity flip isgiven at one-loop by P ( γ (cid:48) | γ ) = | T | , where the amplitude T may be expressed, using lightfront field theory [53–55],as a double integral [21, 56] over two lightfront times φ and θ . Explicitly [21], T = − απ b ∞ (cid:90) −∞ d φ ∞ (cid:90) d θ θ (cid:18) K (cid:0) θµb (cid:1)(cid:0) ¯ a θ a θ − ¯ a φ a φ (cid:1) + K (cid:0) θµb (cid:1)(cid:0) ¯ a φ a θ − ¯ a θ a φ (cid:1)(cid:19) , (16)in which we have defined a φ = ∂ φ (cid:15). (cid:104) a (cid:105) , ¯ a φ = ∂ φ (cid:15) (cid:48) (cid:63) . (cid:104) a (cid:105) , a θ = ∂ θ (cid:15). (cid:104) a (cid:105) , ¯ a θ = ∂ θ (cid:15) (cid:48) (cid:63) . (cid:104) a (cid:105) , (17)and the two K -functions are combinations of modified Bessel functions arising from the integral over the virtuallightfront momentum fraction in the loop [21] K ( x ) = ixe − ix (cid:0) K ( ix ) − K ( ix ) (cid:1) = (cid:90) d s exp (cid:20) − ix s (1 − s ) (cid:21) , K ( x ) = (1 − i∂ x ) K ( x ) . (18) A. High intensity
Consider a linearly polarised plane wave. (The same results are found for other polarisation choices, but theintermediate expressions are not as clear.) Then the high intensity, LCFA, approximation to the flip amplitude is [21] T ∼ αb (cid:90) d φ (cid:90) d s z (cid:0) Ai (cid:48) ( z ) − i Gi (cid:48) ( z ) (cid:1) , where z = (cid:18) χ γ ( φ ) s (1 − s ) (cid:19) / , (19)and Gi is the Scorer function. As before, we simply replace the φ -integral by a volume factor for the constant fieldcase. The individual dependence on a and b seen in the integrand of (16) [21] is again replaced by a dependence onlyon their product, χ . The asymptotic behaviour of (19) is easily extracted from the known expansion of the Airy andScorer functions as T ∼ αb (cid:90) d φ χ / γ ( φ ) ∼ αχ / b = ⇒ P ( γ (cid:48) | γ ) ∼ α χ / b . (20)This is also the scaling in the constant field case [2, 10, 57, 58] (for which the numerical value given by the φ -integralagain becomes a length factor). B. High energy
We turn to the high-energy limit of helicity flip. Consider the K -term in (16). The high energy limit may becalculated simply by expanding in powers of the small parameter /b , replacing K → . By writing the averages interms of Fourier integrals it is then easy to perform the φ and θ integrals. Turning to the K -term, we first integrateby parts in θ (the boundary terms are zero), such that we take the derivative of θ K . (The integration of the averagesis easily performed in Fourier space.) This has the effect of improving convergence under the integral, such that wecan again expand in powers of /b . Expanding − ∂ θ ( θ K ) in this way the leading order terms are log b − − γ E − log( iθµ/
2) + O (cid:18) b (cid:19) + O (cid:18) log bb (cid:19) . (21)This multiplies a function which falls like /θ for large θ . Hence for large enough b the log b term will dominate, andwe can, for pulsed fields, neglect the log θµ term. (This can be confirmed by changing variables θ → θb as before andexpanding again; one finds only higher order terms in /b .) Carrying out the integrals, the final result may be writtenin terms of the Fourier transform ˜ a µ ( ν ) as T ∼ απb (cid:90) d ν π ¯ (cid:15) (cid:48) . ˜ a ( ν ) (cid:15). ˜ a (cid:63) ( ν ) (cid:18) − iπ sign ( ν ) log b (cid:19) + . . . (22)For generic pulses the Fourier integral converges. The first term in (22) is nonzero for helicity states, but whetherthe second term survive depends on the polarisation of the field. Generically, then, for a helicity- and polarisation–dependent constant σ we have, T ∼ αa b (1 + σ b ) = ⇒ P ( γ (cid:48) | γ ) ∼ α a b (cid:0) σ log b (cid:1) . (23) C. Discussion
In order to compare the high- χ (intensity) limit (20) of helicity flip with the high- χ (energy) limit (23) we againrewrite the latter in terms of χ , increasing with b at fixed a . This gives, to leading order, P ( γ (cid:48) | γ ) χ ∼ a →∞ ∼ α χ / b , P ( γ (cid:48) | γ ) χ ∼ b →∞ ∼ α a χ (cid:0) σ log χ (cid:1) . (24)These two high- χ limits exhibit the same trends as for NLC, above. In the high- χ (energy) limit, the probabilityof helicity flip is manifestly decreasing with χ , unlike in the high intensity limit. This agrees with the numericalevaluation of the exact result (16) given in [21]. The LCFA approximation again fails to capture the logarithmicbehaviour at high energy, and for the same reason as above; the high intensity limit precludes taking the high energylimit. Note again that the high- χ (energy) limit coincides with that of the lowest-order perturbative calculation ofthe process. There are, as for the other processes considered here, higher-order corrections depending on a /b . V. CONCLUSIONS
We have considered the behaviour of QED scattering probabilities in background plane wave fields, in the high- χ limit. This was in the context of the conjectured breakdown of (Furry picture) perturbation theory in the regime αχ / ∼ . We have shown, though, that there is no unique high χ behaviour. The high intensity and high energylimits, which both give high χ , yield different scalings. At high intensity, both constant field results and locallyconstant field approximations of pulsed field results show the same power law scaling with αχ / . The high energylimit, on the other hand, shows for general pulsed fields a logarithmic dependence on χ which is typical of QED.Further, we have seen that observables tend to fall with χ at high energy, rather than rise, as at high intensity.This suggests that the high- χ limit of scattering processes in strong fields, reached via high energy, may not exhibitthe perturbative breakdown attributed to the high- χ limit reached via high intensity. More work is however neededto better judge this, especially since only a few diagrams have been calculated in the level of detailed needed.We have also seen that the high-energy behaviour of the considered rates cannot be recovered once the locallyconstant crossed field approximation (LCFA) has been made: making that approximation on the parameter spaceprecludes being able to take the high-energy limit. An immediate consequence is that LCFA-based particle-in-cellsimulations used to model laser-matter interactions do not correctly capture high-energy quantum effects.Note that we do not disagree with previously calculated constant field scalings. However, we have stressed that it isimportant to consider observables. It would in particular be interesting to examine inclusive observables, and to seehow the consistent inclusions of both loops [3, 4, 59] and emission affects their behaviour. It may be that this reducesprobabilities, as happens with the exponentiation of infra-red corrections in QED [36] (and in background planewaves [40]), and so brings the perturbative expansion back under control. Alternatively, if we view the breakdownof perturbation theory as a breakdown of the background field approximation, then it may be that including someform of back-reaction [60] at each order of perturbation theory is enough to give a better behaved series. These arechallenging and interesting topics for future study.Finally, we note that the high energy limit of the observables we have considered is equal to that obtained ifthe plane wave background is treated perturbatively, instead of exactly using the Furry expansion. This may offersimplifications for the calculation of more complex processes at high energy. ACKNOWLEDGMENTS
Thanks to Tom Heinzl and Ben King for useful discussions. As this manuscript was being prepared, [23] appeared onthe arXiv, the subject of which is very similar to our own. Encouragingly, while our methods differ, our conclusionsare in agreement. Thanks to Antonino Di Piazza for correspondence on this matter. [1] W. H. Furry, Phys. Rev. (1951) 115.[2] V. Ritus, Sov. Phys. JETP , 1181 (1970).[3] N. Narozhny, Phys. Rev. D , 1313 (1979).[4] N. Narozhny, Phys. Rev. D , 1176 (1980). [5] D. Morozov and N. Narozhny, Sov. Phys. JETP , 23 (1977).[6] D. Morozov and N. Narozhny, Sov. Phys. JETP , 1103 (1981).[7] V. Yakimenko et al. , (2018), arXiv:1807.09271 [physics.plasm-ph].[8] T. G. Blackburn, A. Ilderton, M. Marklund, and C. P. Ridgers, (2018), arXiv:1807.03730 [physics.plasm-ph].[9] C. Baumann, A. Pukhov, E. N. Nerush, and I. Yu. Kostyukov, (2018), arXiv:1811.03990 [physics.plasm-ph].[10] A. M. Fedotov, Proceedings, LPHYS’16: Yerevan, Armenia, July 11-15, 2016 , J. Phys. Conf. Ser. , 012027 (2017),arXiv:1608.02261 [hep-ph].[11] H. Terazawa, Prog. Theor. Phys. (1968) 830. https://doi.org/10.1143/PTP.40.830[12] B. E. Lautrup and E. de Rafael, Nucl. Phys. B (1974) 317 Erratum: [Nucl. Phys. B (1974) 576].https://doi.org/10.1016/0550-3213(74)90600-2, https://doi.org/10.1016/0550-3213(74)90481-7[13] J. P. Miller, E. de Rafael and B. L. Roberts, Rept. Prog. Phys. (2007) 795 https://doi.org/10.1088/0034-4885/70/5/R03arXiv:hep-ph/0703049.[14] V. I. Ritus, J. Russ. Laser Res. , 497 (1985).[15] A. Di Piazza, C. Muller, K. Z. Hatsagortsyan, and C. H. Keitel, Rev. Mod. Phys. , 1177 (2012), arXiv:1111.3886[hep-ph].[16] N. Narozhny and A. Fedotov, Contemporary Physics , 249 (2015), https://doi.org/10.1080/00107514.2015.1028768.[17] B. King and T. Heinzl, High Power Laser Science and Engineering , e5 (2016), arXiv:1510.08456 [hep-ph].[18] D. Seipt, in Proceedings, HQ 2016: Dubna, Russia, July 18-30, 2016 (2017) pp. 24–43, arXiv:1701.03692 [physics.plasm-ph].[19] A. Hartin, Int. J. Mod. Phys.
A33 , 1830011 (2018), arXiv:1804.02934 [hep-ph].[20] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Nucl. Phys.
B563 , 361 (1999), arXiv:hep-ph/9908320 [hep-ph].[21] V. Dinu, T. Heinzl, A. Ilderton, M. Marklund, and G. Torgrimsson, Phys. Rev.
D89 , 125003 (2014), arXiv:1312.6419[hep-ph].[22] H. Gies, F. Karbstein, and R. Shaisultanov, Phys. Rev.
D90 , 033007 (2014), arXiv:1406.2972 [hep-ph].[23] T. Podszus and A. Di Piazza, (2018), arXiv:1812.08673 [hep-ph].[24] V. Dinu, C. Harvey, A. Ilderton, M. Marklund, and G. Torgrimsson, Phys. Rev. Lett. , 044801 (2016), arXiv:1512.04096[hep-ph].[25] G. Sarri et al. , Phys. Rev. Lett. , 224801 (2014), arXiv:1407.6980 [physics.optics].[26] J. M. Cole et al. , Phys. Rev. X8 , 011020 (2018), arXiv:1707.06821 [physics.plasm-ph].[27] K. Poder et al. , Phys. Rev. X8 , 031004 (2018), arXiv:1709.01861 [physics.plasm-ph].[28] H. Euler and B. Kockel, Naturwiss. , 246 (1935).[29] A. Akhiezer, L. Landau, and I. Pomeranchuk, Nature , 206.1 (1936).[30] K. Scharnhorst, (2017), arXiv:1711.05194 [physics.hist-ph].[31] A. I. Nikishov and V. I. Ritus, Sov. Phys. JETP , 529 (1964), [Zh. Eksp. Teor. Fiz.46,776(1964)].[32] M. Boca and V. Florescu, Phys. Rev. A80 , 053403 (2009).[33] T. Heinzl, D. Seipt, and B. Kampfer, Phys. Rev.
A81 , 022125 (2010), arXiv:0911.1622 [hep-ph].[34] D. Seipt and B. Kampfer, Phys. Rev.
A83 , 022101 (2011), arXiv:1010.3301 [hep-ph].[35] F. Mackenroth, A. Di Piazza, and C. H. Keitel, Phys. Rev. Lett. , 063903 (2010), arXiv:1001.3614 [physics.acc-ph].[36] D. R. Yennie, S. C. Frautschi, and H. Suura, Annals Phys. , 379 (1961).[37] M. Lavelle and D. McMullan, JHEP , 026 (2006), arXiv:hep-ph/0511314 [hep-ph].[38] M. Lavelle, D. McMullan, and T. Steele, Adv. High Energy Phys. , 379736 (2012), arXiv:1008.3949 [hep-ph].[39] A. Ilderton and G. Torgrimsson, Phys. Rev. D88 , 025021 (2013), arXiv:1304.6842 [hep-th].[40] A. Ilderton and G. Torgrimsson, Phys. Rev.
D87 , 085040 (2013), arXiv:1210.6840 [hep-th].[41] V. Dinu, Phys. Rev.
A87 , 052101 (2013), arXiv:1302.1513 [hep-ph].[42] T. W. B. Kibble, A. Salam, and J. A. Strathdee, Nucl. Phys.
B96 , 255 (1975).[43] C. Harvey, T. Heinzl, A. Ilderton, and M. Marklund, Phys. Rev. Lett. , 100402 (2012), arXiv:1203.6077 [hep-ph].[44] V. Dinu, T. Heinzl, and A. Ilderton, Phys. Rev.
D86 , 085037 (2012), arXiv:1206.3957 [hep-ph].[45] C. N. Harvey, A. Ilderton, and B. King, Phys. Rev.
A91 , 013822 (2015), arXiv:1409.6187 [physics.plasm-ph].[46] A. Di Piazza, M. Tamburini, S. Meuren, and C. H. Keitel, Phys. Rev.
A98 , 012134 (2018), arXiv:1708.08276 [hep-ph].[47] A. Ilderton, B. King, and D. Seipt, (2018), arXiv:1808.10339 [hep-ph].[48] A. Gonoskov, S. Bastrakov, E. Efimenko, A. Ilderton, M. Marklund, I. Meyerov, A. Muraviev, A. Sergeev, I. Surmin, andE. Wallin, Phys. Rev.
E92 , 023305 (2015), arXiv:1412.6426 [physics.plasm-ph].[49] F. Hebenstreit, A. Ilderton, M. Marklund, and J. Zamanian, Phys. Rev.
D83 , 065007 (2011), arXiv:1011.1923 [hep-ph].[50] T. G. Blackburn, D. Seipt, S. S. Bulanov, and M. Marklund, Phys. Plasmas , 083108 (2018), arXiv:1804.11085[physics.plasm-ph].[51] M. Kaku, Quantum field theory: A Modern introduction (New York, USA: Oxford Univ. Pr. (1993) 785 p, 1993).[52] A. K. Khokonov and M. K. Khokonov, Tech. Phys. Lett. , 154 (2005).[53] S. J. Brodsky, H.-C. Pauli, and S. S. Pinsky, Phys. Rept. , 299 (1998), arXiv:hep-ph/9705477 [hep-ph].[54] T. Heinzl, Lect. Notes Phys. , 55 (2001), arXiv:hep-th/0008096 [hep-th].[55] B. L. G. Bakker et al. , Nucl. Phys. Proc. Suppl. , 165 (2014), arXiv:1309.6333 [hep-ph].[56] S. Meuren, K. Z. Hatsagortsyan, C. H. Keitel, and A. Di Piazza, Phys. Rev. D91 , 013009 (2015), arXiv:1406.7235 [hep-ph].[57] N. Narozhnyi, Sov. Phys. JETP , 371 (1969).[58] F. Karbstein, Phys. Rev. D88 , 085033 (2013), arXiv:1308.6184 [hep-th].[59] D. Morozov and V. Ritus, Nuc. Phys. B , 309 (1975).[60] D. Seipt, T. Heinzl, M. Marklund, and S. S. Bulanov, Phys. Rev. Lett.118