A scalar product for computing fundamental quantities in matter

We introduce a systematic way to obtain expressions for computing the amount of fundamental quantities such as helicity and angular momentum contained in static matter, given its charge and magnetization densities. The method is based on a scalar product that we put forward, which is invariant under the ten-parameter conformal group in three-dimensional Euclidean space. Such group is obtained as the static restriction (frequency $\omega=0$) of the symmetry group of Maxwell equations: The fifteen-parameter conformal group in 3+1 Minkowski spacetime. In an exemplary application, we compute the helicity and angular momentum squared stored in a magnetic Hopfion.


I. INTRODUCTION AND SUMMARY
Research in light-matter interactions benefits from theoretical tools whose generality allows one to treat light and matter in similar ways.Symmetries and conservation laws are prominent examples of such generic tools [1].For example, linear and angular momentum are physical quantities that are tied to symmetry transformations and that apply to both light and matter.These quantities have in common that they are the generators of transformations in symmetry groups that are relevant in physics, such as the Poincaré group of special relativity [2].
In classical electrodynamics, light and matter are treated in a rather similar way, namely, as continuous fields, such as the electric and magnetic fields E(r, t) and B(r, t) representing radiation, and the densities of charge ρ(r, t), current J(r, t), polarization P(r, t), and magnetization M(r, t), representing matter.Light can be treated by using the tools of Hilbert spaces [3][4][5][6][7], facilitating the consideration of material symmetries and their consequences for light upon light-matter interaction [7].The conformal invariance of Maxwell equations [8][9][10] and the corresponding conformally invariant scalar product between free electromagnetic fields [3,5,6] play an important role in such algebraic approach to electrodynamics.In particular, the scalar product allows one to obtain expressions for computing the amount of quantities such as angular momentum, energy, and helicity contained in a given radiation field [5,6,11].For the radiation field, helicity is essentially its polarization handedness.
In this article, we extend the algebraic approach towards the matter side and introduce a way to obtain expressions for computing the amount of fundamental quantities contained in static material objects from their charge and magnetization densities, ρ(r) and M(r), respectively.The densities are assumed to be spatially confined and real-valued.The expressions are obtained as a scalar product ⟨Φ|Γ|Φ⟩, where |Φ⟩ represents * ivan.fernandez-corbaton@kit.edu {ρ(r), M(r)}, and Γ is the self-adjoint (Hermitian) operator representing the particular quantity of interest.The expression of the scalar product for static matter is derived for the static fields that are bijectively connected to ρ(r) and M(r).To this end, the largest group of transformations that leave Maxwell equations with sources invariant, that is, the fifteen-parameter conformal group in 3+1 Minkowski spacetime [8,9] C 15 (3,1) is considered first, then the transformations that would not preserve the frequency ω = 0 condition of static fields are excluded.The remaining transformations also form a group, the ten-parameter conformal group in threedimensional Euclidean space C 10 (3) [12,Chapter 24].We put forward an expression for a scalar product and show that it is invariant under C 10 (3).In particular, it allows for explicit numerical computations of ⟨Φ|Γ|Φ⟩.We discuss the importance of invariant scalar products for the consistent definition of measurements and quantities such as ⟨Φ|Γ|Φ⟩, which provides the motivation for finding a relevant group of transformations and an invariant scalar product for the static case.
In an exemplary application, we set ρ(r) = 0 to focus on magnetization textures and compute the helicity and total angular momentum squared stored in a Hopfion [13][14][15][16][17] hosted inside a FeGe disk under zero external field.The definition of helicity that we use in this article [Equation (5)] is the one commonly used in optics and field theory, and can be understood as the sense of screw.This is not the same as the definition often used for magnetic solitons [18].One important difference is that the latter can be defined for two-dimensional objects, while the former needs three spatial dimensions.With the definition in Equation (5), the helicity of the Hopfion is equal to -129.1h.This number is a lower bound for the number of circularly polarized photons of positive helicity that would be needed in a helicity-dependent alloptical switch of the Hopfion to its mirror image of opposite handedness: ⌈129.1 × 2⌉ = 259.The number -129.1h also bounds the helicity that can be radiated from the Hopfion as it loses its chirality, for example by the action of a large magnetic bias aligning its magnetization density vector along the same direction at all points in the magnet.
More generally, we find that the net angular momentum and linear momentum of any {ρ(r), M(r)} along any given axis vanishes, which is consistent with the assumption of static matter.
The rest of this article is organized as follows.Section II introduces the setting and context of the work.The relevance of invariant scalar products for the consistent definition of measurements and quantities such as ⟨Φ|Γ|Φ⟩ is discussed in Section III.The scalar product for matter in static equilibrium is presented in Section IV, where its conformal invariance under C 10 (3) is shown, providing the basis for a new approach to computing properties of static material objects of finite volume from their charge and magnetization densities.In Section V we provide explicit formulas for computing the helicity, angular momentum, and angular momentum squared of a static magnetization texture M(r) and show that the net angular momentum vanishes along any given axis.In Section V B, we use the presented formalism for a quantitative study of a Hopfion in a FeGe disk.Section VI concludes the article.
The approach presented here makes the computation of fundamental quantities in matter very similar to corresponding computations for the electromagnetic field, and can be readily applied to analytically derived [13], numerically obtained [19], or experimentally measured three-dimensional charge [20] and magnetization density distributions [21].We foresee that the methodology will particularly be useful for the design and analysis of experiments involving the switching between stable states of a material system, such as, for example, when using circularly polarized light to switch the magnetization direction in magnetic films [22], which indicates a path towards much faster and more energy efficient computer memories, and whose underlying mechanisms are under intense scrutiny.

