Today there is no effective and universal mathematical model to calculate magnetic fields of the Earth, not to mention other planets of the Solar System. Even the most elaborated mechanism of electromagnetic dynamo still needs full mathematical description performed, and it also needs to be confirmed by successful laboratory experiments.

The physically justified method of calculation magnetic fields of planets, presented in this article, is a solution as a proper model of the magnetic fields of Earth and other planets of the Solar system is still not elaborated. Primary elements of the said method are described in the work of A.N. Dmitriev [1] .

The purpose of the present article is to consider the method of calculation magnetic fields of “hot” planets. Its effectiveness and advantages in comparison to the model of hydromagnetic dynamo were proved by finding out the reasons why the Jupiter’s magnetic field has such a complicated structure. The method is based on the quantum theory of electroconductivity of metals and semiconductors, which is oriented to calculate mainly thermal processes, always taking place inside the bowels of “hot” planets. Main elements of those processes are planetary geotemperature gradients, thermo EMF and radially directed thermoelectrical currents, which are associated with Seebeck effect. These currents are always in the bowels of “hot” planets, possess an enormous amount of energy, and they generate and maintain magnetic fields for a long time, mostly dipole magnetic field.

Today Jupiter is not studied as good as Earth is. The data that has been being obtained since 1972 using spacecrafts and the Hubble Space Telescope let us describe the inner structure of the planet just in broad outlines. But it is definitely that Jupiter has the most intensive magnetic field in the Solar System. Its intensity on the equator is about 4.2 G (4.2 ´ 10⁻⁴ T) on the visible clouds surface level, and it is 14 G (1.4 ´ 10⁻³ T) on the Jupiter’s pole [2] - [8] that is 12 times more than the magnetic field on the Earth’s equator; furthermore, it is even 5 - 10 times more than the polar (dipole) magnetic field of the Earth. There are some theories of the Jupiter’s magnetic field nature. The main one is that there must be a dynamo, similar to the Earth’s one, but, as scientists believe, the Jupiter’s electricity conductor is a mostly hydrogen layer, unlike the Earth’s Fe-Ni core [6] .

But can we believe this point of view to be the only one right, considering recently achieved modern data? Anyway, to find the answer we must know the inner structure of the planet. It is possible to make a more accurate and modern model of the inner structure of the Jupiter (Figure 1(а)) with the available information, including space exploration, computer modeling of the Jupiter’s magnetic field [9] [10] , mathematical modeling of helium physical properties at high р and Т [11] and physical experiments to study electrical properties of H and He at high р and Т [12] . In particular, it has a gas layer (atmosphere), a liquid molecular hydrogen layer, a layer of liquid hydrogen with metal properties (hereinafter Me_H), stone core of ice and solid rocks [13] . As for the latter, there is

an assumption that the core’s diameter can be 2 times larger than it was considered to be and be 29 ´ 10³ km [7] , according to the results of a recent computer modeling [9] [10] . Besides, in the center of that core there may be another small inner core of heavy metals, probably iron and nickel [9] [10] . In this case the model of the Jupiter’s inner structure can be viewed as a spherical five-layer model.

That five-layer model contains levels which consist of mostly conductive substance and contain electrons of different concentrations. The third layer of liquid metal hydrogen and the metal core should be referred to those levels. At certain favorable conditions those levels are the levels where big thermoelectrical currents can flow, which can both create and maintain the Jupiter’s magnetic field for a long time. In particular, planetary temperature gradient, that can maintain difference of temperatures ΔТ = (20000 − 10000) К = 10000 К between the inner and the outer shells, is such a favorable condition for the Me_H layer [14] [13] [8] [10] . That layer is also characterized by a notable gradient of pressure from 200 - 300 GPa (the upper shell) to 4000 - 4500 GPa (the lower shell) [7] .

