On electromagnetic and colour memory in even dimensions

We explore memory effects associated to both Abelian and non-Abelian radiation getting to null infinity, in arbitrary even spacetime dimensions. Together with classical memories, linear and non-linear, amounting to permanent kicks in the velocity of the probes, we also discuss the higher-dimensional counterparts of quantum memory effects, manifesting themselves in modifications of the relative phases describing a configuration of several probes. In addition, we analyse the structure of the asymptotic symmetries of Maxwell's theory in any dimension, both even and odd, in the Lorenz gauge.


Introduction and Outlook
Memory effects are permanent changes in physical quantities pertaining to probes facing the passage of a burst of radiation close to or at null infinity. In this paper we illustrate various types of memories concerning matter probes in interaction with electromagnetic or Yang-Mills radiation, in arbitrary even-dimensional spacetimes.
Following earlier works concerning gravitational memories [1,2], Bieri and Garfinkle [3] first pointed out the existence of analogous observables for electrically charged matter in D = 4. For charged test particles close to null infinity, they identified the relevant observable in the form of a "velocity kick" that the probes would experience as a consequence of the passage of a radiative perturbation. They also observed, in the same context, that the possibility to source those effects either with matter that does not reach null infinity or with massless sources allows one to identify two types of effects, termed "linear" (or "ordinary") in the former case and "nonlinear" (or "null") in the latter.
Memory effects caused by electromagnetic fields were further discussed in [4] and later connected to asymptotic symmetries of four-dimensional Maxwell's theory, and thus identified as the classical counterparts of soft theorems, in [5,6]. See also [7,8]. Moreover, the possible existence of analogous effects in dimension four or higher was investigated from various perspectives, with different (and sometimes not fully compatible) outcomes in [9][10][11][12][13][14][15]. Those works adopted various approaches, focussed either on the classical solutions to the corresponding wave equations or on the perturbative expansion of the gauge potentials with suitably assigned boundary behaviour.
One of the main goals of this work is to describe ordinary and null memory effects for both Maxwell and Yang-Mills theories in even-dimensional spacetimes. To this end, we inspect the perturbative behaviour of the corresponding potentials in the Lorenz gauge and connect it to underlying residual gauge symmetries.
Following [9,10], we begin in Section 2 by analysing the classical fields generated by specific background currents, so as to construct a number of concrete examples with which to compare our general results of the ensuing sections. In particular, we compute the null memory induced by a charged particle moving at the speed of light on a test particle initially at rest in even dimensions. In Section 3 we discuss various types of memories for Maxwell's theory in even-dimensional spacetimes. Upon employing a recursive gauge-fixing procedure, analogous to the one adopted in [16] in the gravitational setup, we clarify the interpretation of both ordinary and null memory in any even D as a residual gauge transformation acting at order r 4−D . We generalise this procedure to the non-Abelian case in Section 4, where we obtain similar results and interpretation for linear and non-linear memory effects concerning matter in interaction with Yang-Mills radiation.
Both ordinary and null memory effects concern the permanent modifications affecting kinematical properties of the probes (the velocity in the case of spin-one radiation, the relative displacement in the gravitational case) due to the interaction of the latter with a transient radiation field close to null infinity. On the other hand, the interaction itself takes place in force of the matter probes possessing the required quantum numbers, electric charge or colour charge, whose overall configuration, whenever several particles are involved, can only be defined by providing the value of the corresponding connection field on the celestial sphere. In this sense, one should expect that the configuration itself be sensitive to the passage of radiation whenever the latter entails a vacuum transition between connections that differ by a (large) gauge transformation, even in the idealised case of unperturbed kinematical variables (i.e. charges kept fixed at specified angles on the sphere). We refer to these types of memory effects as "phase memory" for the electromagnetic case and "colour memory" in the non-Abelian setting, to distinguish them from the memory effects affecting kinematical properties of the probes, that we collectively term simply as "electromagnetic" and "Yang-Mills" memories, respectively. To our knowledge, in D = 4 phase memory was first discussed in [17] while colour memory was analysed in [18]. Here we explore these phenomena in the higher, even-dimensional case, in Sections 3.4 and 4.4, respectively. Let us observe that, differently from the ordinary and null memories that are purely classical effects, phase memory and colour memory naturally arise at the quantum level. It would be interesting to explore possible gravitational counterparts of the colour memory effect.
Memory effects in D > 4 appear at a subleading order with respect to radiation. Correspondingly, their counterparts at the level of gauge transformations do not comprise properly-defined infinite-dimensional asymptotic symmetries. The very existence of the latter, on the other hand, was sometimes doubted in the literature because of the need for too-slow falloffs on the gauge potentials that could in principle lead to divergences in physically sensible quantities. The asymptotic structure of spin-one gauge theories, in particular in connection to soft theorems was widely investigated in the literature [19][20][21][22][23][24][25][26][27][28][29][30][31][32]. More recently, the issue concerning the higher-dimensional extensions of those results was further explored in a number of works [33][34][35][36]. On our side, in Section 5 we provide an analysis of the asymptotic symmetries of Maxwell's theory in any D, even and odd, stressing analogies and differences between the two cases, always working in the Lorenz gauge. In addition, we describe an explicit procedure for defining the surface charges on I .
In the appendices we collect a summary of our conventions (Appendix A) together with a number of ancillary results. In Appendix B we provide a detailed discussion of the classical solutions to wave equations and memories for scalar fields in any D, comprising among other things the calculation of null memory induced by a charged particle moving at the speed of light on a test particle initially at rest in odd dimensions. In Appendix C we provide a closed-form solution of ǫ = 0 in any even dimension with boundary condition ǫ(u, r, n) → λ(n) for r → ∞. In Appendix D we complete our exploration of the infrared triangle for electromagnetism [37] by describing the link between Weinberg's soft theorem and the asymptotic symmetry in the Lorenz gauge, in even dimensions.

