Conserved charges for rational electromagnetic knots

We revisit a newfound construction of rational electromagnetic knots based on the conformal correspondence between Minkowski space and a finite S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S^3$$\end{document}-cylinder. We present here a more direct approach for this conformal correspondence based on Carter–Penrose transformation that avoids a detour to de Sitter space. The Maxwell equations can be analytically solved on the cylinder in terms of S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S^3$$\end{document} harmonics Yj;m,n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_{j;m,n}$$\end{document}, which can then be transformed into Minkowski coordinates using the conformal map. The resultant “knot basis” electromagnetic field configurations have non-trivial topology in that their field lines form closed knots. We consider finite, complex linear combinations of these knot-basis solutions for a fixed spin j and compute all the 15 conserved Noether charges associated with the conformal group. We find that the scalar charges either vanish or are proportional to the energy. For the non-vanishing vector charges, we find a nice geometric structure that facilitates computation of their spherical components as well. We present analytic results for all charges for up to j=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j=1$$\end{document}. We demonstrate possible applications of our findings through some known previous results.


Introduction and summary
Topologically non-trivial solutions to source-free Maxwell's equations on Minkowski space R 1,3 were first discovered by Rañada [1] based on the Hopf fibration. Other methods for constructing such linked and knotted rational solutions involving Bateman's complex Euler potentials, conformal inversion and Penrose twistors have been found since then (see [2] for a review). Key features characterizing these solutions are helicity (conserved under conformal transformations) and Noether charges associated with the conformal group SO (2,4) of R 1,3 [3,4]. Quite interestingly, some simple EM knot configurations, including the one with figure-8 topology, have been produced in the laboratory using laser beams with knotted polarization singularities [5]. With such a discovery, this line of research no longer remains just a theoretical endeavor; more complicated knot configurations discussed here, or elsewhere, could be realized experimentally in the near future. Trajectories of charged particles in the background of these knotted electromagnetic fields have also been studied recently [6].
It was shown in [7] that a complete family of such rational knot solutions on R 1,3 , labelled by the so-called left-right S 3 harmonics Y j;m,n , can be constructed; any finite-action rational Maxwell solution can be expanded in terms of these basis solutions. The construction relies on the conformal correspondence between the de Sitter space dS 4 and R 1,3 via a finite Lorentzian S 3 -cylinder over conformal time τ . Here, we establish this conformal correspondence between the cylinder and R 1,3 more directly using Carter-Penrose transformation. These "knot basis" solutions are first obtained on the cylinder by employing the SO(4) isometry of S 3 coupled with the right-multiplication action of SU (2) on S 3 (being the group manifold of SU (2)) and are later transferred to R 1,3 via the conformal map. In another paper [8], co-authored by one of us, general linear combination of such "knot basis" solutions in terms of arbitrary complex coefficients j,m,n was considered and, for a fixed j, its helicity was found to be proportional to the conserved energy.
Here, we reconsider these generic solutions for a fixed j and proceed to compute the remaining charges associated with the conformal group, namely momenta, angular momenta, boosts, dilatation and special conformal transformation (SCT) charges. The charge densities are evaluated on de Sitter space at τ = 0 (without loss of generality, since they are conserved) where considerable simplifications occur demonstrating the usefulness of this "de Sitter method". It is found that the dilatation charge vanishes while the scalar SCT charge V 0 is proportional to the energy E. Furthermore, the boost charges K i vanish and the vector SCT charges V i are proportional to the momentum charges P i . Interestingly, for the vector charge densities viz. momenta p i , angular momenta l i and vector SCT v i we find that a one-form constructed with them on the spatial slice R 3 → R 1,3 , e.g. p i dx i is proportional to a similar one on de Sitter space. This correspondence allows us to compute additional charges ( p r , p θ , p φ ) by the action of such one-forms on spherical vector fields (∂ r , ∂ θ , ∂ φ ). At j = 0 it turns out that there are only four independent nonzero charges: the energy E and the momenta P i . The situation for higher spin j is more complicated, but some of the components of the charges in spherical coordinates are found to vanish for arbitrary j. The action of the so(3) generators D a on the indices of j,m,n can easily be a e-mail: kaushlendra.kumar@itp.uni-hannover.de (corresponding author) obtained for a fixed j owing to the SO(4) isometry. This allows for an action of these generators on the charges. For (the Cartesian components of) the vector charges this action is found to inherit the original so(3) Lie algebraic structure, as expected. We also compute the correct coefficients 0;0,n corresponding to two very interesting generalizations of the Hopfian solution obtained via Bateman's construction in [4], which allow us to validate our generic formulae of these charges. Further works in this direction could be to relate these new knot solutions and their conserved charges to other known field configurations in [4], such as plane waves or to study Fourier analysis of general electro-magnetic knot configurations to understand single mode behaviour.
The paper is organized as follows. In Section 2, we first establish the conformal map between Minkowski space and the Lorentzian S 3 -cylinder using Carter-Penrose transformation and review the necessary computational tools related to S 3 , including its harmonics. We then construct the rational Maxwell solutions on the cylinder in terms of the S 3 harmonics and illustrate some knotted field configurations using figures. Section 3, at first, reviews the Noether current and associated charges arising from the action of the conformal group on free Maxwell theory and presents a computational strategy. The 15 conformal charges are then computed in four subsections and the action of so(3) generators D a on some of the vector charges is analyzed. Finally, in Section 4 we validate the formulae of the charges computed in Section 3 using other known results.

