Articles | Volume 41, issue 1
https://doi.org/10.5194/angeo-41-253-2023
https://doi.org/10.5194/angeo-41-253-2023
Regular paper
 | 
31 May 2023
Regular paper |  | 31 May 2023

Concerning the detection of electromagnetic knot structures in space plasmas using the wave telescope technique

Simon Toepfer, Karl-Heinz Glassmeier, and Uwe Motschmann
Abstract

The wave telescope technique is broadly established in the analysis of spacecraft data and serves as a bridge between local measurements and the global picture of spatial structures. The technique is originally based on plane waves and has been extended to spherical waves, phase-shifted waves and planetary magnetic field representation. The goal of the present study is the extension of the wave telescope technique using electromagnetic knot structures as a basis. As the knots are an exact solution of Maxwell's equations they open the door for a new modeling and interpretation of magnetospheric structures, such as plasmoids.

Dates
1 Introduction

The classification and mathematical modeling of spatial structures are among the major missions of theoretical physics. Our extraterrestrial space environment in particular provides a diversity of spatial structures with different characteristics. For example, oscillating structures can be classified into plane waves (e.g., MHD waves), spherical waves generated at the bow shock, surface waves triggered by instabilities at the magnetopause and phase-shifted waves caused by field line resonances (Plaschke et al.2008; Narita et al.2022). On the other hand, global planetary magnetic fields can be interpreted in terms of a multi-pole series based on spherical harmonics (Gauss1839; Glassmeier and Tsurutani2014; Toepfer et al.2020a, b, 2021). For the characterization of such structures, empirical models, such as magnetospheric models or models based on a set of specific basis functions spanning the solution space of differential equations, are required.

In general, any spatial structure can be expanded into a set of mathematical basis functions, such as plane waves or spherical harmonics. Plane waves are the simplest spatial structures forming a basis for the representation of spatial fields. The contribution of any plane wave with its characteristic spatial scale to the total field is described by the spectrum of the field. However, in the worst case, infinitely many elements forming the basis have to be incorporated to describe the structure, resulting in an infinite set of expansion coefficients that have to be determined from the measurements. In this case, it is desirable to choose a new representation based on a new set of basis functions that are well-adjusted to the symmetry of the structure with fewer unknown parameters.

Electromagnetic knots, proposed by Cameron (2018), are a special superposition of infinitely many plane waves, forming such a new basis set for localized, divergence-free structures, namely the electromagnetic ring and the electromagnetic globule. The geometry of these basis elements is depicted in Fig. 1a and b. A variety of electromagnetic field topologies can be constructed by spatially distributing and superposing several rings and globules as illustrated in Fig. 1c. The complexity of the emerging field geometries prompts the naming electromagnetic knots (Cameron2018).

https://angeo.copernicus.org/articles/41/253/2023/angeo-41-253-2023-f01

Figure 1Vector representation of the electromagnetic ring (a), the electromagnetic globule (b) and spatially distributed, superposed globules (c) after Cameron (2018).

The electromagnetic ring and the electromagnetic globule are an exact solution of Maxwell's equations and provide a new tool in the context of plasma physical and electrodynamical modeling. Based on the elaboration of Cameron (2018), the mathematical foundations of electromagnetic knots are revisited in the present study. Within this context, the formalism is reformulated in terms of the classical wave telescope technique (Motschmann et al.1996). Additionally, the applicability of describing and interpreting spatial structures in planetary magnetospheres via knots is discussed. The wave telescope technique enables the classification of spatial structures in planetary magnetospheres from a limited number of satellite positions and has successfully been applied to several problems in space physics (Glassmeier et al.2001; Narita et al.2003, 2009, 2013, 2022). Originally, the method was based on a plane wave representation and was later extended to spherical waves (Contantinescu et al.2006), phase-shifted waves (Plaschke et al.2008) and planetary magnetic fields (Narita2019; Toepfer et al.2020a, b). The goal of the present study is the extension of the variety of spatial structures that can be analyzed from a limited set of measurement positions by considering the electromagnetic knots a new basis set for the wave telescope. The method is tested against synthetically generated magnetic field data describing a plasmoid as a two-dimensional magnetic ring structure.

2 The classical wave telescope

Maxwell's equations represent a set of coupled partial differential equations for the magnetic field B(x,t) and the electric field E(x,t). These equations can be transformed into a set of algebraic equations via the Fourier transform. In the following discussion we will focus on the magnetic field.

The measurement position x and the measured field B(x,t) are known from a set of magnetometer measurements. Due to the high temporal resolution of the magnetometer, the temporal Fourier transform can be applied to the data, delivering the spectral amplitude B(x,ω) (Motschmann et al.1996). In general, this spectral amplitude is a continuous function of ω. However, in the practical application outstanding points of the spectrum, for example sharp maxima, are of major interest. Thus, the data are evaluated at a peak, where ω=ω0, with the corresponding amplitude B(x,ω0). So far, the magnetic field can be written as

(1) B x , ω 0 = B ^ 0 k , ω 0 e i k x d 3 k ,

where B^0(k,ω0) is the spectral amplitude of the magnetic field with respect to the wave vector k. As the magnetic field measurements are solely available at a limited number of measurement points, the spatial Fourier transform is not applicable. Thus, the spectral amplitudes B^0(k,ω0) and the corresponding wave vectors k are to be determined by the data fitting procedure. Although a variety of inversion techniques are available (Haykin2014, e.g.,), we will focus on the wave telescope technique (Motschmann et al.1996).

Suppose that the magnetic field vector B(x,ω0) is measured at N positions xi (i=1,,N), summarized into the 3N-dimensional vector B(ω0). Thus, the determination of the spectral amplitude B^0(k,ω0) results in an overdetermined inversion problem. Following Motschmann et al. (1996), Narita (2019) and Toepfer et al. (2020b), the magnetic field model can be rewritten as

(2) B ω 0 = H ( k ) B ^ 0 k , ω 0 d 3 k ,

where

(3) H ( k ) = I e i k x 1 I e i k x N R 3 N × 3

is the shape matrix and IR3×3 denotes the identity matrix. The magnetic field measurements can be arranged into the data covariance matrix

