Asymptotic behavior of the QED perturbation series

In this talk, I will summarize the present state of a long-term effort to obtain information on the high-order asymptotic behaviour of the QED perturbation series through the effective action. Starting with the constant-field case, I will discuss the Euler-Heisenberg Lagrangian in various dimensions, and up to the three-loop level. This Lagrangian holds the information on the N-photon amplitudes in the low-energy limit, and combining it with spinor helicity methods explicit all-N results can be obtained at the one-loop and, for the"all +"amplitudes, also at the two-loop level. For the imaginary part of the Euler-Heisenberg Lagrangian, an all-loop formula has been conjectured independently by Affleck, Alvarez and Manton for Scalar QED, and by Lebedev and Ritus for Spinor QED. This formula can be related through a Borel dispersion relation to the leading large-N behaviour of the N-photon amplitudes. It is analytic in the fine structure constant, which is puzzling and suggests a diagrammatic investigation of the large-N limit in perturbation theory. Preliminary results of such a study for the 1+1 dimensional case throw doubt on the validity of the conjecture.

behaviour of the N-photon amplitudes. It is analytic in the fine structure constant, which is puzzling and suggests a diagrammatic investigation of the large-N limit in perturbation theory. Preliminary results of such a study for the 1+1 dimensional case throw doubt on the validity of the conjecture.

Motivation
In 1952 Dyson [1] shocked the high energy physics community by declaring that, quite generally, the QED perturbation series cannot converge. Writing the series as F (e 2 ) = c 0 + c 2 e 2 + c 4 e 4 + . . . , (1.1) Dyson argues: "Suppose, if possible, that the series converges for some positive value of e 2 ; this implies that F (e 2 ) is an analytic function of e at e = 0. Then for sufficiently small values of e, F (−e 2 ) will also be a wellbehaved analytic function with a convergent power-series expansion". He then argues that, on physical grounds, this cannot be the case, since for e 2 < 0 the QED vacuum will be unstable due to a runaway production of e + e − pairs which coalesce into like-charge groups.
Only shortly later C.A. Hurst [2] already provided a mathematical proof of this fact for scalar λφ 3 theory. The proof is essentially based on the following three elements: 1. The use of the inequality to establish lower bounds for arbitrary Feynman diagrams (in the Euclidean).
2. Proof that the number of distinct Feynman diagrams at n th loop order grows like ( n 2 )!n! 3. Absence of sign cancellations between graphs.
In 1979, 't Hooft [3] found another very general, but very different, argument against convergence of the perturbation series based on renormalon chains. Thus today it is believed that the perturbation series in nontrivial quantum field theories generically is asymptotical rather than convergent, so that summation methods must be used. Of those by far the most important one is Borel summation, since it is ideally suited to the typical factorial growth of perturbation theory coefficients. Let me remind you that, for a factorially divergent series, one defines the Borel transform as If B(t) has no singularities on the positive real axis and does not increase too rapidly at infinity, one can also define the Borel integral F is the Borel sum of the original series F . F is asymptotic toF by construction, although the true physical quantity represented by the series F might still differ fromF by nonperturbative terms. The Borel transform remains a useful concept even when it leads to singularities, since those usually contain information on the large-order structure of the theory. In many cases they can be traced either to instantons, renormalons or Euclidean bounces. Until recently, there was a dearth of nontrivial examples for field theory models where sufficient information would be available to decide the question of Borel summability in a definite manner. Fortunately, this has changed through the advent of supersymmetry; in recent years Borel summability (or Borel non-summability) has been rigorously demonstrated in a number of supersymmetric models [4].
Even when Borel summability does not apply, Borel analysis can still be very useful through the use of Borel dispersion relations. This goes as follows. Assume that a function F (g) has an asymptotic series expansion where the expansion coefficients c n have the leading-order large n behaviour c n ∼ ρ n Γ(µn + ν) with some real constants ρ > 0, µ > 0 and ν. It is easy to see that such a series is not Borel-summable, since the Borel integral (1.5) can never converge (for example, in the textbook case µ = ν = 1 it has a pole at t = 1/ρ) . Nevertheless, applying a dispersion relation to this integral one can show that the leading contribution to its imaginary part for small g is given by Coming back to the case of QED, given the arguments by Dyson and 't Hooft it is certainly save to exclude a nonzero convergence radius of the full QED perturbation series. However, despite of the immense work that has gone into low-order perturbative QED computations, presently still little is known about the precise large-order behavior of the coefficients. Contrary to the case of scalar field theories mentioned above, straightforward estimates based on lower bounds for individual diagrams cannot be used in gauge theory, since here Feynman diagrams come with different signs, and gauge invariance is known to lead to cancellations between them. And these cancellations are particularly extensive in the abelian case, where there are no obstructing color factors. Thus QED in this respect is more difficult than QCD, which is made worse by the absence of (spacetime) instantons in QED, which in the nonabelian case can provide some large-order information. In 1977 Cvitanovic [5] suggested, based on an analysis of the calculation of the three-loop anomalous magnetic momentum g − 2 which he had done with Kinoshita [6], that these cancellations should be taken into account by counting the number of classes of gauge-invariant diagrams, rather than the number of individual diagrams. He also conjectured that, for the case of g − 2, they reduce the coefficients of the perturbation series sufficiently to make it convergent in the quenched approximation. This conjecture, although nowadays forgotten, is actually still standing, since neither Dyson's nor 't Hooft's arguments work in the absence of fermionic bubbles.
Here I will summarize the state of a long-term effort [9,10,11,12,13,14,15,16,17] to get information on the high-order behaviour of the QED perturbation series using the Euler-Heisenberg Lagrangian and its higherloop radiative corrections.