Construction of rational electromagnetic knots
We present below the construction of rational electromagnetic knots by completely skipping the detour to de Sitter space in contrast with [8] , where these solutions were first presented. For a clear presentation, we split the construction in three parts.

Minkowski to cylinder: Carter-Penrose transformation
The metric on the Minkowski space R 1,3 in polar coordinates is given by We first employ the lightcone coordinates (u, v) to transform the metric in the following way In the second step, we compactify the spacetime with the help of the coordinate (U,V) as follows: Finally, we rotate the coordinate system back using (τ, χ) to obtain the desired form of the metric: We realize that the Minkowski metric is conformally equivalent to the metric on a Lorentzian cylinder I × S 3 ; I = (−π, π) with a conformal factor that can be recasted in terms of Minkowski coordinates using the above transformations (2.2-2.4): where we have made γ dimensionless using a global scale factor , chosen to be the de Sitter radius so as to make the connection with previous works apparent. We depict in Figure 1 the causal structure of the Minkowski spacetime (preserved under a conformal transformation) using the Penrose diagram. Note that the spatial infinity i 0 (where t = 0 and r → ∞) and the future/past temporal infinities i ± (where r = 0 and t → ±∞) are not included in the figure. Furthermore, the boundaries of the Penrose diagram are (a) r = 0 line, (b) past null infinity I − corresponding to r − t → ∞, r + t = constant and (c) future null infinity I + corresponding to r + t → ∞, r − t = constant. The natural embedding coordinates ω A with A = 1, 2, 3, 4 of S 3 → R 4 satisfying ω 2 A = 1 has a natural parametrization in terms of hyperspherical coordinates (χ, θ, φ) of (2.4): ω a = sin χω a , ω 4 = cos χ, with ω 1 ,ω 2 ,ω 3 = (sin θ cos φ, sin θ sin φ, cos θ ) .
It is not difficult to obtain the following useful relations γ t = sin τ and γ r = sin χ, (2.7) which can be used to find the Jacobian of transformation between the coordinates y m ∈ {τ, χ, θ, φ} and x μ ∈ {t, r, θ, φ}: (2.8)

