On the dynamic toroidal multipoles from localized electric current distributions.

We analyze the dynamic toroidal multipoles and prove that they do not have an independent physical meaning with respect to their interaction with electromagnetic waves. We analytically show how the split into electric and toroidal parts causes the appearance of non-radiative components in each of the two parts. These non-radiative components, which cancel each other when both parts are summed, preclude the separate determination of each part by means of measurements of the radiation from the source or of its coupling to external electromagnetic waves. In other words, there is no toroidal radiation or independent toroidal electromagnetic coupling. The formal meaning of the toroidal multipoles is clear in our derivations. They are the higher order terms of an expansion of the multipolar coefficients of electric parity with respect to the electromagnetic size of the source.

Introduction

After the introduction of toroidal multipoles in the static case[1], the dynamic toroidal multipoles were presented as a new independent multipole family that had been previously ignored[2, 3]. The new family was derived from the split of the transverse multipoles of electric parity into two parts. These parts are often referred to as electric and toroidal. The dynamic toroidal multipoles have also been analyzed in refs 4 and 5, and are currently being considered in the areas of metamaterials, plasmonics, and nanophotonics[6-19].

In this article, we analyze the split between the dynamic electric and toroidal multipoles of a localized source distribution. We find that the two parts cannot be separately determined by measuring the electromagnetic fields produced by the source outside the source region. The two parts can also not be separately determined by measuring the coupling between the source and externally incident electromagnetic waves. The result applies to both near and far field situations, and implies that there is no independent coupling of toroidal character between the sources and the electromagnetic field, and that there is no radiation of pure toroidal character. The toroidal multipoles can hence not be considered an independent family. Rather, our analysis makes clear that the electric and toroidal parts correspond, respectively, to the lowest and the higher order terms of the expansion of the exact multipolar coefficients of electric parity with respect to the electromagnetic size of the source. In this respect, our article contains analytical expressions that inherently contain all correction orders: They are exact for any source size.

Our results do not question the usefulness of considering the higher order terms. For example, some experimental results can only be explained by adding the toroidal dipole contributions to the lowest order terms (see Fig. 3c in ref. 7, and Fig. 4 in ref. 16). The same is true for non-radiating spherical configurations (see Fig. 2 in ref. 20). These cases are examples where the toroidal terms are the dominant dipolar terms.

In the following, we first outline the mathematical setting in which we carry out the analysis. We then highlight the direct connections between our momentum space approach and the physics of multipolar couplings between sources and fields. On the one hand, three kinds of degrees of freedom are needed to describe the sources. Two of them are transverse and have opposite (electric vs. magnetic) parity, and the other one is longitudinal. On the other hand, only the two transverse degrees of freedom are needed to describe the electromagnetic fields produced by the sources outside the source region, where the field produced by the longitudinal degrees of freedom is identically zero.

Then, we analyze the physical meaning of the toroidal multipoles. Our analysis starts from recently obtained exact analytical expressions for the transverse multipolar coefficients of electric parity. These expressions are valid for any source size[21]. We take the dipolar case as the guiding example, make the small source approximation keeping the two lowest orders, and show that: The lowest order term is the well known approximate expression for the electric dipole of electromagnetically small sources; the term of second lowest order is the approximate expression for the toroidal dipole of electromagnetically small sources [Eq. (2.11) in ref. 3]. We then analyze the original split[2, 3] in the general multipolar case, and establish that both electric and toroidal parts contain non-radiative components, which appear due to the separation of terms of different order in the size of the source. These non-radiative components, which cancel each other when both parts are summed, preclude the separate determination of each part by measuring the radiation from the source outside the source region, or by measuring the coupling between the source and externally incident electromagnetic waves. Afterwards, we analyze a different class of splits where the two parts are free of this kind of non-radiative components. We find that, in this case, what precludes their independent measurement is the absence of longitudinal radiation. We show in Supp. Inf. Sec. IV that this is the same reason that precludes the independent measurement of the recently proposed re-definition of toroidal multipoles (Box 2 in ref. 19).

Finally, we use the obtained insights to clarify some statements that are often found in the literature.

Mathematical setting

We start with an electric current density distribution J(r, t) embedded in an isotropic homogeneous medium characterized by real valued permittivity and permeability μ. We assume that J(r, t) is confined in space so that J(r, t) = 0 for |r| > R. We consider its energy-momentum Fourier representationand treat each frequency ω separately. The lower limit of the integral in dω excludes the static case ω = 0, which we do not treat in this paper. This formal setting is independent of the cause of the current. For example, in a scattering situation the current is induced in the scatterer by an external field.

