Asymptotic structure of electrodynamics revisited

We point out that recently published analyses of null and timelike infinity and long-range structures in electrodynamics to large extent rediscover results present in the literature. At the same time, some of the conclusions these recent works put forward may prove controversial. In view of these facts, we find it desirable to revisit the analysis taken up more than two decades ago, starting from earlier works on null infinity by other authors.


Introduction
In a series of recent articles-see [24] and [8] and works cited there-a group of authors report on new conservation laws and symmetries both in classical and quantum electrodynamics. These structures are supposed to be encoded in the asymptotic properties of electrodynamics. Later, an analog of this structure in gravitation theory was used to analyze anew the black hole information paradox [17]. Also, in the context set by these investigations there is a noticeable interest in the literature in the so-called memory effects, both in gravitation, as well as electrodynamics (see [24,28] and works cited there).
B Andrzej Herdegen herdegen@th.if.uj.edu.pl 1 Institute of Physics, Jagiellonian University, ul. S. Łojasiewicza 11,Poland In this article we want to comment on the electromagnetic part of this programme in view of results known from the literature. We shall point out that the structure of asymptotic electrodynamics, along lines similar to those proposed by the authors mentioned above, has received more extensive attention in the past. On the other hand, we shall also point out that some claims in their programme may prove controversial. We use this opportunity to summarize the programme taken up by the present author more than two decades ago.
We sketch our main points and the plan of the article.
(i) On the classical level, a more extensive analysis of the asymptotic structure of electrodynamics is proposed in the article [18], which extends many ideas present earlier in the literature, in particular the null infinity analysis by Bramson [4] and Ashtekar and Streubel [2]. The null infinity 'matching conditions' (so called by Strominger et al. [24]) were known before, and in [18] they are derived in full generality. After recalling an important tool in the field of the wave equation in Sect. 2, we discuss the classical electrodynamics, and the asymptotics of fields in Sects. 3, 4 and 5. Invariant structures in the asymptotic regions are introduced in Sect. 6. (ii) Our discussion consequently uses Lorenz gauge potentials in electrodynamics.
These potentials have well-defined dynamics and null asymptotes appropriate for the asymptotic structure. In Sect. 7, still on the classical level, we ask what other gauges are admitted by this structure. We point out that although the gauge transformation called by the authors of [24] and [8] (and others) 'large gauge transformation'-LGT in the following-produces potentials with asymptotes admitted by the structure, it is not a symmetry of this structure. Therefore, the changes induced by LGT may be incorporated within this structure, but with interpretation wholly different from a gauge transformation: they modify the electromagnetic fields of the state. (iii) Commenting in Sect. 8 on the relations between long-range structure and the so-called memory effects, we recall that the earliest effect of this type in electrodynamics was discovered long ago by Staruszkiewicz [31]. What concerns us here, Staruszkiewicz's effect is a clear argument in favor of the observability of the long-range degrees of freedom in electrodynamics. (iv) The classical asymptotic structure was quantized and analyzed extensively in a series of papers [19][20][21][22]. 'Asymptotic quantization' by Ashtekar [1] may be viewed as the reference starting point of this analysis, but the discussion goes much further toward asymptotic algebraic formulation of matter-radiation system. The authors mentioned at the beginning, as it seems, may have motivations similar to my own: they treat variables at spatial infinity as genuine quantum observables. We make general remarks on specific infrared problems of quantum electrodynamics in Sect. 9 and then summarize the construction of our asymptotic quantum theory in Sects. 10, 11 and 12. (v) We construct a unitary operator built from the long-range quantum variables at spacelike infinity in Sect. 13 . We indicate that the construction of such an operator demands an explicit use of a concrete representation of the theory and, within that representation, a subtle limiting. We show that this unitary operator induces a nontrivial transformation of the basic quantum variables of the asymptotic algebra. This transformation contains elements corresponding to the classical effect induced by LGT. As in the classical case, the changes modify the physics of the state, and their interpretation within the limits of the asymptotic theory is far from a gauge symmetry. We show in Appendix D that our rigorously constructed operator corresponds, in loose terms, to the 'LGT generator' of [24] and [8] (which is not rigorously defined). (vi) However, if a scattering matrix may be defined, then it is true that the operators constructed by smearing of long-range variables should be conserved. This is a quantum version of a 'continuous conservation' law for long-range tails in classical electrodynamics. We mention earlier statements on this problem in Sect. 14. (vii) Appendix contains some more technical material needed in the main text.
The article is based mainly on the articles [18][19][20][21][22]. We give a specific reference only in some special cases.

Invariant measure on the set of lightlike directions
The first fact worth wider knowledge is the existence of a Lorentz-invariant measure on the set of lightlike directions. The existence of this measure is a classical result [14,29,31,35], but-strangely enough-ignored by most physicists despite its crucial role in the field of the wave equation and massless fields.
We consider the flat spacetime with a fixed origin, thus described by the Minkowski vector space M with the signature (+, −, −, −). We shall write l for any future-pointing lightlike (nonzero) vector and denote We shall write t for any future-pointing timelike, unit vector and denote (a unit sphere in the hyperplane t · x = 1). It is well known that the measure dμ 0 (l) = d 3 l/l 0 on the lightcone is Lorentz invariant. However, equally important, but less known, is the following construction. Let f (l) be a measurable function on C + , homogeneous of degree −2: where dΩ t (l) is the angle measure on the unit sphere, does not depend on the choice of the vector t (Lorentz invariance 1 ). We give a simple proof of this fundamental fact in Appendix A. It follows that if we denote (a, b, etc. are spacetime indices 2 ) -the generators of the Lorentz transformations, intrinsic differential operators on the lightcone, then