Classical EM Solutions and Memory Effects
In this section we calculate the fields generated by specific background charged particles, both massive and massless. In particular, we evaluate in any even D the memory effects induced on a test particle at a large distance r from the origin, much in the spirit of [9,10,12], while also including the effects induced by massless charged particles, thus extending the four-dimensional analysis of [3]. Throughout the paper we shall work in the Lorenz gauge, which in Cartesian coordinates reads ∂ µ A µ = 0, so that the equations of motion reduce to a set of scalar wave equations Thus, much of the calculations are essentially the same as those of the scalar case, that we detail in Appendix B. First, we consider the unphysical but instructive case of a static point-like source with charge q created at the origin, 1 where u µ = (1, 0, . . . , 0). The retarded solution is where ϕ denotes the corresponding solution (B.7) for the scalar field that we report here for completeness Here θ(u) denotes the Heaviside distribution, while the coefficients c D,k are specified in (B.5). The coordinates u and r are part of the retarded Bondi coordinates, in which the Minkowski metric takes the form The solution describing a source created at the origin with velocity v is obtained by boosting the previous one according to (B.12) and Moving to retarded coordinates gives where n = n(x i ) is the unit vector identifying points on the (D − 2)-sphere and v i = v · ∂ i n.
Focussing for simplicity on the case of D = 6, let us now consider the radiation field generated by a massive particle with charge q sitting at rest in the origin for t < 0, that starts moving with velocity v at t = 0. Such a field is obtained by matching the solution for a charge destroyed at the origin to that of a charge there created with velocity v. Proceeding as described in Appendix B, one obtains, for large values of r, where we kept the orders relevant to the calculation of the memory effect. The change in the angular components of the momentum P µ of a test charge Q, initially at rest, gives rise to a linear (or ordinary) memory effect that, to leading order, reads where F µν is the Maxwell field strength.
With hindsight, having in mind in particular the results of [10], in order to interpret this leading memory effect in terms of a symmetry, it is useful to perform a further gauge fixing in order to get rid of the residual symmetries affecting the radiation order. Choosing a gauge parameter allows us to cancel the leading term of A u , compatibly with the condition ǫ = 0. The resulting field after performing this gauge transformation satisfies 2 so that, in particular, A i is fully responsible for the memory formula (2.8). The effect is proportional to the variation of the latter field component between u > 0 and u < 0 and takes the form of a total derivative on the sphere, close to null infinity: .
This result provides an explicit connection between the memory effect and a residual symmetry acting, for large values of r, at Coulombic order. 3 Let us turn our attention to the case of null memory. The four-dimensional case was illustrated in [3]. Let us consider, in any even dimension D ≥ 6, the field generated by a charge moving in thex direction at the speed of light: in Cartesian coordinates, with v µ = (1,x) and |x| = 1. Taking into account the corresponding retarded solution (B.24) for the scalar field ϕ, we then have A µ = −v µ ϕ and, moving to retarded coordinates, where (∆ − D + 4) −1 denotes the Green function for the operator ∆ − D + 4 and we have omitted terms of the form (see (B.24)), which would be leading with respect to those displayed in (2.14), for suitable values of k, but which do not contribute to the memory effect because they integrate to zero. The null memory formula then reads 16) where this result is exact for the solution considered, while for more general solutions of the wave equation in the form (B.23) it would provide the leading-order contribution.

Electromagnetic Memory
In the previous section we saw how to establish a connection between local symmetries acting at large r and memory effects for specific matter configurations. Now we would like to investigate this connection beyond those examples by studying the general structure of the solution space in the Lorenz gauge. To this end, we will analyse the asymptotic behaviour of the fields in a power-series expansion in the radial coordinate in D > 6, which is sufficient to the description of the memory effects. More general possibilities, also encompassing the four-dimensional case, will be discussed in Section 5, where we will explore the full structure of the asymptotic symmetries in the Lorenz gauge.

Electromagnetism in the Lorenz gauge
The Lorenz gauge condition reads where D i is the covariant derivative on the Euclidean unit (D − 2)-sphere with metric γ ij . The residual symmetry parameters then satisfy ǫ = 0, namely The equations of motion reduce to A µ = 0, which, component by component, read

Asymptotic expansion
We assume an asymptotic expansion of the gauge field and of the gauge parameters in powers of 1/r where the summation ranges are, for the moment, unspecified. Equations (3.1) and (3.2) then give Eqs. (3.5)-(3.9) appear in particular to order r −(k+1) in the asymptotic expansions of the original equations.
As we shall see, to the purpose of analysing the electromagnetic memory effects the leading falloffs can be chosen so as to match the corresponding radiation falloffs: More general options are possible and influence the structure of the asymptotic symmetry group. We will be concerned with these more general aspects of the discussion in Section 5. The significance of the choice (3.10) lies in the fact that the derivatives with respect to u of the field components are unconstrained to leading order: these conditions for the asymptotic expansion are well-suited to identifying the boundary data for a radiation solution with an arbitrary wave-form. Such components also provide the energy flux at a given retarded time, according to where S u is the section of I + at fixed u and dΩ is the measure element on the unit (D − 2)-sphere.
The asymptotic behaviour of radiation differs in higher dimensions with respect to the characteristic Coulombic falloff r 3−D that, in its turn, can be identified as the leading falloff for stationary solutions. (See also the discussion in Section 3.3 on this point.) As we shall see, Coulomb fields give non-vanishing contributions to the surface integral associated with the electric charge as well as to the memory effects.