For the radially directed movement of electrons, i.e. thermoelectrical currents, there must be electrical potential gradient of the same orientation. But the latter can happen if there is an external so-called potential-forming factor or, in other words, an external EMF of non-electric origin. In a planet’s body function of the latter is performed by the difference o temperatures ΔТ, that is radially orientated, according to the temperature gradient. At the same time, there is the Seebeck effect, which happens when temperature difference ΔТ is attached to the both of the points of an electronic conductor, what causes potential difference Δ φ T between those points because “hot” electrons move to the “cold” point of the conductor [15] [16] :

Δ φ T = ∫ T 1 T 2 β d T (1)

where β is temperature coefficient, which is 0.0001 V/deg for pure metals and 0.0015 V/deg for some semiconductors [16] .

Hence, unidirectional effect of the Jupiter’s planetary temperature gradient, in the Me H layer between its inner and outer spherical borders in particular, makes electric intensity E Me H Т stay lasting between those borders:

E Me H Т = Δ φ T / l Me H ,   V / m ; (2)

where l Me H is the spherical layer thickness of metal hydrogen, m.

The Me_H layer thickness corresponds to the distance between the solid core surface and the outer core of the layer, it is 44500 km. So, if ΔТ is 10000 K, it defines E Me H Т as 2.25 ´ 10⁻⁸ V/deg.

Then E Me H Т makes free electrons forcibly move in radial direction, that generates flows of spherically directed thermoelectrical currents I Me H moving from the lower warmer part of the layer to its cold part. At the same time, both gravitation and inertial Coriolis forces begin to affect electrons. Those forces make electrons change radial direction into the side opposite to rotation, the more intensively the closer to the upper border of the Me H layer electrons are (Figure 1(b)). As for gravitational forces, centrifugal force F centr = m w 2 R i in particular, the latter always diverts electrons both to the axis of planet’s rotation and oppositely (Figure 2(a)). The quality F_centr is proportional to the radius-vector R_i, which is represented by different lengths of horizontal arrows on Figure 2.

Directed currents of thermoelectrical currents I Me H cause toroidal magnetic field inside a spherically enclosed conductor, which is, in its turn, generates a poloidal field, observable beyond that enclosed conductor. But is the I Me H energy enough to generate the Jupiter’s magnetic field, which is 4 ´ 10⁻⁴ T intensive on the equator and 1.4 ´ 10⁻³ T intensive on the pole?

To answer this question one should firstly calculate density of the thermoelectrical current j Me H that appears just on 1 m² on the surface of the shell of the Me_H layer. According to the quantum theory of metal electrical conductivity, current j density value is defined by [16] - [21] :

j Me H = e 2 n E Me H Т l p F (3)

where e is effective electron charge, C (1.6 ´ 10⁻¹⁹); n is quantity (concentration) of electrons in a volume unit of a metal conductor, m⁻³, assuming that metal hydrogen has all the properties of that conductor.

Then n is defined by:

n = ( 5 m е ℏ 2 ( 3 π 2 ) − 2 / 3 ) 3 / 5 p 3 / 5 (4)

where l is the average distance electrons run to the E T electric field between concussions, it is defined by:

l = σ ℏ ( 3 π 2 ) 1 / 3 e 2 n 2 / 3 ,   m (5)

m e is mass of electron, 9.11 ´ 10⁻³¹ kg;

ℏ is Planck constant, 1.05 ´ 10⁻³⁴ J×s;

p is pressure in the metal hydrogen layer, P;

σ is electrical conductivity of the metal hydrogen layer Me_H, S/m;

p F is Fermi momentum, which is defined as:

p F = ℏ ( 3 π 2 n ) 1 3 (6)

Formula (4) defines concentration of electrons in metals, generally, of great density, whereas density of liquid hydrogen is 70.8 kg/m³ [22] . Therefore if we take metal hydrogen as liquid metal in (4) we should make a correction for low density of metal hydrogen b = ρ Me H / ρ Fe empirically (Fe is used as most common in cores of planets). Then:

n = b ( 5 m е ℏ 2 ( 3 π 2 ) − 2 / 3 ) 3 / 5 p 3 / 5 (7)