M=Bω0Bω0R3N×3N,

where the angular brackets denote the statistical average of the data. The spectrum of the wave can be estimated via

(4) P ( k ) = tr H M - 1 H - 1 ,

where the dagger denotes the Hermitian conjugate and trHM-1H-1 is the trace of the matrix HM-1H-1. The maximum values of P(k) may be interpreted as the spectrum of the field. If only a finite number of sharp peaks emerges, the magnetic field may be interpreted as a superposition of plane waves with discrete k values. As P(k) is a nonlinear function of the vector k, the whole three-dimensional k space needs to be scanned to identify the peaks (Motschmann et al.1996).

3 Electromagnetic knots

The classical wave telescope technique does not assume any symmetry or relation between different k vectors of the spectrum. However, to be able to use electromagnetic knots as a system of basis structures, the geometry of the k space needs to be specialized. In this respect, the classical wave telescope technique differs from its extension presented here. The following mathematical derivation of electromagnetic knots is based on Cameron (2018).

3.1 Construction of the knots

For the specific evaluation of the integral in Eq. (1), spherical coordinates (k,φ,θ) in the k space are introduced:

(5) k = k sin θ cos φ sin θ sin φ cos θ = : k e k ,

where the corresponding unit vectors are given by

ek=sinθcosφex+sinθsinφey+cosθez,eφ=-sinφex+cosφey,eθ=cosθcosφex+cosθsinφey-sinθez.

The vectors ex, ey and ez denote the unit vectors of the Cartesian coordinate system.

In this case, the magnetic field in Eq. (1) can be rewritten as

(6) B x , ω 0 = B ^ 0 k , ω 0 e i k x d 3 k = B ^ 0 k , θ , φ , ω 0 e i k e k x d 3 k = 0 0 2 π 0 π B ^ 0 k , θ , φ , ω 0 e i k e k x k 2 sin θ d θ d φ d k = 0 0 2 π 0 π B ^ 0 k , θ , φ , ω 0 e i k e k x sin θ d θ d φ k 2 d k = 0 b ̃ k , ω 0 , x k 2 d k ,

where

(7) b ̃ k , ω 0 , x = 0 2 π 0 π B ^ 0 k , θ , φ , ω 0 e i k e k x sin θ d θ d φ

is the spectral amount of the field corresponding to k.

Due to Maxwell's equations, the magnetic field (as well as the electric field in the absence of free charge carriers) is solenoidal,

(8) x B ( x , t ) = 0 ,

such that

(9) B ^ 0 k , θ , φ , ω 0 e k = 0 .

To guarantee the solenoidality of the magnetic field, the ansatz

(10) B ^ 0 k , θ , φ , ω 0 = α k , φ , θ , ω 0 e φ + β k , φ , θ , ω 0 e θ

is chosen, which results in

(11) b ̃ k , ω 0 , x = 0 2 π 0 π [ α k , φ , θ , ω 0 e φ + β k , φ , θ , ω 0 e θ ] e i k e k x sin θ d θ d φ

and

(12) B x , ω 0 = 0 0 2 π 0 π [ α k , φ , θ , ω 0 e φ + β k , φ , θ , ω 0 e θ ] e i k e k x sin θ d θ d φ k 2 d k ,

where α(k,φ,θ,ω0) and β(k,φ,θ,ω0) are complex functions of (k,φ,θ) and ω0. In the following this ansatz is specified by constraining the geometry of the three-dimensional k space.

Equation (11) represents the spectral amplitude of the magnetic field for a fixed value of k. Thus, it is useful to separate the angular dependency (φ,θ) of the spectral amplitude from the k dependency by choosing the functions α and β as

(13)αk,φ,θ,ω0=αφ,θ,ω0B02πK(k)k2,(14)βk,φ,θ,ω0=βφ,θ,ω0B02πK(k)k2,

where 𝒦(k) is a function of k alone, α(φ,θ,ω0) and β(φ,θ,ω0) are complex functions of (φ,θ,ω0), and B0 is a real constant.

In this respect, the spectral amplitude (Eq. 11) can be rewritten as

(15) b ̃ k , ω 0 , x = B 0 2 π 0 2 π 0 π [ α ( φ , θ ) e φ + β ( φ , θ ) e θ ] e i k e k x sin θ d θ d φ K ( k ) k 2 ,

where the functions α(φ,θ) and β(φ,θ) weight the summation over the k space with respect to the angulars φ and θ. Introducing the abbreviation

(16) B ̃ k , ω 0 , x = B 0 2 π 0 2 π 0 π [ α ( φ , θ ) e φ + β ( φ , θ ) e θ ] e i k e k x sin θ d θ d φ

provides

(17) b ̃ k , ω 0 , x = B ̃ k , ω 0 , x K ( k ) k 2

such that

(18) B x , ω 0 = 0 B ̃ k , ω 0 , x K ( k ) d k .

In the following, the functions α(φ,θ) and β(φ,θ) are specified to evaluate the spectral amplitude B̃(k,ω0,x) with regard to electromagnetic knots (Cameron2018).

Each spectral amount (corresponding to a fixed k value) of the field may be characterized by a superposition of plane waves with the same amplitude propagating in every direction (independent of φ and θ) such that

(19) α ( φ , θ ) = α 0 = const . C , β ( φ , θ ) = β 0 = const . C .

In this case, the spectral amplitude results in

(20) B ̃ ( k , ω 0 , x ) = B 0 2 π 0 2 π 0 π α 0 e φ + β 0 e θ e i k e k x sin θ d θ d φ ,

representing a superposition of infinitely many plane waves of the same amplitude with the spectrum

(21) S k = k R 3 | k x 2 + k y 2 + k z 2 = k 2 .

Therefore, the distribution in k space is completely characterized by the value k.

Using the definitions of the unit vectors eφ and eθ, the magnetic field can be further expanded into the form