The 1-loop Euler-Heisenberg Lagrangian.
The Euler-Heisenberg Lagrangian ("EHL") is the one-loop QED effective Lagrangian for a constant external field. Euler and Heisenberg [7] obtained for it in 1936 the following well-known proper-time representation: Here a, b are the two invariants of the Maxwell field, related to E, B by The analogous result for Scalar QED was obtained by Weisskopf [8], but will be called "Scalar Euler-Heisenberg Lagrangian" in the following. 3) The Euler-Heisenberg Lagrangian ("EHL") holds information on the oneloop N photon amplitudes, but only in the low energy limit (since a constant field can emit only zero energy photons). Thus diagrammatically L (1) (F ) is equivalent to the sum of the Feynman graphs shown in Fig. 1, where all photon energies are small compared to the electron mass, ω i m.
This formula (called 'AAM formula' in the following) is highly remarkable for various reasons. Despite of its simplicity it is a true all-loop result; the rhs receives contributions from an infinite set of Feynman diagrams of arbitrary loop order, as sketched in fig. 1   In [12] it was shown how to carry out explicitly the construction of these limiting low-energy amplitudes starting from the weak field expansion of the EHL, It turned out that, if one fixes the number of photons, their momenta k 1 , . . . , k N and a helicity assignment for each photon, then in this limit all the dependence on momenta and polarizations is carried by a unique invariant. Thus the magnitude of the amplitude can be specified by a single number, which will be essential for our whole approach. Except for the purely magnetic case, the EHL has also an imaginary part related to vacuum pair creation by the electric field component (to be called "Sauter-Schwinger pair creation" in the following) [18,19]. In the purely electric case one finds, from the poles in the T -integration, the following decomposition due to Schwinger [19], . In the following we will focus on the weak field limit β 1 where only the first of these "Schwinger exponentials" is relevant.
The nonperturbative dependence of the Schwinger exponentials on the field supports the interpretation of field-induced pair creation as a vacuum tunneling effect, as proposed by Sauter as early as 1931 [18].
As usual in quantum field theory, the real and imaginary parts of the EHL are related by a dispersion relation. For the N -photon amplitudes at full momentum, this would be a standard dispersion relation performed diagram-by-diagram, relating the diagrams of Fig. 1 to the "cut diagrams" shown in Fig. 2, involving on-shell electrons. However, in the zero-energy limit the cut diagrams all vanish, since a finite number of zero-energy photons cannot create a pair on-shell. Thus what counts here is only the asymptotic behavior for a large number of photons, and instead of an ordinary dispersion relation we have to use a Borel dispersion relation. This works in the following way [9].
Consider the purely magnetic EHL. Expanding it out in powers of the field yields with an effective expansion parameter g = eB m 2