Homogeneous Maxwell equations
Let A(x) be a Lorenz gauge potential (we often suppress spacetime indices) of a free electromagnetic field where δ is the Dirac measure, ε(ω) is the sign of ω, and where the relations reflect reality of A, the Lorenz condition and a possible gauge transformation respectively, defines a very wide class (also possibly distributional) of such fields, and one needs a physically motivated selection criterion. An obvious restriction is to consider fields with finite energy-momentum, which is given by This is still a very large class and its often considered subclass-having mathematical advantages-consists of the so-called regular wave packets-solutions with initial data (on a Cauchy hyperplane) having smooth Fourier transforms with compact support not containing the zero vector. However, this class does not include fields of the Coulomb-like decay at spatial infinity-such fields are produced in scattering of charged particles. We describe further selection in three steps.

Step 1: Coulomb-like spacial decay; conservation of spacelike tails
The selection criterion which admits fields produced in scattering events, but at the same time leaves more infrared singular solutions outside, is this: for y 2 < 0 the potential A has a well-defined scaling limit A as (y) := lim λ→∞ λA(λy), y 2 < 0.
In terms of the Fourier transform this is equivalent to the existence of the limit note that a as (l) is homogeneous of degree −1. This type of singularity of a(k) is consistent with the finite energy demand, 3 and we further consider fields satisfying both restrictions. For fields satisfying our selection criterion not only the scaling limit centered at zero is well defined, but also for each x ∈ M and y ∈ M, y 2 < 0, there is This may be viewed, if one wishes, as a 'continuous conservation law': the form of the limit does not depend on the choice of the reference point x in Minkowski space. We note that: This 'conservation law' has nothing to do with any gauge symmetry-the field A may be of any algebraic type-and is only a property of the wave equation together with the selection criterion (3).
Note also that the two parts on the rhs have definite parities with respect to the reflection y → −y: the first part is odd, while the second is even. 4 Now, all standard free fields taking part in scattering processes in electrodynamics are of the second type: the spacelike tails of their Lorenz potentials are even (and the tails of the fields themselves are correspondingly odd). For example, the potential of the radiation field (i.e. A ret − A adv ) of a point particle with charge q scattered instantaneously in x = 0 is given by where v and u are the initial and the final velocity, respectively (the singularity at the cone x 2 = 0 is an artifact due to the instantaneous nature of the scattering event). The same is true for radiation fields produced by sources due to massive fields (as in the Maxwell-Dirac system). This absence of an odd part may be used as a further restricting criterion, but before doing this we want to give another form to the above formulae.
Let l ∈ C + and denoteV where overdot denotes differentiation with respect to the real variable s, and the function V will be obtained by integration. It is easy to see thatV is a real function homogeneous of degree −2, and-in the case of a vector function representing a Lorenz potential-transversal: Using this function in the Fourier representation of A(x), one obtains another useful integral representation of the solution of the wave equation whose advantage is that the integration goes over a compact set. If A is a Lorenz potential of an electromagnetic field, then the field is expressed by Now, our selection criterion-the existence of a as (k), Eq. (3)-implies that ω Re a(ωl) is continuous in ω at zero, while ω Im a(ωl) has a jump at zero of the magnitude 2 Im a as (l). Correspondingly, the behavior ofV (s, l) for large |s| is governed by 5

Step 2: null asymptotes
The null asymptotic behavior of the field, which we now want to investigate, depends crucially on this estimate. If Im a as (l) = 0, then R A(x ± Rl) is of order log R for large R. Moreover, one can show that also angular momentum radiated over finite time into any solid angle is not well defined for R → ∞. This is another reason to apply the more restrictive selection criterion mentioned before. Thus, from now on we assume that Im a as (l) = 0. (10) In that case the estimate (9) implies that the functionV (s, l) is absolutely integrable over s ∈ R and we fix V (s, l) by demanding that it vanishes for s → ∞: then for s → −∞ it has a well-defined limit obtained from the inversion of formula (6): We also define and then With these definitions, the null asymptotic behavior of the potentials is easily expressed [18]: and the null asymptotes of the electromagnetic fields are then: It is evident that the relations between V and V are gauge independent, and quantities l ∧ V (s, l) and l ∧ V (s, l) are gauge invariant. The real parameter s appearing in our formulae has a twofold interpretation. On the one hand, if in the above asymptotic formulas one sets x = st (t as defined in Sect. 2) and scales l's to t ·l = 1, then s is the retarded or advanced time, according to the case. On the other hand, if one compactifies the spacetime a la Penrose, then s is an affine parameter along generators of future or past null infinity. In particular, V (+∞, l) and V (−∞, l) are the values of the respective functions on the future (past) edge of the past (future) null infinity respectively, and their equality (14) is a free field version of a general 'matching property' to be discussed below. Fields with V (−∞, l) = 0 are called infrared-regular, and for them the spacelike tail vanishes, as seen from formula (4). Fields with V (−∞, l) = 0 are called infrared-singular.
At this point we want to stress once more a crucial fact concerning the definition of V (s, l). As indicated above, this function is obtained by integratingV (s, l) over s, so it may appear one has a freedom of transformation V (s, l) → V (s, l) + V 0 (l), with V 0 constant in s. However, A(x) is uniquely determined byV (s, l) alone according to formula (7), so its null asymptotes (past and future) are also unique. It is only with the specification (11) [and (12) for the past null asymptote, with (13) as a consequence] that the asymptotic formulae (15) are true.
We exploit this knowledge to characterize more fully possible functions of gauge transformations. Suppose that two Lorenz potentials differ by a gauge transformation, (1) (+∞, l) = 0, thus also α(+∞, l) = 0. According to formula (7), we have where ε + is a constant number. We set here x = st ± Rl and take the limit R → ∞.
As l · l > 0 almost everywhere, we find We summarize: Among Lorenz gauges there can be no gauge transformation ∂Λ whose asymptote lα(s, l) remains constant along null infinity generators (i.e., is independent of s). For any admissible gauge function Λ, the null asymptotic limits (16) are two, in general different, constant numbers. 6 To close this section, we would like to illustrate the null asymptotic behavior and comment on the possibility of asymptotic expansions of the form 7 A(st + Rl) = ∞ n=1 R −n A (n) (s, l), in the simplest possible case of a scalar field A, spherically symmetric in a chosen frame. Let A be given by (7) withV (s, l) = (t ·l) −2Ẇ ((t ·l) −1 s), where t is a fixed time axis vector andẆ (τ ) is any regular function decaying at least like |τ | −1− for |τ | → ∞. This field satisfies all our selection criteria, and the explicit integration in formula (7) gives One can now immediately explicitly confirm all our statements on the null asymptotes, Eqs. (12)-(15) and subsequent comments. Let us choose two examples: where α, β and τ 0 are constants. With these choices, both fields A 1 and A 2 fall into the class specified in this section; the first one is infrared-regular, while the second is infrared-singular. However, for < 1 the asymptotic expansion of these fields does not exist in the form of the series in 1/R beyond the leading order. In general, this problem is even more involved, but for the electromagnetic fields due to the potentials defined in the present section the two leading orders do exist.