(22) B ̃ ( k , ω 0 , x ) = B 0 2 π { α 0 0 2 π 0 π [ - sin φ e x + cos φ e y ] e i k e k x sin θ d θ d φ + β 0 0 2 π 0 π [ cos θ cos φ e x + cos θ sin φ e y - sin θ e z ] e i k e k x sin θ d θ d φ } = B 0 2 π { e x 0 2 π 0 π ( - α 0 sin φ + β 0 cos θ cos φ ) e i k e k x sin θ d θ d φ + e y 0 2 π 0 π ( α 0 cos φ + β 0 cos θ sin φ ) e i k e k x sin θ d θ d φ - β 0 e z 0 2 π 0 π sin 2 θ e i k e k x d θ d φ } .

For the evaluation of the integrals in Eq. (22) it is useful to introduce a cylindrical coordinate system (ρ,ϕ,z) in the position space:

(23) x = ρ cos ϕ ρ sin ϕ z = ρ e ρ + z e z ,

where ρ=x2+y2. The corresponding unit vectors are given by

eρ=cosϕex+sinϕey,eϕ=-sinϕex+cosϕey,ez=ez.

The scalar product of the k vector and the position vector results in

(24) e k x = ρ cos ϕ cos φ sin θ + ρ sin ϕ sin φ sin θ + z cos θ = ρ sin θ ( cos ϕ cos φ + sin ϕ sin φ ) + z cos θ .

Using x=ρcos ϕ and y=ρsin ϕ provides

(25) e k x = sin θ ( x cos φ + y sin φ ) + z cos θ .

For the further evaluation of the integrals in each component of Eq. (22), the abbreviations

η1(θ):=kxsinθandη2(θ):=kysinθ

are introduced. By means of these preparations, the φ integration can be solved analytically, delivering the Bessel functions of the first kind:

Jn(x)=12π-ππei(nτ-xsinτ)dτ.

The detailed evaluation of the integrals can be found in the Appendix, resulting in

(26) B ̃ k , ω 0 , x = Re B 0 i α 0 f e ϕ + β 0 g e ρ + h e z

and

(27) B ( x , ω 0 ) = 0 B ̃ k , ω 0 , x K ( k ) d k = 0 Re { B 0 [ i α 0 f ( x , k ) e ϕ + β 0 g ( x , k ) e ρ + h ( x , k ) e z ] } K ( k ) d k ,

where

f(x,k):=0πsinθcos(kzcosθ)J1(kρsinθ)dθg(x,k):=-0πsinθcosθsin(kzcosθ)J1(kρsinθ)dθ

and

h(x,k):=-0πsin2θcos(kzcosθ)J0(kρsinθ)dθ.

The complex constants α0 and β0 are the free parameters of the magnetic field in Eq. (26) and can be chosen independently of each other. The first part of the field,

(28) B ̃ r k , ω 0 , x := Re B 0 i α 0 f e ϕ or B r x , ω 0 = 0 Re B 0 i α 0 f e ϕ K ( k ) d k ,

that corresponds to the expansion coefficient α0 is called the magnetic ring (see Fig. 1a). The second part,

(29) B ̃ g ( k , ω 0 , x ) := Re B 0 β 0 g e ρ + h e z or B g ( x , ω 0 ) = 0 Re B 0 β 0 g e ρ + h e z K ( k ) d k ,

corresponding to the expansion coefficient β0, is the magnetic globule (see Fig. 1b).

It should be noted that the electromagnetic knot structures do not form an entire set of mathematical basis functions. Regarding the derivation presented here, electromagnetic knots can be written as a superposition of infinitely many plane waves, as plane waves represent an entire set of basis functions. However, the inverse is not true. The functions α(φ,θ) and β(φ,θ) in Eq. (16) control the angular dependency in the k space. By choosing α(φ,θ)=const. and β(φ,θ)=const., infinitely many plane waves propagating in every direction contribute to the field. The resulting field structures are solenoidal and spatially localized. Thus, the magnetic ring and the magnetic globule can be interpreted as a set of basis functions for isotropically localized, divergence-free structures. Choosing different shapes for the functions α(φ,θ) and β(φ,θ) enables the modeling of structures beyond electromagnetic knots.

3.2 Electric field

The electric field and the magnetic field are connected via Ampère's law. Under the absence of ohmic currents, Ampère's law reduces to

(30) x × B ( x , t ) = 1 c ph 2 t E ( x , t ) ,

where cph is the phase velocity. Fourier transformation provides

(31) i k × B ( x , ω ) = - i ω c ph 2 E ( x , ω ) .

Using k=kek yields

(32) e k × B ( x , ω ) = - ω k c ph 2 E ( x , ω )

such that

(33) E ( x , ω ) = - k c ph 2 ω e k × B ( x , ω ) .

Ampère's law is valid for every k vector that contributes to the spectrum of the field, yielding the ansatz

(34) E ̃ ( k , ω , x ) = - k c ph 2 ω e k × B ̃ ( k , ω , x ) .

Using

ek×eφ=-eθandek×eθ=eφ

delivers

(35) E ̃ ( k , ω , x ) = - k c ph 2 ω B 0 2 π 0 2 π 0 π [ β 0 e φ - α 0 e θ ] e i k e k x sin θ d θ d φ

such that the real part can be expressed as

(36) E ̃ ( k , ω , x ) = - k c ph 2 ω Re { B 0 [ i β 0 f ( x , k ) e ϕ - α 0 g ( x , k ) e ρ + h ( x , k ) e z ] } .

Thus, the electric field is given by

(37) E ( x , ω 0 ) = - 0 k c ph 2 ω 0 Re { B 0 [ i β 0 f ( x , k ) e ϕ - α 0 g ( x , k ) e ρ + h ( x , k ) e z ] } K ( k ) d k .

3.3 Electric current density

When ohmic currents j(x,t)0 are present, Ampère's law can be written as

(38) x × B ( x , t ) = μ 0 j ( x , t )

under the assumption of stationarity or if the displacement current is negligible. Again, Fourier transformation provides

(39) i k × B ( x , ω ) = μ 0 j ( x , ω )

such that

(40) j ̃ ( k , ω , x ) = i k μ 0 e k × B ̃ ( k , ω , x ) .

In analogy to the electric field, the current density can be calculated via