Recursive gauge fixing
The gauge variations imply a number of restrictions on the allowed gauge parameters, in order to keep the corresponding falloffs. From δA r = O(r − D−2 2 ) we read off ǫ (k) = 0 for k < (D − 4)/2 and k = 0, while δA does not give rise to further constraints at this stage. This leads to a gauge parameter of the following form ǫ(u, r, where we conventionally set the global part of ǫ to 1. However, by (3.6), which implies ǫ ( D−4 2 ) = 0 since the Laplacian on the sphere is negative semidefinite. 4 Thus, we need to search further down in the asymptotic expansion of ǫ, employing a recursive on-shell gauge-fixing procedure. The residual symmetry is parametrised as follows ǫ(u, r, (3.15) Equation (3.6) leaves the u-dependence of ǫ ( D−2 2 ) unconstrained, and therefore we may use it to set A leaving a residual parameter ǫ ( D−2 2 ) (x i ) arbitrary. Setting k = D/2, (3.7) and (3.6) reduce to 4 Interestingly, however, let us observe that the corresponding putative chargẽ is finite and non-vanishing as r → ∞.
respectively. Thus, A ( D 2 ) u is a function of the angles x i only, while δA can be expressed in terms of ǫ ( D−2 2 ) , which can be used to set A ( D 2 ) u = 0, while still leaving ǫ ( D 2 ) (x i ) arbitrary. We proceed recursively, assuming for some q > D/2. Then, for k = q − 1, (3.7) and (3.6) give Therefore, we may employ u to zero provided that the differential operator on the right-hand side is invertible, which is true for any q < D − 2.
We shall now consider two options. We may first choose to truncate the recursive gauge fixing right after the step labelled by q = D − 4, which leaves us with the asymptotic expansions where A ( D−2 2 ) r = 0 on shell. The residual symmetry is then given by (3.21) whose corresponding charge, evaluated in the absence of radiation close to the past boundary I + − of I + , reads 5 Alternatively, we may also perform the recursive gauge-fixing until the very last allowed step, q = D − 3, thus obtaining The latter choice highlights the possibility of making the components gauge-invariant, and hence in principle responsible for any observable effect due to radiation impinging on a test charge placed at a large distance r from a source. Indeed, consistently with the examples of Section 2, in the next section we shall see that electromagnetic memory effects appear at the Coulombic order A

Electromagnetic memory
A test particle with charge Q, initially at rest at a large distance r from the origin, will experience a leading-order momentum kick due to the presence of an electric field according to We are neglecting the contribution from the magnetic field, assuming the interaction time between the test particle and the radiation is sufficiently small. In particular, (3.26) holds to leading order for an ideally sharp wave front, as we checked in Section 2. Now we shall consider solutions that are stationary before u = u 0 and after u = u 1 . For such solutions, the Maxwell tensor does not vanish outside this interval, in general, because static forces are present. However, it does not contain radiation, and thus all the components of the gauge potential associated to radiation are to vanish, or, more generally, are to be pure gauge. In particular, the radiation field components before u 0 and after u 1 are to be identical or are to differ by a gauge transformation.
Let us combine the information of the previous recursive gauge-fixing with the requirement that the solution be stationary before u 0 and after u 1 . The coefficients F (k) iu of the asymptotic expansion of the electric field u are zero, thanks to the gauge-fixing. Thus, (3.28) With respect to our previous observation, let us notice that the components of A i that enter the subleading terms O(r 3−D ) are those connected with non-vanishing stationary properties of the field (we assume no long-range magnetic forces), while all the leading components explicitly written enter the radiation behaviour and thus their difference after u 1 and before u 0 can be at most the angular gradient of given functions. In addition, we shall immediately see that, combining this information with the equations of motion will allow us to appreciate that they all vanish with the exception of the last one at order r 4−D .
Indeed, let us note that for a stationary solution, in our gauge, equations (3.5), (3.7), (3.8) and (3.9) read, for k < D − 2, and = 0 consistently with the recursive gauge-fixing. Equations (3.29) and (3.31) then give, for 1 < k < D − 2, (3.36) By equation (3.28), the condition on A i implies that the leading-order memory effect appears at O(r 4−D ), In view of the discussion below (3.28), we conclude that the momentum shift must take the form where g(x i ) is a u-independent function, which depends on the shape of the radiation train and in particular on u 0 and u 1 (see for instance the example (2.11)). Let us note that, as it must be, this difference is not affected by the action of the residual gauge transformation (3.21), because the latter is u-independent and thus cannot alter the difference A In this sense, whether or not one performs the last step of the recursive gauge fixing is irrelevant to the extent of calculating the electromagnetic memory.
To conclude, we have established a formula that, for any even D, exhibits a momentum kick characterising the transition between the initial and final vacuum configurations, parametrised by the gauge transformation g(x i ), induced by the exposure to electromagnetic radiation reaching null infinity. In particular, the norm of this effect scales as r 3−D .
Up to this point we have only been dealing with an ordinary memory effect; in order to encompass null memory, we must modify the equations of motion (3.3) by adding a suitable source term to the right-hand side, namely a current density J µ allowing for the outflow to future null infinity of charged massless particles. The falloff conditions on such a current can be taken as follows This is clearly displayed by the current for a single massless charge q moving in thex direction, whose components read in Minkowski coordinates 40) and in retarded coordinates (for t = u + r > 0) To the purposes of the recursive gauge fixing, the only modification is thus the introduction of a source term J (D−2) u in the right-hand side of (3.7) when k = D − 3, which now actually forces us to stop the gauge fixing after the use of the said equation for k = D − 4 (the step labelled by q = D − 4 in the previous section) and leaves us with the falloff (3.20). On the contrary, reaching (3.23) is not allowed, and thus (3.20) comprises a complete gauge fixing. However, also in view of the above considerations, the discussion of the memory effect and its relation to the symmetry acting at Coulombic order remain unaltered.

Phase memory
Let us consider a pair of electric charges q and −q that are pinned in the positions (r, n 1 ) and (r, n 2 ), for large r. We will now derive the expression for an imprint that the passage of a radiation train leaves on the properties of these particles that is encoded in the phase of their states. What follows is inspired by the four-dimensional discussion of [17]. For a quantum treatment of electromagnetic kick memory see instead [14].
To this purpose, let us assume that, as in the previous section, radiation impinges on the charges only during the interval between two given retarded times u 0 and u 1 . As we have seen, this means that for a suitable angular function g(n). We will assume for simplicity that the gauge field before the onset of radiation is the trivial one. Let |ψ 1 = |q and |ψ 2 = | − q be the initial states in which the charged particles are prepared, which are uniquely labelled by their charges since translational degrees of freedom have been suppressed. Before u 0 , the state |ψ 2 , n 1 obtained by the parallel transport of the second state |ψ 2 to the position n 1 of the first the charge is because A µ = 0, so that the corresponding tensor state evaluated in n 1 is given by At any fixed u ≥ u 1 , instead, the same operation of parallel transport must be performed by calculating where we have employed (3.42). Therefore, after the passage of radiation, which displays how the transition between two different radiative vacua, already experimentally detectable by the occurrence of a non-trivial velocity kick for a test charge, is also signalled by the variation of the relative phases in the states obtained by parallel transport of charged particles. Such a phase can be non-trivial provided that the function g is non-constant, namely when there is a non-trivial memory kick (3.38).
A point that should be stressed is that, in this setup, the states |ψ 1 and |ψ 2 do not evolve, since each particle is kept fixed in its position (its translational quantum numbers are frozen) while electromagnetic radiation cannot change its charge. The relative phase difference occurs entirely as an effect of the evolution of A µ , which undergoes a transition between two underlying radiative vacua. As we shall see in the following, this aspect is qualitatively different in a non-Abelian theory such as Yang-Mills, where radiation can alter the colour charge.

Yang-Mills Memory
In this section we extend the previous analysis to the non-Abelian case.

Yang-Mills theory in the Lorenz gauge
We consider pure Yang-Mills theory of an anti-Hermitian gauge field We impose the Lorenz gauge condition ∇ µ A µ = 0, which leaves as residual gauge parameters those that satisfy ǫ + [A µ , ∂ µ ǫ] = 0. Furthermore, the equations of motion Adopting retarded Bondi coordinates, the Lorenz gauge condition is tantamount to while the constraint on residual transformations is The equations of motion are instead

Asymptotic expansion and colour flux
Performing the usual asymptotic expansion in inverse powers of the radial coordinate r, one obtains the following set of equations. Equations (4.1) and (4.2) give and respectively, while from (4.3) one obtains for the corresponding components of the equations of motion. Eqs. (4.4)-(4.8) appear to order r −(k+1) in the asymptotic expansions. We choose to adopt the same "radiation" falloff conditions (3.10) that we imposed in the linearised theory: the asymptotic expansions of A u and A r start at order r −(D−2)/2 and that of A i starts at order r −(D−4)/2 .
For completeness, and as a cross-check of the choice of falloffs, let us now discuss the definition of the total colour charge Q(u) = Q A (u)T A and its dependence on the retarded time u. In the Lorenz gauge, this quantity is given by the surface integral in the limit r → ∞. Combining the Lorenz condition (4.4) and the r component of the equations of motion (4.7), one obtains ] . (4.10) We see that (4.4) implies ∂ u A ( D−2 2 ) r = 0 and that (4.10) reduces to for k < D − 2. The final expression for the colour charge may be cast in the form by means of the Lorenz condition (4.4).
Concerning the dependence of Q A on retarded time, recalling that ∂ u A ( D−2 2 ) r = 0 by the Lorenz condition, and employing (4.6), (4.10) for k = D − 3, we get (4.13) This formula provides the flux of the total colour charge due to non-linearities of the theory. Note in particular that the right-hand side involves the radiation components, representing the flux of classical gluons across null infinity.