Step 3: fields of electric type
There is one final step needed to complete our selection criterion. The fact that there are no magnetic monopoles and free fields are thus produced by scattering electric charges is encoded in the property of the spatial tail of the field being of electrical type. This is equivalent, as it turns out, to the condition Smooth functions, homogeneous of degree −1, transversal: l · V (−∞, l) = 0, and satisfying the above condition form a linear space, whose properties are described in Appendix C. It follows that where Φ(l) is a homogeneous function given by Further properties of the relation of V with Φ, which will be used later, are to be found in Appendix C. We end this section with a simple illustration of the null structures for the potential of the scattered particle (5). In this case V (s, l) = θ(−s)V (−∞, l) (θ is the Heaviside step function; the singularity at s = 0 is again due to unrealistic instantaneous

Matter-radiation system: null and spacelike infinity
Let now the inhomogeneous Maxwell equations be a part of a closed electrodynamic system. We again assume Lorenz gauge, thus in Gaussian units, which we use in this paper. 8 We consider a scattering situation; thus the conserved current J (x) describes sources which stabilize in remote past and future. This stabilization, more specifically, means that for v on the future unit hyperboloid, the leading behavior of the current is with some scalar functions ρ ± (v) on the hyperboloid. This is precisely true for asymptotically freely moving classical particles, but also, up to oscillating terms not contributing to the leading behavior of A, for matter fields like massive Dirac field. Also, we assume that J (x) vanishes sufficiently fast in spacelike directions. It follows then that the function defined by where δ is the Dirac measure, is well defined, bounded and for s tending to ±∞ has finite limits the integral over the unit hyperboloid. These limits are fully determined by the asymptotic forms of the current in the future/past. The conserved charge carried by the current is Moreover, due to (19) one has The reason to define this function is that the null asymptotes of the retarded and advanced fields are then given by It also follows that the radiation field A rad = A ret − A adv is the free field given by Formula (7) in which one should substituteV J forV , and it satisfies all our selection criteria.
The total potential is decomposed in two ways appropriate for incoming and outgoing characterization, A = A ret + A in = A adv + A out , respectively, where A in and A out are free fields. It is now evident that if the incoming field A in satisfies our selection criteria, then also the outgoing field A out does. We again denote by V out and V in the future asymptotes, and by V out and V in the past asymptotes of A out and A in , respectively. Then it follows that the null asymptotes of the total potential and the total field are given by where Putting here s = ±∞, one finds that so these characteristics are totally due to the outgoing/incoming currents, respectively, while so these characteristics are the sums of the source and the free infrared singular contributions. Moreover, adding the sides of equations (21) The second equality in this formula constitutes a matching property of the values at the past edge of the future null infinity, and the future edge of the past null infinity (see [18], Eq. (2.26) and the following discussion). We stress once more that to this matching equality contribute both retarded/advanced as well as free 'in'/'out' fields. 9 The spacelike asymptotic behavior is also governed by this characteristic, for y 2 < 0: where δ is the derivative of the Dirac delta distribution. 9 It is not true that 'finite energy wave packets die off at i 0 ' as the authors of [24] put it.
At this point, an observation on gauge of the total potential is due. As we work in a scheme which assumes Lorenz gauge, the gauge transformation A → A + ∂Λ demands Λ = 0. Therefore, ∂Λ is a free potential satisfying all our selection criteria. In particular, if lα(s, l) is its future null asymptote, then α(+∞, l) = 0, as in the free case with no sources (and similarly α (−∞, l) = 0 for past asymptote). Thus, the limit values V (+∞, l) and V (−∞, l) are gauge independent (in the Lorenz gauge class). 10 We end this section with an illustration of the 'matching property' (22). Consider a scattering event in which there are n incoming particles and n outgoing particles, with charges and four-velocities given by q i , v i and q j , v j , respectively. It is easy to show that for a single free particle with charge q and velocity v, one has V J (s, l) = q v/v ·l. Therefore, in our scattering event we have Then, Eq. (22) gives This relation is a clear announcement of the problems to come in quantum theory: suppose the space of incoming states does not contain infrared-singular 'in' fields, so the 'in' fields are represented in the usual Fock space of photons. However, this is only possible if the profiles of incoming photon states satisfy lim ω→0 ωa in (ωl) = 0. On the other hand, the rhs of Eq. (25) cannot vanish for all scattering amplitudes (depending on velocities v i and v i ) in the given representation. The relation tells us then that lim ω→0 ωa out (ωl) = 0, which contradicts the possibility of representing the 'out' profiles in Fock space. We shall come later to existing strategies for the resolution of this difficulty.

Matter-radiation system: timelike infinity
We start by recalling the Fourier representation of the solution of the free Dirac equation. A convenient version of this integral representation is: where dμ(v) = d 3 v/v 0 is the standard invariant measure on the unit future hyperboloid and f (v) is a Dirac spinor-valued function on this hyperboloid. A more standardlooking form is obtained by noting that . While the free electromagnetic field is fully encoded in its null asymptote, the free Dirac field is fully encoded in its timelike asymptote, which is Remember that the free field has no gauge freedom of the second kind, so ψ and f are unique up to a constant phase.
Suppose now that the current J appearing in the last section is produced by the Dirac electron-positron field, and the fields form a closed theory with the Dirac equation added. It is well known that the full Dirac field in standard Lorenz gauge does not approach then in simple manner any free field for time tending to plus/minus infinity, its leading asymptotic term containing electromagnetic contributions. The standard technique usually employed for handling this problem was proposed by Dollard [12] and in quantum electrodynamics was applied by Kulish and Faddeev [26]. However, we find it more convenient, for reasons to become clear below, to use another procedure.
We introduce a new gauge of the total electromagnetic potential: if A is the Lorenz potential discussed in the preceding sections, then where S is any scalar function assumed to satisfy For x = ±λv, the leading asymptotic terms of A(±λv) for λ → ∞ are of order 1/λ. It follows then that in this limit v · A tr (±λv) = O(λ −1− ). Now we can use one of the results of [18], which says the following. In the gauge defined above, the full Dirac field again has the asymptotic behavior given by (27), with some spinor functions f ∓ for 'in' and 'out' asymptotes, respectively. There are no electromagnetic field modifications of these asymptotic terms, and one could define free incoming/outgoing fields by plugging f ∓ into the integral representation (26). We indeed intend to perform such operation, but with one additional modification explained in the next section.
We note an important fact that there is no gauge freedom of the second kind in the definition of asymptotes f ∓ . Indeed, one can show that f ∓ do not depend either on the choice of particular Lorenz gauge of A, or of the specific function S satisfying (28).

Invariant structures and conserved quantities
The causal 'in' and 'out' infinities are equipped with two Poincaré invariant structures: symplectic form at the null infinity and scalar product in the timelike infinity. In the following we restrict attention to the 'out' infinity; similar structures exist in the 'in' case. Namely, the Poincaré transformations acting by leave invariant the symplectic form and the pre-Hilbert scalar product Here For a pair of free, infrared-regular electromagnetic fields, one has the standard symplectic form defined by the integral over any Cauchy surface by {F 1 , . This form, when expressed in terms of null asymptotes, is equal to the form (29). As can be easily seen, the integral (29) extends to all free fields in the class considered in this article (it is absolutely convergent also for infrared singular fields 11 ). Note that for free fields the form (29) is gauge invariant: it does not change under the transformation The symplectic form (29) extends also to the asymptotes of total fields, with the integral staying absolutely convergent. In this case the gauge transformation causes the following change of this quantity: l) is the charge of the field V i . This addition will play a role in the quantum case.
Let us return to the Poincaré transformations of the asymptotic structures. The generators of these transformations defined by , l), acts intrinsically in the hyperboloid. Consider now the energy-momentum and 4-angular momentum going out into timelike and null infinity, carried by massive and electromagnetic fields, respectively. One finds that these quantities are elegantly expressed by where V (s, l) is the asymptote of the total electromagnetic potential, and f + (v) is the asymptotic variable of the Dirac field as defined at the end of Sect. 5. All four expressions on the rhs of these equations are gauge invariant: the scalar product parts because of the second-kind gauge invariance of f + , and the invariance of the symplectic form parts is proved by an easy calculation using the form of generators and identity (31). If analogous incoming quantities P in a and M in ab are formed, then one finds the expected conservation laws P in = P out and M in = M out . In case of the energymomentum, this is also equal to the integral of the total energy-momentum density over any Cauchy surface, but in case of the angular momentum, a similar statement needs a comment. The 1/R 2 spacelike tail of the electromagnetic field has the consequence that the density of the angular momentum is not absolutely integrable over a Cauchy surface. However, the oddness of the tail of electromagnetic field implies that this density is also asymptotically odd, which allows one to regularize the integral over spacelike hyperplane by taking a limit of integrals over balls of finite radius. The quantity so regularized is equal to the incoming and the outgoing angular momentum. As an aside, we note the following interesting fact. In a theory in which both electric and magnetic charges are present, the analogous statements are not true, namely M in = M out -the angular momentum leaks into the spacelike infinity [18].
If we now decompose V (s, l) = V (+∞, l) + V out (s, l), then {V, r a V } = {V out , r a V out }-the electromagnetic energy-momentum is expressed in terms of the free 'out' field. However, the angular momentum turns out to contain a mixed adv-out term: 1 2 {V, n ab V } = 1 2 V out , n ab V out + V out , n ab V (+∞, .) .
Using the asymptotic form of the Dirac field, one finds where is the asymptotic density of particles moving with velocity v and V e a (v, l) = ev a /v · l is the null asymptote of the Lorenz potential of the Coulomb field surrounding the particle with charge e moving with constant velocity v.
The latter form of this term allows its absorption into the matter part by a phase transformation. We introduce (the second form following from {V e (v, .), V (+∞, .)} = 0) and then the conserved quantities take the form These expressions look formally as sums of independent free fields contributions. Therefore, we identify the Dirac out-field by placing g + in place of f in the integral representation (26). However, this Dirac field must then describe massive particles together with their Coulomb fields. This will be confirmed in quantization. We end this section by giving another form to the f → g transformation: if one uses the definition (18), then the identity (53) applied to the pair V out (−∞, l), Φ out (l) implies