The electromagnetic fields radiated by the source at each frequency ω are solely determined by the part of J (p) in the domain that satisfies . This theorem has been proven in refs 22 and 23. The authors showed that this part of the current completely determines the electromagnetic field generated by the source at any point r such that |r| > R, which, in particular, includes near fields. This result has been extended beyond spheres to confining volumes of any smooth enough shape (see Chap. 9 in ref. 24): The part of J (p) in the domain that satisfies determines the electromagnetic field generated by the source at any point outside the source region, for a wide class of shapes of such source region. Figure 1 depicts a spatially confined monochromatic source distribution and the fields that it generates outside the source region.

Figure 1

The monochromatic electric current density distribution J (r) is confined to a sphere of radius R: The source region. In this article, we consider the electromagnetic fields produced by the source outside its region, i.e. [E (r), H (r)] for |r| > R. These fields are completely determined by the Fourier components of the current density which meet |p| = ω/c (see Fig. 2). The other components do not produce fields outside the source region. This is a general result that is also valid when the confining volume defining the source region is not spherical (see text).

Figure 2

Only the transverse components of the current density J (p) with produce electromagnetic fields outside the source region (see Fig. 1). (a1) Spherical shell of radius in momentum space. The (a2) region depicts two different orthonormal bases for vector functions defined on the shell. The and the . The (j, m) indexes are omitted in the figure, and the subscript in the Y’s is the difference between their second and first indexes. The two bases are related to each other as written in Eq. (7). With respect to their polarization vectors: , , and are orthogonal (transverse) to the momentum vector p, is parallel (longitudinal) to p, and and lay in the plane and are neither transverse nor longitudinal, as seen in (a3). The transverse vectors are represented by dashed blue arrows, the longitudinal ones by solid black arrows and the ones of mixed character by dotted red arrows. For the figure, we took j = 1 in Eq. (7) for the relationships between and .

We denote by the components of J (p) in the spherical shell of radius |p| = ω/c as a function of the solid angle [see region (a1) in Fig. 2].

In this domain, the relevant scalar product is , where † denotes conjugate transpose, and runs over the entire spherical shell.

We expand in an orthonormal basis for vector functions defined in a spherical shell:The basis is composed by the three families of vector multipolar functions in momentum space (Sec. B.3 in ref. 25)where the are the scalar spherical harmonics and the three components of the vector L are the angular momentum operators for scalar functions.

Each of the functions in Eq. (3) is an eigenstate of the angular momentum squared J 2, the angular momentum along one axis , for which we choose , and the parity operator Π. As depicted in Fig. 2, the polarization of and is orthogonal (transverse) to , and the polarization of is parallel (longitudinal) to [see region (a2) in Fig. 2]. In coordinate space (r), this distinction corresponds to the distinction between divergence free (transverse) and curl free (longitudinal) vector fields. Table 1 contains the eigenvalues for the aforementioned three operators and the transverse (⊥) or longitudinal () character of each vector multipolar function.

Table 1

Vector multipolar functions: Polarization character and eigenvalues of J 2, J and parity Π.

	J 2	J z	Π	Polarization
X_jm(pˆ)	j(j + 1)	m	(−1)j+1	⊥
Z_jm(pˆ)	j(j + 1)	m	(−1)j	⊥
W_jm(pˆ)	j(j + 1)	m	(−1)j	∥

The symbol ⊥ means transverse polarization. The symbol means longitudinal polarization. Both j and m are integers, and . For and , j > 0. For , j ≥ 0.

We provide expressions for the in Supp. Info. Sec. I, which we derived in ref. 21.

Multipolar coupling between sources and fields in momentum space

The decomposition in Eq. (2) is useful for discussing the physics of multipolar couplings between the sources and the fields. The complex scalars determine the coupling of the source to the transverse electromagnetic field, a.k.a radiation field, a.k.a electromagnetic waves. The one-to-one relationship between the coefficients of the source and the multipolar fields of frequency ω radiated by it is exposed in App. B of ref. 26, Chap. 9.10 in ref. 27, and ref. 23, for example, where the analysis is made in the r domain. The represent independent degrees of freedom in the following sense. Each individual coefficient is connected to the radiation (absorption) of a particular kind of electromagnetic wave: A transverse multipolar field with well defined frequency, angular momentum, and parity. The are connected to the multipolar waves of electric parity and the to those of magnetic parity (see Table 1). The multipolar waves connected to different coefficients are orthogonal to each other. In principle, each individual can be independently measured through the interaction of the source with the radiation field.

