Solar photosphere magnetization

Context. A recent review shows that observations performed with different telescopes, spectral lines, and interpretation methods all agree about a vertical magnetic field gradient in solar active regions on the order of 3 G km−1, when a horizontal magnetic field gradient of only 0.3 G km−1 is found. This represents an inexplicable discrepancy with respect to the divB = 0 law. Aims. The objective of this paper is to explain these observations through the law B = μ0 (H + M) in magnetized media. Methods. Magnetization is due to plasma diamagnetism, which results from the spiral motion of free electrons or charges about the magnetic field. Their usual photospheric densities lead to very weak magnetization M, four orders of magnitude lower than H. It is then assumed that electrons escape from the solar interior, where their thermal velocity is much higher than the escape velocity, in spite of the effect of protons. They escape from lower layers in a quasi-static spreading, and accumulate in the photosphere. By evaluating the magnetic energy of an elementary atom embedded in the magnetized medium obeying the macroscopic law B = μ0 (H + M), it is shown that the Zeeman Hamiltonian is due to the effect of H. Thus, what is measured is H. Results. The decrease in density with height is responsible for non-zero divergence of M, which is compensated for by the divergence of H, in order to ensure div B = 0. The behavior of the observed quantities is recovered. Conclusions. The problem of the divergence of the observed magnetic field in solar active regions finally reveals evidence of electron accumulation in the solar photosphere. This is not the case of the heavier protons, which remain in lower layers. An electric field would thus be present in the solar interior, but as the total charge remains negligible, no electric field or effect would result outside the star.


Introduction
The problem of the large magnitude difference between the observed horizontal and vertical magnetic field gradients in and around sunspots has been known for a long time. Balthasar (2018) wrote a detailed review of observations, where it is shown that typical values of 3 G km −1 and 0.3 G km −1 are obtained for the vertical and horizontal gradients of the magnetic field, respectively, regardless of the telescope, spectral line(s), and measurement interpretation method used. Using these values would surprisingly lead to a non-zero divergence of the observed magnetic field, which is a priori not acceptable. In his review "Sunspots: An overview", Solanki (2003, p. 184) expressed the problem as follows: No satisfactory solution has been found as yet for the unexpectedly small vertical gradients obtained by applying the div B = 0 condition [to the observed horizontal gradients]. Balthasar (2018) tries to explain the discrepancy by simulating unresolved magnetic structures. This does not manage to fully explain the observations. In Sect. 2 we rule out the effect of measurement inaccuracies. In Sect. 3, we rule out a pure effect of spatial resolution on the basis of two mathematical theorems. We discuss the difference between derivatives and finite differences (the observations), and we show that these mathematical theorems prove that the observed non-zero divergence reveals the existence of at least one non-zero contribution from an element of the averaged region.
In this paper we explain the observations by applying the Maxwell relation in magnetized media B = µ 0 (H + M), provided that it can be proved that what is measured by the Zeeman effect is H and not B. This is the object of Sect. 5. To establish this point it is necessary to go back to the microscopic scale of an atom embedded in the medium magnetized at the macroscopic scale. In addition, we study the atom potential magnetic energy. As introduced in Sect. 4, in this paper we denote as B the magnetic induction and as H the magnetic field, which are related by the law B = µ 0 (H + M), where M is the magnetization.
However, the photosphere magnetization is negligible when evaluated via the usual models. In Sect. 6 we propose a model that is able to increase magnetization in the photosphere. It is considered that in the solar interior, at 0.5 R , the electron thermal velocity of 12 Mm s −1 largely surpasses the escape velocity of 850 km s −1 ; however, this is not the case for protons whose thermal velocity is 290 km s −1 , only due to their much higher mass. A similar effect occurs in the solar corona (Meyer-Vernet 2007). Following Allen (1973), the mass inside the sphere of 0.5 R radius is 0.94 of the total solar mass, and the temperature is 3.4 × 10 6 K there. This leads to the above velocity values. Thus, gravity separates the charges, in the sense that electrons escape when protons do not. More precisely, the underlying protons do not completely prevent electrons from escaping. The result is a free electron accumulation and a space charge at the Sun's surface. We show that the electron migration is a quasi-static process. As a result, the photosphere magnetization due to the plasma diamagnetism, which is itself due to these charges, could become non-negligible. Under the usual Fig. 1. Histograms of inversion results of noised theoretical polarization profiles. The introduced noise level is 1.5 × 10 −3 , in terms of polarization degree. The method used is a Milne-Eddington inversion, but including a magnetic filling factor α (see text). As a result, the inaccuracy can be estimated to be 20 G for the x and y components (left and middle plots), whereas it is 10 G for the z component (right plot). The Oz-axis lies along the line of sight. The assumed line is Fe i 6302.5 Å. These histograms were published (in spherical coordinates) in Fig. 5 of Bommier et al. (2007). photospheric conditions (the VALC model; Vernazza et al. 1981), the magnetization would be about 10 −4 lower than B/µ 0 , due to the low electron and ion densities (Bommier 2015). The accumulation of escaping electrons could increase the magnetization up to a non-negligible level.
In Sect. 7 we present the most direct electron density measurements in the solar photosphere to our knowledge. We propose that these measurements show evidence of an electron overdensity in the solar photosphere. We conclude in Sect. 8. Balthasar (2018) published a review of observations of the magnetic field gradients in and around sunspots. The measurement inaccuracy obviously depends on the method and the instrument used. However, we find it interesting to publish here the histograms of magnetic inaccuracy obtained within the UNNOFIT inversion method that we recently developed (Bommier et al. 2007), in magnetic field cartesian coordinates. These histograms can be found in Fig. 5 of the above-mentioned paper, but in spherical coordinates, which prevents us from clearly discriminating between longitudinal and transverse field components. In Fig. 1 of the present paper, we publish the histograms for each of these components. As described in Bommier et al. (2007), these histograms were obtained from inversion of theoretical noised polarization profiles. The noise level was assumed to be 1.5 × 10 −3 , in terms of polarization degree. The magnetic field strength ranged from 100 to 3000 G. The inversion method is of the Milne-Eddington type, but assuming a two-component atmosphere represented with a magnetic filling factor α. The originality of our method lies in the fact that the magnetic filling factor is entered as the ninth free parameter in the Levenberg-Marquardt algorithm. Figure 1 shows that the magnetic inaccuracy is on the order of 20 G for the transverse component, while it is on the order of 10 G for the longitudinal component. We applied our method to observations of the Fe i 6301.5 and 6302.5 Å lines, which are formed with a height difference of about 100 km in active regions, as shown in Fig. 4 of Khomenko & Collados (2007). As the field strength there is on the order of one thousand gauss, with a difference of about 300 G between the two line measurements, leading to the vertical gradient value of 3 G km −1 , as reported by Balthasar (2018), the 10 G inaccuracy on the longitudinal field is much smaller than the finite difference of the measurements. The horizontal gradient is found to be only about 0.3 G km −1 , which corresponds to the reversal of a 1500 G typical horizontal component from one side to the other of a sunspot with a typical diameter of 10 000 km. Assuming a 0.3 arcsec pixel, which is 250 km wide, the observed horizontal field variation is about 75 G, which is not very large, yet larger than the 20 G inaccuracy on the transverse field. As a result, the difference between the observed vertical and horizontal gradients cannot be simply assigned to measurement inaccuracies.