Other gauges; LGT is not a gauge symmetry of the asymptotic structure
We have thus completed the discussion of the asymptotic structure of the classical electrodynamics. With this picture before our eyes, let us now once more return to the question of the gauge symmetry of this structure. In our constructions we have assumed the class of the Lorenz gauges. This class is distinguished by a well-defined dynamics and allows for gauge transformations within the class, which are encoded in the null asymptotes by V 2 (s, l) = V 1 (s, l) + lα(s, l). The question we want to discuss now is the following: are there other/wider classes of gauges for which this structure is preserved?
The class of electromagnetic fields F ab (x) considered here is characterized by the existence of the Lorenz gauge potentials A b (x) with well-defined null asymptotes V b (s, l). We now want to admit other gaugesÂ = A + ∂Λ, but we keep the restriction demanding that the potentials have well-defined null asymptotes, i.e. for eachÂ a (x) there is lim -the denominator t·l guarantees that the asymptoteV b (s, l) has the same homogeneity property as V b (s, l):V (λs, λl) = λ −1V (s, l), λ > 0 (see the discussion of the variables R, s, l in Appendix B). It follows that for the gauge function Λ(x) the limit has to exist. This condition restricts the gauge functions to those with the asymptotic expansion with the sufficiently regular rest, so that the leading term of ∂Λ is obtained by differentiation of the two leading explicitly written terms; note that ε + (l) is homogeneous of degree 0, and β t (s, l) is homogeneous of degree −1.
Let V + b (l) be a vector function entering with ε + (l) into the relation described in Appendix C for the pair V and Φ V , formulae (50)-(53). Using these relations, and applying the formula for differentiation, Eq. (49) in Appendix B, we obtain The form of the expansions (36) and (37) remains valid if the time axis vector t is replaced by another such vector t . However, while the function ε + (l) remains unchanged (as the notation suggests), the function β t undergoes the following transformation: -this is shown by writing (s /t · l ) + R l = (s/t · l) + Rl and expressing R , s , l in terms of R, s, l in leading orders in 1/R. This shows that β t (s, l) − st · V + (l)/t · l = γ (s, l) does not depend on t, and with this notation the most general transformation of the null asymptote of the potential is given byV b (s, l) = V b (s, l) + V + b (l) + l bγ (s, l). However, if |γ (s, l)| = O(|s| −1− ) for |s| → ∞, which is needed forV (s, l) to satisfy the demands of the asymptotic structure, then the term l bγ (s, l) is just the gauge transformation previously considered in the Lorenz class, 12 so this addition brings nothing new. Therefore, we set γ (s, l) = 0 and we are left with the transformation This is, in a somewhat different guise, the 'large gauge transformation' considered by the authors of [8,24], and in numerous other recent articles.
We now look at the transformation (38) from the intrinsic point of view of the asymptotic structure as described in preceding sections. Does the new asymptotê V (s, l) fit into the picture? It does, but with interpretation wholly different from the gauge transformation which led to the relation (38).
First of all, the addition V + (l) to the asymptote of the potential, if it results from the gauge transformation as discussed above, does not have the meaning of the initial data at infinity as before, as no equation is given which propagates these data into the bulk of the Minkowski space. Any rigid choice of the bulk functions is rather a transformation into a different sector, than a symmetry.
Secondly, and even more importantly, the transformation (38) is not a symmetry of the structure. The symplectic structure is changed under this transformation: if l). This addition to the symplectic form does not vanish even in the free field case (characterized by V i (+∞, l) = 0). Moreover, the electromagnetic field contribution to the angular momentum (32) changes: 1 2 V , n abV = 1 2 V, n ab V + V, n ab V + .
On the other hand, the transformation (38) may be interpreted within the asymptotic structure of the preceding sections, but it then changes the physics of the system. Namely, if we decompose V (s, l) into the sum of the advanced part of the outgoing current (see the discussion in Sect. 4) and the free outgoing field V out (s, s, l). The transformation (38) does not change the asymptote of the potential of the free 'out' field. On the other hand, we show in Appendix C that for almost all V + (the precise sense of this statement is to be found there) there exists the representation where ρ(v) is a smooth function of compact support on the unit hyperboloid, and such that ρ(v) dμ(v) = 0. This uniquely induces the interpretation of the transformation (38) within the limits of the asymptotic structure: The transformation (38) adds the advanced electromagnetic field of an external, asymptotic, charge-free current J + (x) = xρ(x), ρ(x) homogeneous of degree −4.
The current J + is not unique: the function ρ(v) depends on three real variables, while the asymptote V + (l) is a function of two independent real variables. As we shall see, the above interpretation will also be confirmed in the quantum case.