Structure of three-sphere and its harmonics
The presence of S 3 in the metric (2.5) is an added advantage, which can be exploited to write a SO(4)-invariant vector potential and obtain the corresponding electromagnetic fields by solving Maxwell's equations on the Lorentzian cylinder. These quantities can be later exported to the Minkowski spacetime owing to the conformal invariance of the vacuum Maxwell equations in four dimensions.
To this end, we start with the group SO(4), which is isomorphic to two copies of SU (2) (up to a Z 2 grading). Each of these SU (2) has a group action that generates a left (right) multiplication (a.k.a. translation) on S 3 . This can easily be checked using the map This parametrization of g ensures that the identity element e = 1 2 of the group SU (2) can be obtained from (0, 0, 0, 1), i.e. the North pole of S 3 . It is well known that S 3 is the group manifold of SU (2). Keeping this in mind, we consider the Cartan one-form being the SU (2) generators. The left-invariant one-form 1 e a can, alternatively, be expressed using the so-called self-dual 't Hooft symbol η a BC : They satisfy the following useful identities δ ab e a e b = d 2 3 and de a + ε abc e b ∧ e c = 0. (2.12) The left-invariant vector fields L a generating the right translations are dual to e a and are given by In a similar way, the right-invariant vector fields R a generating the left translations are given by where the anti-self-dual 't Hooft symbolsη a BC are obtained from (2.11) by flipping the B, C = 4 sign. The stability subgroup

SO(3) ⊂ SO(4) under the group action on the coset space SO(4)/SO(3) S 3 is generated by
(2.15) Furthermore, the vector fields L a and R a act on the one-forms e a via their Lie derivative 2 yielding L a e b = 2 ε abc e c and R a e b = 0. (2.16) We can now write the differential d of the functions f ∈ C ∞ I × S 3 using L a as Furthermore, the restriction of such functions on S 3 has a natural expansion in terms of left-right harmonics Y j;m,n with 2 j ∈ N 0 and m, n ∈ {− j, − j + 1, ..., j} that are eigenfunctions of Casimirs 3 L 2 := L a L a , R 2 := R a R a , and operators L 3 , R 3 : (2.18) The action of the corresponding ladder operators on the harmonics Y j;m,n is, therefore, given by The normalized harmonics Y j;m,n can be expanded in terms of functions α, β and their complex conjugates as They satisfy the orthonormality condition These adjoint harmonics can be written in terms of the standard S 2 spherical harmonics Y l,M and the associated Legendre polynomials of the first kind P b a : We can write down the cylinder one-forms in Minkowski coordinates using relations (2.6-2.8) in (2.11) (2.25) We further note down the expressions for the vector fields L a in terms of S 3 angles by using (2.6) in (2.13) for later purposes: (2.26)

Rational solutions to vacuum Maxwell equations
We start with the following ansatz [8] for a real-valued Maxwell gauge potential A on the Lorentzian cylinder which can be extended to the Minkowski space since the Maxwell's theory is conformally invariant here. We further impose the following "temporal+Coulomb" gauge on this space Note that the first part in (2.28) is not the usual temporal gauge A t = 0. In fact, we can make use of the inverse Jacobian (2.8) while promoting the gauge potential to the Minkowski The field strength F is easily obtained with this gauge potential using (2.12): The corresponding vacuum Maxwell equations d * F = 0 read Solutions to these coupled linear wave equations on S 3 were obtained in [7] as an infinite tower of basis solutions labelled by j and divided into two categories: type-I with frequency j = 2( j+1) and type-II with frequency j = 2 j. We shall confine our discussions in this paper to type-I solutions. The corresponding results for type-II solutions can easily be obtained using a duality transformation (see [7]) where the basis solutions of type-I at level j are related with the basis solutions of type-II at level j+1 via the exchange m ↔ n. We take a general complex linear combination of these basis solutions and obtain (see [8]) the real-valued and (2.38) It may be noted here that, due to the compactness of the Lorentzian cylinder, the sphere-frame energy and action are always finite. One can obtain, by a cumbersome calculation, the electric/magnetic fields in Minkowski space by using the results (2.25) in (2.36) and realizing that The actual expression for electric/magnetic fields are rather involved and it bodes well to combine them as the Riemann-Silberstein vector S: The expression of S becomes more and more complicated with increasing j. 4 The (self) knotted field lines corresponding to such electromagnetic configurations also increase in complexity with j. We illustrate this point by producing field lines for j = 0, 1/2 and 1 in Figure 2.