The represent the longitudinal degrees of freedom of , but the longitudinal electric field with |p| = ω/c is zero outside the source region defined in Fig. 1. This can be explained in the following way. The longitudinal components of the current () generate a longitudinal electric field outside the source region defined in Fig. 1, but such field is canceled by the longitudinal electric field produced by the associated charge density. This statement follows from the continuity equation, as we show in Supp. Info. Sec. II. The cancellation can also be recognized in §13.3 p1875–1877 in ref. 28, and in the expression of the electric field in Eq. (2.39) in ref. 5, where the two contributions to the longitudinal field cancel exactly. In ref. 5 the longitudinal wave functions are .

The bottom line is that the contain the necessary and sufficient information to determine both the fields of frequency ω produced by the sources outside the source region, and the coupling of the sources to external electromagnetic waves of such frequency.

The lack of longitudinal radiation means that, even though there are three kinds of degrees of freedom in the current, , only two of them are observable from outside the source region. The powerful theory developed by Gustav Mie for spherical scatterers[29] is probably the simplest example of this. In this case, the current density is induced in the sphere by an external incident field. The of this induced current are determined by the Mie coefficients, which can be obtained analytically thanks to the symmetry of the problem. Using the Mie coefficients, the field produced by the induced current can be determined exactly everywhere outside the sphere. This allows the exact calculation of near fields, which is crucial in many applications like for example: Determination of the transition rates of atoms placed near spheres[30], determination of the forces and torques caused by evanescent fields onto a sphere near a substrate[31], modeling the tips of scanning tunnelling optical microscopes[32], and general optical trapping analysis[33].

The dynamic toroidal multipoles

We will now study the split of the multipolar coefficients of electric parity () which lead to the appearance of the dynamic toroidal multipoles[2, 3]. We use the dipolar case j = 1 as our guiding example and then generalize the findings to any j ≥ 1.

We first show that the well known expressions for the electric and toroidal dipoles of small sources are, respectively, the lowest and second lowest order terms in a series expansion of the exact expression of .

Formal meaning

We start from the recently obtained exact expression for the dipolar vector of electric parity in the spherical vector basis (Eq. (21) in ref. 21):where k = ω/c, and † means conjugate transpose. Let us now make the small argument approximation to the spherical Bessel functions with terms up to second order [ and ], and then group the contributions with the same power of k.

We find that the (k)0 term, in Eq. (5), is the well known small source approximation of the electric dipole moment of a current density distribution [Eq. (9.14) in ref. 27]. The k 2 term, in Eq. (6), is the small source approximation of the toroidal dipole moment [Eq. (2.11) in ref. 3]. The integral in Eq. (6) contains contributions from the two integrals in Eq. (4): It is the sum of the terms of order k 2 coming from both j 0(kr) in the first integral of Eq. (4) and j 2(kr) in the second.

It is now clear that is nothing but the next to leading order term in the electromagnetically small source approximation of the exact expression of . All higher order terms are easily obtained from the Taylor series of the spherical Bessel functions in Eq. (4).

The same structure revealed in Eqs (4)–(6) underlies the general case of the electric/toroidal split for j ≥ 1. In refs 2 and 3, the lowest order term in , of order (k), is isolated as in Eq. (5). The rest, in our example, is called toroidal multipole form factor. Then, the lowest order term of the toroidal multipole form factor, of order (k) [like our Eq. (6)], is called toroidal multipole moment. For simplicity, we refer to the toroidal multipole form factor as the toroidal part. It corresponds to the toroidal multipolar family as defined in refs 2 and 3. In the dipolar example, the toroidal part () is equal to plus all other higher order corrections coming from the terms (k) of the expansion of the spherical Bessel functions in Eq. (4).

The formal origin of the toroidal multipoles is recognized in the literature[4, 15, 18, 20]. We now investigate their physical meaning.

Physical meaning

The splitting of into two parts[2, 3] has a physically relevant consequence. Both parts, often called electric and toroidal, contain Fourier components of the current from outside the |p| = ω/c shell. The |p| ≠ ω/c components do not couple to electromagnetic waves[23]. For example, it follows from the properties of the Fourier transform that the result of the integral in the expression of the electric part in Eq. (5) is proportional to J (p = 0). Since |p| = 0 ≠ ω/c, it follows that the toroidal part, , must also contain a contribution proportional to J (p = 0). This is so because does not have any out of shell contribution [Eq. (4.12a) in ref. 23]: The out of shell contributions in the electric part must be canceled by those in the toroidal part when the two parts are summed. Since the |p| ≠ ω/c components do not couple to the transverse electromagnetic field, it is impossible to determine the toroidal (electric) parts of by measuring the fields produced by the source outside its region, or by externally exciting the source with electromagnetic waves. As previously discussed, this conclusion does not depend on whether the measurement(excitation) occurs in(from) the near, mid, or far field zones, as long as they occur outside the source region. For the sake of discussion, let us now consider the hypothetical case where the first term in the Taylor series of a given component of is zero at a particular frequency. We make no judgment about whether this hypothesis is physically realizable. In this case, the toroidal part would be equal to the said component of and, as such, it would couple to the transverse electromagnetic field. Nevertheless, in this hypothetical case the statement that the toroidal part couples to the field means nothing else than that couples to the field, and it would not amount to separate coupling of the electric and toroidal parts.