Electromagnetic memory
The long-range electromagnetic variables play a substantial role, both in the classical structure described above and in the quantum theory to be developed below. We stop here to ask what is the observational status of such degrees of freedom.
In this connection, it is appropriate to invoke certain analogy between electromagnetic and gravitational waves. Within the limits of linearized gravity, it has been pointed out in 1974 by Zel'dovitch and Polnarev [34] that if a pair of test bodies is exposed to a gravitational wave burst, then their relative separation undergoes a finite and permanent shift between remote past, when the burst has not arrived yet, and remote future, when it has passed away. This effect was discussed again and named the 'memory' of the gravitational wave burst in 1985 [3]. The size of this effect depends on the initial and final masses and velocities of the bodies which produce gravitational wave burst (see, for example, [33]). In 1991, it was shown by Christodoulou [11] that a substantial, and potentially more prominent, contribution to the displacement of the test bodies is caused by nonlinear effects of gravitation.
In the meantime, in 1981, Staruszkiewicz has discussed an electromagnetic effect which may also be interpreted as a kind of memory of wave packets [31]. To my best knowledge, this was the first observation of this kind in electrodynamics. Moreover, Staruszkiewicz's effect shows in the most clear way the connection between 'memory' and infrared degrees of freedom. We describe the effect in our notation. Let A(x) be a free Lorenz potential with the asymptote V (s, l), thus given by formula (7). The scaled potential A λ (x) = λ −1 A(x/λ), λ > 0, with the asymptote V λ (s, l) = V (s/λ, l), has the same spacelike tail (as it depends on V λ (−∞, l) = V (−∞, l), see (23)), but with λ → ∞ the energy content of the scaled electromagnetic field tends to zero; we call this large λ regime the infrared limit. In this limit the field is too weak to change the momentum of a charged test particle. However, as shown by Staruszkiewicz in the quasi-classical approximation, the phase of the wave function of the particle undergoes, between times minus and plus infinity, the momentum-dependent phase shift where q is the charge of the particle and p its momentum. For particle wave packets this produces interference effects, which result in an adiabatic finite shift of the wave packet of the particle. The shift is solely due to the long-range degrees of freedom of the electromagnetic field; the field itself may have a negligible energy content. Later this effect was confirmed in more elevated settings in [18] and [23] (in the latter reference a quantum Dirac field in the external infrared limit field was considered). 13 The Staruszkiewicz effect may be simply explained for a classical test particle with charge q and mass m placed in the free external electromagnetic field F ab (x). The equation of motion for the trajectory z a (τ ): where τ is the proper time and the overdot denotes differentiation with respect to τ , may be written in the equivalent form For τ tending to ±∞, the integrals Δ a ± = −(q/m) ±∞ 0 F a b (z(σ ))ż b (σ ) σ dσ give the vectors separating the 'out'/'in' asymptotic free trajectories (straight lines) from the point z a (0). For scaled F, in the infrared limit we haveż(τ ) v, a constant four-velocity, and then the 'out' asymptotic trajectory is parallel to the 'in' asymptotic trajectory, but shifted by Δ a = a + − Δ a − . Substituting in these integrals z(σ ) z(0) + vσ , using the representation (8), and integrating by parts in σ we find in the infrared limit 14 Now we use the relation (17), followed by (2), to obtain All steps in this sketchy derivation could be made watertight, what we are not going to discuss here.