If we take boundary pressures in the Me_H layer from 200 GPa to 4500 GPa considering [7] [13] [23] and use them as substitutes in (7), we will get average concentration of electrons in a unit of layer volume n = 5.64 × 10 27 m − 3 . Constant b in the calculations was taken as 0.009 (70.8/7870 kg/m³). However despite using the b correction the derived quantity n cannot be objectively used for the Me_H layer, because average density of Jupiter generally is just 1.326 g/cm³ [14] [23] [7] . That is why it may be better to use classical electronic theory of metals where moving of electrons obeys the laws of statistical methods and n is determined by formula [16] :

n = ρ N A A , (8)

where ρ is density of metal hydrogen, 70.8 kg/m³;

N A is Avogadro constant, 6.02 ´ 10²³ mol⁻¹;

А is molar mass of hydrogen, 0.0010079 kg/mol.

If there is one free electron for an atom, which corresponds to monoatomic hydrogen H, then (8) results in n = 4.23 × 10 22 m − 3 .

At the same time, to calculate l , j and I Me H , we must know the electrical conductivity σ of the Me H layer. However this parameter is still unknown. Physical experiments are held nowadays [24] [25] [26] [27] , but specific and repeatable results cannot be gained because they cannot be stored for a long time. There are first efforts in mathematical modeling of metal helium, results of which inspire the authors that there may be successful prospering of electrical parameters [11] .

That is why we need to use a logical reasoning. The layer cannot be superconducting, because it would require temperature close to absolute zero. At the same time that layer cannot be high conductive as usual metals (10⁵ - 10⁷ S/m), because of its liquid state with extremely low density and, particularly, very high temperature in the planet’s body. But the layer can have heightened electrical conductivity (20 - 10³ S/m) even compared to semiconductors (0.1 - 10 S/m). So, the start point of electrical conductivity for the Me_H layer can be, for instance, σ = 470   S / m , low value of which corresponds to the tendency of H and He to have quite low electrical conductivity at low р, Т ( σ = 1.94   S / m ), that is proved by a physical experiment [12] . If we put the above value of electrical conductivity into formula (5) and substitute the calculated values of E Me H Т , n and p F into (3), we obtain the magnitude of thermoelectrical current density j Me H = 1.06 × 10 − 5 A / m 2 .

Then thermoelectrical current #Math_48#, that goes through the outer shell S Me H of the third Me H layer, is defined by:

I Me H = j ⋅ S Me H , (9)

where S Me H = 4 π R Me H 2 and it is 4.15 ´ 10¹⁶ m² at R Me H = 5.75 × 10 7   m .

Then we can find the final total current magnitude I Me H = 4.39 × 10 11   A from (9).

If we know I j u and geometrical parameters of the planet, we should use the Biot-Savart formula to calculate the magnetic field B Me H of the Me_H layer. However the said formula, that describes magnetic induction for a spherical conductor with electric current moving from the sphere center to its outer surface, is unknown to the author. To obtain it we can use the following model: a small spherical conductor with radius R Σ is put inside a spherical conductor with radius d R i . The currents, moving from the center of the sphere, make a magnetic field inside the volume of the small sphere. This magnetic field is always directed to the small sphere’s volume surface tangentially. Then, according to circuital theorem, we get

d B j u 4 3 π d R i 3 = μ 0 4 π I i 4 π d R i 2

where the I i currents are currents that move through the surface of the small sphere 4 π d R i 2 . If current spreads over the sphere section evenly we have

I i = I Me H ( d R i R ) 2 . Hence, shortening and substituting we get Biot-Savart formula

for calculating magnetic field induction inside a sphere (metal core). Here it is, rationalized:

#Math_63# or B j u = 3 μ 0 I Me H R Me H 4 π R 2 sin α , (10)

where R is radius of Jupiter, R Me H is radius of outer spherical surface of the Me_H layer, d B j u = d B Me H . Because of spherical symmetry angle α = 90˚ at any point of outer spherical surface when we calculate B Me H .

When calculating magnetic field on the planet surface, the spherical Me_H layer, in its turn, depends not on the radius of that layer R Me H , but on its thickness R Me H . Then (10) changes as follows:

B Me H = 3 μ 0 I j u Δ R Me H 4 π R 2 sin α (11)

Putting the appropriate meanings into (11), we can see that magnetic induction B Me H e on the equator is 1.15 ´ 10⁻³ T at 1Rj = R e = 7 .15 × 10 7 m , while the magnetic induction on the pole B Me H p is 1.23 ´ 10⁻³ T at R = R e = 6.68 × 10 7   m .

According to the data obtained by spacecrafts via measuring magnetic field, magnetic field changes at the equator from 3.2 - 4.3 G to 13 - 14.3 G on the poles (let us take B j u p as 1.4 ´ 10⁻³ T) [2] [8] [28] [29] [30] [31] . However, it is not always clear from the quoted literature to what point of the Jupiter’s equator the measured values 3.2 - 4.3 G should be referred, whether they should be referred to the Jupiter’s surface for 1Rj or to the remote points for nRj. As it appears, the most reliable data can be derived from the connection graph B n = f ( n Rj ) , where n is measured from 2 to 20 [2] . In particular, the value of 4.2 - 4.3 G on this graph corresponds to 3Rj distance. That is why hereinafter we use B j u e which is 4.2 ´ 10⁻⁴ T. It is noteworthy that magnetic field at poles is not homogeneous, it is 1.4 ´ 10⁻³ T on the North pole and 1.07 ´ 10⁻³ T on the South pole of the Jupiter [32] .

To watch the Jupiter’s magnetic field calculations and its results comfortably here is Table 1.

It shows that the calculated value of magnetic field on the Jupiter’s equator B Me H e is almost 3 times more than the measured one B j u e , while on the pole the measured magnetic induction B Me H p differs from the calculated one B j u p by a

little bit more than 12%. However this quite coincident result cannot be taken as the final one, because the planet has additional resources of magnetic fields, metal core in particular.

To calculate the magnetic field produced by the core its main physical parameters should be defined: a) difference of temperatures between the planet’s center and the core’s outer shell ΔT_С, predicted by the known data [9] [10] [13] [14] and b) electroconductivity σ ranging from 0.1 to 1 million S/m [33] [34] [35] .

According to p. a), let us take Т₁ = 20000 K at the outer border of the core and Т₂ = 23000 K in the center of the corner, that will make ΔT_С be 3000 K. Comparing to thermal conductivity of pure metals this value can be real, if we agree that the planet’s core can consist predominantly of iron carbide Fe₇C₃ just like on Earth [36] with available data of abnormally high р, Т [37] . According to p. b) let us take the electroconductivity of the Jupiter’s core σ = 1.1 × 10 5   S / m , which is close to the minimal Earth one and can be explained due to a great amount of carbon.

Now let us calculate electric and magnetic parameters of the solid core, which has temperature dynamics. To do this, let us use the mentioned formulas (2 - 3, 5 - 7, 9, 11).

In particular, we must use the initial data when calculating:

Т of the metal core layer: outer shell/core center - 20000/23000 K;

Metal core radius, R C = 1.0 × 10 7 m ;

Radius of the core with the rock shell, R C + S = 1.3 × 10 7 m ;

Area of the core spherical surface with radius S C = 4 π R C 2 = 1.26 × 10 15 m ; Equatorial radius of Jupiter, R e = 7 .15 × 10 7 m ;

Polar radius of Jupiter, R p = 6.68 × 10 7 m ;

Core electrical conductivity, σ = 110000 S/m;

Then

Thermofield intensity E C T in the core when ΔT_С = 3000 К, can be found from the following formula:

E C T = β * Δ T С / R C (12)

where β = 0.000 1   V / deg .

From the formula (7) is сoncentration of electrons in the core when pressure is 4.5 ´ 10¹² GPa, n = 1.48 × 10 30 m − 3 ; from the formulas (6, 5) is Fermi momentum and electron mean free path, respectively p F = 3.72 × 10 − 24 m ⋅ kg / s , l = 0.0108   nm . To calculate the current density, use the formula:

j C = e 2 n E C T l p F , (13)

which is defined by a magnitude of 3.3 ´ 10⁻³ A/m². In its turn, total current value through its outer surface S C , I C calculated by the formula:

I C = j C S C (14)

Then find the magnetic induction on the Jupiter’s pole B C p :

B C p = 3 μ 0 I C R C 4 π R p 2 sin α (15)

and magnetic induction on the Jupiter’s equator B C e :

B C e = 3 μ 0 I C R C 4 π R e 2 sin α , (16)

which are defined by 2.79 ´ 10⁻³ T and 2.43 ´ 10⁻³ T respectively.

For the readers’ convenience, all the calculations are presented in the compact Table 2.

Table 2 shows that the calculated Jupiter’s core magnetic field on the pole B C p is 2 times more that the measured one B j u p . At the same time that magnetic

field on the equator B C e is now 6.1 times more than the measured one B j u e . If we suppose that both thermoelectrical currents of the core and of the liquid metal hydrogen layer are directed to the surface of the planet, then the calculated total magnetic fields on the pole and on the equator are B j u p p = B Me H p + B C p = 4.01 × 10 − 3   T and B j u e e = B Me H e + B C e = 3.58 × 10 − 3   T respectively. However these values are way greater than the magnetic fields, measured on the Jupiter’s pole and equator. That may be because there are some mistakes in the calculation or because there are thermoelectrical currents moving in the opposite directions.

The latter is possible if the outer rock shell of the Jupiter’s core is not just a transitive zone but a radiation zone, where the temperature of rocks is maximum towards the inner core and the layer, just like on Earth. In this case thermoelectrical current of the metal core must be directed to its center (Figure 2(b)), and the total Jupiter’s magnetic field can be defined by the residual B j u p p = B Me H p + ( − B C p ) = − 1.56 × 10 − 3   T , which is close to the measured B j u p by modulus. In its turn, B j u e e = B Me H e + B C e is −1.29 ´ 10⁻³ T.

As we can see, the calculated total magnetic fields are negative because | − B C p | > B Me H p and | − B C e | > B Me H e . Hence thermoelectrical current of the core, depending on its direction, determines locations of the N and S poles of the Jupiter’s magnetic field. In particular, as it is shown on the Figure 2(b), the latter coincide with the planet’s geographical poles of the same names. It is known that at the present time the orientation of the poles of the Jupiter’s magnetic field is opposite to the orientation of the poles the Earth’s magnetic field, that was previously found by spacecrafts Pioneer 10 and Pioneer 11 and other [7] .

Earth is now in the state of planetary warming, caused by flows of the core and low mantel’s thermoelectrical currents moving to the same direction of the Earth’s surface. Because of the oppositely directed thermocurrents of the core and of the metal hydrogen layer, Jupiter now is in the state of planetary cooling. This is proved by research showing that the Earth’s inner core now tends to cool down at 50 K/10⁹ years [38] , i.e. 7.2˚С/10⁹ years. And this is possible when thermoelectrical currents of the core leave it and begin moving to the lower mantel.

Let us look at the calculated magnetic induction on the equator B j u e e again, it looks like a glaring error in calculations that its value differs almost 6 times from the measured one #Math_161#. However, there is data from Pioneer-10 showing that magnetic field density on the equator at a distance of 3Rj is about 4.2 G [2] . Therefore magnetic induction of the Me_H layer B Me H e and magnetic induction of the metal core B C e should be recalculated at a distance of three radii of Jupiter. Using formula (11) again we get B Me H e = 1.27 × 10 − 4   T and B C e = − 2.7 × 10 − 4   T . If we sum up these components considering the sign of B C e , we get a rough magnitude B j u e = − 1.43 × 10 − 4   T at a distance of 3Rj, which is almost 3 times less than the calculated one B j u e . This difference is quite insignificant, because the existing data of almost all the Jupiter’s parameters is approximate. That is why the value B j u e e = − 1.29 × 10 − 3   T at 1Rj should be regarded as quite admissible and close to the real one. But summing up B Me H e and B C e by modulus gives a value of magnetic induction at a distance of 3Rj as 3.98 ´ 10⁻⁴ T, which is quite close to the calculated one 4.2 ´ 10⁻⁴ T. The latter should be regarded as a coinsidence.