The same conclusion applies to any value of j ≥ 1, as we show in Supp. Inf. Sec. III. The root cause is the breaking up of the spherical Bessel functions in the split[2, 3]. Their lowest order term, proportional to (kr), goes to the electric part and the rest of the terms go to the toroidal part. When unsplit, the spherical Bessel functions in exact integral expressions like Eq. (4) or Supp. Inf. Eq. (2) completely reject the |p| ≠ ω/c components of J (r). The Fourier transforms of spherical Bessel functions correspond to momentum space radial deltas δ(|p| − ω/c) (see end of Sec. III in ref. 21). After the split, the term (kr) by itself does not provide such rejection. This allows |p| ≠ ω/c components to “leak into” the electric part, and to be forcefully present in the toroidal part as well since the sum of the two parts does not contain any out of shell contribution.

We have hence established that the electric and toroidal parts cannot be separately determined by measuring the electromagnetic fields produced by the source outside its region, or by measuring its coupling to external electromagnetic waves. Our result implies that, as opposed to electric and magnetic multipoles whose distinct character is given by their parity, there is no electromagnetic field of pure toroidal character. It also implies that toroidal multipoles cannot be independently excited, i.e. there is no independent toroidal coupling. This precludes the existence of selection rules of toroidal character. The two different parities (electric and magnetic), and their linear combinations are necessary and sufficient to characterize electromagnetic transitions (see e.g. Chap. 9.8 p. 436 in ref. 27, and Chap. XII 2. B in ref. 26). The lack of longitudinal coupling together with the binary character of the parity operator leave no room for a third option.

While the have an independent physical meaning, the individual electric and toroidal parts do not have it. The toroidal multipoles cannot be considered an independent family. The same can be said about the isolated electric part, or about any split of the terms in the Taylor series of the spherical Bessel functions in Supp. Info. Eq. (2) or (4). Any such split is bound to introduce non-radiative out of shell components in each of its parts.

Since the toroidal multipoles are higher order terms in the approximation of physically meaningful quantities, they can be used to improve analytical models that use only the lowest order terms. The inclusion of the next to leading order term improves the accuracy and prediction ability of such models, as can be seen in refs 7, 11, 12, 14–18. In ref. 20, and in the context of non-radiating source configurations, the inclusion of the second order term explains a zero in the induced , which cannot be explained by the first order term alone (see Fig. 2 in ref. 20). The improvement of the models will be maximized by using expressions like Eq. (4), Eq. (20) in ref. 21, and Supp. Inf. Eq. (2). These expressions inherently contain all correction orders because they are exact for any size of the source.

Splits without out of shell components

We have seen that the out of shell non-radiative components appear when splitting by separating terms of different orders in the size of the source. Nevertheless, we observe that the in Supp. Inf. Eq. (2) [ in Eq. (4)] have two formally distinct contributions, and that each of them contains a spherical Bessel function which guarantees the restriction to |p| = ω/c, thereby avoiding the intrusion of the out of shell components. However, we will now show that, due to the absence of longitudinal radiation, the two on-shell contributions cannot be separately determined. We consider another orthonormal basis for vector functions on the momentum shell: The vector spherical harmonics . They are related to the asThe first and third lines of Eq. (7) can be found in Eq. (16.91) and Exercise 16.4.4 of ref. 34, respectively. The second line can be obtained from the equation . The result is reached by using the expression of the orbital angular momentum operator in momentum space , the formula , and the gradient formula Eq. (16.94) in ref. 34.

Figure 2 illustrates Eq. (7). Table 2 contains information about the . Notably, we deduce from Eq. (7) and see in Fig. 2 that, while is transverse, and are of mixed character, i.e., neither transverse nor longitudinal.

Table 2

Vector spherical harmonics: Polarization character and eigenvalues of J 2, J , Π, and L 2.

	J 2	J z	Π	L 2	Polarization