(41) j ( x , ω 0 ) = 0 k μ 0 Re { i B 0 [ i β 0 f ( x , k ) e ϕ - α 0 g ( x , k ) e ρ + h ( x , k ) e z ] } K ( k ) d k .

Thus, the current density of the magnetic ring follows the topology of a globule and vice versa.

3.4 Spatially distributed knot structures

Within the derivation of the knot structures, the magnetic ring and the magnetic globule are defined with respect to the same origin of the cylindrical coordinate system (ρ,ϕ,z). The resulting structures are also known as (electro)magnetic disturbances of the first kind (Cameron2018). However, in general the structures can be defined with respect to different (local) coordinate systems, spanned by the local unit vectors (eρq,eϕq,ezq), where q=1,,Q, with different origins 𝓞q. The resulting structures,

(42) B x , ω 0 = 0 Re { B 0 q = 1 Q [ i α 0 q f x q , k e ϕ q + β 0 q ( g x q , k e ρ q + h ( x q , k ) e z q ) ] } K ( k ) d k ,

where

x=Oq+xq,xq=ρqeρq+zqezq,

are a superposition of Q translated and/or rotated (electro)magnetic disturbances of the first kind (see Fig. 1c) and are also called (electro)magnetic disturbances of the second kind (Cameron2018). The field is characterized by 8Q free parameters, i.e., the expansion coefficients α0q and β0q, the origins 𝓞q, and the orientation of the local coordinate system that can be described, for example, via Euler angles (Cameron2018).

3.5 Discussion of the knot structures

Within the derivation presented above, the spectral distribution of the field with respect to k is controlled by the function 𝒦(k). Electromagnetic knots, as originally described by Cameron (2018), are superpositions of infinitely many monochromatic plane waves, i.e., K(k)=δ(k-k0), with the same amplitude, propagating in every direction with the spectrum

(43) S k 0 = k R 3 | k x 2 + k y 2 + k z 2 = k 0 2 .

In contrast to single plane waves, knots are localized structures, similar to wave packages. The localization of the structures results from the spatial distribution of the wave phases:

(44) F ( θ , φ ) := e k x = sin θ x cos φ + y sin φ + z cos θ .

Thus, the knots are a superposition of plane waves with different phases ℱ(θ,φ) at all points in space despite its central point. At the origin of the structure (x=y=z=0) the phases of the waves are all equal: F(θ,φ)=0, resulting in a constructive interference with a maximum amplitude at the central point. The scale size of the knot is determined by k0, representing a set of infinitely many k vectors with the same length. The superposition of the plane waves is schematically illustrated in Fig. 2.

https://angeo.copernicus.org/articles/41/253/2023/angeo-41-253-2023-f02

Figure 2Illustration of superposed, monochromatic plane wave fronts (gray lines) with the wave length λ0=2π/k0. The knots are localized in the origin of the red coordinate system spanned by the vectors eρ, eϕ and ez.

Download

Equation (27) represents the magnetic field with respect to the position vector x an the frequency ω0. However, the spatial structure of the field can also directly be analyzed from the measurement data B(x,t) evaluated at different time steps t, and thus no Fourier transform with respect to time is required.

4 Extension of the wave telescope

Following this short derivation and discussion of the electromagnetic knots, the knot model needs to be reformulated in terms of the wave telescope technique to estimate the spectrum of the knots.

4.1 Reformulation of the model

After performing the temporal Fourier transform, the magnetic field (Eq. 27), measured at the position xi, i=1,,N, can be rewritten as

(45) B x i , ω 0 = 0 B ̃ ( k , ω 0 , x i ) K ( k ) d k = 0 Re H ( k ) B ^ 0 ( k , ω 0 ) i α 0 B ^ 0 ( k , ω 0 ) β 0 K ( k ) d k ,

where

(46) H i ( k ) = - f ( x i , k ) sin ϕ g ( x i , k ) cos ϕ f ( x i , k ) cos ϕ g ( x i , k ) sin ϕ 0 h ( x i , k )

is the corresponding shape matrix of the position xi. Summarizing the measurements into a 3N-dimensional vector B(ω0), the magnetic field can be rearranged as

(47) B ( ω 0 ) = 0 Re H ( k ) B ^ 0 ( k , ω 0 ) i α 0 B ^ 0 ( k , ω 0 ) β 0 K ( k ) d k ,

where

(48) H ( k ) := H 1 ( k ) H N ( k ) R 3 N × 2 .

Again, the determination of the amplitudes α0B^0(k,ω0) and β0B^0(k,ω0) results in an overdetermined inversion problem. In analogy to the classical wave telescope technique, the spectrum of the ring can be estimated via

(49) P ( k ) = tr H M - 1 H - 1 .

Since P(k) is a nonlinear function of k, the whole k space has to be scanned to estimate the spectrum of the field (Motschmann et al.1996).

Solely considering the magnetic ring (Eq. 28), the shape matrix transfers onto the shape vector (Narita2019)

(50) h r ( k ) := f ( x , k ) e ϕ = - f ( x , k ) sin ϕ f ( x , k ) cos ϕ 0 .

In this case, the spectrum of the ring can be estimated via

(51) P r ( k ) = 1 h r ( k ) M - 1 h r ( k ) .
https://angeo.copernicus.org/articles/41/253/2023/angeo-41-253-2023-f03

Figure 3Reconstructed spectrum Pr(k) resulting from different measurement positions (red dots) with respect to the origin of the plasmoid. The length scale of the plasmoid is chosen to be k04.19RE-1.

Download

4.2 Application to plasmoids

For the first application of electromagnetic knots in the context of magnetospheric structures, we consider the modeling of plasmoids via a magnetic ring (Zhang et al.2013). Plasmoids are a consequence of magnetic reconnection in the far-tail region of a planetary magnetosphere triggered by the Dungey cycle (McPherron1995, e.g.,). The structures are characterized by a magnetic ring along the neutral sheet line with a length scale of the order of the solar wind's obstacle (e.g., McPherron1995; Zong et al.2004).

We model the magnetic field in the tail region by superposing a stationary magnetic ring (α0=-i, Eq. 28),