Aside from the studied autonomous sources of the Jupiter’s core and Me_H layer own magnetic fields, there is at least one more source that can be related to the central part of the planet―equatorial band, which rotates with increased velocity (rotation period is 9 hours 50 minutes and 30 seconds) in comparison to other sectors of the planet, which rotation period is 9 hours 55 minutes and 41 seconds [39] .

The increased velocity of the gas clouds is obviously transmitted to the Ме_Н layer and to some part of the liquid metal hydrogen layer because of viscosity of gases and liquids. If we suppose that the Ме_Н layer viscosity is not great, then three spherical layers can be involved into the equatorial zone as a sector, which geometrical dimensions are shown on the Figure 1(a). Electrical current occurs only in the sectoral part of the Ме_Н layer because of the spinning free electrons, that is why it is possible to find current density j e z and amperage I e z for this part of the Ме_Н layer [16] [40] :

j e z = ρ e z v , (12)

where ρ e z is volume electric charge density, C/m³, it can be defined by:

ρ e z = e n (13)

When we know the electron charge e and the previously calculated by the (8) formula electrons concentration per unit volume for the Ме_Н layer and calculate linear velocity of the sector for its medial part, we get j e z = 4.23 × 10 7   A / m 2 from the (12) formula. In its turn, current value can be found from:

I e z = j e z S e z , (14)

where S e z is the section of the spherical sectoral circle in the φ angle alignment. Can be found from:

S e z = R e z L S sin ( π 180 φ ) , (15)

where R e z is radius of the equatorial band sector with maximum rotation. It is 4.45 ´ 10⁷ m, what is equal to the Ме_Н layer’s thickness;

L S is the chord of the L_S sector (see Figure 1(a)), which is 18.5 ´ 10⁷ m;

φ is the central angle of the Jupiter’s equatorial band visibility, deg.

Having substituted the above numbers into (15) we can derive the value of the current in the Ме_Н layer sectoral circle from (14), I e z = 6.16 × 10 19   A . In its turn, magnetic induction on the Jupiter’s equator B e z e can be defined from [15] [16] :

B e z e = μ 0 e n S e z 4 π ( R e − R e z ) 2 sin α (16)

Its value turned out to be 3.14 ´ 10⁻² T, that is more than 22 times greater than the calculated B j u e on the planet’s equator. It definitely indicates that the Ме_Н layer cannot be concatenated to the upper part of the planet’s equatorial band to the entire depth. That is why we can suppose that because of the harder condition of the Ме_Н layer (increased density) concatenation can be accepted only within its thickness d e z not more than 250 m. Then, having calculated j ′ e z , S ′ e z and I ′ e z parameters according to the above formulas we can find the value of the magnetic field for the sectoral plate circle from:

B d e = μ 0 e n S ′ e z 4 π ( R e − R Me H ) 2 sin α (17)

This magnitude is 1.42 ´ 10⁻⁶ T on the equator, at the same time B d p on the pole is 3.58 ´ 10⁻⁸ T, if ( R e − R Me H ) 2 would be changed into R d 2 = R e 2 + R p 2 in (17). The gathered constituents of magnetic field B d e and B d p for the spherical circle plate of the Ме_Н layer are way less than the calculated values of the Jupiter’s magnetic fields, so they can be related as contingently real. First, this magnetic field is narrowly located along the circle, because at such distances as 1Rj the constituents B d e and B d p tend to zero. At the same time, the value of the current I ′ e z is 10³ - 10⁴ times greater than the currents I Me H of the Ме_Н layer and the I C core. It clarifies the reason why there is such colossal energy of lightning discharges, that can be observed in the middle part of Jupiter and are several times more powerful than the Earth’s ones [41] . Actually, the calculations show that the thermocurrents of the Earth’s core are not greater than 3.37 ´ 10⁹ А [1] . Besides, characteristic instability of the way the equatorial band rotates toward the neighbor zones is an obvious cause of equatorial jerks.