Y_j,j,m(pˆ)	j(j + 1)	m	(−1)j + 1	j(j + 1)	⊥
Y_j,j−1,m(pˆ)	j(j + 1)	m	(−1)j	(j−1)j	mixed
Y_j,j+1,m(pˆ)	j(j + 1)	m	(−1)j	(j + 1) (j + 2)	mixed

L 2 is the orbital angular momentum squared. Both j ≥ 0 and are integers. For j = 0, , Y 0,0,0 = 0, and Y 0,−1,0 is not defined.

We can expand as:It follows from Eq. (7) that:We invert Eq. (9) to obtain:The current dependent integrals in r space that define contain only the spherical Bessel function of order l. This can be shown using Eq. (10) and the expressions for in Supp. Inf. Eq. (2). For example, in the dipole case of Eq. (4), the integrand that contains j 0(kr) is only due to , and the integrand that contains j 2(kr) is only due to . Both terms have zero out of shell components.

Let us assume that we want to determine and by measuring the electromagnetic radiation from the source at points r such that |r| > R. As is clear from Eq. (10), and can be visually appreciated in Fig. 2, the contain both transverse and longitudinal components. Their measurement involves the detection of both transverse and longitudinal fields generated by the current distribution. But, as previously discussed the longitudinal field due to the current distribution outside the source region is exactly canceled by that due to the charge distribution. Without access to the longitudinal degrees of freedom of , it is impossible to separate the two contributions to .

We see that there are splits of that, while both parts remain free of out of shell components, they contain longitudinal components which render them non-separable. This is the case for a definition of the toroidal multipoles that has been recently given in Box 2 of ref. 19. In Supp. Info. Sec. IV, we first show that this definition is different from the original definition in refs 2 and 3. In the split of the proposed in Box 2 of ref. 19, the integrands contain entire spherical Bessel functions which prevent the appearance of out of shell components. We then show that each of the two parts contains longitudinal components that cannot be detected from outside the source region.

Discussion and conclusion

To finalize, we discuss some statements that are often found in the toroidal literature.

It is often stated that the toroidal multipole family is a third family, independent of the electric and magnetic ones[5, 7, 13, 19]. Originally, the authors of ref. 2 suggested in Sec. 4 to replace the by , where are the coefficients of the toroidal part of . We have shown that the cannot be separately determined, and can hence not constitute a third independent multipolar family.

It is also often stated that three multipolar families, including the toroidal one, are needed in order to expand a general charge-current distribution, while the radiation that they emit can be described with only two multipolar families[5, 14, 35]. Starting from Eq. (2), our analysis makes clear that, indeed, three multipolar families are needed to describe the source, and only two are needed to describe the fields radiated by it. It also makes clear that the third family is constituted by the longitudinal multipoles, and not the toroidal multipoles. The fact that the longitudinal field is zero outside the sources explains why the number of necessary multipolar families decreases by one when going from the sources to the fields.

Finally, it is also stated in refs 5 and 19 that the toroidal multipoles can be separated using spectroscopic techniques. The separation is to be enabled by exploiting the extra frequency dependent weighting of the toroidal multipoles [k 2 = (ω/c)2] with respect to the electric multipoles. It is straightforward to prove that, should this be actually possible, there would exist a doubly infinite number of independent multipolar families, and not just a third one. This unphysical outcome shows that the conjectured separability is not possible. In order to show this, let us assume that it is possible. The steps going from Eqs (4)–(6) reveal that the extra k 2 scaling comes from different terms in the series expansion of the spherical Bessel functions. We stopped at the second lowest order, but each of the infinite number of terms in the complete expansion has an extra k 2 weight with respect to its predecessor. Additionally, one can make a similar expansion of the magnetic multipoles (see Sec. 4 in ref. 21). If the separate measurement of the terms with additional k 2 factors where possible, as conjectured in Box 2 of ref. 19, there would be not just an extra independent multipolar family, but a doubly infinite number of them. The physical reality is that each of this doubly infinite number of formally obtained terms is not an independent degree of freedom. The multipolar electromagnetic transitions that determine the spectroscopic response of a system of charges and currents are completely characterized by using the two possible parities, electric and magnetic. The spectroscopic measurements of already contain all the available information. On the other hand, should the analysis of the spectroscopic data make use of models limited to the lowest order, the inclusion of the second lowest order terms can only improve it, and the use of exact expressions like Eq. (4), Supp. Inf. Eq. (2), and Eq. (20) in ref. 21 will optimize it.

In conclusion: Our analysis proves that the dynamic toroidal multipoles do not have an independent physical meaning with respect to their interaction with electromagnetic waves. Their formal meaning is clear, however: They are higher order terms of an expansion of the transverse multipolar coefficients of electric parity in the electromagnetic size of the source.