Yang-Mills memory
Starting from the radiation falloffs (3.10), we may employ the residual gauge symmetry of the theory to perform a further, recursive gauge fixing. Actually, since the non-linear corrections to equations (3.19) appear to order q = D − 3 or higher, the discussion of this gauge fixing proceeds as in Section 3.2.1 up to the next to last step. Indeed, in accordance with the fact that Yang-Mills theory must encompass both ordinary/linear and null/non-linear memory, the gauge fixing stops at q = D − 4 and cannot be performed up until q = D − 3, as already observed in the case of null electromagnetic memory.
The resulting falloffs after this procedure are thus with residual symmetry parameter where c A are constant coefficients. A coloured test particle with charge Q = Q A T A interacts with the background Yang-Mills field by the Wong equations [38] P µ = tr(QF µν )ẋ ν ,Q = −ẋ µ [A µ , Q] . (4.16) Focussing on the region near null infinity, a test "quark", initially at rest, invested by radiation between the retarded times u 0 and u 1 , will therefore experience the following leading-order momentum kick (4.17) where we have chosen the vacuum configuration at u = u 0 to be A µ = 0. On the other hand, the colour of the test quark will change, to leading order, according to where the u-dependence of A u is governed, according to equation (4.7) with k = D − 3, by (4.19) Equation (4.19) characterises the evolution of the charge of the quark in terms of the leading outgoing radiation terms.
In order to better understand the dependence of this momentum kick and of colour evolution on the incoming radiation near null infinity, let us further analyse the equations of motion. Combining equation (4.4) and (4.7) we see that and for k = D − 3 in dimensions D > 4, where we have defined and, taking the derivative of the previous equation with respect to u, Starting from eq. (4.22) and employing (4.20) recursively, we find that and, by equation (4.24), (4.28)