Quantum theory
Construction of interacting quantum field theories faces numerous problems. Not only perturbative solutions of nonlinear equations of physically realistic models involving quantum fields need elaborate techniques in order to avoid the well-known ultraviolet 13 There is, recently, a considerable activity in the field of 'electromagnetic memory' and its connection with the long-range structure of electrodynamics, apparently not aware of Staruszkiewicz's work, and the articles just mentioned; see [28] and the literature cited there. 14 We assume that the field F satisfies all our selection criteria, in particular the condition (10).
divergencies, but the very existence of these equations is not clear beyond perturbative formulation. Whether this state of art is due only to technical problems, or is some better physical input needed, in our opinion is still not clear. However, on the perturbative level the Epstein-Glaser technique [13,30] of proper time splitting of distributions supplies the most fundamental, conceptually clear, even if calculationally formidable, solution of the problem. However, here we are interested rather in the opposite, in terms of momentum transfer, end of scale, the infrared problems. In perturbative calculations these difficulties manifest themselves in the appearance of divergencies in Feynmann diagrams for small momenta. But this perspective only touches the tip of an iceberg: in this case the source of the difficulties is definitely not merely technical, but related to the long-range character of the interaction. To understand the source of the problem, one has to consider the theory from its algebraic structure and representation perspective.
One of the paradigms on which the standard quantum field theory is build is relativistic locality. This is best expressed in algebraic terms [16] and consists of two assumptions: (i) basic quantum observables are localized in compact spacetime regions, and all other observables are limits of such quantities; (ii) observables localized in regions spacelike separated commute. For electrodynamics this has the following consequence [5,6,15]: Suppose that in the representation space one can define the spacelike asymptotic field of the type indicated by the classical expression (24). As the localization of this field becomes spacelike to any compact set in the limit, this asymptotic field commutes with all observables. In any irreducible representation this implies that the field is proportional to the unit operator and its specific value is a superselection label of the representation. In the language of previous sections this means that in such setting V (−∞, l) = V (+∞, l) is a classical function characterizing the choice of representation. In view of this knowledge, we look once more at relations (22) and (25). If quantum charged free 'in' and 'out' particles are ascribed their Coulomb fields, then incoming and outgoing free electromagnetic fields cannot be both represented in unitarily equivalent way (e.g., as mentioned before, if a in is Fock, then a out cannot be Fock). Thus, the description of scattering falls into troubles.
Standard way to deal with this problem is to 'dress' charged particles, in addition to their Coulomb fields, with some permanent 'radiation clouds.' These clouds must be chosen so as to decouple the long-range tail of the field attached to the charged particle from its specific momentum and make the rhs of the relation (25) vanish. 15 Another potentially possible strategy is to use representations which do not satisfy the above assumptions and are so singular in the infrared/long-range regime as to be able to mask radiative soft additions (e.g., so-called KPR infravacua [25]).
In what follows we want to avoid any of the above two paths. We insist that free charged particles carry their Coulomb fields, and these fields alone form their structural parts. The only way not to come in conflict with the above discussion is to admit certain degree of nonlocality in the theory: the long-range variables will acquire genuinely quantum character. Having agreed to this, we proceed as follows.