(52) B r ( x , ω 0 = 0 ) = 0 Re B 0 i α 0 f ( x , k ) e ϕ δ ( k - k 0 ) d k = Re B 0 i α 0 f ( x , k 0 ) e ϕ = B 0 h r ( k 0 ) ,

composed of monochromatic plane waves, representing the plasmoid, with the field generated by the neutral sheet current (Harris neutral sheet, Harris1962) such that

(53) B ( x ) = B r x , ω 0 = 0 - B s tanh y L e x = B 0 h r x , k 0 - B s tanh y L e x ,

where the x axis points towards the night side magnetosphere, the y axis points from the southern geographic pole to the northern geographic pole and the z axis completes the right-handed system. Thus, we model the plasmoid as a two-dimensional structure in the xy plane (Zhang et al.2013). The value B0 represents an arbitrarily chosen background amplitude, Bs=0.3 B0, and the length scale of the current sheet is chosen to be L=10-3RE, where RE is the planetary radius, e.g., the terrestrial radius. The characteristic length scale of the plasmoid is chosen to be λ0=1.5 RE, corresponding to k0=2π/λ04.19RE-1.

The resulting magnetic field data are evaluated at N=7 synthetically generated spacecraft positions, representing a HelioSwarm-like configuration (Klein and Spence2021). As plasmoids are highly dynamical, traveling structures, the measurement positions are shifted along the x axis with respect to the origin of the plasmoid (left, mean, right), representing different time steps. The length scale λ0 (or equivalently k0) of the plasmoid is estimated from the virtual spacecraft data via Eq. (51). The resulting field geometry (blue arrows) and the measurement positions (red dots) as well as the corresponding spectra are illustrated in Fig. 3.

When the measurement positions are distributed around the origin of the plasmoid (mean), the implemented value of k0 can be reconstructed with high precision from the data. In the other cases, the spatial length scale is slightly overestimated and the relative error results in about 6 % (left) and 4 % (right). Thus, the wave telescope technique is capable of (1) separating the plasmoid from the neutral sheet part and (2) estimating the characteristic length scale of the plasmoid from a limited number of measurement positions.

In analogy to the classical wave telescope technique, the accuracy of the reconstruction depends of the relation between the plasmoid's length scale λ0 and the mean distance d between the spacecraft positions (Narita et al.2022, e.g.,). For example, if dλ0, the measurement positions do not properly cover the spatial extend of the plasmoid, resulting in ambiguities within the reconstruction procedure. In the case of dλ0, the magnetic field structure of the plasmoid is not detectable. Thus, the mean distance between the spacecraft positions has to be of the order of the plasmoid's spatial scale dλ0, which will be realized by the configuration of the planned HelioSwarm multiscale mission.

Furthermore, the amplitude of the ring B0 has to be of the same order as or larger than the sheath field Bs to guarantee a precise reconstruction result. For example, in the case of Bs=10 B0 no peak occurs within the spectrum Pr(k) and the ring cannot be discerned from the background field. On the other hand, the peak within the spectrum becomes sharper in the case of Bs=0.1 B0.

4.3 Further applications

The application presented above of electromagnetic knots indicates the potential of the representation. Spatially distributed electromagnetic knots as described by Cameron (2018) enable the modeling of more complex structures, provide generalized spectral information and open the door for further applications, delivering an alternative interpretation of magnetospheric structures. For example, the magnetic field configuration resulting from a field-aligned current can be modeled as a superposition of magnetic rings stacked on top of each other. Due to Ampère's law, the corresponding current density is given as a superposition of globules. Thus, the inner structure of field-aligned currents can be analyzed directly from the magnetic field measurements (Toepfer et al.2021). Also, the current system of Alfvén wings can be described as a superposition of rings (Vernisse et al.2018, e.g.,) so that the corresponding magnetic field topology follows the structure of superposed globules. Furthermore, field line resonances (Glassmeier et al.1999; Plaschke et al.2008) may be described as a special superposition of magnetic rings.

5 Conclusions

Electromagnetic knots are a superposition of infinitely many monochromatic plane waves with a spherical symmetric spectrum and represent an exact solution of Maxwell's equation. The resulting basis elements, i.e., the electromagnetic ring and the globule, form a basis set for localized, divergence-free spatial structures. For this reason, the concept of electromagnetic knots opens the door for a completely new description and interpretation of spatial structures in planetary magnetospheres.

The classification of spatial structures evaluated at a limited number of measurement points describes an overdetermined inversion problem. The wave telescope technique serves as a robust data analysis tool for the global interpretation of spacecraft measurements in terms of expected physical structures. By reformulating the formalism of electromagnetic knots in terms of the wave telescope technique, we extended the zoo of spatial structures that can be analyzed by the method. In this sense, the present study can be interpreted as a generalization of the wave telescope technique to a structure telescope technique.

For a first validation, the concept of electromagnetic knots has been applied to the modeling of a plasmoid. Using a HelioSwarm-like satellite configuration, the wave telescope technique is capable of separating the plasmoid, modeled as a magnetic ring, from the field generated by the neutral sheet current and enables the estimation of the length scale of the ring. Thus, the presented extension of the wave telescope technique serves as a new data analysis tool for multi-spacecraft missions, such as the planned HelioSwarm mission. However, the application of electromagnetic knots for characterizing further structures, such as field-aligned currents or Alfvén wings, should be analyzed in future studies. In general, we conclude that the modified wave telescope technique outlined here bears the potential for a new representation and physical description of complex spatial structures existing in space plasmas.

Appendix A: Evaluation of the integrals

The x component of the magnetic field in Eq. (22) can be rewritten as

(A1) 0 2 π 0 π - α 0 sin φ + β 0 cos θ cos φ e i k e k x sin θ d θ d φ = 0 π sin θ e i k z cos θ 0 2 π - α 0 sin φ + β 0 cos θ cos φ e i ( k x cos φ + k y sin φ ) sin θ d φ d θ = - α 0 0 π sin θ e i k z cos θ 0 2 π sin φ e i η 1 cos φ e i η 2 sin φ d φ d θ + β 0 0 π sin θ cos θ e i k z cos θ 0 2 π cos φ e i η 1 cos φ e i η 2 sin φ d φ d θ = - α 0 0 π sin θ e i k z cos θ I 1 ( θ ) d θ + β 0 0 π sin θ cos θ e i k z cos θ I 2 ( θ ) d θ ,