The calculations of the magnetic field for the sectoral plate circle are brought together into Table 3 for the results to be easily accessible.

As for the direction of the vector of a rotating circle magnetic induction, according to the right-hand rule, it is the same as the direction of the planet’s core magnetic field, i.e. lines of magnetic induction comes out of the North pole and go to the South pole (Figure 3).

There are still studied only three autonomic Jupiter’s magnetic fields with different directions of their field lines, and it is quite interesting to predict the structure of total magnetic field B j u . The latter is determined by the orientation of resultant vectors B₁ and B₂ of the magnetic induction when summing up the vectors of the Ме_Н layer magnetic field D₁ and of the core magnetic field C₁ and, respectively, D₂C₂ in the North half of Jupiter (Figure 3).

It is clear that B₁ and B₂ vectors extend field lines oppositely to the North pole. Otherwise, the resultant vectors of magnetic induction B₃ and B₄ at summing up the D₃C₃ vectors and the D₄C₄ respectively in the South half of the planet draw the field lines to the South pole. As a result, the external outline of the resulting

Jupiter’s magnetic field B j u should have a trapezoidal or a “sauser” form.

Previously there was an opinion on the Jupiter’s magnetic field structure that the reason why the field lines are stretched to the 20Rj distance, forming the magnetic field structure as a thin pancake is that both the centrifugal force and the thermal pressure of the plasma affect the plasma [42] . Our calculations show another reason, though it can be an addition to the said opinion.

Besides, if we look at the vector diagrams D₁B₁C₁ and D₂B₂C₂ it gets clear that the B₁ и B₂ vectors cause “pulling out” field lines with higher energy in the North pole area. At the same time the B₃ and B₄ vectors cause “dragging in” field lines

with lower energy through the South pole. That is why the value of magnetic induction at the Jupiter’s North pole is expected to be higher while the one at the South pole must be lower, i.e. there must be pseudo anisotropy of the planet’s magnetic field. The latter coincides with the known measurements of the magnetic field’s densities at the planet’s poles, which are 14 G at the North pole and 10.7 G at the South pole [32] .

Due to lack of accurate knowledge of main forming parameters’ values ΔT_С, R_C, σ, oscillation range of the Jupiter’s magnetic field intensity in dependence of those parameters’ values changes should be carried on. Those very parameters determine the value of density j and amperage I, which form the needed intensity of the planet’s magnetic field B. When there is no temperature difference in the core there are no thermo e.m.f and, accordingly, there are no thermocurrents, even if the core is metal. Both current density j and amperage I cannot be enough to maintain the needed level of the magnetic field В at low electrical conductivity σ, even if ΔT_С and R_C are of great value. At the same time importance of the metal core radius R_C increases significantly only when the quantity of electrons in the core (or in other conductive shells) is enough to generate the needed flow of thermocurrents, which generate the corresponding magnetic fields.

The calculations of the studied parameters are brought together to the Table 4. The algorithm of the calculations was as follows: when, for instance, the first parameter of the discretization value changed, the magnetic induction B C p was fully recalculated. That allowed to take into consideration effects of every discrete on the whole system, making the evaluation of each parameter of the system objective.

Table 4 shows that B C p increases proportionally at uniform increase of temperature difference ΔT_С with 1000˚С step. If discretization step decrease to, for instance, 100 K (deviation 4) then proportionality for B C p stays the same and has mutual linear connection B C p = 0.000 3 Δ T С + 0.00 25 at every discretization of ΔT_С. As the radius R_С increases linearly, the magnetic induction increases according to quadratic polynomial B C p = 3 × 10 − 5 × R C 2 + 0.000 5 R С + 0.00 23 , i.e. at every new magnitude of R_С the value of B C p always gets more than two times larger (or smaller). Here is also a connection between electroconductivity of metals and magnetic induction according to the dependence: B C p = 0.000 9 σ +   0.00 19 .