2
, and coefficients c (1) n that can be written in terms of the Bernoulli numbers B n : .
Here c (1) n holds information on the N = 2n photon amplitudes. The asymptotic behavior of the coefficients can be easily studied using well-known properties of the Bernoulli numbers. One finds Thanks to the factor (−1) n , the individual terms on the right hand side of (2.8) all give convergent Borel integrals. This is one (rather roundabout) way of seeing that the purely magnetic EHL has no imaginary part, and does not pair create. The analogous expansion for the purely electric field case is almost the same: , but with the same c (1) n . However, the additional factor (−1) n makes the series non-alternating, which is crucial, because now the termwise use of the expansion (2.8) leads to divergent Borel integrals. These divergent integrals do, however, all possess well-defined imaginary parts, by a (now ordinary) dispersion relation. One finds a perfect match between the expansion (2.8) and Schwinger's expansion (2.5): The Euler-Heisenberg Lagrangian at higher loops.
Proceeding to higher loop orders, the two-loop EHL L (2) is generated by the diagrams shown in Fig. 3 (here and in the following it is understood that internal photon corrections are put in all possible ways).
This formula (called 'AAM formula' in the following) is highly remarkable for various reasons. Despite of its simplicity it is a true all-loop result; the rhs receives contributions from an infinite set of Feynman diagrams of arbitrary loop order, as sketched in fig. 1. As in the one-loop case, L (2) will have an imaginary part iff the field is not purely magnetic. This imaginary part corresponds to the one-loop "cut diagrams" depicted in Fig. 4. The two-loop EHL was first studied by V.I. Ritus, both for Spinor [20] and Scalar QED [21]. These calculations, as well as later recalculations [22,23,14], resulted in a type of rather intractable two-parameter integrals. However, the first few coefficients of the weak-field expansions of the twoloop EHLs have been computed [23,9,14]. As to the imaginary parts, the Schwinger formulas (2.5) generalize to the two-loop level as follows [24]:

Number of external legs Number of loops
(these coefficients, also called c n , should not be confused with the weakfield expansion coefficients introduced above). Thus at two-loop the kth Schwinger-exponential appears with a prefactor which is still a function of the field strength, and of which presently only the lowest order terms in the weak-field expansion are known. Still, things become very simple at leading order in this expansion: Adding the one-loop and two-loop EHL's, one finds, e.g. for the spinor QED case [24], and this result is spin-independent (but for the normalization). In ( [24] )it was further noted that, if one assumes that in this weak-field approximation higher order corrections just lead to an exponentiation, then the e απ factor can be absorbed into the Schwinger factor e − π β by the following mass-shift, Moreover, the existence of this mass-shift can be independently confirmed in two different ways. First, the same mass shift had been found by Ritus already before in the crossed process of electron propagation in the electric field [25]. Second, in the tunneling picture it can be interpreted as the correction to the Schwinger pair creation rate due to the pair being created with a negative Coulomb interaction energy at a definite distance, taking the Coulomb interaction into account at the one-photon exchange level [24]. And, although this was not known to the authors of [24], an analogous exponentiation had already been conjectured two years before for the Scalar QED case by Affleck et al. [26]. However, those authors used a totally different approach based on worldline instantons. To explain this concept, we first have to discuss Feynman's worldline path integral representation of the QED effective action.