Conformal group and Noether charges
It is well known [9] that free Maxwell theory on R 1,3 arising from the action is invariant under the conformal group SO (2,4). Furthermore, the above action is also invariant under the gauge transformations: where the vector fields ξ μ obey the conformal Killing equations: The conserved Noether current J μ is obtained by equating the "on-shell variation" of the action where the variations δ A μ are arbitrary and the fields A μ satisfies the Euler-Lagrange equations with its "symmetry variation" where the fields A μ are arbitrary but variations δ A μ satisfy the symmetry condition. The correct variation δ A μ is obtained by imposing the gauge invariance on the Lie derivative of A μ w.r.t. the vector field ξ μ : Finally, the conserved current is obtained as which satisfies the continuity equation and gives the conserved (in time) charge Q 5 : Before proceeding further, a couple of remarks pertaining to the subsequent calculations are in order: • All the charges Q are computed at t = 0 owing to the simple J 0 expressions on this time-slice. To that end, we record the following useful identities at t = τ = 0: (3.6) Moreover, the electromagnetic fields at t = 0 are given in terms of tetrads e τ = e τ μ dx μ , e a = e a μ dx μ : using the two-form (2.36) and its Minkowski counterpart (2.39). • Furthermore, these charges Q are computed for a fixed spin j and, thus, we will suppress the index j from now onwards, unless necessary. Note that the sphere-frame EM fields for fixed j can be obtained by using the expansion (2.32) in (2.38) as withZ a denoting the complex conjugate of Z a . • The simplifications below for the charge density J 0 are carried out using the harmonic expansion (2.35) while using the results (3.6) in (3.7). • We frequently use below the well-known fact that an odd {ω A } integral over S 3 vanishes because of the opposite contributions coming from the antipodal points on the sphere. In particular, it can be checked that the following integral vanishes 6 : Having made these remarks, we now proceed to compute the charges Q for various conformal transformations ξ μ obeying (3.2) in following four categories.

Translations
An easily seen solution to (3.2) is the set of four constant translations which also partly generates the Poincaré group and give rise to the usual stress-energy tensor of electrodynamics corresponding to the μ-component of the translation for an arbitrary ν . The corresponding charges are the energy E and the momentum P.
Energy. The expression of the energy density e := T 00 simplifies to which, in turn, simplifies the expression for the energy E to The resultant expression for E in terms of complex parameters m,n and its conjugate was already computed in [8]. The ω 4 part of the above integral vanishes as discussed in an earlier remark and we obtain (3.14) Table 1 Expressions of P 3 , P φ and L 3 for j = 1/2 and 1 Momentum. For the momentum densities p i = T 0i = (E × B) i we obtain an interesting correspondence relating the one-form p := p i dx i on Minkowski space with a similar one on de Sitter space: A straightforward calculation then yields the expression of momenta P i : with e a i given by (3.6). The results for j = 0 are (3.17) As a consistency requirement, we check that the vector charges P i are rotated according to the algebra of D a (2.15) (see appendix A): We also note down P 3 for j = 1/2 and 1 in Table 1. One can compute the corresponding P 1 and P 2 for j = 1/2 and j = 1 by employing the action of an appropriate D a of Table 3. We can additionally compute the spherical components of the momentum (P r , P θ , P φ ) by letting the one-form e a in (3.15) act on the vector fields (∂ r , ∂ θ , ∂ φ ). In practice, we first write ∂ r = − γ ∂ χ using (2.8) at t = 0 and then invert the vector fields (∂ r , ∂ θ , ∂ φ ) in terms of the left-invariant vector fields (L 1 , L 2 , L 3 ) using (2.26). Finally, using the duality relation e a (L b ) = δ a b we obtain From the expression of P r we see that the integrand over S 2 , i.e.ω a P a (2.6) is an odd function, 7 which makes P r vanish. We also find with explicit calculations (verified for up to j = 1) that P θ vanishes. For j = 0 we find that P φ is proportional to P 3 : (3.20) The expressions of P φ for j = 1/2 and 1 are recorded in Table 1. 7 The terms of P a using (2.24) are proportional to Y l,M Y l ,M , which is even inω a for a fixed j.