The defined connections show that the magnetic field intensity is closely connected to the main forming parameters―even a slight change of them enhance (or reduce) the Jupiter’s magnetic induction B C p instantly. So, if you need the calculated magnitude of magnetic induction B C p coincide accurately with the one measured by spacecrafts B j u p , you have to sort out all the forming parameters ΔT_С , R_C , σ when there is no accurate data of their values. Calculations show that sorting out helps to exclude formal (deliberate) results with changes of one or another parameter. When the quality B C p = B j u p for the forming parameters is achieved, their optimal magnitudes can be defined.

The conclusion is that in the future if there are magnitudes of a planet’s magnetic field elements of any exactness, then solving the inverse problem lets one predict the inner structure of that planet with the same exactness, but only for its electrically conductive elements (shells, layers, zones).

The approbation of the suggested method of calculation magnetic fields of planets, of Jupiter in particular, results in the following:

- Autonomous sources of the Jupiter’s magnetic field are the spherical liquid metal hydrogen layer, the metal core and the spherical circle of liquid metal hydrogen external part plate rotating in the alignment of Jupiter’s equatorial band with maximal speed.

- The Jupiter’s main magnetic field is the magnetic field of the metal core B C , the main additional one is the magnetic field of liquid metal hydrogen layer B Me H , the accompanying one is the local magnetic field of the Ме_Н layer’s circular band (plate) B d .

- Orientation of main magnetic field B C poles fits Jupiter’s geographical poles.

Besides, the performed calculations made it possible to evaluate the most possible magnitude range of electroconductivity of both liquid metal hydrogen layer and the Jupiter’s metal core for the first time. The reason why the Jupiter’s magnetic fields orientation is opposite to the Earth’s one is also explained. For the first time there is a supposition that cores of electronically conducting planets, and even stars, have universal structures, because there is always an intermediary radiation band with most heated and melted masses of matter between the “core” and the “mantel”.

Even though there are significant differences between the existing models of the Jupiter’s inner structure, the introduced method of calculating magnetic fields of planets turned out to be not only working for Jupiter, but able to give a reasonable explanation to a lot of obscure questions regarding the structure of its complicated magnetic field. The results of the performed mathematical modeling are expected to be confirmed by the Jupiter Polar Orbiter Juno soon. The inner structure of the planet, regarding its electrical conductivity elements (the metal core and the liquid metal hydrogen layer), is going to be specified significantly with help of the expected detailed data and the introduced method.

In conclusion, it could be considered that the method of calculating planets’ magnetic fields generally corresponds the meaning of the thermoelectrical model of those planets’ magnetic fields. The latter, in comparison to the magnetodynamo model, has undoubted advantages, because of its simpler physical and mathematical description. The thermoelectrical model can also explain some peculiarities of the Earth and Jupiter’s magnetic fields conduct.

Thanks to Abramovich S.А. for checking mathematical calculations and Proskurnina T.V. for efforts to translate the manuscript to English.

The manuscript “To the question about the nature and structure of Jupiter’s magnetic field” was not published in any magazine and it is not being edited in any editorial office.

An article close by meaning was published in 2016:

Dmitriev A.N. A New Look on the Nature of the Earth’s Magnetic Field//IOP Conf. Series: Earth and Environmental Science 44 (2016) 022001. https://doi.org/10.1088/1755-1315/44/2/022001

I made the manuscript on my own in the Tyumen Industrial University, where I am working. Here is no competing financial or any other interest of any organization as I did not get any subsidies, grants or ordinary shares. I, A.N. Dmitriev, am the only author of the manuscript and there is no legal claim to me from anyone.

The manuscript author, A.N. Dmitriev.

Dmitriev, A.N. (2017) To the Question about the Nature and Structure of Jupiter’s Magnetic Field. International Journal of Geosciences, 8, 1091-1110. https://doi.org/10.4236/ijg.2017.89062