Some mathematical theorems
In this section we determine whether the lack of spatial resolution is able to explain the observed discrepancy. First of all we note that if the lack of spatial resolution were able to explain the effect, we would expect the non-zero value observed for the divergence to depend on the pixel size. As far can be seen in the review by Balthasar (2018), this is not the case. The review deals with various instruments, and therefore various pixel sizes, but the results are homogeneous, even though obtained using different active regions. In order to investigate this question in a more quantitative manner, we considered the spectropolarimetric data of NOAA 10953 acquired by the SOT/SP experiment on board the HINODE satellite on 1 May 2007 at 10:46 UT. The pixel size was 0.16 . We then prepared artificially degraded resolution data by averaging the Stokes parameters on 2 2 , 4 2 , and 8 2 pixels. We then submitted these artificially degraded resolution data together with the original data to our inversion code UNNOFIT (Bommier et al. 2007). We did not detect any significant variation in the difference between vertical and horizontal magnetic field gradients as a function of the pixel size.
It can then be objected that the typical variation length for the magnetic field may be very small with respect to these pixel sizes (or typical pixel sizes or current instrument pixel sizes). In this case no variation with the pixel size could be observed because of the scale difference between the pixel size and the magnetic field typical variation length. However, such a case can be treated by the theorems we introduce below.