where

(A2) I 1 ( θ ) := 0 2 π sin φ e i η 1 cos φ e i η 2 sin φ d φ

and

(A3) I 2 ( θ ) := 0 2 π cos φ e i η 1 cos φ e i η 2 sin φ d φ .

Analogously, the y component in Eq. (22) results in

(A4) 0 2 π 0 π α 0 cos φ + β 0 cos θ sin φ e i k e k x sin θ d θ d φ = 0 π sin θ e i k z cos θ 0 2 π α 0 cos φ + β 0 cos θ sin φ e i η 1 cos φ e i η 2 sin φ d φ d θ = α 0 0 π sin θ e i k z cos θ 0 2 π cos φ e i η 1 cos φ e i η 2 sin φ d φ d θ + β 0 0 π sin θ cos θ e i k z cos θ 0 2 π sin φ e i η 1 cos φ e i η 2 sin φ d φ d θ = α 0 0 π sin θ e i k z cos θ I 2 ( θ ) d θ + β 0 0 π sin θ cos θ e i k z cos θ I 1 ( θ ) d θ

and

(A5) - β 0 0 2 π 0 π sin 2 θ e i k e k x d θ d φ = - β 0 0 π sin 2 θ e i k z cos θ 0 2 π e i η 1 cos φ e i η 2 sin φ d φ d θ = - β 0 0 π sin 2 θ e i k z cos θ I 3 ( θ ) d θ

for the z component in Eq. (22), where

(A6) I 3 ( θ ) := 0 2 π e i η 1 cos φ e i η 2 sin φ d φ .

Here, the magnetic field is given by

(A7) B ̃ ( k , ω 0 , x ) = B 0 2 π { [ - α 0 0 π sin θ e i k z cos θ I 1 ( θ ) d θ + β 0 0 π sin θ cos θ e i k z cos θ I 2 ( θ ) d θ ] e x + [ α 0 0 π sin θ e i k z cos θ I 2 ( θ ) d θ + β 0 0 π sin θ cos θ e i k z cos θ I 1 ( θ ) d θ ] e y - β 0 0 π sin 2 θ e i k z cos θ I 3 ( θ ) d θ e z } = B 0 2 π { α 0 [ - 0 π sin θ e i k z cos θ I 1 ( θ ) d θ e x + 0 π sin θ e i k z cos θ I 2 ( θ ) d θ e y ] + β 0 [ 0 π sin θ cos θ e i k z cos θ I 2 ( θ ) d θ e x + 0 π sin θ cos θ e i k z cos θ I 1 ( θ ) d θ e y ] - β 0 0 π sin 2 θ e i k z cos θ I 3 ( θ ) d θ e z } .

At least, the remaining integrals I1(θ), I2(θ) and I3(θ) have to be evaluated.

For the evaluation of the integral

(A8) I 3 ( θ ) = 0 2 π e i η 1 cos φ e i η 2 sin φ d φ = 0 2 π e i η 1 cos φ + η 2 sin φ d φ

we define

tanγ0:=η2η1

such that

sinγ0=η2η12+η22

and

cosγ0=η1η12+η22.

Thus, the argument of the complex exponential can be rewritten as

(A9) η 1 cos φ + η 2 sin φ = η 1 2 + η 2 2 ( η 1 η 1 2 + η 2 2 cos φ + η 2 η 1 2 + η 2 2 sin φ ) = η 1 2 + η 2 2 ( cos γ 0 cos φ + sin γ 0 sin φ ) = η 1 2 + η 2 2 sin ( φ + γ 0 ) .

Substituting τ:=φ+γ0+π and using sin(φ+γ0)=sin(τ-π)=-sinτ delivers

(A10) I 3 ( θ ) = - π + γ 0 π + γ 0 e - i η 1 2 + η 2 2 sin τ d τ .

As the integrand is a 2π-periodic function, the integral is independent of γ0 so that

(A11) I 3 ( θ ) = - π + γ 0 π + γ 0 e - i η 1 2 + η 2 2 sin τ d τ = - π π e - i η 1 2 + η 2 2 sin τ d τ .

Making use of the definition of the Bessel functions of the first kind,

Jn(x)=12π-ππei(nτ-xsinτ)dτ,

yields

(A12) I 3 ( θ ) = 2 π J 0 ( η 1 2 + η 2 2 ) = 2 π J 0 ( k ρ sin θ ) .

The integral

(A13) I 1 ( θ ) = 0 2 π sin φ e i η 1 cos φ e i η 2 sin φ d φ

can be evaluated using the identity

(A14) η 2 e i η 2 sin φ = i sin φ e i η 2 sin φ

so that

sinφeiη2sinφ=-iη2eiη2sinφ

and results in

(A15) I 1 ( θ ) = - i η 2 0 2 π e i η 1 cos φ e i η 2 sin φ d φ = - 2 π i η 2 J 0 η 1 2 + η 2 2 = - 2 π i η 2 η 1 2 + η 2 2 J 0 η 1 2 + η 2 2 = - 2 π i sin ϕ J 0 η 1 2 + η 2 2

by means of Eq. (A12). Using

(A16) J 0 ( x ) = - J 1 ( x ) ,

delivers

(A17) I 1 ( θ ) = 2 π i sin ϕ J 1 ( k ρ sin θ ) .

Analogously, the integral

(A18) I 2 ( θ ) = 0 2 π cos φ e i η 1 cos φ e i η 2 sin φ d φ

can be evaluated using

(A19) η 1 e i η 1 cos φ = i cos φ e i η 1 sin φ

such that

(A20) cos φ e i η 1 cos φ = - i η 1 e i η 1 cos φ

and results in

(A21) I 2 ( θ ) = - i η 1 0 2 π e i η 1 cos φ e i η 2 sin φ d φ = - 2 π i η 1 J 0 η 1 2 + η 2 2 = - 2 π i η 1 η 1 2 + η 2 2 J 0 η 1 2 + η 2 2 = 2 π i η 1 η 1 2 + η 2 2 J 1 η 1 2 + η 2 2 = 2 π i cos ϕ J 1 ( k ρ sin θ ) .