Worldline representation of the QED effective action
In 1950 Feynman presented, in an appendix to one of his groundbreaking papers on the modern, manifestly relativistic formalism of QED [27], also an alternative first-quantized formulation of scalar QED, "for its own interest as an alternative to the formulation of second quantization". There he provides a simple rule for constructing the complete scalar QED S-matrix by representing the scalar particles in terms of relativistic particle path integrals, and coupling them through photons in all possible ways. Upon restriction to the purely photonic part of the S-matrix (no external scalars), and to the "quenched" contribution (only one virtual scalar), this "worldline representation" can be stated very compactly in terms of the (quenched) effective action Γ[A]: (4.1) Here T denotes the proper-time of the scalar particle in the loop, m its mass, and x(T )=x(0) Dx(τ ) a path integral over all closed loops in spacetime with fixed periodicity in the proper-time. The worldline action S[x(τ )] has three parts, They are given by The kinetic term S 0 describes the free propagation of the scalar, S ext its interaction with the external field, and S int generates the corrections due to internal photon exchanges in the loop. Expanding out the two interaction exponentials leads back to Feynman diagrams, however with the important difference that no particular ordering of the photon legs along the loop needs to be fixed. Thus the term S ext alone upon expansion yields the diagrams of Fig. 1 (where each leg now stands for an interaction with the arbitrary field A(x)).
The "worldline instanton" of Affleck et al. [26] is an extremal trajectory of the worldline path integral for a stationary phase approximation. For the case of a constant electric field in the z direction this extremal action trajectory is given by a circle in the (euclidean) t − z plane: It can be shown that in the weak field (= large mass) limit the imaginary (although not the real) part of the effective Lagrangian can be wellapproximated by replacing the path integral with this single trajectory: This is easily evaluated to be (4.6) Thus the contribution of S 0 + S ext just reproduces the leading (one loop) Schwinger exponential of (2.5), and the one of S int the e απ factor. Thus Affleck et al. arrive, with very little effort, at the same exponentiation for Scalar QED that Lebedev and Ritus find in Spinor QED: Their argument assumes the field to be weak, but there is no restriction on the strength of the coupling α. We note: • Formula (4.7), if true, constitutes a rare case of an all-loop summation of an infinite series of graphs of arbitrary loop order. Those graphs are shown in Fig. 5.
• According to [26], the contribution of all non-quenched diagrams gets suppressed in the weak-field limit.
• Perhaps most surprisingly, the scalar mass appearing in (4.7) is already the physically renormalized one, implying that the worldline instanton approach automatically takes all mass renormalization counterdiagrams into account. This is remarkable considering that the determination of the physical mass parameter for the EHL becomes a rather nontrivial issue already at two-loops [20,22,23].
Thus according to Affleck et al. a true all-loop summation has produced the factor e απ , which is not only unreasonably simple, but also perfectly analytical in the fine structure constant α! According to what has been said above, this would seem to point towards extensive cancellations between Number of external legs Number of loops Moreover, the mass appearing in (1.15) is argued to be still the physical renormalized mass, which means that the above figure should strictly speaking include also the mass renormalization counter diagrams which appear in EHL calculations starting from two loops.
The derivation given in [33] is very simple, if formal. Based on a stationary path approximation of Feynman's worldline path integral representation [34] of L scal (E), it actually uses only a one-loop semiclassical trajectory, and arguments that this trajectory remains valid in the presence of virtual photon insertions. This also implies that non-quenched diagrams do not contribute in the limit (1.15), which is why we have shown only the quenched ones in fig. 1.
Although the derivation of (1.15) in [33] cannot be considered rigorous, an independent heuristic derivation of (1.15), as well as extension to the spinor QED case (with the same factor of e απ ) was given by Lebedev and Ritus [31] through the consideration of higher-order corrections to the pair creation energy in the vacuum tunneling picture. At the two-loop level, (1.15) and its spinor QED extension state that Feynman diagrams. However, neither [26] nor [24] made this point; perhaps, because the Schwinger pair creation rate is a rather peculiar quantity. Thus our next goal will be to transfer the exponential factor from the imaginary to the real part of the EHL, by a Borel dispersion relation, and from there to the low-energy photon S-matrix through the same procedure as at one-loop: First we need to see whether our one-loop Borel dispersion relation can be extended to the multiloop level. For this it will be useful to consider the simplest possible nonzero constant field background, which is the self-dual one.