Construction of the algebra
We want to use the heuristic method of 'relativistic quantization,' in which one demands that the classical expressions for energy-momentum and 4-angular momentum should become the Poincaré generators of the quantum theory. This is not in conflict with the knowledge that Lorentz symmetry may be broken in electrodynamics: the heuristic idea is applied only on the algebraic level, while breaking of symmetry happens on the level of representations. The form of conserved quantities (32) suggests that f and V should be quantized independently. We denote the corresponding quantum variables f q and V q and then the quantization conditions are easily obtained: Here f i and V i are test variables. The variable V q has a well-defined physical meaning as the total electromagnetic field at null infinity. On the other hand, as argued above, the physical field equipped with its Coulomb field should be defined by (34), thus on the quantum level g q (v) = exp i{V e (v, .), V q } f q (v). We now use the pair g q , V q as generating a closed algebra. Denoting one obtains [Ψ (g 1 ), Ψ (g 2 )] + = 0, Ψ (g 1 ), Ψ (g 2 ) * + = (g 1 , g 2 )1, where One notes that , when acting on a state, adds to this state the field V 2 . However, the Coulomb fields carried by particles are now attached to the field Ψ (g), so the test fields V i should be free fields test functions; therefore, we demand V i (+∞, l) = 0. Function Φ(l) is related to V (−∞, l) as described in Sect. 3.3. Relations (40)-(43) may be now considered detached from the heuristic considerations which led to them. One shows that they generate a C * -algebra, which simply means that they can indeed be represented by bounded operators in a Hilbert space. Physically meaningful representations must satisfy some selection criteria, and we shall briefly comment on the form of possible representations below. But first we want to clarify the question of gauge (in)dependence of the algebra.

Gauge invariance
First, we recall that the test functions V i (s, l) are free field functions, in particular l · V i (s, l) = 0. It follows that any gauge transformation of the quantum field V q (s, l) → V q (s, l) + lα(s, l) leaves the quantity {V i , V q } unchanged 16 . Thus W (V i ) should be interpreted as exponentiations of the total electromagnetic field, and not merely potential. Second, the field f was obtained in a gauge-independent way, and the transition from f to g is given by (35). Thus, on the (heuristic) quantum level the gauge transformation causes the transformation g q → e ieσ g q , where σ = 1 4π α(−∞, l)d 2 l. In terms of algebraic elements: Ψ (g i ) → e ieσ Ψ (g i ). Summarizing, gauge transformations lead only to transformations of the first kind of the elements of the algebra.
Another question related to gauge symmetry is whether the elements of the algebra depend on the gauge of test fields V i (s, l). To answer this, it is sufficient to consider W (lα)-element with pure gauge test field. Then {lα, V q } = −σ Q q , where Q q = l · V q should have the interpretation of the total charge (see Eq. (31)). Therefore, one should have W (lα)Ψ (g)W (lα) * = e −ieσ Ψ (g), what agrees with the relations (42), (43). This also shows that the only dependence on gauge is through the factor appearing in the definition (43) of S Φ . Therefore, we interpret the additive gauge constant in eΦ(l) as a phase variable and identify eΦ ∼ eΦ mod 2π .

Representations
Physically meaningful representations of the algebra must satisfy some selection criteria. We impose two standard conditions: regularity of representation of W (V ) (which means that these unitary operators are in fact exponentiations of self-adjoint generators) and positivity of energy. The latter condition means that translations are represented by unitary transformations and their generators-energy-momentum operators-have joint spectrum in the forward lightcone. It was shown in [19] that such representations have the following form. The representation space is H = H F ⊗ H r , where H F is the usual Fock space for Dirac fermions, on which Ψ (g) = (g, g q ) act in the usual way. The space H r carries some regular, positive energy representation W r (V ) of commutation relations (40). The representation of W (V ) is then given by where V q (+∞, l) is the quantum version of (33): -normal ordering in n q .
Representation W r (V ) is not uniquely defined. However, it is important to realize that the algebra contains in a substantial way variables at infinity. One of the consequences is that W r (V ), the electromagnetic part of the representation, cannot be the usual Fock representation or any coherent state representation. More than that, there can be no vector state with zero energy. However, a class of representations which are close to standard structures may be formed by a direct integral over coherent state representations [19]. In the following, we shall use properties of these representations.