Therefore, the magnetic field is given by

(A22) B ̃ ( k , ω 0 , x ) = B 0 { i α 0 [ - sin ϕ 0 π sin θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e x + cos ϕ 0 π sin θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e y ] + β 0 [ i cos ϕ 0 π sin θ cos θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e x + i sin ϕ 0 π sin θ cos θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e y ] - β 0 0 π sin 2 θ e i k z cos θ J 0 ( k ρ sin θ ) d θ e z } = B 0 { i α 0 0 π sin θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e ϕ + β 0 i 0 π sin θ cos θ e i k z cos θ J 1 ( k ρ sin θ ) d θ e ρ - β 0 0 π sin 2 θ e i k z cos θ J 0 ( k ρ sin θ ) d θ e z } .

The remaining integrals can be expanded into the form

(A23)Φ:=0πsinθeikzcosθJ1(kρsinθ)dθ=0πsinθcos(kzcosθ)J1(kρsinθ)dθ+i0πsinθsin(kzcosθ)J1(kρsinθ)dθ(A24)Γ:=i0πsinθcosθeikzcosθJ1(kρsinθ)dθ=-0πsinθcosθsin(kzcosθ)J1(kρsinθ)dθ+i0πsinθcosθcos(kzcosθ)J1(kρsinθ)dθ

and

(A25) Ξ : = 0 π sin 2 θ e i k z cos θ J 0 ( k ρ sin θ ) d θ = 0 π sin 2 θ cos ( k z cos θ ) J 0 ( k ρ sin θ ) d θ + i 0 π sin 2 θ sin ( k z cos θ ) J 0 ( k ρ sin θ ) d θ .

Introducing the transformation

θ̃=θ-π2orequivalentlyθ=θ̃+π2

such that

sinθ=sinθ̃+π2=cosθ̃

and

cosθ=cosθ̃+π2=-sinθ̃

shows that the integrands of the imaginary parts are symmetric functions with respect to the value θ̃=0 in the interval θ̃-π2,π2 so that

(A26)ImΦ=-π/2π/2cosθ̃sin(-kzsinθ̃)J1(kρcosθ̃)dθ̃=0,(A27)ImΓ=-π/2π/2cosθ̃(-sinθ̃)cos(kz(-sinθ̃))J1(kρcosθ̃)dθ̃=0,(A28)ImΞ=-π/2π/2cos2θ̃sin(kz(-sinθ̃))J0(kρcosθ̃)dθ̃=0,

whereas the real parts do not vanish in general, as illustrated in Fig. A1. Thus, introducing the abbreviations

(A29)f(x,k):=Re0πsinθeikzcosθJ1(kρsinθ)dθ=0πsinθcos(kzcosθ)J1(kρsinθ)dθ,(A30)g(x,k):=Rei0πsinθcosθeikzcosθJ1(kρsinθ)dθ=-0πsinθcosθsin(kzcosθ)J1(kρsinθ)dθ,

and

(A31) h ( x , k ) : = - Re 0 π sin 2 θ e i k z cos θ J 0 ( k ρ sin θ ) d θ = - 0 π sin 2 θ cos ( k z cos θ ) J 0 ( k ρ sin θ ) d θ ,

the measurable part (i.e., the real part) of the magnetic field can finally be written as

(A32) B ̃ ( k , ω 0 , x ) = Re B 0 i α 0 f e ϕ + β 0 g e ρ + h e z .
https://angeo.copernicus.org/articles/41/253/2023/angeo-41-253-2023-f04

Figure A1Real (blue) and imaginary part (orange) of the integrand of Φ in the interval θ0,π.

Download

Data availability

Data can be provided upon request.

Author contributions

ST performed the mathematical derivations. KHG came up with the idea of extending the wave telescope technique. ST and UM wrote the first draft of the manuscript. All authors contributed equally to the conception and design of the study; all authors read and approved the submitted version.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

We are grateful to Norbert Fürstenau (DLR Braunschweig) for pointing our interest to electromagnetic knots. We acknowledge support by the German Research Foundation and the Open Access Publication Fund of the Technische Universität Braunschweig.

Financial support

This work has been supported by the German Research Foundation and the Open Access Publication Fund of the Technische Universität Braunschweig. The contribution by Karl-Heinz Glassmeier was financially supported by the German Bundesministerium für Wirtschaft und Klimaschutz and the Deutsches Zentrum für Luft- und Raumfahrt under 50OC1803.

This open-access publication was funded by Technische Universität Braunschweig.

Review statement

This paper was edited by Elena Kronberg and reviewed by two anonymous referees.

References

Cameron, R. P.: Monochromatic knots and other unusual electromagnetic disturbances: light localised in 3D, J. Phys. Commun., 2, 015024, https://doi.org/10.1088/2399-6528/aa9761, 2018. a, b, c, d, e, f, g, h, i, j, k

Capon, J.: High resolution frequency-wavenumber spectrum analysis, Proc. IEEE, 57, 1408–1418, https://doi.org/10.1109/PROC.1969.7278, 1969. 

Constantinescu, O. D., Glassmeier, K.-H., Motschmann, U., Treumann, R. A., Fornaçon, K.-H., and Fränz, M.: Plasma wave source location using CLUSTER as a spherical wave telescope, J. Geophys. Res., 111, A09221, https://doi.org/10.1029/2005JA011550, 2006. a

Gauss, C. F.: Allgemeine Theorie des Erdmagnetismus: Resultate aus den Beobachtungen des magnetischen Vereins im Jahre 1838, edited by: Gauss, C. F. and Weber, W., 1–57, Weidmannsche Buchhandlung, Leipzig, 1839. a

Glassmeier, K.-H. and Tsurutani, B. T.: Carl Friedrich Gauss – General Theory of Terrestrial Magnetism – a revised translation of the German text, Hist. Geo Space. Sci., 5, 11–62, https://doi.org/10.5194/hgss-5-11-2014, 2014. a