The self-dual case
As we mentioned above, the two-loop correction to the EHL for a purely electric or purely magnetic field are known only in terms of intractable integrals, and only the first few coefficients have been calculated so far. However, this case is not the simplest one that one can consider; mathematically much better behaved is the one of a (Euclidean) constant self-dual field, defined by F µν = 1 2 ε µναβ F αβ . The field strength tensor can be written as At the one-loop level, the self-dual ("SD") EHLs for Scalar and Spinor QED are special cases of (2.1), (2.3). Anticipating the result of the two-loop calculation below, it will be useful here to eliminate the proper-time integral, and perform a change of variables from f to κ ≡ m 2 2ef . This leads to [10] L (1)(SD) scal where the function Ξ(x) is defined as follows: The spinor EHL in this SD case after renormalization differs from the scalar one only by a trivial global factor of −2 (the reason for the independence of spin is that the Dirac equation in such a background possesses a hidden supersymmetry [28]). Remarkably, for the SD case it is possible to do all integrals in closed form not only at one-loop, but even at two-loop, in both Scalar and Spinor QED. The results can be written compactly in terms of the digamma function ψ(x) ≡ Γ (x)/Γ(x):

(5.4)
Here κ ≡ m 2 2ef and (note that ξ(x) = Ξ (x)). Using the well-known expansion of the digamma function at x = ∞ in terms of the Bernoulli numbers, one finds the following closed-form expressions for the one-and two-loop weak-field expansion coefficients c (1,2)(SD) n (we write them down for the spinor case): , .
Further, in this self-dual case there is also an analogue of the distinction between a purely magnetic and a purely electric field. For f real the SD EHL turns out to have a weak-field expansion with alternating coefficients, so that it is Borel summable, and there is no imaginary part. Thus we call this case "magnetic-like". Taking f imaginary removes the alternating sign and creates a pole in the Borel integral, which implies an imaginary part for the EHL. Thus we call this case "electric-like". This imaginary part of the self-dual EHL with complex f is obtained from (5.4) simply by using the analytic continuation of the digamma function, and thus also known in closed form. Studying the self-dual case turned out to be useful in three ways: • For this case we could verify that the Borel dispersion relation (1.8) can be used to construct the imaginary part of the EHL from the weak-field expansion even at the two-loop level. That is, the asymptotic threeparameter matching (1.7) works, and (more nontrivially) the Borel summation procedure does not miss any non-perturbative terms (see [29] for a case where such a thing actually occurred, even at one-loop).
• The effective action for a self-dual field is unphysical, since such a field cannot be real in Minkowski space. Nevertheless, it still carries information on the physical photon amplitudes; the self-duality condition corresponds precisely to a projection on the 'all +' (or 'all −') photon amplitudes [30]. Thanks to the closed-form expressions (5.4), even at the two-loop level we are still able to write down a closed-form all-N expression for this particular polarization choice: 6 Synthesis: a conjecture for the photon S matrix We are now ready to state a conjecture for the N photon amplitudes at arbitrary loop level l [13]. Consider the l -loop correction to the purely electric EHL, and define its weak-field expansion coefficients by (note the change of convention with respect to (2.6)). Assuming that the AAM formula (4.7) holds, and that the Borel dispersion relation (1.8) works at each loop order, we can conclude that the leading asymptotic factorial growth rate must be the same at each loop order, namely ∼ Γ(2n − 2): And here c (l) ∞ relates to the leading Schwinger exponential at l loops, At two loops, the numerical calculations of [9] confirm this, but only if physical mass renormalization is used! For generic mass renormalization one finds instead a leading factorial behavior of Γ(2n), and it is only through a cancellation of this leading order term between the unrenormalized EHL and its mass renormalization counter term that this leading factorial behavior gets reduced to the same Γ(2n − 2) behavior as at one loop. At the l -loop level, it is still not difficult to establish the leading factorial growth of the weak field expansion coefficients before renormalization, which is Thus at higher loop orders the correctness of the AAM conjecture requires increasingly extensive cancellations in the mass renormalization process to cut the leading factorial growth all the way down to Γ(2n − 2). Now let us consider the ratio of the l -loop to the one-loop coefficients. Combining (6.2), (6.3) with the AAM conjecture (4.7), we find at any fixed loop order At this stage, let us switch to the self-dual case. This is not essential for our argumentation, but we prefer it for two reasons: First, more is known explicitly about the SD EHL; second, as mentioned above the SD EHL directly translates into one particular helicity component of the N -photon amplitude, the one with all helicities equal ("all +" or "all -"). Using our above rule for the conversion of the self-dual weak-field expansion coefficients into the "all +" photon amplitudes, which is independent of the loop order, we get the following statement for the "all +" amplitudes in the limit of large photon number N = 2n, Summing this relation over l we get Assuming sufficient uniformity in l of the convergence for N → ∞, one could now conclude that the amplitude must be analytic in α for some sufficiently large N . But analyticity of the complete amplitude can certainly be safely excluded by renormalons and other arguments. Therefore uniformity must fail, and it is easy to see how this comes about diagramatically. In Fig. 6 we show the diagrams contributing to the EHL up to four loops, not showing the external legs.
In the worldline-instanton based derivation of the AAM conjecture (4.7) , only quenched diagrams contribute to the weak-field limit of the imaginary part of the electric EHL, thus this must also be true for the leading asymptotic terms in the large N expansion of the weak-field expansion coefficients; non-quenched diagrams must get suppressed for N → ∞. However, the number of such diagrams is strongly growing with the loop order, so that the process of the suppression of the non-quenched contributions by increasing N should slow down with increasing l. This provides a good reason for uniformity to fail for the whole amplitude, but there is no obvious reason to expect such a non-uniformity if one stays inside the class of quenched diagrams from the beginning. This led G.V. Dunne and one of the authors in 2004 [13] to conjecture that perturbation theory converges for the QED photon amplitudes in the one electron-loop approximation. If true, this would imply enormous cancellations between Feynman diagrams, presumably due to gauge invariance.
Only afterwards we learned that, as mentioned in the beginning, Cvitanovic [5] had conjectured the analogous statement for the electron g − 2 factor.