Colour memory
Let us now consider a pair of colour charges that are pinned in the positions (r, n 1 ) and (r, n 2 ), for large r. Analysing the effect of the passage of Yang-Mills radiation will provide the non-Abelian analogue of the phase memory effect highlighted in Section 3.4, namely a relative rotation between their states in colour space. The content of this section is basically inspired by [18], where colour memory was first discussed and whose four-dimensional results we extend to arbitrary even-dimensional spaces. We consider coloured particles in the fundamental representation of the gauge group that, for definiteness, we first take to be SU (2). Assuming that radiation is nontrivial only between two given retarded times u 0 and u 1 , as we have seen, the angular components of the gauge field satisfy for a suitable α(n) = α A (n)T A , where T A are the su(2) generators. We have assumed for simplicity that the pure gauge configuration before u 0 is the trivial one. Let |ψ 1 and |ψ 2 be the initial states in which the coloured particles are prepared. Before u 0 , the state |ψ 2 , n 1 that results from the parallel transport of the second state |ψ 2 to the position n 1 of the first the charge is because A µ = 0. In particular, one can take it to be one of the eigenstates of T 3 , i.e. |+ or |− . In order to build a colour singlet state in n 1 it then suffices to prepare the superposition We remark that, in order for the concept of singlet to be well-defined when SU(2) transformations can depend on the position, it is crucial to build a singlet out of states that have been parallel-transported to the same point. A given state |ψ evolves according to the covariant conservation equation [39] 6 where we recall that A u = O(r 3−D ) for large r. Thus, the state after the passage of radiation differs with respect to its initial value by Now, at u 1 one must also take into account the effect of parallel transport from n 2 to n 1 on the sphere, given by  where P denotes path ordering and we have employed (4.29) to establish the second equality. We conclude that the actual state after radiation has passed is no longer the singlet (4.31), but rather: that manifests itself to order O(r 4−D ) in the final state. In other words, a pair of particles initially prepared in a singlet will no longer be in a singlet after the passage of radiation. This colour memory can be non-trivial provided that the function α is non-constant, namely whenever there is also potentially a non-trivial memory kick (3.38).
From a technical point of view, it should be noted how the effect of radiation on each colour state is to induce a time evolution to order O(r 3−D ), according to (4.33), which allows one to disentangle it from the leading effect due to the vacuum transition to (4.36), due to the non-trivial configuration (4.29) attained by A i .

