The Mie representation for Mercury’s magnetospheric currents

Poloidal–toroidal magnetic field decomposition is a useful application of the Mie representation and the decomposition method enables us to determine the current density observationally and unambiguously in the local region of magnetic field measurement. The application and the limits of the decomposition method are tested against the Mercury magnetic field simulation in view of BepiColombo’s arrival at Mercury in 2025. The simulated magnetic field data are evaluated along the planned Mercury Planetary Orbiter (MPO) trajectories and the current system that is crossed by the spacecraft is extracted from the magnetic field measurements. Afterwards, the resulting currents are classified in terms of the established current system in the vicinity of Mercury.


Introduction
The fluxgate magnetometer Heyner et al. 2021) on board the Mercury Planetary Orbiter (MPO) of the BepiColombo mission (Benkhoff et al. 2010) will provide precise measurements of the magnetic field in the vicinity of Mercury which enable the detailed characterization of Mercury's internal magnetic field. Due to the plasma physical interaction of Mercury with the solar wind, the measurements are composed of the desired internal contributions as well as external fields resulting from currents flowing within the magnetosphere. For the reconstruction of the internal field, each contribution has to be parametrized properly. Besides the reconstruction of the internal magnetic field, the characterization of the magnetic field contributions within the Open Access *Correspondence: s.toepfer@tu-braunschweig.de 1 Institut für Theoretische Physik, Technische Universität Braunschweig, Mendelssohnstraße 3, 38106 Braunschweig, Germany Full list of author information is available at the end of the article magnetosphere allows the analysis of the external field and the corresponding currents that are crossed by the spacecraft. Since the conventionally used Gauss representation (Gauss 1839;Glassmeier and Tsurutani 2014) is based on the irrotational structure of the magnetic field, the Gauss potential theory does not allow to accurately reconstruct the planetary magnetic field if electric currents are flowing outside the planet and actively generate rotational field structures. Typically, magnetospheric models are used to estimate the magnetic field in current-carrying regions (e.g., Korth et al. 2004Korth et al. , 2015Korth et al. , 2017Alexeev et al. 2008). These models approximate the magnetic field using a system of modules for the internal field and the fields generated by the tail and the magnetopause currents. As an alternative for the application of global magnetospheric models a combination of the Gauss representation with the Mie representation (toroidal-poloidal decomposition) (Backus 1986(Backus , 1996Olsen 1997), called the Gauss-Mie representation (Toepfer et al. 2021a), enables the parametrization of the magnetic field in current-carrying regions. This parametrization generalizes the international geomagnetic reference field (IGRF) model for the Earth's core field (Thébault et al. 2015) and has successfully been applied in order to reconstruct Mercury's internal magnetic field from simulated magnetic field data (Toepfer et al. 2021a) as well as for the description of the Earth's magnetosphere (Andreeva and Tsyganenko 2016). Especially, the Gauss-Mie representation delivers an analytical parametrization for the magnetic field and therefore enables the analytical calculation of the (local) current density that is crossed by the spacecraft. Within the Earth's magnetosphere the Gauss-Mie representation has successfully been applied to reconstruct ionospheric F-layer currents from MAG-SAT and Champ data (Olsen 1997;Engels and Olsen 1999;Bayer 2001;Mayer and Maier 2006). Tsyganenko et al. (2021) reconstructed magnetospheric storm-time dynamics within the terrestrial magnetosphere using the Mie representation. Furthermore, the magnetic field can be parametrized via Euler potentials (e.g., Stern 1967Stern , 1970Cheng and Zaharia 2003;Romashets and Vandas 2020).
The goal of our study is the determination of the poloidal and toroidal current density flowing in Mercury's magnetosphere in order to validate the parametrization of the magnetic field via the Gauss-Mie representation for the reconstruction of Mercury's internal magnetic field (Toepfer et al. 2021a). In preparation for the observations soon to be obtained by the BepiColombo mission, the plasma interaction of Mercury with the solar wind is simulated numerically using the 3D hybrid model AIKEF (adaptive ion kinetic electron fluid) (Müller et al. 2011). First of all, the established current system at Mercury in comparison with the Earth's magnetosphere is summarized. Then the mathematical foundations of the Gauss-Mie representation are revisited and the parametrization is applied in order to reconstruct poloidal and toroidal currents in Mercury's magnetosphere flowing in the vicinity of the trajectories of the Mercury Planetary Orbiter (MPO) from the simulated magnetic field data. After the reconstruction of the currents, the question arises how these currents are related to a current system around Mercury and how this current system is closed. In analogy to the Earth the structure of Mercury's magnetosphere depends on the orientation of the interplanetary magnetic field (IMF) (e.g., Grygorov et al. (2017)). Thus, the resulting poloidal and toroidal current systems are presented for six different IMF-orientations which may be regarded as mathematical basis orientations.