Theorem 1: filtering and derivation commute
We define the filtered quantity through a filter not necessarily isotropic as a local average. The filtering operation is a convolution product of the quantity with the filter function. Theorem 1 is that this filtering (or convolution) operation commutes with the derivation operation. The result is that the derivative or divergence of the local average is the local average of the derivatives or divergences respectively. We recall below our demonstration, presented in Bommier (2013Bommier ( , 2014.
Due to the linearity of the Zeeman effect, as previously discussed, the line-of-sight integration can be modeled as where ϕ is the contribution function, and where we have assigned a "height of formation" z to the final result. This height of formation is generally close to the maximum of the contribution function, as seen in the examples computed by Bruls et al. (1991). The contribution function acts as a filter in the above equation.
Analogously, the pixel integration in x or y can be modeled with a convolution by a crenel function. If now we compute the divergence of the observed magnetic field, we can apply the derivation of a convolution product as recalled below. Given the convolution product defined as we evaluate F (x + dx) in order to evaluate the derivative. We obtain by changing the variable x = x − dx. We then have the derivative F (x) where the derivative of the convolution product is the convolution (by the same function) of the derivative.
Returning to the case of the line-of-sight integration for the divergence, which is a combination of derivatives: An analogous derivation can be made for the filtering with a crenel function for the pixel integration. It results from the derivation property of the convolution product that any filter would lead to the same result, which is that the divergence of the filtered quantity is the filtering applied to the divergence of the local quantity. This result can be more easily derived in the Fourier space, where convolution products are transformed into simple products. If we denote asĤ x (k) (resp. y, z) the spatial Fourier transform of the magnetic field component H x (r) (resp. y, z), the Fourier transform of the field divergence is We denote as ϕ (r) the 3D spatial filter to be applied to model the observations, and we accordingly denote asφ (k) its Fourier transform. Becausê we have ϕ [div H] = div ϕ [H], which is the above-mentioned result. This does not assume any isotropy of the filter. Different filter sizes and types may be assumed in x, y, z, as they are in the case of the observations, and the result will be maintained.

Derivatives and finite differences
In practice, the divergence of the magnetic field vector in the observations is evaluated by means of finite differences and not mathematical derivatives. In order to compare derivatives and finite differences, we have to discriminate between typical lengths of field variation shorter than, similar to, or longer than the typical pixel size used as the typical length for local averaging.
The field variations shorter than the pixel size may have local derivatives larger or smaller than the derivative of the averaged function, which is also the average of the derivatives over the pixel. The averaged derivative over the pixel corresponds to larger variations with respect to the pixel size. By integrating over the pixel, the average function variations yield typical variation lengths that are longer than the pixel size. Only those remain, but the smaller values are included in the average.
When the typical variation lengths are longer than the pixel size, function derivatives and finite differences at the pixel size are close together. The derivatives computed by finite differences at the pixel size are then very close to the averaged derivatives over the pixel.
The case of typical variation lengths similar to the pixel size is a limit case. In this case, the effect of a difference between derivative and finite difference, would be sensitive to the pixel size. This is not what is observed, as discussed at the beginning of this section.
The typical height difference between the observed lines is on the order of 63 km in the quiet photosphere, as directly observed by Faurobert et al. (2009) for the two usual lines, Fe i 6301.5 and 6302.5 Å, and slightly more in active regions following the simulation by Khomenko & Collados (2007, Fig. 4). We also note that when such a height difference is combined with the observed vertical gradient of 3 G km −1 , this leads to a field decrease of about 250 G in this interval, as observed, when the field itself is about 2000 G. Therefore, the ∆z used to determine the vertical gradient is significantly smaller than the magnetic field strength height scale, which is the height difference for the magnetic field strength being divided by two, which increases the significance of the result.
For the observed horizontal gradient, its order of magnitude of 0.3 G km −1 is fully compatible with the typical sunspot diameter (10 000 km) and the typical horizontal field component in the penumbra (1500 G), which reverses from one side of the sunspot to the opposite side. It is probably not underestimated.
Alternatively, if it is assumed that the vertical gradient is overestimated, and if consequently it were instead ten times smaller (on the order of 0.3 G km −1 also), the height formation difference between the two lines would then be ten times larger, from 700 to 1000 km, which is inconsistent with the photosphere visible thickness.

Theorem 2: the content of the non-zero average
When the averaged quantities are all zero, the average is accordingly zero. By taking the opposite of this statement, we obtain the second theorem: if the averaged quantity is non-zero, there is at least one quantity among the averaged quantities that is nonzero.
Applying this second theorem to our problem, we obtain that if a non-zero divergence of the magnetic field is obtained from the averaged quantities, there is at least one point within the pixel where the divergence is non-zero. A non-zero divergence cannot result from the lack of spatial resolution alone because it is made up of averaging operations. The observed divergence is probably true.
It should be recalled that the Zeeman effect, which is responsible for the effect of the magnetic field on the spectral lines, is an intrinsically linear effect because the sublevel energy variation is linear as a function of the magnetic field strength.

Magnetic induction, magnetic field, and magnetization
We recall that the three quantities of magnetic induction, magnetic field, and magnetization are related by the Maxwell equation in magnetized media where B is the magnetic induction, H is the magnetic field, and M the magnetization. For the definition of these quantities, we refer to Jackson (1975), in particular Sect. 6.7, "Derivation of the Equations of Macroscopic Electromagnetism", where Eq. (9) is derived by averaging from the microscopic scale, with an averaging length that is large with respect to the microscopic scale, but small with respect to the macroscopic scale. Equation (9) is then macroscopic. At the microscopic scale, the electric and magnetic fields are denoted e and b, respectively, in lowercase characters. They obey the Maxwell equations in a vacuum, where the charge and current densities include all the charges present in the medium. By separating the effects of free and bound charges, the four macroscopic quantities D, H, P, and M are found at the macroscopic scale, where the polarization P and the magnetization M are due to the contribution of the bound charges. Jackson (1975) tells us (p. 233) that when the medium is a plasma that also contains free charges, and if these free charges also possess intrinsic magnetic moments, these magnetic moments can be simply included in the definition of M. Delcroix & Bers (1994, vol. 1, p. 89) evaluate this intrinsic magnetic moment of the plasma free charges in their spiral motion about the magnetic induction due to the Lorentz force qu× B. They show for the ensemble of orbital magnetic moments of the individual particles where β is the plasma usual parameter, but computed with the free charge density. This magnetic moment has the opposite sign with respect to the magnetic induction. The effect is therefore called plasma diamagnetism. The photospheric VALC model (Vernazza et al. 1981) was derived from spectroscopic observations, statistical equilibrium, radiative transfer, and the Saha law for deriving the electron density within the hypothesis of medium electric neutrality. At the formation height of the Fe i lines used for the magnetic field diagnostic (about 250 km) the temperature is 4780 K, the neutral hydrogen density is 2.3 × 10 16 , and the electron density is 2.7 × 10 12 cm −3 . The plasma β expressed in terms of the electron density is β = 4.4 × 10 −5 (for B = 1000 G) and M is negligible with respect to H and B/µ 0 . The main second source of magnetization is the neutral hydrogen paramagnetism. Its susceptibility is χ = 3 × 10 −5 at 250 km (correcting an error present in Bommier 2013Bommier , 2014, with M = χH. The neutral iron paramagnetism susceptibility is weaker χ = 2×10 −8 . We have in this case B µ 0 H, but if there were free electrons coming "from below" and accumulating in the photosphere as we propose, and if their density were finally comparable to the neutral hydrogen density, we would have β = 0.38 (for B = 1000 G). M would then be comparable to H and B/µ 0 , and they would all be different.
In the Maltby-M sunspot umbra model (Maltby et al. 1986), the formation height of the Fe i lines is about 100 km. The β expressed in terms of the neutral hydrogen density is β = 1.27 at this height and for B = 1000 G, which is β = 0.20 for the magnetic field B = 2500 G typical of sunspot umbra. We assume for a moment an oversimplified model of sunspot umbra with a purely vertical field or induction along z. The Mawxell law div B = 0 then results in ∂B z /∂z = 0. We then have where β is computed with the electron density. If we assume that this density is equal to the neutral hydrogen density, following our proposal, we obtain B z ∂β/2∂z ≈ −3.4 G km −1 for B = 1000 G (and −1.4 G km −1 for B = 2500 G) at the formation height of the Fe i lines used for the measurements and for the neutral hydrogen density of the Maltby-M sunspot umbra model, which is comparable to the observed value ∂µ 0 H z /∂z ≈ −3 G km −1 . This value is −1.2 G km −1 for the quiet Sun VALC model and for B = 1000 G. The case of laboratory plasmas is very different. From Eq. (10), when M is comparable to H and B/µ 0 , the plasma β is on the order of unity (we temporarily assume here that plasma is made of free charged particles). When β is on the order of unity, the gas pressure is comparable to the magnetic pressure. Therefore, particles can escape under the effect of the gas pressure, and the magnetic field will be insufficient to keep then. M is an obstacle to plasma confinement (Delcroix & Bers 1994). Consequently, in usually observed laboratory plasmas, M probably remains weak with respect to H and B/µ 0 . Therefore, we have B ≈ µ 0 H; B and µ 0 H are equivalent in these plasmas. In this paper we consider a plasma where M, H, and B/µ 0 are all comparable, and we study the effect of M on this plasma. This plasma is the solar photosphere and is naturally formed and maintained, but it is different from the laboratory confined plasmas.

Magnetic field produced by the conduction currents
The contribution H c to the magnetic field due to the (macroscopic) conduction currents J is such that H c can be computed by applying the Biot & Savart law to J. As a consequence, H c is a divergence-free field Generally, however, we have in magnetized media where H d is a curl-free field as we introduce in the following section.

Magnetic masses
There is another contribution to the magnetic field that is due to the magnetization resulting from the plasma diamagnetism. This magnetization M is a magnetic moment density, which is continuously spread outside the atom. In the solar photosphere and following our proposal of electron accumulation in the solar photosphere, we have because of the vertical density gradient. We can then introduce where the curl-free field H d has the mathematical form of an electric or gravitation field due to the magnetic mass density ρ M , which may be positive or negative. The magnetic mass density is defined by the leftmost member of the above Eq. (17). We recall that the magnetic mass is a mathematical tool only.
In the case of a finite magnetized volume, surface magnetic masses are formed on the volume surface. Denoting by n the unit vector perpendicular to the surface and oriented outward, the surface magnetic mass density is We consider again the oversimplified model of a purely vertical field or induction along z. We assume that B z is positive and constant. The magnetization resulting from the plasma diamagnetism in Eq. (10) is then also vertical, but with M z negative and with modulus decreasing with height. As a result ρ M > 0. However, if ρ M is spatially constant, the resulting field H d is zero, due to reasons of symmetry and Eq. (17) cannot be satisfied. But this is not the case: the β gradient itself decreases with height in such a way that ρ M > 0 decreases with height. As a result H d is oriented upward, and decreases with height. The total field H of Eq. (14) then decreases with height and has a non-zero divergence as observed, provided that it is proven that what is measured is the magnetic field H and not the magnetic induction B, which is the object of the next section.
The H d field is called the demagnetising field because it is opposite to M.

Contribution of the magnetization current
A magnetization current appears when which is the case when there are density variations such that the small loops (made by the charged particles about the magnetic induction and responsible for the magnetization) do not counterbalance. We then have so that B can be computed by applying the Biot & Savart law to J and J M .
In the case of a finite magnetized volume, surface magnetization currents are formed on the volume surface. Denoting as n the unit vector perpendicular to the surface and oriented outward, the surface magnetization current density is The effect of M on H or B can be evaluated either via the magnetic masses and their demagnetising field, or via the Biot & Savart law applied to the magnetization current. The two approaches lead to the same result, as can be easily verified in the case of a sphere with constant magnetization (also considering the surface magnetic masses and magnetization current).

The Zeeman Hamiltonian
The Zeeman Hamiltonian is the interaction energy between the atom having an elementary magnetic momentum m and the magnetic field.
In atomic physics, this Hamiltonian is obtained in particular by developing the atom impulse in the presence of a magnetic potential (P − q A) 2 , which leads to the well-known Zeeman Hamiltonian −m · B, but the atom here considered lies in a vacuum. In the present paper we study the case of an atom embedded in a magnetized material, and we study in particular the effect of the surrounding material on the atom.
In the following, we present four demonstrations all agreeing that the Zeeman Hamiltonian for the atom embedded in the magnetized material is −µ 0 m·H, where the magnetic field H includes the demagnetising field H d created by the surrounding magnetization, whose effect on the atom is thus taken into account. This is a result currently used in magnetized materials physics (see, e.g., Gignoux & Schmitt 1993;Garnier et al. 1998;Zhang et al. 1994).

Demonstrations 1 and 2: evacuating an elementary matter cylinder to determine its magnetic energy
The aim of this demonstration is to determine the energy required to evacuate to infinity or, alternatively, to fill from infinity a cylinder of matter containing a single atom of magnetic momentum m. At the microscopic scale, the medium is made of particles (free electrons, ions, and neutral atoms, themselves made of the nucleus and electrons), and the Maxwell law Eq. (9) is the result of a macroscopic averaging as described in Jackson (1975), Sect. 6.7, "Derivation of the Equations of Macroscopic Electromagnetism". This averaging is performed on lengths that are large with respect to the microscopic scales. As a result the magnetization is represented by the quantity M, which is a density of magnetic moments (i.e., a spatially smoothed quantity in cm −3 ; see Eq. (26) below), whereas the medium is made of particles at the microscopic scale.
In order to determine at what field the atom is submitted in the medium, it is then necessary to go back from the macroscopic smoothed scale to the microscopic scale where the atom is in a vacuum, with the other atoms or ions or free electrons far from it. This atom is however submitted to the field created by the other particles even far from it. Considering that in sunspot umbra the charged particles move along the magnetic field and drag the neutral particles along with them in such a way that there is no displacement perpendicular to the field in a first approximation, the vacuum created around the atom takes the form of a very long cylinder along the magnetic field, as represented in Fig. 2. In this figure, we have also represented the magnetic field H c due to the conduction currents (see previous section), and the demagnetising magnetic field H d due to the non-zero divM, which is the way the magnetization affects the atom.
The magnetic field variations in induction units increases up to 3 G km −1 following the observations, which corresponds to a variation of 3 × 10 −10 G along a scale of 1000 Å. This variation can be neglected at the cylinder diameter scale, so that B, H, and M can be considered constant close to the atom and to the cylinder center.
The cylinder dimensions are assumed to be large with respect to the atom dimensions, but small with respect to the macroscopic characteristic lengths, and such that the cylinder contains only one atom to be evacuated to infinity in order to determine its magnetic energy. The cylinder length is assumed to be very large with respect to the cylinder width. Magnetic masses appear at the cylinder extremities where the separation surface is not parallel to M. However, the second demagnetising field created by these masses is negligible at the cylinder center because the cylinder is very long. This may be the second reason for considering the cylinder form to describe the vacuum volume made around the atom. For any other form, for instance a sphere, such magnetic masses would also be present and would induce a non-negligible second demagnetising field, which would also be present outside the volume and would add a perturbation to the local field, induction, and magnetization. The long cylinder form (or at least an elongated form parallel to the field) is the only possible form considered that does not add any perturbation to the local fields. A third argument in favor of the cylinder form is the symmetry of the medium, which is cylindrical along the field, induction, and magnetization, therefore not spherical.

Demonstration 1: via magnetic masses
Outside the cylinder, the magnetic induction is B, the magnetization is M and the magnetic field is H = H c + H d . Inside the cylinder, close to its center, the tangential component of H and the normal component of B are transmitted through the separation surface. The cylinder axis is parallel to B and to H, so that H is transmitted inside the cylinder. There is no matter inside the cylinder, which is assumed to be empty for single atom energy study purposes; therefore, no magnetization is present. Therefore, the magnetic induction inside the cylinder is which is different from the induction B outside the cylinder, where there is some magnetization and where B = µ 0 (H + M).
We assume constant H at the atom scale, so at this scale, which is the one of our present problem, divB i = 0 is satisfied. The cylinder is parallel to B and to H. As a consequence, the free electrons, which move along the magnetic field, do not enter the cylinder, which contains a single atom. We note, however, that the electron Larmor radius in typical photospheric fields is about 25 microns, which is large with respect to the atomic or microscopic typical dimensions. However, the magnetization due to these electrons was averaged as M, which affects the atom via H d , which is transmitted inside the cylinder.

Demonstration 2: via magnetization currents
As described in Sect. 1.2.3 of Chap. 2 of du Trémolet de Lacheisserie et al. (2002), the magnetic induction B i inside the empty cylinder can also be evaluated by considering the magnetic induction ∆B created by the matter, which is evacuated in order to create the vacuum in the cylinder. The aim of this operation is to determine the matter cylinder energy in the field by bringing back the matter cylinder from infinity to the empty cylinder. This matter cylinder is made of constant magnetization M, and the magnetic induction inside a cylinder of constant magnetization M is ∆B = µ 0 M (see, e.g., du Trémolet de Lacheisserie et al. 2002, Sect. 1.1.6). This magnetic induction can be computed from the magnetization currents at the surface of the cylinder. Therefore inside the vacuum cylinder, the induction is B i = B − ∆B, and Eq. (22) is recovered.

Agreement of the two demonstrations
These two determinations of B i , one based on the magnetic masses evaluation, and the second one based on the magnetization currents evaluation, fully agree, as expected.
The energy of the atomic magnetic momentum denoted as m is then

Demonstration 3: from the magnetic potential energy
The total magnetic energy of a medium with magnetic field H and induction B, which is either paramagnetic or diamagnetic, as ours is, can be written (Jackson 1975, Sect. 6.2, "Energy in the Magnetic Field") By applying Eq. (9), this integral can be expanded as A40, page 6 of 11 where the first member on the right-hand side is the field energy and the second is the magnetization, or magnetized matter, energy in the field. The magnetization M is the spatial average of local elementary magnetic moments (Jackson 1975, Eq. (6.98)) where δ is the Dirac function and where is the spatial average. By considering an elementary volume about the atom of magnetic moment m, with only the atom inside the volume, the atom energy is then, in terms of the total energy, However, the system is not totally isolated because there are external current sources. When a small matter piece is introduced, or a small displacement is performed, energy is transferred from the sources to the system. In his Sect. 6.2, Jackson (1975) writes that "we can show that for a small displacement the work done against the induced electromotive forces is twice as large as, and of the opposite sign to, the potential-energy change of the body". The energy of the atom embedded in the magnetized matter is therefore where the effect of the surrounding dipoles on the atom is taken into account in the demagnetising field contribution H d included in H.

Concluding remarks about the Zeeman Hamiltonian
All the methods lead to that conclusion that the energy of the atom embedded in the magnetized matter is W = −µ 0 m · H. Accordingly, the Zeeman Hamiltonian H M describing the interaction between the atom (paramagnetic of momentum m) and the magnetic field is When LS -coupling is valid for describing the atomic states, this may be simply rewritten as where J is the atomic total kinetic momentum J = L + S, g J the level Landé factor and µ B the Bohr magneton µ B = q /2m e , where q and m e are the electron charge and mass, respectively. As a result, the magnetic field H, which includes the demagnetising field H d contribution, governs the Zeeman Hamiltonian describing the interaction between the atom embedded in magnetized matter, and the magnetic field. This can be read, including the contribution of the demagnetising field, in Gignoux & Schmitt (1993), Garnier et al. (1998), andZhang et al. (1994). The result is that what is measured by interpretation of the Zeeman effect in spectral lines in a magnetized medium like the solar photosphere, is H, whose divergence may be non-zero at the macroscopic scale, as observed.
Otherwise, it can be argued that −m · B is inappropriate to express the energy of the magnetic momentum m embedded in the magnetized matter and the magnetic field. When expanded following Eq. (9), products of the type −m · M appear, which express, following Eq. (26), that the matter acts on itself, even if the different values of m are associated with different particles of it. As remarked by Landau & Lifshitz (1973, Sect. 35), when examining the case of a conductor carrying a current and placed in a magnetic field, "the field of the conductor itself cannot, by the law of conservation of momentum, contribute to the total force acting on the conductor". This does not mean that the effect of the magnetic moments surrounding the atom is not taken into account; as previously stated, their effect is at the origin of the demagnetising field H d included in the magnetic field H, which enters the energy instead of B.

Demonstration 4: considering the mean atom submitted to the macroscopic mean magnetic field
From the atomic point of view, the Zeeman Hamiltonian depends on the magnetic field external to the atom, which is the magnetic field at the atom position, with the exclusion of the possible field generated by the atom itself (the self-energy term must be excluded). Let us consider the mean atom submitted to the macroscopic mean field (see the averaging process described at the beginning of Sect. 4). As we consider the mean atom, M is made of the repeated atomic magnetic dipole, whose contribution to the Zeeman Hamiltonian should then be ignored as an internal atomic contribution. As M is part of the magnetic induction B (see Eq. (9)), this indicates that the Zeeman Hamiltonian is not governed by the magnetic induction B, but by the magnetic field H. The magnetic field generated by the surrounding dipoles is the demagnetising field H d included in H; the effect of the surrounding dipoles is therefore fully taken into account. Langevin (1905) demonstrated that the magnetic energy of a moment M placed in an external magnetic field H is W = −µ 0 M · H. Although he determined divH = 0 for the external field in the exterior where B = µ 0 H, the analysis of his experiment is unambiguous about the actual roles of M, B, and H. In the case of non-negligible M, M is contained in a magnetized barrel introduced in a coil such that it is H and not B that is transmitted from the exterior to the interior of the barrel, because M, B, and H are all parallel to the barrel surface inside the coil, similarly to Fig. 2. Weiss (1907) showed that the effect of the surrounding dipoles is the demagnetising field, which contributes to H and is then fully taken into account in the magnetic energy.

Electron accumulation at the surface of the Sun
Our proposal is based on the consideration that electrons, which are much lighter than protons, have a thermal velocity inside the Sun much higher than the escape velocity, when protons remain submitted to the Sun's gravity. As a result, internal electrons tend to escape. The cloud of free electrons spreads out toward the star exterior. However, as we explain below, the escaping electrons come to be finally partly retained by the protons, the protons themselves are retained by the star gravity. The combined result is an accumulation of electrons in the star's more external layers, in particular the surface layer. It should be noted that the physical conditions in terms of temperature and density are specific in the surface layer (see Fig. 3). We present and discuss this mechanism in the following.

Escape velocity from gravity
The fact that the electron thermal velocity surpasses the escape velocity is well known in the solar corona. As Meyer-Vernet (2007) writes on p. 251, "At a temperature of 10 6 K, their thermal A40, page 7 of 11 A&A 634, A40 (2020) Fig. 3. Temperature (top) and electron density (bottom) inside the Sun, from Allen (1973) and assuming local electric neutrality. These figures show the different behaviors between the interior and the surface layer.
speed 5.5 × 10 6 m s −1 is nearly 10 times greater than the escape speed". Instead, protons do not escape because "electrons and protons have opposite charges, but their masses differ by the factor m p /m e 1837, so that electrons have a thermal speed greater than protons by a factor of order of magnitude m p /m e 43 (because their temperatures have generally the same order of magnitude)".
Inside the Sun, at 0.5 R , the temperature is analogous, being 3.4 × 10 6 K, and 94% of the solar mass is contained inside the sphere of 0.5 R radius (Allen 1973). If we recall that the escape velocity from a star is given by where G is the gravitation constant, M < is the mass inside the sphere of radius R, where R is the distance from the star's center where the escape velocity is evaluated. In these conditions, the electron thermal velocity is found to be 14 times higher than the escape velocity at 0.5 R , when the proton velocity is only 0.34 times the escape velocity.

Thermal escape
Thermal escape is a well-known process in planet atmosphere studies (see, e.g., Chamberlain 1963, for a review). It is a very efficient mechanism for evaporating light elements like hydrogen from the Earth atmosphere, but the conditions in our case are very different for two aspects. First, in these atmospheres the escape velocity remains significantly higher than the quadratic average thermal velocity defined as where m is the mass of the particle under study (electron, proton), k B is the Boltzmann constant, and T the temperature. Then, in the case of planet atmospheres, only the tail of the velocity Maxwell distribution is affected by escaping. Second, escaping takes place in very low density regions, where the collision mean free path is larger than a typical length, such as the pressure scale height, which means that this atmosphere layer can be considered collision free. In the lower lying neighbor layer, however, the few collisions reestablish the thermal equilibrium and repopulate the escaping velocities. This is the source of the efficiency of the mechanism in planet atmospheres. Our case of the solar interior is different for these two points. First, the escape velocity is much lower than the electron quadratic average thermal velocity, by a factor of 14 at 0.5 R . As a consequence, nearly all the electron velocity Maxwell distribution is affected by escaping. Second, however, the mean free path for the electron-proton collision is about 7 × 10 −9 m at 0.5 R , when the density scale height is on the order of 55 000 km in the solar interior, according to the Allen (1973) data. The electronelectron collision mean free path is about twice the electronproton collision mean free path, as explained by Beck & Meyer-Vernet (2008) in their note 5: "For collisions between electrons the calculation must be done in the center-of-mass frame, and the mean free path for change in the velocity direction is greater by a factor of about 2". As the proton and electron densities are similar in the solar interior, the proton-proton collision mean free path is similar to that of the electron-electron collision. In terms of collisions, the conditions in the solar interior are completely different from those of a planet atmosphere.

Escape velocity from one proton
However, even in the presence of collisions, the electron velocity values are nearly all higher than the escape velocity as it results from gravity, as defined in Eq. (31). We also need to consider the attracting effect exerted by protons on electrons, when the proton velocity remains generally lower than the escape velocity and the protons submitted to the star gravity. We propose the following model to analyze this effect.
We first assume local electric neutrality and the same number of mixed electrons and protons. Then we assume that their electric fields roughly balance two by two in an electron-proton pair, even if they are not linked to each other. Thus, an escaping electron is submitted to the electric field of a single proton. The calculation of this electric field requires knowing the distance between the electron and the proton, which can be taken as the mean distance between particles, which is the −1/3 power of their density, which is roughly the same for electrons and protons. This electric field derives from a potential, as does the gravitation field. Thus, there is also an escape velocity, as there is for gravitation. The escape velocity for freeing the electron from the corresponding single proton is then v esc = 2q 2 4πε 0 m e n −1/3 , where q is the elementary charge, m e is the electron mass, and n the density, which is nearly the same for electrons and protons.
In the solar interior conditions given by Allen (1973), we then obtain that this escape velocity is nevertheless six times lower than the quadratic average thermal velocity for electrons, which then also escape from the protons, at least initially. As the escape velocity is obtained by equating the particle kinetic and potential energies, the total escape velocity from gravity and proton for an electron is v esc = 2GM < R + 2q 2 4πε 0 m e n −1/3 , and we obtain that for electrons it is 5.5 times lower than the quadratic average thermal velocity at 0.5 R . The electrons then initially escape from gravity and protons.

Escape velocity from several protons
As the electron cloud is free at the beginning, in a first approximation it begins to expand. The result is a decrease in the electron density, while the proton density remains unchanged. As a consequence, the number of electrons and protons is not yet the same in a given volume. As a result, an escaping electron is now submitted to the attraction of more than one mean proton, depending on the density ratio. If we denote as ρ this density ratio where n e and n p are the electron and proton density, respectively, the total escape velocity of the electron submitted to ρ > 1 protons and to gravity becomes because the mean distance between the electron and the other particles is n −1/3 e . If we denote as x the ratio between the electron escape and thermal velocities we obtain from the preceding equations which permits us to compute the proton-to-electron density ratio ρ from the desired velocity ratio x. The electron escaping will certainly be suppressed when x = 3, which means that the escape velocity will be at least three standard deviations higher than the thermal velocity for each velocity component. In this case, ρ is found to be 6 × 10 3 at 0.5 R in the Allen (1973) conditions. This corresponds to a mean distance between electrons 18 times larger than that between protons. These quantities are found to be quite constant along the solar radius.

Spreading time
One important question is where exactly do the initially free electrons escape, and the answer is to the distance where the electron cloud density is divided by the ratio ρ defined above.
We propose that the density scale height is the typical dimension of the cloud. We denote as h this density scale height. Then the cloud extends at most up to 18/2 = 9 times h from its initial position. The division by two comes from the fact that the electron cloud may extend in two directions from its initial position. The density scale height h is on the order of 55 000 km inside the Sun, but sharply decreases at the surface, up to 110 km (pressure scale height) in the photosphere. The time for the initially free electron cloud to reach these stopping distances can be evaluated from the Fick's laws of diffusion. If we denote as D the diffusion coefficient, which is where the quadratic average particle velocity is the thermal velocity v th and is the shortest electron collision mean free path discussed above, the time t necessary to reach the distance d is The diffusion velocity v diff can then be defined as v diff = d/t and it is found that this diffusion velocity is very weak. The cloud slowly spreads out. The time needed to reach the 9h distance is found to be greater than the age of the Universe for electrons of the more internal layers where h is about 55 000 km. Thus, these electrons do not leave the star. Only the electrons lying in the layers close to the surface can reach their limit distance, and this distance 9h is also found close to the surface because the density scale height is much smaller there. As a result, the electrons do not leave the star, and there is an accumulation of them in the surface layer. The surface electrons reach their limit distance 9h at a velocity of 1.3 × 10 −6 m s −1 , which is 1.7 m year −1 . Thus, the spreading is quasi-static.

Surface electron density
We finally propose an attempt to evaluate the surface electron density. If we consider that the density scale height is roughly the length along which the density is divided by three, and if we also neglect for a moment the height variation of the density scale height h and assume a constant h, the number of electrons per volume unit will be obtained by summing the number of electrons in each of the nine internal layers, each expanded up to 9h, and each being three times more populated than its more Fig. 4. Linear polarization degree with polarization direction parallel to the solar limb, observed 25 arcsec inside the limb in the continuum and far from active regions, as a function of the wavelength. Observations: squares Leroy (1972Leroy ( , 1974, triangles Wiehr (1975), circles • Mickey & Orrall (1974), crosses × Goutail (1978). Line: model result based on the HSRA atmosphere model (Dumont & Pecker 1971). From Leroy et al. (1977), used with permission. Red line: starting from 8000 Å, the polarization degree ceases to decrease following Rayleigh scattering on neutral hydrogen. It remains constant, as in Thomson scattering on electrons; HSRA modeling would give Thomson scattering dominant only after 10 000 Å.
However, on the surface, a part of these electrons associate with protons to form neutral hydrogen atoms. The initial proton density was ρ times the electron density of the first expanded layer, with ρ about 6 × 10 3 at 0.5 R , but similarly 5 × 10 3 in the more external layer in the Allen (1973) conditions. Thus, the remaining free electron density is comparable to the neutral hydrogen density. In Sect. 4, we show that this is the order of magnitude where the resulting magnetization divergence can explain the observed magnetic field divergence.

Electron density measurements in the solar photosphere
To our knowledge, the electron density measurements in the photosphere are only indirect. Results are derived from spectroscopic analysis of the observed spectra (e.g., the models by Vernazza et al. 1981). The hypothesis and equation of electric neutrality are the basis of the analysis. The electron density then results, within this hypothesis, from the ionization equilibria of the different species present in the atmosphere. The most direct measurement we know is the investigation of the continuum linear polarization observed close to the solar limb and far from active regions, as done by Jean-Louis Leroy at the Pic-du-Midi Observatory in the 1970s (Leroy 1972(Leroy , 1974. In Leroy et al. (1977, Fig. 3, p. 167), he presented the measurement and modeling synthesis, which we reproduce here in Fig. 4 with his permission. New measurements were made later on by Wiehr (1978; who improved the results of Wiehr 1975 included in Fig. 4) and Wiehr & Bianda (2003; in excellent agreement with Wiehr 1978), but in such good agreement with the measurements and with the theoretical model presented in Fig. 4 that we did not find it necessary to redo the figure. These more recent measurements fall in the visible range and not in the infrared range, which is the one covered in our discussion.
The continuum limb polarization results from Rayleigh scattering by the neutral hydrogen atoms. Rayleigh scattering is wavelength-dependent, which explains the decrease in the polarization degree as a function of increasing wavelength. The line in the figure is the result of a model by Debarbat et al. (1970) and improved by Dumont & Pecker (1971), based on the HSRA model atmosphere. More recently, for the purpose of stellar studies, Kostogryz & Berdyugina (2015) redid the calculations but with a series of model atmospheres including those of Vernazza et al. (1981) and the HSRA model. Figure 5 of Kostogryz & Berdyugina (2015) shows that all models lead to very close results for the theoretical limb polarization.
However, the observed polarization ceases to decrease with increasing wavelength at about 8000 Å. This suggests that Thomson scattering on free electrons then becomes dominant. Thomson scattering is wavelength independent. This is indicated by the horizontal red line in Fig. 4. However, the theoretical results based on the HSRA model would obtain a dominant Thomson scattering only after 10 000 Å. Thus, observations seem to reveal an electron overdensity, with respect to the HSRA model, which is based on electric neutrality. Even if the orders of magnitude do not seem to exactly correspond to our hypothesis of comparable free electron and neutral hydrogen atom densities, this is a confirmation (to be investigated by new observations and modeling) of our hypothesis of increased free electron density that can explain the difference between vertical and horizontal magnetic field gradients in the photosphere in and around sunspots via the B = µ 0 (H + M) law.

Conclusion
We have proposed an explanation of the non-zero divergence of the observed magnetic field (see the observation review by Balthasar 2018) by investigating the law B = µ 0 (H + M) in magnetized media. In order to obtain a non-negligible magnetization M at the surface of the Sun where the measurements are performed, we have invoked electron thermal escape from the more internal layers, where the electron thermal velocity noticeably exceeds their escape velocity. We have paid attention to the retaining role of protons and show that it does not totally prevent electrons from escaping. The result is an electron accumulation in the surface layer.
Positive free charges could contribute in an additive manner to the magnetization. The positive charges gyrate about the induction or field in the opposite sense with respect to the negative charges (see Fig. 2.8 of Meyer-Vernet 2007), but as the charges are also opposite themselves, the elementary corresponding loop currents are similar for the two charge types, leading to their additive contribution to the magnetization. However, the usual photosphere models like VALC and HRSA are built on the hypothesis of electric neutrality. Then, if there is electric neutrality in the layer, the charge densities and then M remain low, following the spectroscopic analysis results. We cannot keep both non-negligible M and electric neutrality in the layer.
The escaping electrons would accumulate at the Sun's surface, when the protons would remain lower. Electric fields would then appear inside the star, but the electric effects would remain inside. As the global charge would remain zero or very weak, no electric effects would follow outside the star, in a first approximation.
The final result is that what is measured by Zeeman effect analysis is the magnetic field H and not the magnetic induction B. However, the magnetohydrodynamical modeling requires the knowledge of B. It can be deduced from the measurements H by applying the law B = µ 0 (H + M), provided M is also known. M would have to be reconstructed from its divergence div M, which is opposite to that of the magnetic field H, because div B = 0, which implies div M = −div H, which is measured. The conduction currents J, however, can be directly derived from the measurements H because curl H = J, as detailed in Eqs. (12)-(15).