Variables at spacelike infinity and quantum transformations induced by them
Let us discuss the spacelike asymptotic structure of the electromagnetic field in some detail. As was seen in (23) and (24), the spacelike tail is determined by V (−∞, l).
In quantum case this has to be smeared with some test function, so let us consider the quantity where V + (l) is a smooth test function related to a homogeneous function ε + (l) as described in Sect. 7. In spite of the smearing in l, the quantity (45) is not yet completely defined, as it involves a formal value at s = −∞. Now, we split The first part involves only the Coulomb fields of outgoing particles, and in our representation has a well-defined meaning (44), giving the third equality by an analogue of (53). The second part V q (−∞, l) − V q (+∞, l) involves only the free quantum outgoing field, thus it should be possible to obtain the corresponding part of U by some limiting process U free (V + ) = 1 ⊗ lim β 0 W r (V + β ), with appropriately chosen V + β (s, l). This indeed is the case, and we sketch the solution. For functions of s we denote the Fourier transform: Let h(ω, l) be a smooth function, fast vanishing for |ω| → ∞, such that h(λ −1 ω, λl) = h(ω, l) (λ > 0), and h(0, l) = 1. Denote 17 Then V + β (s, l) is homogeneous of degree −1 and vanishes for |s| → ∞. In consequence, as one can show, its symplectic form with any classical function V (s, l) with asymptotes V (±∞, l) may be written as Now, using the definition (46) one finds This is a classical calculation, but remembering that W r (V + β ) = e −i{V + β ,V q } one can ask whether a similar limit exists for this operator-valued expression. It turns out that in the representations mentioned above it does: for β 0 the operators W r (V + β ) converge weakly to a unitary operator, which implies that the convergence is in fact in the strong operator sense. Thus, we define the unitary operators: The transformation induced by these unitaries on the basic variables of the algebra is easily calculated and yields: where S ε + is the transformation defined in (43). For the interpretation of this result, we note what follows.
1. The asymptotic algebra is generated by the second-kind gauge-invariant quantities W (V ) and Ψ (g), which are nontrivially transformed under the action (47) and (48).
2. The action of the free part U free (V + ) transforms W (V 1 ), while leaving Ψ (g 1 ) unchanged. Recall the interpretation of the elements of the algebra supplied by (39), in particular that W (V 1 ) = exp[−i{V 1 , V q }], and also note that Therefore, the rhs of (47) should be interpreted as exp[−i{V 1 , V q + V + }], which means that this transformation adds a classical constant function V + (l) to V q (s, l). As discussed in concluding paragraphs of Sect. 7, this addition is naturally interpreted as the asymptote of the Lorenz potential of the Coulomb field of some classical external, charge-free current added by this transformation. Such addition does not change the quantum field Ψ (g 1 ).
3. The action of the Coulomb part U Coul (V + ) transforms Ψ (g 1 ), while leaving W (V 1 ) unchanged. To interpret this transformation we choose any free Lorenz potential asymptote V (s, l) (thus V (+∞, l) = 0) such thatV (s, l) is even in s, and V (−∞, l) = V + (l). This asymptote defines an even potential A(x) according to (7). The scaled function V λ (s, l) = V (s/λ, l), λ > 0, produces the scaled potential A λ (x) = λ −1 A(x/λ), which for λ → ∞ gives the electromagnetic field with energy content tending to zero, but with unchanged spacelike tail. The operator W (V λ ) creates this electromagnetic field as discussed after (43). Using the algebraic relations (40) and (41), one finds that the transformation of Integrating by parts, we obtain V (s/λ, l) ·V 1 (s, l) dsd 2 l, which vanishes in the limit λ → ∞, as V (0, l) = V + (l)/2 due to evenness oḟ V (s, l). Now, the operators W (V λ ) do not converge in this limit, but the transformation induced by them on the basic variables agrees in the limit with that induced by U (V + ), which is therefore interpreted as the infrared limit transformation. 4. The above discussion shows that the variables at spacelike infinity have a kind of duality structure: the Coulomb tail variables are conjugate to the infrared free tail variables.

Scattering
The constructions started at the beginning of Sect. 6 and further developed up to now refer to the causal future-outgoing fields. Similar constructions of the asymptotic algebra and its representation may be performed for causal past-incoming fields. Smearing functions g in Dirac fields Ψ in (g) are then as in outgoing elements Ψ out (g). A slight modification occurs for electromagnetic fields. Recall that the smearing test functions in W out (V ) are free field future null asymptotes V (s, l) vanishing for s → +∞. Recall, also, that the past null asymptote of this field is then given by V (s, l) related to V (s, l) by (12). Therefore, we take as test functions for 'in' case these past asymptotes: W in (V ). The mapping Ψ in (g) → Ψ out (g), W in (V ) → W out (V ) is then a canonical isomorphism of the two algebras. If a complete theory could indeed be developed along the lines sketched here, then there should exist the scattering operator S in the representation space such that on the level of representations W out (V ) = S * W in (V )S, Ψ out (g) = S * Ψ in (g)S. 18 By an appropriate limit described in the last section, one obtains the long-range variables U in (V + ) and U out (V + ), which are thus also related by this adjoint transformation. However, in fact the variables at spacelike infinity are conserved, so U in = U out ≡ U . Therefore, one obtains This is convergent with an observation made by the authors of [24].

Concluding remarks
The problem of understanding, from fundamental point of view, the long-rangeinfrared structure of quantum electrodynamics is still open. There are many ideas how to approach its solution, some of them mentioned above. The persistent resistance of the field to be fully understood provokes radical attempts, including a recent proposal [7] to deny the spacelike infinity quantities an experimentally accessible status. This article summarizes a proposal going in the opposite direction (as compared to [7]). The scheme summarized above goes outside the strict local paradigm and admits into the theory some nonlocal observables at spacelike infinity. Whether it may be developed further to become indeed the asymptotic description of interacting theory is still to be seen. The scheme seems to be a rather direct quantization of the classical causal asymptotic structure, if free charged particles are to carry their Coulomb fields, not extended by some 'clouds'-a theoretical construct, rather not having experimental motivation. As we have seen, the scheme necessarily leads outside traditional structures of popular quantum field theory, like Fock space, and needs a substantial mathematical care in defining and handling quantum variables. The scheme has been further developed in other directions, including spacetime localization properties and scattering with classical currents [21,22].
As the Dirac distribution δ(k 2 ) is also homogeneous of degree −2, we have ∂ · k f (k) δ(k 2 ) = 0 outside any neighborhood of k = 0. Thus, using the Gauss theorem we obtain for vectors t, t (θ is the Heaviside step function): which shows that I t in fact does not depend on t (note that the expression in braces has a compact support outside zero, so there are no boundary terms).
In fact, as mentioned in a footnote in Sect. 2, the invariance of the integral is much larger. Let ρ(k) be a homogeneous function of degree +1, such that ρ(k) = 1 is a Cauchy surface cutting the forward lightcone (any such surface can be defined in this way in the neighborhood of the lightcone). Replacing in the last equation t · k by ρ(k), one shows that I t is also equal to

Appendix B: Variables R, s, l
Let B(x) be any differentiable field on Minkowski space. We choose a unit, futurepointing timelike vector t, and for R ≥ 0, s ∈ R and l ∈ C + write Differentiation: as the R-dependence of b(R, s, l) is determined by the s, ldependence and the scaling property, it is sufficient to apply ∂ s ≡ ∂/∂s and L ab . We find Contracting the second equality with t a and using the first equality, one finds ∂ b B(x) = l bḃ (R, s, l) + t a R t · l L ab b(R, s, l), whereḃ ≡ ∂ s b.