Glassmeier, K.-H., Othmer, C., Cramm, R., Stellmacher, M., and Engebretson, M.: Magnetospheric field line resonances: A comparative planetology approach, Surv. Geophys., 20, 61–109, https://doi.org/10.1023/A:1006659717963, 1999. a

Glassmeier, K.-H., Motschmann, U., Dunlop, M., Balogh, A., Acuña, M. H., Carr, C., Musmann, G., Fornaçon, K.-H., Schweda, K., Vogt, J., Georgescu, E., and Buchert, S.: Cluster as a wave telescope – first results from the fluxgate magnetometer, Ann. Geophys., 19, 1439–1447, https://doi.org/10.5194/angeo-19-1439-2001, 2001. a

Harris, E. G.: On a plasma sheath separating regions of oppositely directed magnetic fields, Nuovo Cim. 23, 115–121, https://doi.org/10.1007/BF02733547, 1962. a

Haykin, S.: Adaptive Filter Theory, 5th Edn., International edition, Pearson, ISBN 10: 0-273-76408-X, 2014. a

Klein, K. and Spence, H. and the HelioSwarm Science Team: HelioSwarm: Leveraging Multi-Point, Multi-Scale Spacecraft Observations to Characterize Turbulence, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6812, https://doi.org/10.5194/egusphere-egu21-6812, 2021. a

McPherron, R. L.: Magnetospheric dynamics, in: Introduction to space physics, Cambridge University Press, Cambridge, 400–458, edited by: Kivelson, M. G. and Russell, C. T, ISBN-10: 0-521-45104-3, 1995. a, b

Motschmann, U., Woodward, T. I., Glassmeier, K.-H., Southwood, D. J., and Pinçon, J.-L.: Wavelength and direction filtering by magnetic measurements at satellite arrays: Generalized minimum variance analysis, J. Geophys. Res., 101, 4961–4966, https://doi.org/10.1029/95JA03471, 1996. a, b, c, d, e, f

Narita, Y.: A note on Capon's minimum variance projection for multi-spacecraft data analysis, Front. Phys., 7, 8, https://doi.org/10.3398/fphy.2019.00008, 2019. a, b, c

Narita, Y., Glassmeier, K.-H., Schäfer, S., Motschmann, U., Sauer, K., Dandouras, I., Fornaçon, K. H., Georgescu, E., and Rème, H.: Dispersion analysis of ULF waves in the foreshock using cluster data and the wave telescope technique, Geophys. Res. Lett., 30, 1710, https://doi.org/10.1029/2003GL017432, 2003. a

Narita, Y., Kleindienst, G., and Glassmeier, K.-H.: Evaluation of magnetic helicity density in the wave number domain using multi-point measurements in space, Ann. Geophys., 27, 3967–3976, https://doi.org/10.5194/angeo-27-3967-2009, 2009. a

Narita, Y., Nakamura, R., and Baumjohann, W.: Cluster as current sheet surveyor in the magnetotail, Ann. Geophys., 31, 1605–1610, https://doi.org/10.5194/angeo-31-1605-2013, 2013. a

Narita, Y., Glassmeier, K.-H., and Motschmann, U.: The wave telescope technique, J. Geophys. Res., 127, e2021JA030165, https://doi.org/10.1029/2021JA030165, 2022. a, b, c

Plaschke, F., Glassmeier, K.-H., Constantinescu, O. D., Mann, I. R., Milling, D. K., Motschmann, U., and Rae, I. J.: Statistical analysis of ground based magnetic field measurements with the field line resonance detector, Ann. Geophys., 26, 3477–3489, https://doi.org/10.5194/angeo-26-3477-2008, 2008. a, b, c

Toepfer, S., Narita, Y., Heyner, D., and Motschmann, U.: The Capon method for Mercury's magnetic field analysis, Front. Phys., 8, 249, https://doi.org/10.3389/fphy.2020.00249, 2020. a, b

Toepfer, S., Narita, Y., Heyner, D., Kolhey, P., and Motschmann, U.: Mathematical foundation of Capon's method for planetary magnetic field analysis, Geosci. Instrum. Method. Data Syst., 9, 471–481, https://doi.org/10.5194/gi-9-471-2020, 2020. a, b, c

Toepfer, S., Narita, Y., Exner, W., Heyner, D., Kolhey, P., Glassmeier, K.-H., and Motschmann, U.: The Mie representation for Mercury’s magnetospheric currents, Earth Planet. Space, 73, 65, https://doi.org/10.1186/s40623-021-01536-8, 2021. a, b

Vernisse, Y., Riousset, J. A., Motschmann, U., and Glassmeier, K.-H.: Simulations of stellar winds and planetary bodies: Magnetized obstacles in a super-Alfvènic flow with southward IMF, Planet. Space Sci., 152, 18–30, https://doi.org/10.1016/j.pss.2018.01.010, 2018. a

Zhang, Y. C., Shen, C., Liu, Z. X., Rong, Z. J., Zhang, T. L., Marchaudon, A., Zhang, H., Duan, S. P., Ma, Y. H., Dunlop, M. W., Yang, Y. Y., Carr, C. M., and Dandouras, I.: Two different types of plasmoids in the plasma sheet: Cluster multisatellite analysis application, J. Geophys. Res., 118, 5437–5444, https://doi.org/10.1002/jgra.50542, 2013.  a, b

Zong, Q.-G., Fritz, T. A., Pu, Z. Y., Fu, S. Y., Baker, D. N., Vogiatzis, I., Glassmeier, K.-H., Korth, A., Daly, P. W., Balogh, A., and Reme, H.: Cluster observations of earthward flowing in the tail, Geophys. Res. Let., 31, L18803, https://doi.org/10.1029/2004GL020692, 2004. a

Download
Short summary
The present study discusses the modeling and interpretation of magnetospheric structures via electromagnetic knots for the first time. The mathematical foundations of electromagnetic knots are presented, and the formalism is reformulated in terms of the classical wave telescope technique. The method is tested against synthetically generated magnetic field data describing a plasmoid as a two-dimensional magnetic ring structure.