Three predictions for the three-loop EHL
To either disprove or further corroborate this conjecture, a calculation of the EHL at the three-loop level is called for. We would like to see the following three things happen: 1. We should see the next term of the exponentiation: 2. At three loops there is already a non-quenched contribution, and it should be suppressed in the large N limit.
3. The convergence of c (3) (n) c (1) (n) should not be slower than the one of c (2) (n) when only quenched diagrams are taken.
However, a calculation of the three-loop EHL in D = 4 seems presently technically out of reach.

QED in 1+1 dimensions
The proper-time representation of the one-loop EHL is essentially independent of dimension. In 2006 M. Krasnansky [31] studied the Scalar EHL in various dimensions also at two loops, and found, in particular, the rather surprising fact that the Scalar EHL in 1 + 1 dimensions even at two-loop has a structure almost identical to the one of the self-dual Scalar EHL in 3+1 dimensions. Let us contrast the two cases: above we wrote down the self-dual field strength tensor for D = 4, We also gave the Scalar EHL for this background, In 2D the field strength tensor is F = 0 f −f 0 , and the two-loop Scalar QED EHL comes out as [31] Since higher-loop calculations are somewhat easier in two dimensions, this suggests to use the 2D case as a toy model for studying the AAM conjecture. An effort along these lines was started in [15], however switching from Scalar to Spinor QED. Here we derived an analogue of the AAM conjecture in 2D, also using the worldline instanton approach, and established the correspondences between the 4D and 2D cases shown in Table 1:  There are two essential differences. First, in 2D the fine structure constantα is not dimensionless. Thus the exponent of the AAM formula (rhs of third line) here involves also a factor of κ 2 , which in the formula for the asymptotic behavior of the weak-field expansion coefficients (rhs of fourth line) leads to a shift in the argument between the l -loop and the one-loop coefficients. Thus in 2D the leading asymptotic growth of the coefficients increases with increasing loop order, as it does in the 4D case before mass renormalization, and correspondingly it can be shown that the contributions to the EHL from mass renormalization are asymptotically subleading, and thus irrelevant for our purposes (although the fermion propagator in 2D does not have UV divergences, mass renormalization is still a quite nontrivial issue, see [32] and refs. therein). Presumably this relates to the fact that QED in 2D is confining.
In any case, all three of our three-loop predictions above have an analogue in the 2D case. In [15] we also obtained the following formulas for the oneand two-loop EHLs in 2D Spinor QED: where we have now abbreviated Comparing with the Scalar QED result (8.1), we see that the spinor QED one is significantly simpler, as it involves the digamma function only linearly. This is another surprise, since in 4D the Scalar and Spinor EHLs do not show structural differences.
Remarkably, the two-loop EHL can (up to an irrelevant constant) even be written in terms of derivatives of the one-loop EHL: From (8.2) we find for the one-and two-loop weak-field expansion coefficients From this we can, using properties of the Bernoulli numbers, easily show that This verifies the 2D AAM-like formula of Table 1 at the linearized level. to the AAM prediction (normalized such that the limit is unity) .
In Fig. 7 we show the convergence to the asymptotic limit, which is rather rapid.
Even in the 2D case, the calculation of the three-loop EHL turned out be a formidable task, and it is only very recently that we were able to obtain it in a form suitable for computing a sufficient number of the weak-field expansion coefficients [16,17].
At three loops, we have the three Feynman diagrams shown in Fig. (8). Here the solid lines denote the electron propagator in the constant field. A and B are quenched, C is non-quenched.