More on Asymptotic Symmetries for the Maxwell Theory
So far, we focussed on the symmetries directly related to memory effects. As we saw, these symmetries act at Coulombic order so that the corresponding charge (3.22) displays a falloff O(r 4−D ). Thus, they do not comprise asymptotic symmetries, stricto sensu, since the associated charges vanish at I + . 7 In this section we aim to complement our discussion by performing a more general analysis of the asymptotic symmetry algebra for the Maxwell theory in higher dimensions. 8 7 In addition, one should check whether the symmetry is actually canonically generated by those charges. This type of analysis, however, lies beyond the scope of our work. 8 Let us mention, however, that contrary to the calculation of memory effects, which appear to be tightly related to the properties of the wave equation in even-dimensional spacetimes, these further investigations of asymptotic symmetries admit an extension to spacetimes of not only even but also odd dimension. While focussing on the former case in the main discussion, we will comment along the way on the differences arising in the latter.
In the radial gauge, in D = 4, the standard power-like asymptotic expansion easily allows one to identify the existence of infinite-dimensional asymptotic symmetries. The same result, in the same gauge, can be obtained in higher dimensions upon choosing a power-like ansatz with leading-order terms as weak as 1/r [23]. (Showing the consistency of this choice in the gravitational case allowed to identify the BMS group as the relevant group of asymptotic symmetries for asymptotically flat spacetimes of arbitrary even dimension [40].) We would like to perform the same analysis in the Lorenz gauge.
Differently from the radial gauge, in the Lorenz gauge one first observes that a mere power-like ansatz with leading-order falloffs as weak as 1/r does not allow to recover a non-trivial structure for the asymptotic symmetry group. In order to display the latter, one finds the need for generalising the asymptotic expansion by allowing for logarithmic terms. In this section we will provide an explicit realisation to this effect, while also discussing the non-obvious issue of the finiteness of the corresponding soft surface charges.
The observed difference between radial and Lorenz gauge is one concrete facet of a general issue concerning the possible gauge-dependence of the asymptotic analysis [35,[41][42][43]. While in some instances, such as the present one, this problem has been clarified, the general systematics may still deserve further attention.

Polyhomogeneous expansion for D ≥ 4
In the Lorenz gauge, in order to investigate the possible existence of large gauge symmetries acting at I + we look for solutions to equation (3.2) with the boundary condition lim r→∞ ǫ(u, r, n) = ǫ (0) (u, n) , (5.1) for some non-constant function ǫ (0) (u, n). As it turns out, a power-law ansatz effectively selects only the global symmetry in any even dimension D ≥ 4. Indeed, considering (3.6) with k = 0, 1, . . . , D−2 2 , we have The last equation sets ǫ ( D−4 2 ) to zero for any D ≥ 6, and to a constant for D = 4. Then, the above system recursively sets to zero the other components ǫ (k) all the way up to ǫ (0) since the corresponding operators are all invertible. Eventually, ǫ (0) must be a constant by ∂ u ǫ (0) = 0 and ∆ǫ (0) = 0.
One concludes that the power-law ansatz does not allow for an enhanced asymptotic symmetry sitting at the order O(1). In order to retrieve it, one needs a more general ansatz. We find that it is sufficient, to this purpose, to consider the following type of asymptotic expansion 9 involving also a logarithmic dependence on r [29,[33][34][35]: In this fashion, the last equation in (5.3) becomes modified by the presence of the logarithmic branch and yields an equation that determines the u-dependence of λ ( D−2 2 ) , thus still allowing for an arbitrary ǫ (0) (n).
Indeed, the recursion relations expressing the equation ǫ = 0 read 6) and hence can be solved by direct integration with respect to u. Explicitly, for k < D−2 2 , they reduce to the familiar system so that the solution is given by a polynomial in u with angle-dependent coefficients C j,k (n), with 0 ≤ j ≤ k, and is uniquely determined by specifying the integration functionsǫ (k) (n): where D l was defined in (4.21). If k = D−2 2 , equation (5.6) reduces to (5.5), so that with a suitable integration function, while ǫ ( D−2 2 ) (u, n) stays unconstrained. For k > D−2 2 , we can recast (5.6) in the slightly more suggestive form (5.11) In this way, it becomes clear that these two equations can always be solved by first finding the integral of the first equation directly, and then substituting it in the second one where it acts as a "source" term on the right-hand side.
The situation is simpler in odd dimensions, where the recursion relation (3.6) can be solved compatibly with the boundary condition (5.1) without the need of introducing logarithmic terms. The technical reason is that, when D is odd, the left-hand side of (3.6) never vanishes for integer k.
By means of the transformations parametrised by the gauge parameters thus obtained, we can then act on solutions to Maxwell's equations in the Lorenz gauge characterised by the radiation falloffs (3.10) and generate a wider solution space. In particular, in view of the expansion (5.4) of ǫ, this type of solutions will generically exhibit the following asymptotic behaviour as r → ∞: in any D ≥ 5, while in dimension D = 4.
It is important to remark that, since such solutions differ from those characterised by (3.10) by a pure gradient, they will still retain all the corresponding physical properties, such as the finiteness of the energy flux at any given retarded time.
As far as (5.12) alone is concerned one could in principle consider parameters ǫ that obey arbitrary asymptotics near I + . However, the ones of interest to us are those that define asymptotic symmetries, i.e. parameters that give rise to finite and nonvanishing asymptotic surface charges. In the next section we shall see that this is indeed the case for our class of solutions identified by (5.1).

Asymptotic charges
As it is usually the case in the presence of radiation, the naive surface charge associated to the symmetry (5.4), namely is formally ill-defined, because the right-hand side contains terms of the type e.g.
which do not vanish, even after imposing the equations of motion, precisely due to the presence of an arbitrary parameter ǫ (0) (n). Such difficulties are absent in the case of the global charges because the equations of motion always ensure that these potentially dangerous terms are actually zero. Indeed, ∇ µ F µν = 0 away from electric charges, so that Su,r F µν dx µν cannot depend on r, nor in fact on u, because it must be independent of the specific (D − 2)-surface under consideration.
It should be noted however that, for the general variation with parameter (5.4), the above difficulties arise only if one attempts to calculate the surface charge by integrating over a sphere at a given retarded time u and radius r and then lets r tend to infinity. On the other hand, the calculation of the charge on a Cauchy surface still gives a well-defined result [35]. For instance, under the simplifying assumption that the electromagnetic field due to radiation vanish for u < u 0 in a neighbourhood of future null infinity, the calculation of the surface charge indeed yields for fixed u < u 0 , since all radiation components F For u < u 0 , the quantity Q ǫ (u) must match the analogous surface integral calculated at spatial infinity because in both cases the Noether two-form is integrated over the boundary of a Cauchy surface, in view of the requirement that no radiation be present in a neighbourhood of I + for u < u 0 . In force of the considerations above, we now define the charge Q ǫ (u) for all the values of u as follows. For u < u 0 we define it as in (5.17), while for u ≥ u 0 , even in the presence of radiation, we define it as the evolution of (5.17) under the equations of motion. Indeed, the Maxwell equation ∇ · F r = 0 gives while ∇ · F u = 0 reads and hence d du Analogous considerations allow one to introduce well-defined surface charges evaluated at I − . From the perspective of the analysis performed in Section 2, it is possible to explicitly calculate the soft charges according to the above strategy. Restricting to the case of a massive charge in dimension 6 that starts moving at t = 0, the integral of F ur ǫ (0) on a sphere at fixed retarded time u and radius r yields Except for the case of the electric charge ǫ (0) = 1, where the first integral vanishes identically, the limit of this surface charge as r → ∞ is ill-defined in the presence of radiation, namely on the forward light-cone u = 0, due to the linear divergence appearing in the first line. However, the charge is well-defined on I + before and after the passage of radiation, u = 0, and reads For ǫ (0) = 0, this quantity reduces to the (constant) electric charge Q = 1, while for more general parameters ǫ (0) , the soft charge exhibits a jump discontinuity at u = 0, which measures the fact that the particle is no longer static for u > 0 in a manner akin to the memory effect itself. Performing instead the limit r → ∞ at fixed time t, it is also possible to verify the matching between the surface charge evaluated at null infinity before the onset of radiation, for u < 0 (or, equivalently, at I + − ), and the Hamiltonian charge H ǫ (t), obtained by integrating on a slice at fixed time t. Indeed, taking (B.18) into account and writing the result in terms of polar coordinates t, r and n, we have, for the scalar field, The corresponding electromagnetic potential is given by for t > r, and A µ = (ϕ, 0), for t < r. The radial component of the electric field then yields F tr = 3r −4 as r → ∞ for fixed t and hence Similar arguments showing the finiteness of the Hamiltonian charge in higherdimensions have been given, in the case of linearised spin-two case in retarded Bondi gauge, in [44], while a renormalisation procedure has been recently proposed, for the Maxwell theory in the radial gauge, in [45] (for the general definition of surface charges see [46]). The relation between asymptotic charges and the soft photon theorem in any D was recently clarified in [33,34]. In Appendix D we present the details of the computation for the case of interest in our work.
he also acknowledges the hospitality of ESI. A.C. and C.H. gratefully acknowledge the University Roma Tre for hospitality.

A Notation
Retarded Bondi coordinates are a retarded time u = t − r, a radial coordinate r, and angular coordinates x i on the Euclidean unit (D − 2)−sphere, with metric γ ij . The Minkowski metric, in such coordinates, reads while the (nonvanishing) Christoffel symbols are With D i we denote the covariant derivative associated to γ ij and ∆ = D i D i is the corresponding Laplace-Beltrami operator. In particular, the d'Alembert operator acting on a scalar ϕ takes the explicit form We find it useful to also employ the notation n = n(x i ) for unit vectors identifying points on the sphere, in terms of which γ ij = ∂ i n · ∂ j n.

B Classical Scalar Solutions and Memory Effects
In this appendix, we derive explicit solutions of the scalar wave equation with different types of background sources and calculate the associated memory effects. This provides a useful warming-up for our discussion of electromagnetic memory effects in the Lorenz gauge in even dimensions, while also allowing us to shed some light on the nature of memory effects in odd dimensions.

B.1 Scalar fields in even D
Let us first consider a particle with charge q under the scalar field ϕ that is created in the origin at t = 0. The field generated by this process is obtained by solving the wave equation . The solution is given by the convolution of the source on the right-hand side with the (D-dimensional) retarded wave propagator G ret D (x), i.e. in this case The field generated by a particle that is destroyed in the origin can then be obtained by time reversal of the above solution, while the field for a moving particle can be calculated by applying a Lorentz boost. Let us recall that [47], for even D ≥ 2, the retarded propagator is where θ is the Heaviside distribution. Restricting to even D ≥ 4, and using the chain rule for the distribution δ ( D−4 2 ) (x 2 ), the propagator (B.3) can be recast as where u = t − r and r = |x|, while c D,k are the coefficients In particular note that where Ω D−2 is the area of the (D −2)-dimensional Euclidean unit sphere. The resulting scalar field is thus Notice that only the term associated to k = D/2 − 2 gives rise to a persistent field for fixed r, while the other terms have support localised at u = 0, namely on the futuredirected light-cone with vertex at the particle's creation. This is a general consequence of the recursion relation (3.6), namely obeyed by ϕ = k ϕ (k) r −k near future null infinity, which requires [∆ + k(k − D + 3)]ϕ (k) = 0 for 0 < k < D − 3, and hence ϕ (k) = 0, for any stationary solution. Now, a test particle with charge Q, held in place at a distance r from the origin, will be subject to a force f µ = Q ∂ µ ϕ(u, r) at a given retarded time u due to the presence of the scalar field. Hence, its D-momentum P µ will in general be subject to the leading-order variation For this very simple example, this quantity can be calculated explicitly for any even D. The variations of P u and P r in particular yield and (B.11) Equation (B.10) simply expresses the fact that the test particle will start feeling the Coulombic interaction energy with the newly created particle in the origin as soon as it crosses the light-cone subtended by the origin of spacetime. On the other hand, (B.11) tells us that the test particle will feel an instantaneous, radial momentum kick, for even dimensions greater than four. Since this process is spherically symmetric, the variations of the angular components P i vanish identically. The field emitted by a particle destroyed in the origin at t = 0 is obtained by sending u → −u in (B.7). The case of a particle moving with velocity v can be instead obtained by boosting (B.7): which gives, for large r, denoting n = x/r, We can then cast the boosted solution in the following form: whereφ is a sum of terms of the type namely, whose support is localised on the light-cone. Let us stress that the terms inφ formally dominate the asymptotic expansion of ϕ as r → ∞. However, these terms will not contribute to the leading u-component of the momentum kick due to the presence of δ(u) and its derivatives. We can therefore conclude that in any even D. This is, not surprisingly, just the analogue of equation (B.10) for the Coulombic energy in which one needs to account for the relativistic length contraction. For a more general scattering process involving a number of "in" and "out" particles destroyed or created in the origin, the result is obtained by linearly superposing solutions and therefore reads (η a = −1 for an incoming particle and η a = +1 for an outgoing one) (B.17) Calculating radial and angular components of P µ requires more effort, since they arise instead from the terms proportional to δ(u) whose number increases with the spacetime dimension. They have been given for any even dimension in [10] in terms of derivatives of a generating function. For our present, illustrative, purposes, it suffices to consider the first relevant case D = 6, where the exact solution in the case of the particle created in the origin with velocity v is given by The corresponding radial and angular memory effects in D = 6 are then, to leading order, where v i = ∂ i n·v is the component of the particle's velocity in the i-th angular direction.
While the above examples illustrate the phenomenon of ordinary memory, associated with the field emitted to massive charges that move in the bulk of the spacetime, we can also consider the wave equation with a source term characterising the presence of a massless charged particle, moving along a given directionx: with |x| = 1. This equation can be conveniently solved for any even D ≥ 6 by going to retarded coordinates, where it reads and performing the usual asymptotic expansion ϕ(u, r, n) = ϕ (k) (u, n)r −k , which gives The latter equation is solved by setting ϕ (k) = 0 for k ≤ D 2 − 2 and for k ≥ D − 3, while, for where the functions C k (n) are determined recursively by Here, (∆ − D + 4) −1 is the Green's function for the operator ∆ − D + 4, which is unique for D > 4. As a consequence, the field gives rise to the null memory effect (note that only the term with k = D − 4 contributes), consisting in a kick along a direction tangent to the celestial sphere.

B.2 Comments on the odd-dimensional case
In odd dimensions D ≥ 3 the retarded propagator is given by [47] Considering instead the limit of ϕ as t → +∞ for fixed r, one sees that only the last term in (B.31) survives and yields (B.35) This means that the Coulombic energy due to the newly created particle is felt by the test charge only at i + , namely after one has waited (for an infinite time) at a fixed distance r that the perturbations due to the dispersion occurring in odd spacetime dimensions have died out. To some extent, this is to be regarded as a smeared-out memory effect, as opposed to memory effects occurring sharply at I + near u = 0 in even dimensions (see also [15]).
The situation does not improve if one considers a particle that is created with a non-zero velocity v. Indeed, boosting the exact solution (B.31) by means of (B.12), one sees that ϕ goes to zero for fixed r as t → +∞. The reason is that, while one waits for the dispersion to die out, the source, moving at a constant velocity, has travelled infinitely far from the test charge.
Shifting our attention to the case of null memory, we see that it is possible to provide the following formal solution to the recursion relations (B.23), which hold in any dimension. We consider ϕ = with the functions C k (n) determined recursively by Thus, although the field is highly singular at u = 0, the resulting null memory effect will be formally identical to the one occurring in even dimensions.
C Exact solution of ǫ = 0 The closed-form solution of the wave equation (3.2) with the boundary condition lim r→∞ ǫ(u, r, n) = ǫ (0) (u, n) (C.1) introduced in (5.1) in any even dimension D is given by where q = (1, q) and the limit ε → 0 + is understood. The introduction of this small imaginary part is needed in order to avoid the singularities occurring in the angular integration for |t| < |x|, namely outside the light-cone. Indeed, it is straightforward to verify that, for any value of D even or odd, while, aligning n along the (D − 1)th direction, we have where dΩ ′ (q ′ ) denotes the integral measure on the (D − 3)-sphere. Letting τ = r(1 − cos θ)/u, for u = 0, the previous expression becomes Re 2r/u 0 u dτ r (C.5) which, as r → ∞, tends to The solution (C.2) which, to the best of our knowledge, was not previously exhibited in closed form in the literature, is thus compatible with the asymptotic expansion (5.4) and generalises the expression given in [29]. In addition, being non-perturbative, it allows to explicitly verify the antipodal matching condition: lim r→∞ ǫ(u = v − 2r, r, n) = ǫ (0) (−n) (C.7) for any fixed advanced time v.

D Soft photon theorem in even D
We would like to show that the surface charges defined in Section 5.2 enter Weinberg's soft theorem [48,49]. More precisely, we will see how the Weinberg theorem implies the validity of the Ward identities associated to such charges in D ≥ 4 [33,34].
The surface charge associated to (5.4), evaluated at I + − reads where we have taken into account the absence of radiation terms for u → −∞. Recasting (D.1) as an integral over the whole of I + , and assuming that no contribution arises at I + + , which is the case in particular if there are no stable massive charge, we find where we have used the fact that, after recursive gauge fixing, A (D−3) u is independent of u on shell. We would like to express (D.2) in terms of the leading radiation field, which, as we shall see below, indeed contains the creation and annihilation operator of asymptotic photons. To this end, we first combine (3.8) and (3.5) and obtain employing (3.7) as well, where D k is given in (4.21). Employing this relation recursively, we find 2 ) . In the above writing, we adopt the convention that for D = 4 the product l D l (which in this case has a formally ill-defined range) reduces to the identity. We can then use (D.5) to recast (D.2) as (D.6) On the other hand, the asymptotic expansion of the free electromagnetic field operator, expressed in terms of creation and annihilation operators, yields, to leading order, wherex µ = x µ /r, while ǫ σ are polarisation tensors for the D − 2 propagating helicities.
This formula provides thus an explicit expression for A ( D−4 2 ) i and hence allows us to make explicit the relation between the charge Q ǫ and the soft photon creation and annihilation operators as follows (we employ the prescription (D.8) Assuming that the charge Q ǫ , together with its counterpart at I − , generates the residual symmetry δψ(u, r,x) = iǫ (0) (x)+O(r −1 ) in a canonical way, and employing suitable antipodal matching and crossing symmetry conditions, we have the Ward identity where the sum on the right-hand side extends to all charged external particles in the amplitude and e n is the electric charge of the nth particle (taking into account with a suitable sign whether the particle is outgoing, resp. incoming). Notably, the left-hand side contains exactly the combination P[ · ] = lim ω→0 + [ω · ] that selects the pole in the amplitude with the soft insertion.
On the other hand, the Weinberg theorem for an amplitude involving external massless particles with momenta p n = E n (1,x n ) and a soft photon emitted with helicity σ pointing along then direction on the celestial sphere reads lim ω→0 + out|ωa σ (ωn)S|in = n e n ǫ σ (n) * · p n p n · (1,n) out|S|in . (D.10) Multiplying this relation by ǫ σ i (n) and summing over σ, we see that this is equivalent to ǫ σ i (n) lim on equation (D.11), multiplying by an arbitrary ǫ (D−4) (x) and integrating over the unit sphere allows one to retrieve the Ward identity (D.9), thanks to the relation (D. 13). This proves that the Weinberg factorisation implies the existence of the asymptotic symmetry Ward identities.
Remarkably, the charge (3.22) associated to the symmetry (3.21), which is responsible for the memory effect, formally differs from (D.1) only by a factor of 1/r D−4 (other than by the substitution ǫ (0) ↔ ǫ (D−4) ), which makes it vanish on I + . However, the corresponding symmetry transformation of the matter fields would be δψ(u, r,x) = iǫ (D−4) (x)/r D−4 + O(r 3−D ), and hence would give rise to Ward identities completely equivalent to (D.9), with the factors of 1/r D−4 cancelling each other on the two sides. This indicates that the large gauge symmetry (5.6) and the residual symmetry (3.21), acting at Coulombic order, can both be seen as consequences of Weinberg's soft theorem. This is also reflected in the observation that the Fourier transform of the soft factor occurring in Weinberg's theorem is strictly related to the memory formulae [10].