Current system at Mercury
Before reconstructing the currents in the vicinity of the planned MPO trajectories, we shortly summarize the established current system at Mercury. As the solar wind is a highly conducting medium, Mercury's internal magnetic field is unable to penetrate the solar wind plasma. Therefore, the internal field is compressed at the dayside and stretched into a tail along the nightside. This deformation of the irrotational internal magnetic field is associated with currents flowing in the magnetosphere. In comparison to the Earth's magnetic field, Mercury's internal magnetic field is much weaker so that Mercury's magnetosphere is much smaller than the Earth's magnetosphere. Since the internal dipole moments of both the planets are orientated nearly antiparallel to the planetary rotation axis, the shape of the magnetospheres is comparable. Therefore, the Chapman-Ferraro currents j cf and especially the equatorial magnetopause current j mp isolating the internal magnetic field from the solar wind plasma as well as the neutral sheet current j ns in the tail are qualitatively the same (e.g., Glassmeier (2000)).
Within the Earth's magnetosphere the existence of a dense ionosphere is of major importance for the current closing of field-aligned currents (Iijima and Potemra 1976;Glassmeier 2000). In contrast to the Earth, Mercury does not possess an ionosphere consisting of a dense plasma with a high electrical conductivity, where these currents can be closed. Since the physical and chemical constitution of Mercury's surface is not fully resolved yet, the question whether the field-aligned currents at Mercury are able to penetrate the surface is under debate for decades (e.g., Glassmeier 2000;Janhunen 2004). The analysis of spectrometer data from the MESSENGER mission suggests the existence of an exosphere which is mainly composed of sodium ions Na + Raines et al. 2014Raines et al. , 2015James et al. 2019). However, the absolute value of the sodium ion density, ranging from 5.1 × 10 −3 cm −3 up to 22 cm −3 is still under discussion. Hybrid simulations of the plasma interaction of Mercury with the solar wind considering the existence of a dense sodium exosphere show that Region 1 currents j R1 are able to close within the exosphere (Exner et al. 2020). On the other hand, the analysis of MESSENGER magnetic field data indicates the closing of field-aligned currents through Mercury's surface .
Since the significance of Mercury's exosphere is still under debate, it is not clear whether Region 2 currents as well as fully developed ring currents occur within Mercury's magnetosphere (e.g., Glassmeier 2000;Baumjohann et al. 2010;Exner et al. 2020). In Table 1, the established current system at Mercury in comparison with the Earth's magnetosphere is summarized.
Although the magnetospheres of Mercury and the Earth are qualitatively comparable, it should be noted that the physical processes for the current closure within the magnetospheres differ from each other (e.g., Glassmeier (2000)). Due to its high electrical conductivity the ability for the closing of magnetospheric currents even without the existence of an ionosphere may be attributed to the plasma. Therefore, it is debatable, whether the naming of the currents that has been established within the description of the Earth's magnetosphere (e.g., Region 1 current) is directly transferable onto the naming of the currents flowing in Mercury's magnetosphere.