Lorentz transformations
Another solution of (3.2) is given by (3.21) which correspond to the six generators of the Lorentz group SO (1, 3). These six together with the above four translations generates the full Poincaré group. The corresponding six charges for these transformations are grouped into the boost K and the angular momentum L.
Boost. The conserved charge densities arising from J 0 (3.4) corresponding to 0i are the boost densities k = ρ x − p t, which simplify for t = 0 to The corresponding charges K i vanish because of the odd integrand as discussed in an earlier remark: Angular momentum. The other three conserved charge densities J 0 corresponding to i j in (3.4) are the angular momentum densities l = p × x, which take a simple form just like in momentum (3.15): The expressions for the charges L i simplify to Explicit calculations show that for j = 0 the angular momenta L i are proportional to the momenta P i : This is, however, not true for higher spin j. The angular momenta L i has the same rotation behaviour as for the momenta P i (3.18). We therefore note down the results of L 3 for j = 1/2 and 1 in Table 1, from which the corresponding expressions of L 1 and L 2 can be obtained using Table 3. We can again compute the spherical components of the angular momentum (L r , L θ , L φ ) using the relations (3.19) by replacing P with L in it. We realize that the S 2 integrand for L r would only have terms likeω aωb P c , which are all odd functions 8 of φ and, therefore, the φ integral over the domain (0, 2π) would make L r vanish. We also find, with explicit computations, that the charges L θ and L φ for j = 0 are proportional to P 3 : As a non-trivial example, we collect these charges for j = 1/2 below:

Dilatation
It is easy to verify that a constant rescaling by λ: is also a solution of (3.2). The charge density corresponding to this single generator of the conformal group is p · x − e t, which for t = 0 simplifies to The corresponding charge D vanishes because of the odd integrand:

Special conformal transformations
A fairly straightforward calculation shows that the following not so obvious transformation also satisfies (3.2). The four generators corresponding to b μ give rise to four different charges V 0 and V.
The expression for the corresponding charge V 0 takes the following simple form: Here again the ω 4 term of the integral, being odd, vanishes and yields These simplify at t = 0 and take a structure similar to the momentum densities p i (3.15): The expressions for the charges V i then simplify to With explicit computation we observe that the charges V i are proportional to the momenta P i (verified explicitly for up to j = 1) As before, we can compute the spherical components (V r , V θ , V φ ) by using the expressions (3.19) and replacing P with V in it. We notice that the charge V r vanishes owing to the odd S 2 integrand 9 just like in the case of P r . However, unlike P θ here the charge V θ is non-vanishing. Explicit calculations show that for j = 0 the charges V θ and V φ are proportional to the momentum P 3 : Additionally, we record below the charges V θ and V φ for the non-trivial case of j = 1/2 One can compute these charges for higher spin j by following the same strategy. 9 Observe that the terms in V a (3.37) are all even inω a .   1 4 sech θ E R

Applications
The method of constructing rational electromagnetic fields presented in this paper has the added advantage that it produces a complete set labelled by ( j, m, n); any electromagnetic field configuration having finite energy can, in principle, be obtained from an expansion like in (2.32-2.34), albeit with a varying j. The operational difficulty involved in this procedure has to do with the fact that this set is infinite as j ∈ N 2 . There are, however, many important cases where only a finite number of knot-basis solutions (sometimes only with a fixed j) need to be combined to get the desired EM field configuration. One such very important case is that of Hopfian solution given by Ranãda, which was already covered in [7]. Below we analyse two very interesting generalizations of Hopfian solution presented in [4] in the context of present construction. It is imperative to note here that while the scope of construction of a new solution from the known ones as presented in [4] is limited, the same is not true for the method presented in this paper, which by design can produce arbitrary number of new field configurations. Some of these possible new field configurations obtained from the j = 0 sector (possibly from j = 1/2 or 1 as well) could find experimental application like, e.g. in [5].
Acknowledgements KK is grateful to Deutscher Akademischer Austauschdienst (DAAD) for the doctoral research grant 57381412. We thank Olaf Lechtenfeld for several insightful discussions and valuable suggestions.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

A Rotation of indices
By construction the gauge potential A is SO (4)  where D a (e b ) are determined from (2.16) while D a (Y j;m,n ) are determined from (2.18-2.20). By collecting the coefficients of various linearly independent e b and Y j;m,n terms in the above expansion for a fixed D a , one gets a set of coupled linear equations for D a ( m,ñ ), which can be easily solved. The action of the generators D a on m,ñ for j = 0, 1/2 and 1 is given in Table 3.