II. MOTIVATION AND PROBLEM SETTING
Figure 1 depicts a light-matter interaction sequence.A beam of electromagnetic radiation approaches a material object of finite size.Before the start of the interaction the object is in equilibrium with the radiation field.Then, the beam and the object interact for a finite period of time.When the interaction stops and equilibrium is reached again, both the beam and the object may have changed.For example, the energy, momentum, and angular momentum contained in the radiation field before and after the interaction may be different.The same can be said about the material system.
There are well-known expressions for computing the amount of a fundamental quantity such as energy or momentum contained in a given electromagnetic field.For example, in SI units which are used throughout this article, and with ε 0 and µ 0 respectively denoting the per-mittivity and permeability of vacuum, we have (1) as the energy and momentum of the field, respectively.The fields E and B are real-valued.As is explained later, we will use complex fields with positive frequencies to describe the dynamic electromagnetic field, which we denote by E and B, and for which the relation with the real-valued fields is X = X + X * .
Expressions such as those in Equation ( 1) can be derived in several different ways.For example, using conservation laws as in [23,Chapter 3], or integrating the electromagnetic stress tensor as in [24,Chapter 12.10].An alternative approach uses the tools of Hilbert spaces.In such a framework, the fields are vectors in the Hilbert space of free solutions of Maxwell equations, that is, electromagnetic fields that are not interacting with matter.Each particular solution {E(r, t), B(r, t)} corresponds to a ket |Φ⟩.The fundamental quantities are represented by self-adjoint operators that act on the kets.Then, the total amount of a given fundamental quantity Γ contained in a given electromagnetic field |Φ⟩ can be written as the scalar product of Γ|Φ⟩ and |Φ⟩: ⟨Φ|Γ|Φ⟩. ( time FIG. 1.A beam of electromagnetic radiation interacts with a material object of finite size during the grayed-out period.Before and some time after the interaction the object is in static equilibrium, where the time derivatives of all macroscopic quantities vanish.The interaction typically changes fundamental quantities of the field such as its energy or momentum.Well known formulas exist for computing such quantities for the electromagnetic field.In this article, we develop a method to compute them for the material object in static equilibrium. The expression of the scalar product for the free radiation fields, which in particular is used for obtaining explicit expressions from Equation (2), reads [3,6]: where c 0 = 1/ √ ε 0 µ 0 , and the two kets |F⟩ and G⟩ are represented by their plane wave components of well-defined helicity λ = 1 and λ = −1, F + (k) and F − (k), respectively, which correspond to left-and right-handed circular polarizations: with k • F λ (k) = 0, and, importantly, the angular frequency is restricted her to positive values ω = c 0 |k| > 0.
The exclusion of ω < 0 is possible in electromagnetism because both sides of the spectrum contain the same information [6, §3.1] [25].The ω = 0 point is also excluded from the domain of the dynamic fields.The F ± (r, t)[F λ (k)] are eigenstates of the helicity operator Λ. Helicity is the projection of the angular momentum J onto the direction of the linear momentum P. The k-space representation of Λ is particularly simple: ( The defining property of the scalar product in Equation ( 3) is that it is conformally invariant [3].That is, the value of ⟨F|G⟩ is identical to the scalar product between X|F⟩ and X|G⟩, for any transformation X in the conformal group in 3+1 Minkowski spacetime.This group is the largest group of invariance of Maxwell equations including sources as spacetime densities [8,9].In particular, free Maxwell fields transform into free Maxwell fields under the conformal group [3].The group consists of spacetime scalings, four special conformal transformations, and the Poincaré group, which consists of four spacetime translations, three Lorentz boosts, and three spatial rotations [26,27].
The quantity under the integral sign in Equation ( 3) is unitless.This can be verified by the direct computation of its units, where it should be take into account that, according to Equation (4), the units of F λ (k) are equal to the units of F ± (r, t) times m 3 .A conformally invariant scalar product must be unitless, because spacetime scalings and special conformal transformations can be interpreted as changes of units [28][29][30].In the case of the spacetime scalings, the change of units is the same for all spacetime points, while in the case of the special conformal transformations, the change varies with space and time.
Together, Equation (2) and Equation (3) are a general and convenient way of computing the amount of any fundamental quantity in the field.This motivates the following question: Can this algebraic approach be used for the material object?
An affirmative answer to this question requires an appropriate mathematical representation of matter, an appropriate group of transformations, and a scalar product that is invariant under all the group transformations.These latter two requirements are explained in the following section.

III. THE IMPORTANCE OF INVARIANT SCALAR PRODUCTS IN THE CONSISTENT DEFINITION OF MEASUREMENTS
The conformal invariance of ⟨F|G⟩ in Equation ( 3) is crucial for a consistent interpretation of measurements [31,32] which, besides quantum mechanics, can also be used in wave mechanics.In this interpretation, an observable property is represented by a self-adjoint operator Γ, and the value of the property in a ket is computed by the trace rule, which for a pure state such as |F⟩ leads to the "sandwiches" in Equation ( 2) that we consider in this paper The fact that Equation ( 6) is useful for Maxwell fields is clearly seen in e.g.[6].For example, it is there shown that the result of ⟨Φ|H|Φ⟩ for the energy operator H, and of ⟨Φ|P|Φ⟩ for the momentum operator vector P are equivalent to the corresponding integrals in Equation ( 1).Besides reproducing results typically obtained by other means, the algebraic approach for new ones to be derived, such as the content of the generators of Lorentz boosts in a given field [6,Equation (4.16)], as well as for an alternative expression for the optical helicity [11].
In order to better appreciate the importance of an invariant scalar product, let us now examine its role in the consistent definition of projective measurements.Consider a setup such as the one in Fig. 2, where the field |F⟩ is measured.We can imagine, for example, that |F⟩ is the outgoing field in the right hand side of Fig. 1, although the exact way in which the field is produced is not important.At a certain point far away from the material object, that is, in the far field, we place a measurement device consisting of an analyzer that selects a particular frequency ω 0 and polarization σ, followed by a photodetector.The number of clicks in the detector will be equal to |⟨D|F⟩| 2 , where |D⟩ is essentially a σ-polarized plane wave with momentum ω 0 /c 0 d lying in the direction connecting the origin of coordinates (assumed to be in or nearby the object) and the location of the measurement device.
Let us now consider the effects of a global transformation X applied to the whole physical system including both the field and the measurement apparatus: (7) In order for the measurement outcome |⟨D|F⟩| 2 to be meaningful, we must require that |⟨D|F⟩| 2 = |⟨F|X † X|D⟩| 2 for any X in the largest symmetry group of Maxwell equations, i.e., the conformal group C 15 (3,1).Measurement outcomes should not change under the allowed changes of reference frame or changes of units.Similarly, we must also require that quantities such as ⟨F|Γ|F⟩, which can be interpreted as the average value of the measurements of Γ on |F⟩, are also invariant under any X in the conformal group: ⟨F|Γ|F⟩ = ⟨F|X † XΓX −1 X|F⟩.
All such invariance requirements are fulfilled because any X in the conformal group is unitary with respect to the scalar product in Equation ( 3) [3], and hence X † X is the identity.
Thus, it is clear that the extension of the algebraic approach to matter requires a group of transformations and an invariant scalar product for the static case.

IV. CONFORMALLY INVARIANT SCALAR PRODUCT FOR STATIC MATTER A. Representation of matter
With respect to electromagnetism, matter in static equilibrium can be represented by its electric charge density ρ(r) and its magnetization density M(r).We assume that the time derivatives of macroscopic quantities vanish, and that there are no static currents [J(r) = 0].We also assume that ρ(r) and M(r) are contained in a finite volume.
The choice of ρ(r) and M(r) for representing matter in static equilibrium is motivated by the existence of electric charge and magnetic spin as fundamental properties of elementary particles.In sharp contrast, while the search continues [33,34], there is no experimental evidence of isolated magnetic charges or static electric dipole moments.We therefore exclude magnetic charge densities ρ mag (r) and static polarization densities P(r) from the description of matter in static equilibrium.We note that the static electric dipoles that are present in certain molecules can be described by the dipole moment of their electric charge density, and that the dynamic polarization density P(r, t) can be understood as arising from a nonstatic magnetization density.In such model M(r, t) and P(r, t) are the space-space and time-space components of a totally antisymmetric tensor, respectively, and P(r, t) vanishes in static equilibrium (P(r) = 0).Such a totally antisymmetric tensor is a common model for point particles in relativistic electrodynamics [35, Chapter II, Section 4], and has been used in studies of the effect of the electron spin on the atomic nucleus [36], and of the relativistic spin precession [37].
The densities ρ(r) and M(r) are equivalent to the static fields that they generate, as per the equations of electrostatics and magnetostatics.The scalar charge density generates the Coulomb field, which is a longitudinal (zero-curl) electric field: The vectorial magnetization density generates the B(r) and H(r) fields.Outside the material object, both fields are proportional to each other and transverse (i.e., have zero divergence).Inside the object, where M(r) ̸ = 0, the B(r) field is transverse, and the H(r) field is longitudinal: Equation ( 8) and the second line of Equation ( 9) can be obtained as the limit of Maxwell equations in static equilibrium [38, Section III], and the first line of Equation ( 9) is the definition of H(r).The expressions in Equation ( 9) can alternatively be obtained by imposing our assumption J(r) = 0 onto the typical magnetostatic equations found in e.g.[39, Equations (2.37),(2.40),(2.41)]or in [24, Equations (5.80)-(5.82)].Effectively, the equations in (9) are identifications of the transverse and longitudinal parts of M(r) with the other fields, and we may as well use only M(r) instead of both B(r) and H(r) together.Equation (8) implies that E(r) and ρ(r) determine each other bijectively.We choose to use E(r) here because it shares some transformation properties with the magnetization M(r), allowing us to avoid different derivations in some cases, should we use ρ(r) instead.We can therefore represent matter in static equilibrium by means of E(r) and M(r), or alternatively by means of their Fourier transforms E(k) and M(k), which can be obtained by integrals in the finite volume V occupied by the object: where the first equality in the second line of Equation ( 10) follows from Equation (8).We note that both E(r) and M(r) are real-valued.
B. Group of transformations and invariant scalar product for the static case In this section, we define an appropriate scalar product for static matter, which, in particular, allows one to use Equation ( 2) for computing the amount of fundamental quantities stored in matter.With such scalar product, each {ρ(r), M(r)} corresponds to a ket |Φ ω=0 ⟩ in the Hilbert space of static matter.This Hilbert space is different from the one containing the radiation fields.One salient difference is the frequency ω, which is equal to zero for static matter and strictly larger than zero in radiation fields.Another difference is that there are helicity zero (longitudinal) components in static matter, while the radiation fields are always transverse.
Motivated by the requirement of invariant measurements formulated in Section III, we pursue the idea of a group invariance for the scalar product by starting with the full conformal invariance of Maxwell equations with sources shown by Bateman [8], Cunningham [9] and Dirac [10], among others [30].While in Dirac's work the sources were electric charge-current densities, both Bateman and Cunningham additionally included magnetization and polarization densities.
The conformal invariance including sources deserves further discussion.On the one hand, the invariance of Maxwell equations with sources under the 15 parameter conformal group in Minkowski spacetime C 15 (3, 1) is often explicitly recognized in the literature; e.g., see [27, Section 3.1], or the review by Kastrup [30] where the conformal symmetry, its role in theoretical physics, and their historical evolution are comprehensively explained.On the other hand, the fact that fixed mass parameters break the scale invariance, and consequently the conformal invariance, could suggest the incorrect conclusion that Maxwell equations are conformally invariant only in the source-free case.In this respect, it is important to note that Maxwell equations including electric chargecurrent densities [ρ(r, t), J(r, t)], and magnetization and polarization densities [M(r, t), P(r, t)] do not feature any fixed mass parameter.The invariance under conformal transformations follows from the way in which the sources transform.The transformations of [ρ(r, t), J(r, t)] can be found in [27,Equations 3.40ab].A simple illustration is the transformation law of the electric charge density under a scaling r → αr: ρ(r) → ρ(r/α)/α 3 in the static case, which is readily shown to preserve the total charge in a volume.The bottom line is that when both fields and sources are transformed, the form of the dynamic Maxwell equations with sources remains invariant.This is what we use below as the starting point for obtaining the sought after expression of a scalar product for the static case.Incidentally, even equations with a mass parameter, such as the Dirac equation, can be shown to be conformally invariant if one allows for a particular rescaling of the mass [30,Section 5.1].
Let us now return to the idea of invariance for the scalar product.Since we are considering matter in static equilibrium, we need to remove all the transformations of C 15 (3, 1) that do not preserve the ω = 0 condition.Lorentz boosts change ω and mix it with the components of k.The time component of the special conformal transformations four-vector also changes ω [26, Section 3].We remove all these transformations.We also remove timetranslations because, while preserving the ω = 0 condition, any time-translations will, for the static case, just degenerate into the identity operator.We are then left with a ten parameter group consisting of spatial translations, spatial rotations, the spatial scaling (r → αr with α ∈ R, α > 0 ), and three special conformal transformations: The ten parameter conformal group C 10 (3) in three-dimensional Euclidean space [12,Chapter 24].It is interesting to see that the static restriction of the conformal group in 3+1 spacetime dimensions results in the conformal group in three spatial dimensions.Accordingly, along with the discussion in Section III, our sought-after scalar product must be invariant under all the transformations in C 10 (3).
We now proceed by considering the C 15 (3, 1)-invariant scalar product expression in Equation ( 3).We first perform a unitary change of basis, going from the F ± (k) to the { √ ϵ 0 E(k), iB(k)/ √ µ 0 } basis by using Equation ( 4) to express the latter as linear combinations of the former. 1 where I is the 3×3 identity matrix.We now replace B(k) by µ 0 M(k).Both quantities have the same units.This replacement is necessary in order to include the longitudinal degree of freedom that M(k) can contain, which is always absent in B(k): Equations ( 3), (11), and (12) motivate us to write the following scalar product expression for the ω = 0 fields: As is the case in Equation ( 3), the quantity under the integral sign in Equation ( 13) is unitless.This is readily seen with the following equalities between the units of different fields: .
When written in the equivalent form in Equation ( 28), Equation ( 13) coincides with the invariant scalar product for the relevant representations of the scale-Euclidean group in [40,Equation (29)].The extra 1/|k| 2 in Equation (29) of [40] with respect to Equation (28) in this paper is compensated by the extra 1/|k| factor in the plane wave decomposition in [40,Equation (57)] as compared to Equation (10).That N = −2 in [40, Equation (57)] follows from the definition of N in [40, Equation (25a)], and the transformation properties of E(r) and M(r) under spatial scalings.Therefore, the scalar product in Equation ( 13) is invariant under the scale-Euclidean group, which is a seven-parameter sub-group of C 10 (3) composed by translations, rotations, and spatial scalings.These transformations act hence unitarily with respect to the scalar product.The invariance of the scalar product in Equation ( 13) under the rest of C 10 (3), that is, under special conformal transformations, can be seen as follows.
The special conformal transformations act on the coordinate vector as [12, Equation (24.4)]: where c is a real-valued 3-vector: c l ∈ R, l = 1, 2, 3. To prove that Equation ( 13) is invariant under such transformations one must show that the special conformal transformations are unitary.This unitary character is shown by considering that the special conformal transformation C l (c l ) along axis l with parameter c l can be obtained from the spatial translation along such axis by the same parameter T l (c l ), and the inversion operation R [12, below Equation (24.4)]: The action of R on the coordinate vector is Since we already known that the translations act unitarily, Equation (15) implies that and if R is unitary, then C l is unitary as well.The unitary character of R is shown in App.A by extending a result contained in Section II of [41].
Once R is known to be unitary, the conclusion that the scalar product in Equation ( 13) is invariant under C 10 (3) can also be reached directly without using its invariance under the scale-Euclidean group.First, we establish that the action of spatial translations leave Equation (13) invariant, which is obvious when substituting the known action of arbitrary translations by a displacement vector c on the Fourier transforms E(k) and M(k) onto Equation (13), where the acquired phases cancel.
Since the translations are unitary, their generators P l are self-adjoint, and so are the generators of the special conformal transformations K l , because of Equation (15).The other generators of C 10 (3), that is, the angular momentum operators and the generator of dilations D, can be obtained as commutators of K l and P i [41, Equation (20)]: Then, because the adjoint of the commutator of two self-adjoint operators is minus itself, it follows, noting that i → −i when taking the adjoint of the right hand side of Equation (18), that D and M li must be self-adjoint, and hence the transformations that they generate must be unitary.We recall that M ll = 0, M li = −M il , M 12 = J 3 , M 23 = J 1 , and M 31 = J 2 .
We adopt Equation ( 13) as the C 10 (3)-invariant scalar product for ω = 0.This provides a new way of obtaining expressions for computing the amount of fundamental quantities stored in matter.
Conveniently, the very similar functional form of the expressions for the scalar products in the dynamic and static cases, Equations ( 3) and ( 13) respectively, implies that known k-space expressions of fundamental operators for the ω > 0 case [5,6,42] can be also used for the ω = 0 case.The same is true for the corresponding r-space expressions.
We note that, while the scalar products for the dynamic case in Equation ( 3) and for the static case in Equation ( 13) are invariant under C 15 (3, 1) and C 10 (3), respectively, the dynamic fields and static densities are not required to exhibit specific symmetries.
In the rest of the article, we will focus on fundamental quantities stored in the magnetization density, for which we set E(r) = 0 in Equation ( 13): (19) Nevertheless, the methodology that we use applies to the case with E(r) ̸ = 0 as well.For example, the representations of M(r) using three scalar complex functions in the domain of linear momentum f λ=−1,0,1 (k), or angular momentum f jm(λ=−1,0,1) (|k|), where λ is the helicity, have their counterparts for representing E(r), except that E(r) contains only the λ = 0 component.The counterpart of Equation ( 25) is and Equations ( 30),( 32), (33), and (34) apply, with the difference that the helicity takes only the value λ = 0.

V. HELICITY AND ANGULAR MOMENTUM IN MAGNETIZATION
We start with the helicity stored in M(r).At each k point, M(k) can be decomposed into the three eigenstates of the helicity operator for vectorial fields: The longitudinal λ = 0 component corresponds to nonvanishing ∇ • M(r).While a net magnetic charge in an isolated object has not been observed, complicated magnetization textures are expected to contain pairs of monopoles of opposite charge [14], and pairs of singularities have indeed been experimentally imaged [21].The longitudinal λ = 0 term in the magnetization is crucial for their description.In contrast, the free dynamic electromagnetic field is divergenceless, and the λ = 0 component vanishes.
In particular, the e λ ( k) used by Bialynicki-Birula [5,6] are related to the Q λ ( k) vectors used by Moses [40,43] in the following way: Using the expression of the helicity operator in Equation ( 5), the helicity stored in M(r) can then be written as: Equations (22,24) = h We note that the charge density cannot store helicity, because i k × E(k) = i k × −i kρ(k) ε0|k| = 0.As shown in [38], the result of Equation ( 26) is proportional to the definition of the static magnetic helicity [44][45][46]: where ∇ × A(r) = B(r), and hence ik × A(k) = B(k).More generally, Equations ( 21), (24), and ( 25) can be used to write the expression of the scalar product in Equation ( 19) as a function of the f λ (k) (28) which allows one to compute ⟨M|Γ|M⟩ from the f λ (k), as long as the action of Γ on f λ (k) is known.
Let us now turn our attention to angular momentum.Rather than using f λ (k), the computation of ⟨M|J i |M⟩ is more conveniently performed in the angular momentum basis: where |k j m λ⟩ denotes a simultaneous eigenstate of P • P = P 2 , J • J = J 2 , J z , and Λ, with eigenvalues h2 |k| 2 = h2 k 2 , h2 j(j + 1), hm, and hλ, respectively, and where the Wigner D-matrices for spatial rotations [47,Chapter 4] enter as follows As we show in App.B, the scalar product ⟨M 1 |M 2 ⟩ can also be written as: The action of angular momenta operators on f jmλ (|k|) is [49, Equations (2.1)-(2. 3)]: with which we can write: Appendix C shows that the angular momenta ⟨M|J τ ∈{x,y,z} |M⟩ and the linear momenta ⟨M|P τ ∈{x,y,z} |M⟩ both vanish.The origin of this vanishing can be found in the following restrictions, which follow from assuming M(r) ∈ R 3 : For example, for ⟨M|J z |M⟩, the proof is straightforward because the sum over m vanishes since the second line of Equation (35) Derivations similar to those in App.C show that the linear and angular momenta in ρ(r)[E(r)] vanish as well since, for E(r) ∈ R 3 : The vanishing of linear and angular momenta for the whole |Φ ω=0 ⟩ is consistent with the static equilibrium condition under zero external fields.In contrast, the squared linear and angular momenta ⟨M|P 2 |M⟩, and ⟨M|J 2 |M⟩, do not vanish [see e.g. the last line of Equation (34)].We note that the non-zero mean value of the square of a property simultaneously with zero mean value of such property can easily occur, and should not be taken as a sign of a quantum effect.Since the eigenvalues of self-adjoint operators are real numbers, the eigenvalues of squares of self-adjoint operators are non-negative real numbers.Then, the value of a property measured as the average of a self-adjoint operator for a given state can be zero, while, simultaneously, the average of the squared operator is non-zero.For example: The average joint linear momentum of two counterpropagating planewaves is zero in any direction, yet the linear momentum squared (P • P) is not.
The f λ (k) and f jmλ (|k|) for a given M(r) can be readily obtained with the following sequence of computations, M(r) → f jmλ (|k|).(37) which can be applied to analytically derived, numerically obtained, or experimentally measured three-dimensional magnetization textures confined to finite volumes.
We note that, given ρ(r), the f 0 (k) and f jm0 (|k|) corresponding to E(k) can be obtained similarly.

A. Discrepancies with existing results
The computation of the linear and angular momenta of magnetic solitons [50][51][52][53][54][55] remains an active field of research [56][57][58][59].Starting from the Landau-Liftshitz (LL) equation, different integral expressions have been proposed, mostly for predicting the dynamics of topologically non-trivial magnetization textures inside an unbounded magnetic medium under, for example, an external magnetic field.Nevertheless, such integrals are assumed to also be valid for the static case without external fields.Using one of the latest versions of such integrals, even a completely static domain wall has a finite linear momentum [58].With our methodology, as shown in App C 2, a static magnetization density has a vanishing value of average linear momentum and angular momentum along any given axis.The difference in the results can be attributed to the fact that the two approaches are applied to systems meeting quite different assumptions.In this article, we consider a bounded three-dimensional domain containing a given static three-dimensional magnetization.In contrast, the LL-based approaches consider a two or three-dimensional domain of infinite size where the magnetization is allowed to vary, as for example when three-dimensional topological textures are stabilized by endless motion with constant velocity [53].

B. Helicity and angular momentum of a Hopfion
As an exemplary application, we will now compute the helicity and angular momentum squared stored in a Hopfion.
A convenient analytical approximation of the magnetization density of a Hopfion has been provided by P. Sutcliffe in [13], where the Hopfion is hosted in a disk of FeGe of height L=70 nm and diameter 3L, and is numerically stabilized to a cylindrically symmetric configuration using the Dzyaloshinskii-Moriya interaction of the chiral magnet and a strong anisotropy perpendicular to the flat ends of the cylinder.With good approximation to the numerical results, the Cartesian components of unit magnetization density vector m(r) of the Hopfion are given in [13,Equation (3.3)] as a function of the cylindrical coordinates of the position vector r ≡ [ρ, θ, z] = [ x 2 + y 2 , atan2(y, x), z].The magnitude of M(r) is assumed to be constant: M(r) = M s m(r).A value of M s =384 kA/m is assumed here.
For the numerical calculations, the cylindrical r domain of the Hopfion was discretized with [107,65,71] points in cylindrical coordinates, and the k space with ⟨M|Λ|M⟩ ⟨M|J 2 |M⟩ -129.1 h 1.30×10 3 h2 TABLE I. Helicity Λ, and angular momentum squared J 2 stored in a Hopfion computed with the last lines of Equation (26) and Equation (34), respectively.An analytical approximation of the Hopfion in a chiral FeGe magnet of cylindrical shape [13] was used in the calculations.The height of the cylinder is equal to the magnetic helical period L, and the diameter is equal to 3L.A magnetization density saturation value of Ms=384 kA/m is assumed.
950 points for |k| between 0 and 20L −1 , and with 2592 points for k at each |k|.The multipolar orders from j = 1 to j = 9 were considered.
All the expected outcomes such as the conditions in Equation (35), or the vanishing of the momenta and angular momenta, are numerically verified for the Hopfion up to numerical inaccuracies at the level of the fifth significant digit.
Table I contains the values for the helicity and angular momentum squared stored in the Hopfion computed with the last lines of Equation (26) and Equation (34), respectively.The stored helicity is equivalent to ≈ 129 righthanded circularly polarized photons, which is about ten orders of magnitude smaller than the number of photons in a circularly polarized femtosecond laser pulse of 10 mJ cm −2 fluence at a central wavelength of 800 nm.
The value of helicity, equal to -129.1h, implies a lower bound for the number of circularly polarized photons that would be needed in a helicity-dependent all optical switching of the Hopfion onto its mirror image of opposite helicity: ⌈129.1 × 2⌉ = 259.Additionally, -129.1h also bounds the helicity that can be radiated by the Hopfion as it loses its chirality, for example by the action of a large magnetic bias aligning its magnetization density vector along the same direction at all points.
Appendix D contains further analysis of the angular momentum content of the Hopfion.In particular, we show that the Hopfion is an eigenstate of J z with angular momentum zero.

VI. CONCLUSION
We have introduced a new way to obtain expressions for the computation of the fundamental quantities in static matter from its charge and magnetization densities.The method is based on a scalar product obtained from requirements of invariance under the ten-parameter conformal group in three-dimensional Euclidean space C 10 (3).This group is obtained as the static (ω = 0) restriction of the symmetry group of Maxwell equations with sources, namely, the fifteen-parameter conformal group in 3+1 Minkowski spacetime.
In an exemplary application, we have used the for-malism to compute the angular momentum squared and helicity stored in a Hopfion inside a FeGe disk.We foresee that this methodology will in particular be useful for the design and analysis of experiments involving the switching between stable states of a material system.
to expand any φ(k) meeting Equation (A2) as: (A5) The unitary character of R under the scalar product follows then readily from the third line of Equation (A5), and from ⟨e nj m|e njm ⟩ = δ nn δj j δ mm .The result that R is unitary can be extended to the f λ (k) functions used in this article as follows.
The angular momentum basis is again convenient for the task.The properties of the Wigner D-matrices imply that the relation inverse to Equation ( 30) is We see that Equation (A7) is very similar to Equation (A3).Importantly, the angular and radial dependences are separated into is somewhat more involved.We start by using Equation (B1) to write: C5) and manipulate the last term and re-labeling the variable m → −m comparing Equation (C8) with Equation (B1) makes it clear that the condition f λ * in Equation (C3) implies Equation (C4).

Angular momentum
The proof that ⟨M|J z |M⟩ = 0 can be found in the main text, just below Equation (35).
Let us now consider ⟨M|J y − iJ x |M⟩ in Equation (34), and show that The steps for the corresponding proof that ⟨M|J y + iJ x |M⟩=0 in the second line of Equation ( 34) are very similar.

Linear momentum
In k-space, the action of the momentum operator P τ ∈{x,y,z} is P τ M(k) = hk τ M(k).Then, the momentum of M(r) in direction τ can be written:  This appendix contains further analysis of the angular momentum content of the Hopfion.Some insight can be gained even without having the f jmλ (|k|) at hand.For example, we now show that the Hopfion magnetization density is an eigenstate of J z with eigenvalue zero.
We start by slightly rewriting the expression in [13,Equation (3.3)] with the help of a rotation matrix:

(D2)
We now apply apply the J z of [61,Equation (5.43)] to the Hopfion.In cylindrical coordinates for r but Cartesian components for the vector, J z has the following expres-sion: Since M s does not depend on θ, we can just apply J z to R z (θ) mθ=0 : The result of Equation (D4) is faithfully reproduced by the f jmλ (|k|) amplitudes obtained from the sequence of computations indicated in Equation (37). Figure 3 shows the distribution of the different helicities across the values of m.We highlight that the values for m ̸ = 0 are not suppressed from the plot, but they are not visible in this scale.The sign of the helicity of the Hopfion in Tab.I can be deduced from the larger value of the negative helicity component in Fig. 3.

FIG. 2 .
FIG.2.In a projective measurement, the outcome of an apparatus for measuring the field |F⟩ can be modeled as |⟨D|F⟩|2 , that is, the modulus square of the projection of |F⟩ onto an electromagnetic mode |D⟩.
Appendix D: Angular momentum content of the Hopfion

j 2 5. 8 ×FIG. 3 .
FIG.3.Angular momentum content of the Hopfion for each helicity computed with the formula on the top right corner.Such expression gives the partial norm squared of the Hopfion in each (m, λ) subspace, and is therefore a conformally invariant unitless number.While the numerical errors in the calculations produce m ̸ = 0 components, they are not visible in this scale because they are much smaller than the m = 0 components, as indicated by the inequality.For clarity, two of the color bars are horizontally offset from the m = 0 point.