Characterization of the magnetic field and the current density
In the following, the basic ideas of the Gauss-Mie representation for modeling the magnetic field and the current density, based on the works of Backus (1986Backus ( , 1996, Olsen (1997) and Toepfer et al. (2021a) are shortly revisited.

Characterization of the magnetic field
The magnetic field B will be measured around Mercury along elliptic orbits like that of the MPO spacecraft. Conceptually covering the orbits by a spherical shell with inner radius a, outer radius c and mean radius b = (a + c)/2 enables us to distinguish between different magnetic field contributions. Due to the superposition principle, the total measured field can be written as which is called the Gauss-Mie representation of the magnetic field (Toepfer et al. 2021a). The internal field B i results from currents flowing beneath the shell in the region r < a . Thus, B i is purely poloidal and especially irrotational within the shell. The currents flowing above the shell ( r > c ) generate an external field structure B e of the same character. Therefore, these contributions can be parametrized via the Gauss representation ( (1) The currents flowing within the shell in the region a < r < c generate (completely rotational) toroidal B sh T and (potentially rotational) poloidal B sh P magnetic field structures that can be parametrized via the Mie representation (Backus 1986(Backus , 1996Olsen 1997) by making use of scalar functions sh T and sh P , so that and is valid, where r = r e r and e r is the unit vector in radial direction. The functions sh T and sh P can be interpreted as a special case of the Euler potentials (Stern 1967(Stern , 1970, as shown within the Appendix. Because of the underlying spherical geometry it is straightforward to expand the scalar functions sh T and sh P into spherical harmonics in analogy to the scalar potentials. Since the exact radial dependence of the corresponding expansion coefficients is unknown, it is useful to perform a Taylor series expansion for the coefficients with respect to the radius r in the vicinity of the mean radius b of the spherical shell (Toepfer et al. 2021a), providing us and where Theoretically, the functions sh P and sh T as well as the potentials i and e have to be expanded into a series containing an infinite set of basis functions. Within the practical application, the summations have to be truncated at a suitable degree l < ∞ , where it is desirable to incorporate as less terms as possible. On the other hand, the magnetic field has to be parametrized adequately.
Thus, it is useful to choose a manageable ansatz for the series expansion which can be extended gradually.
Inserting the parametrization of the magnetic field contributions into Eq. (1) delivers Rearranging the terms of the series expansions into the shape matrices H i , H e , H sh T and H sh P and summarizing the corresponding expansion coefficients into the vectors g i , g e , g sh T and g sh P , the Gauss-Mie representation of the magnetic field can be rewritten in the linear algebraic form The H s mainly contain well-known information about the positions, whereas the g s are the unknown coefficients which need to be determined from the measurements.
For the reconstruction of Mercury's internal magnetic field, the application of the thin shell approximation (Backus 1986(Backus , 1996Olsen 1997;Toepfer et al. 2021a) is a valid assumption. The thin shell approximation allows the negligence of the poloidal field B sh P generated by toroidal currents flowing within the shell, when the width of the shell is smaller than the length scale on which the currents change in radial direction. For the analysis of the full current density inside the shell this negligence, in turn, causes a problem, as for this purpose the knowledge of B sh P is indispensable. The way out works as follows. After the application of the thin shell approximation, the internal Gauss coefficients g m l and h m l , the external Gauss coefficients q m l and s m l as well as the expansion coefficients a m l , b m l , a ′ m l , b ′ m l for the toroidal magnetic field can be estimated by making use of a suitable inversion method such as Capon's method (Capon 1969;Motschmann et al. 1996;Toepfer et al. 2020bToepfer et al. , 2021b, that serves as a powerful tool for the analysis of planetary magnetic fields (Toepfer et al. 2020a). Afterwards, the expansion coefficients c m l , d m l , c ′ m l , d ′ m l for the poloidal magnetic field can be estimated by obtaining an approximate solution of for the coefficient vector g sh P . It should be noted that the expansion coefficients for the poloidal magnetic field cannot be estimated simultaneously with the internal (10) (11) B = H i g i + H e g e + H sh T g sh T + H sh P g sh P .
(12) B sh P = B − B i + B e + B sh T = H sh P g sh P and external Gauss coefficients. As discussed in Toepfer et al. (2021a), these fields follow the same topological structure, so that the poloidal field B sh P cannot be distinguished from the internal B i and external field B e within the reconstruction procedure.

Characterization of the current density
After the reconstruction of the poloidal g sh P and toroidal coefficients g sh T from the magnetic field data, the corresponding current density j sh flowing within the shell can be determined using Ampère's law so that j sh describes a solenoidal vector field. Inserting the parametrization of the magnetic field (cf. Equation 1) and using the current density within the shell is given by Considering the divergence-free nature of the current density (no charge accumulation), j sh can be decomposed into its toroidal part and poloidal part in analogy to the magnetic field. Thus, and simultaneously hold so that the poloidal and the toroidal current systems are closed independently of each other. (13)

Application to simulated Mercury magnetic field data
For a first validation of the Gauss-Mie representation suggested here to reconstruct the current density j sh in the vicinity of the planned MPO orbits, simulated stationary magnetic field data and electric currents system information based on using the AIKEF hybrid code (electron fluid, kinetic ions) (Müller et al. 2011) are analyzed. This code has successfully been applied to several problems in Mercury's plasma interaction (e.g., Exner et al. 2018Exner et al. , 2020. From the simulated magnetic field data, the current system is derived using the above described new approach and the results are compared with the electric currents directly derived from the simulation. Mercury's internal magnetic field is modeled as a dipole field with an internal dipole moment of −190 nT , which is shifted northward by 0.2 R M . This field can equivalently be described as a multipole field with the internal Gauss coefficients g 0 1 = −190 nT for the dipole field, g 0 2 = −78 nT for the quadrupole field and g 0 3 = −20 nT for the octupole field (Anderson et al. 2012;Wardinski et al. 2019). The internal Gauss coefficients determining the internal stationary magnetic field are implemented in the simulation code and the plasma interaction of Mercury with the solar wind is simulated. The interplanetary magnetic field with a magnitude of B IMF = 20 nT (Winslow et al. 2013) is orientated along the vector x, y, z T = (0, 0, −1) T in the Mercury-Anti-Solar-Orbital coordinate system (MASO), i.e., the x-axis is orientated towards the nightside of Mercury (away from the sun), the z-axis is orientated parallel to the rotation axis (i.e., antiparallel to the internal dipole moment) and the y-axis completes the right hand system. The solar wind velocity of u sw = 400 km/s points along the x-axis and the solar wind proton density was chosen to n sw = 30 cm −3 (cf. Winslow et al. (2013)). Mercury's outer mantle ( 0.7 R M ≤ r ≤ 1 R M ) is modeled with a radially symmetric planetary resistivity profile with a resistivity of η S ≈ 840 k�m (or conductivity of 10 −6 S/m ) at the surface and the core possesses a vanishing resistivity Exner et al. 2018Exner et al. , 2020. The influence of Mercury's exosphere is neglected, since the plasma interaction is only affected by the exosphere in the case of an extremely high ion density and the absolute value of the density is unknown (Exner et al. 2020). The resulting simulated magnetic field data are evaluated along the ellipsoid which describes the envelope of the elliptical MPO orbit rotated in longitudinal direction from −50 • (afternoon/ post-midnight sector) over 0 • (noon/midnight, x-zplane) to 50 • (morning/pre-midnight sector) around the rotation axis (z-axis).
For the reconstruction of the current density j sh from the magnetic field data B , the scalar potentials i , e of the internal and external fields are expanded into spherical harmonics up to the third degree and order representing the internal/external dipole, quadrupole and octupole field. The scalar functions sh P , sh T of the poloidal and toroidal magnetic fields are expanded into spherical harmonics up to the fourth degree and order. Furthermore, the scalar function sh P is expanded into a Taylor series with respect to the radius r in the vicinity of the mean radius b of the spherical shell up to the first order. The function sh T is cut off at the zeroth order for the radius describing the influence of the radial currents (Toepfer et al. 2021a) since it is expectable that the radial currents are the dominating poloidal currents in the vicinity of the surface. Thus, the magnetic field is modeled with 102 expansion coefficients in the data analysis, i.e., 15 internal Gauss coefficients, 15 external Gauss coefficients, 24 toroidal coefficients, 48 poloidal coefficients. As a proof of concept, the maximum orders of the series expansions are chosen to achieve a reasonable qualitative agreement between the simulated and the reconstructed currents. For a detailed quantitative analysis, it is worthwhile to incorporate higher orders of the series expansions. The wanted expansion coefficients are estimated with Capon's method (Capon 1969;Motschmann et al. 1996;Toepfer et al. 2020aToepfer et al. , b, 2021b from the simulated magnetic field data. Afterwards, the estimated expansion coefficients g sh P for the poloidal magnetic field and the coefficients g sh T for the toroidal magnetic field are used to reconstruct the current density j sh (cf. Equation 15) within the shell. The simulated and the reconstructed current density j sh along the ellipsoid E on the nightside and dayside of Mercury are displayed in Figs. 1 and 2, respectively. The simulated current system is dominated by equatorial currents flowing from dawn ( y > 0 ) to dusk ( y < 0 ) at the nightside (red arrows in Fig. 1a) and vice versa at the dayside (red arrows in Fig. 2a). The reconstructed equatorial currents (yellow arrows in Fig. 1b and green arrows in Fig. 2b) follow this geometry. At the nightside, the equatorial currents split into a current flowing towards Mercury at the dawnside and depart from the planet at the duskside. The polar regions are characterized by horizontal currents. At the northern pole, the currents penetrate the surface in the region y > 0 (green arrows in Fig. 2a, yellow arrows in Fig. 2b) and leave the planet at y < 0 (red arrows in Fig. 2a, red arrows in Fig. 2b). It should be noted that the reconstructed currents are calculated analytically and therefore, the reconstruction is smoother than the numerical simulation. However, the reconstruction method reproduces the structure of the simulated currents (cf. Figs. 1b and 2b) and thus, the Mie representation is capable of describing the geometrical nature of the currents flowing in the vicinity of the MPO orbit. Since the scalar functions sh P and sh T are expanded into spherical harmonics up to the fourth degree and order, it is expectable that the extension of the series expansion by degrees l > 4 will close the void between the reconstructed and the simulated magnitude of the current density.  Fig. 3 Decomposition of the reconstructed current density j sh into its poloidal and toroidal parts along the ellipsoid E on the nightside of Mercury. The poloidal current density is presented in multiples of j sh P0 = 25 nA/m 2 and the toroidal current density is given in multiples of j sh T 0 = 15 nA/m 2 Figure 3 displays the poloidal j sh P and toroidal parts j sh T of the reconstructed current density. The toroidal currents are characterized by closed loops with the divergence-free nature of j sh T being immediately visible. Since the radial (poloidal) currents cross the ellipsoid E , from the first point of view it is not obvious how these currents are closed. Moreover, the magnitude of the toroidal current density with a maximum value of j sh T ≈ 15 nA/m 2 at the northern pole is smaller than the magnitude of the poloidal current density ( j sh P ≈ 20 nA/m 2 up to 25 nA/m 2 at the northern pole), but both the magnitudes are of the same order.

Poloidal and toroidal current systems under different IMF-orientations
The simulated current density has reasonably been reconstructed from the simulated magnetic field data. Now, the question arises how the currents are related to the current system around Mercury presented in the section "Current system at Mercury" and how the poloidal currents are closed.
Besides the geometry of the internal magnetic field the current system depends on the direction of the interplanetary magnetic field (IMF) (Ganushkina et al. 2015;Milan et al. 2017;Ganushkina et al. 2018). To investigate the qualitative structure of the current system around Mercury, simulated stationary magnetic field data and current densities resulting from the plasma interaction of Mercury with the solar wind under different IMF-orientations, i.e., B IMF = ±20 nT e x , B IMF = ±20 nT e y , B IMF = ±20 nT e z , where e x , e y and e z are the unit vectors of the corresponding main axes in the MASO system, are analyzed. Although Mercury's magnetosphere is a highly dynamic system, the interplanetary magnetic field in the vicinity of Mercury's orbit can be regarded as stationary for a time period of 20-40 min in times of a calm upstream solar wind (He et al. 2017;James et al. 2017). The resulting poloidal and toroidal current systems on the nightside and on the dayside of Mercury for each IMF-direction are sketched in Figs. 4 and 5. The dayside current system is dominated by the Chapman-Ferraro currents j cf . These currents decompose into poloidal (orange) and toroidal currents (green). The toroidal part of the Chapman-Ferraro currents is characterized by closed loops so that this part of the current system is closed within the magnetospheric plasma. The nightside current system is dominated by the neutral sheet current j ns , simplified by just a single dawn-duskdirected arrow. The magnetopause current j mp as well as the neutral sheet current are connected via poloidal currents at the dawn ( y > 0 ) and duskside ( y < 0 ) so that both the currents are characterized by a poloidal topology. Since the structure of the magnetopause current and the neutral sheet current is determined by the internal magnetic field, these currents remain qualitatively unchanged for all IMF-directions. In the polar regions the poloidal part of the Chapman-Ferraro currents transits into poloidal Region 1 currents j R1 . Due to the high conductivity within Mercury's core, the Region 1 currents flowing in radial direction towards the planet at the dawnside and depart from the planet at the duskside are able to penetrate the surface and partially close via the core-mantle boundary, as proposed by Anderson et al. (2014). The currents flowing depart from the planet are closed within the nightside magnetosphere. Within the simulation presented here, no exosphere has been adopted. Considering the influence of Mercury's exosphere on the current system, a significant portion of the Region 1 currents can be closed within a sodium exosphere of sufficient density, in similarity to Earth's ionosphere, as shown by Exner et al. (2020).
Analogously to the plasma interaction of the Earth's magnetic field with the solar wind, there occur fieldaligned currents j fac at the northern and at the southern pole in the case of a non-vanishing B y -component of the IMF (cf. Figs. 4b and 5b) (Leontyev and Lyatsky 1974;Trondsen et al. 1999;Liou and Mitchell 2019). Due to the frozen-in theorem the motional electric field is given by In the case of B y = 0 ( B x = B z = 0 ), the electric field is orientated (anti-)parallel to the z-axis, since the solar wind velocity is orientated along the x-axis. We find that the field-aligned currents flow towards the planet at the northern pole and depart from the planet at the southern pole in the case of B y > 0 and vice versa in the case of B y < 0.
Furthermore, the IMF-direction determines the symmetry of the magnetosphere. A non-vanishing B y -component results in dawn-dusk asymmetries within the tail, whereas the B x -component influences the north-south symmetry of the magnetosphere.
In the case of B IMF = 20 nT e x as well as B IMF = ±20 nT e y (cf. Figs. 4a and b, 5b), within the reconstruction procedure there occur toroidal currents with an amplitude of about 10 nA/m 2 that are oppositely directed to the (poloidal) magnetopause current at the dayside of Mercury and oppositely directed to the (poloidal) neutral sheet current at the nightside. This behavior is founded on the mathematical decomposition of the current density as a vector field. First of all, the quantities j sh P and j sh T are mathematical vector fields, which do not mandatorily exist as physical quantities. Thus, a vanishing current at a point x 0 can equivalently be described as a superposition of oppositely directed non-vanishing (21) E IMF = −u sw × B IMF .
poloidal and toroidal currents with the same magnitude so that is valid locally as schematically sketched in Fig. 6a and b. Due to the finite spatial extent of the currents, (22) j sh P (x 0 ) + j sh T (x 0 ) = 0 at least partially they can flow independently of each other in regions this side of x 0 or beyond x 0 (cf. Fig. 6c).
Although the total (physical) current j is determined by the superposition of the poloidal and the toroidal current, the mathematical decomposition of the current density enables us to predict the potential origin of the current and to analyze the potential trajectories of the particles carrying the current. Thus, the toroidal current flowing antiparallel to the magnetopause current as well as antiparallel to the neutral sheet current can be interpreted as a ring current which superposes with the poloidal neutral sheet current and the magnetopause current. This ring current should be partially trackable along the planned MPO orbits.
In terms of the poloidal-toroidal decomposition, the current system around Mercury as sketched in Figs. 4 and 5 can be summarized as follows. The poloidal current (orange) flows towards Mercury at the dawnside ( y > 0 ) and splits into the neutral sheet current j ns , the (c) Independent currents.

Fig. 6
Superposition of poloidal (orange) and toroidal (green) currents. A vanishing current density j = 0 on the left side of the grey sphere (a) can equivalently be described as a superposition of non-vanishing poloidal and toroidal currents, so that j P + j T = 0 is valid locally (b). The currents can flow independently of each other in a certain spatial distance (c) Region 1 current j R1 and the dayside magnetopause current j mp . The Region 1 currents are able to penetrate the surface. These currents close via the core-mantle boundary and leave the planet at the duskside ( y < 0 ), where they reconnect with the magnetopause current and the neutral sheet current. This poloidal current system remains qualitatively similar for all IMF-directions. The toroidal Chapman-Ferraro currents j cf are characterized by closed loops (green). The position of these loops varies for different IMF-orientations. In the case of B IMF = 20 nT e x as well as B IMF = ±20 nT e y the equatorial parts of the loops at the dayside are orientated antiparallel to the magnetopause current as well as antiparallel to the neutral sheet current at the nightside. For a non-vanishing B y -component, the (poloidal) fieldaligned currents j fac at the northern and at the southern pole also penetrate the surface and close via the coremantle boundary.
Although the plasma interaction of Mercury with the solar wind does not linearily depend on the IMF-orientation, it is expectable that the resulting current system for any other IMF-orientation can be constructed by superposing the cases presented above.

Summary and outlook
The analysis of current systems is of major importance for the comprehension of Mercury's magnetosphere. Since the Mie representation delivers an analytical parametrization for rotational magnetic fields, it enables the analytical calculation of the current density flowing in regions that are crossed by the Mercury Planetary Orbiter (MPO) from the measured magnetic field data. As a proxy for the not yet available MPO data, the magnetic field in the vicinity of Mercury is simulated with the hybrid code AIKEF and the data are evaluated along the planned trajectories of the MPO. The comparison of the simulated and the reconstructed current density shows, that the reconstructed currents are in reasonable agreement with the simulated currents. Thus, the Mie representation serves as a useful model for the analysis of currents that are crossed by the spacecraft. Furthermore, the Mie representation enables the decomposition of the current density into its poloidal and toroidal parts. The reconstructed toroidal currents are characterized by closed loops, whereas the poloidal currents cross the envelope of the trajectories in radial direction. Especially, the mathematical decomposition of the current density enables us to predict the potential origin of the current and to analyze the potential physical trajectories of the particles carrying the current.
After the reconstruction procedure, the question arises how the poloidal currents are closed and how the reconstructed currents are related to a complete current system around Mercury. Since the structure of Mercury's magnetosphere depends on the orientation of the interplanetary magnetic field (IMF), the plasma interaction of Mercury with the solar wind is simulated for six different IMF-directions along the main axes in the MASO system. It turns out that the poloidal current system remains qualitatively unchanged since, for example, the magnetopause current as well as the neutral sheet current is mainly determined by Mercury's internal magnetic field. In the case of a non-vanishing B y -component of the IMF, there occur field-aligned currents in the polar regions which follow the motional electric field. Furthermore, the non-vanishing B y -component yields dawn-dusk asymmetries, whereas the existence of a non-vanishing B x -component violates the north-south symmetry.
Within the present study, the magnetic field data have been evaluated along a complete ellipsoid to reconstruct the current density for each IMF-orientation. Analyzing in situ magnetic field data, the current density can be reconstructed along segments of the MPO orbit for each IMF-orientation and the results can be classified in terms of the complete current system presented here. In addition, it should be noted that the Mie representation requires the series expansion of the magnetic field. Within the practical application, the series expansions have to be truncated at a suitable degree, where the choice of the maximum degree depends on the distribution and the amount of the available data points.
Furthermore, Mercury's exosphere has been neglected within the simulations so that the field-aligned currents close via the core-mantle boundary. Therefore, the influence of Mercury's exosphere on the poloidal and toroidal current systems should be analyzed in future studies. Moreover, it is worthwhile to investigate the influence of dynamical effects on the current system within Mercury's magnetosphere Slavin et al. 2021).
Concerning the BepiColombo mission, this work establishes the basis for the application of the Mie representation to calculate currents from the magnetic field measurements analytically and to analyze their poloidal and toroidal contributions. so that B = ∂ x × A is valid (Toepfer et al. 2021a). Alternatively, the magnetic field can be parametrized via Euler potentials α and β (Stern 1967(Stern , 1970

resulting in
The corresponding vector potential can be written as The Mie representation is based on the decomposition of the magnetic field with respect to spherical coordinates. Expanding Eq. (27) in spherical coordinates and comparing the coefficients with Eq. (25) delivers as well as Thus, the scalar functions P and T can be interpreted as a special case of the Euler potentials α and β.