A B C
The last one is by far the easiest one, and it is straightforward to obtain for it the following compact integral representation: sinh z sinh z sinhẑ sinh z [sinh(z + z ) sinh(ẑ + z )] 2 × e −2κ(z+z +ẑ+z ) sinh z sinh z sinh(ẑ + z ) + sinhẑ sinh z sinh(z + z ) .

(8.4)
This representation turned out to be quite adequate for a numerical calculation of the first 9 weak-field expansion coefficients c C (n) of diagram C. In Fig. 9 we use these nine coefficients to show that this unquenched contribution is indeed asymptotically subleading. This settles point 2. of our wish list above! Diagrams A + B are much more difficult, but the use of the "traceless gauge" choice λ = 2 led to simplifications, and in particular to manifest IR finiteness term-by-term. We managed to compute the first coefficient analytically, and five more coefficients numerically. Using these to plot the ratio  we get Fig. 10.
Thus we are falling even below the asymptotic prediction!

Conclusions and Outlook
Let us summarize: • We have presented first results of a calculation of the three-loop 2D EHL. This is the first calculation of a three-loop effective Lagrangian in QED.
• Although so far we have been able to compute only six coefficients of the weak-field expansion (we should be able to obtain a few more) it seems already likely that the analogue of the AAM conjecture fails in 2D QED. This would throw also serious doubts on the validity of the 4D AAM conjecture.
• However, since the coefficient ratios fall below, rather than above, the asymptotic prediction, the riddle of the unreasonable smallness of loop corrections remains. Presumably the worldline instanton approach captures some valid information on large-scale cancellations between Feynman diagrams, but needs refinement beyond two loops.
• We have also made an effort to make the point that the QED photon amplitudes in the limit of low energy and large number of photons are very natural objects for a study of the asymptotic properties of the QED perturbation series.
• It should also have become clear that physical mass renormalization is essential for asymptotic estimates in QED! Unless mass renormalization is done physically, QED perturbation theory will break down already at the two-loop level, because the two-loop contribution to any helicity component of the N -photon amplitude will, at least in the low energy limit, dominate over the one-loop one for sufficiently large N . This implies, in particular, that approaches to the study of the high-order behavior of the QED perturbation series that are indifferent to the issue of physical mass renormalization ought to be viewed with great caution.
• As an aside, it would be interesting to study also the QCD N -gluon amplitudes for large N from the point of view of mass renormalization.