Free Access
Issue
A&A
Volume 550, February 2013
Article Number A98
Number of page(s) 12
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201118239
Published online 01 February 2013

© ESO, 2013

1. Introduction

In recent years several attempts have been made to investigate the magnetic field vector distribution in the solar internetwork. Initially, these works studied the magnetic field strength is in these regions. Some favored magnetic fields of about a few hundred Gauss or less (Asensio Ramos et al. 2007; López Ariste et al. 2007; Orozco Suárez et al. 2007a, Orozco Suárez & Bellot Rubio 2012) while others found magnetic fields in the kilo-Gauss range (Domínguez Cerdeña et al. 2003, 2006; Sánchez Almeida 2005). These studies were carried out mostly with low spatial resolution data (1′′). Whenever the spatial resolution increased to better than 1 arcsec, this decreased the signal-to-noise ratio (S/N). With the Hinode satellite (Kosugi et al. 2007) it is now possible to obtain spectropolarimetric data (full Stokes vector) with high spatial resolution (0.3′′) and low noise (σ ≈ 10-3 in units of the continuum intensity). Thanks to these new data, it is now also possible to investigate not only the module but the three components of the magnetic field vector. This has led to a new controversy about the angular distribution of the magnetic field vector in the quiet Sun. While some authors (Orozco Suárez et al. 2007a,b; Lites et al. 2007, 2008) found that the magnetic field is mostly horizontal (γ ≈ 90°; with γ being the inclination of the magnetic field vector with respect to the observer’s line-of-sight), others favor a quasi-isotropic distribution of magnetic fields (Martínez González et al. 2008; Asensio Ramos 2009; Stenflo 2010). With a few exceptions (Harvey et al. 2007; Lites et al. 2008; Martínez González et al. 2008), all previous studies were carried out employing data recorded at disk center only. Therefore, to better constrain the angular distribution of the magnetic field vector in the internetwork, we considered spectropolarimetric data recorded at different positions on the solar disk (Sect. 2).

In addition, Asensio Ramos (2009), Stenflo (2010), and Borrero & Kobel (2011; hereafter Paper I) warned that the highly inclined magnetic fields obtained by some studies could be caused by the noise in the linear polarization profiles. This yields a distribution of B (component of the magnetic field vector that is perpendicular to the observer’s line-of-sight) with a peak at around 50–90 Gauss. To avoid this problem, these authors proposed to include only those profiles in the analysis that have a S/N > 4.5 in the linear polarization (Stokes Q and U). Although this selection criterion allows one to retrieve reliable distributions for the magnetic field vector, it has the disadvantage of excluding most of the Stokes profiles within the field-of-view from the analysis (see Borrero & Kobel 2012; hereafter Paper II; cf. Bellot Rubio & Orozco Suárez 2012). In this paper we adopt an alternative approach based on inverting the histograms of the observed Stokes vector (Sect. 4) over the entire field-of-view instead of inverting the Stokes vector at each pixel over the observed region. Under a number of simplifying assumptions, whose limitations are described in Sect. 5, we were able to reach some important, albeit preliminary, conclusions about the angular distribution of the magnetic field vector in the solar internetwork and its variation across the solar disk (Sect. 6).

2. Observations and datasets

The data employed in this work were recorded with the spectropolarimeter (SP; Ichimoto et al. 2008) attached to the Solar Optical Telescope (SOT, Tsuneta et al. 2008; Suematsu et al. 2008; Shimuzu et al. 2008) onboard the Japanese spacecraft Hinode (Kosugi et al. 2007). The spectropolarimetric data comprise the full Stokes vector (I,Q,U,V) around the pair of magnetically sensitive spectral lines Fe I 6301.5 Å (geff = 1.67) and Fe I 6302.5 Å (geff = 2.5). geff refers to the effective Landé factor calculated under LS coupling. The spectral resolution of these observations is about 21.5 mÅ per pixel, with 112 pixels in the spectral direction. The spatial resolution of Hinode/SP observations is 0.32′′. For this paper we selected three maps at three different heliocentric positions. In all three maps the spectrograph’s slit was kept at the same location on the solar surface for the whole duration of the scan. This means that, while the vertical direction (Y-axis or direction along the slit) contains information about different spatial structures on the solar surface, the horizontal direction (X-axis or direction perpendicular to the spectrograph’s slit) samples the same position at different times. Each spectrum was recorded with a 9.6 s exposure, yielding a noise of about σ = 7.5 × 10-4 in units of the quiet-Sun continuum intensity. Each map records data for a period of time ( > 1 h) that includes several turnovers of the granulation, thus breaking down the temporal coherence and providing spatial information (in a statistical sense) along the X-axis.

In Paper I we have demonstrated that photon noise plays an important role in determining the magnetic field vector from spectropolarimetric observations. To further decrease the level of noise in our observations we averaged every seven slit positions (temporal average of about 67.1 s), which yields a new noise level of about σ = 3 × 10-4 (in units of the quiet-Sun continuum intensity). However, averaging means that the original map is shortened by a factor of seven in the direction that is perpendicular to the slit (X-axis). This decreases the number of points available for statistics. Fortunately, Hinode/SP data have a sufficient number of pixels to ensure good statistics even after averaging (see Sect. 6). In the following we briefly describe each map individually.

2.1. Map A

This map was recorded on February 27, 2007 between 00:20 UT and 02:20 UT. It originally consists of 727 slits positions, of which 103 remain after temporal averaging. The center of slit was located at approximately the following coordinates on the solar surface: X =  −31.7′′ and Y = 7.7′′. This corresponds to a heliocentric position of μ = cosΘ ≈ 1 (Θ is the heliocentric angle) and to a latitude of Λ ≈ 0°. The noise level is σ = 2.8 × 10-4. This map (original and temporally averaged) corresponds to Maps B and C in Paper I, and it was also employed (with and without temporal averaging) by Lites et al. (2008) and Orozco Suárez et al. (2007a).

2.2. Map B

This map was recorded on February 6, 2007 between 11:33 UT and 15:51 UT. It originally consists of 1545 slits positions, of which 222 remain after temporal averaging. The center of slit was located at approximately the following coordinates on the solar surface: X = 493.6′′ and Y = 491.3′′. This corresponds to a heliocentric position of μ = cosΘ ≈ 0.7 and to a latitude of Λ ≈ 30°. The noise in this map is very similar to that in Map B: σ = 3.1 × 10-4.

thumbnail Fig. 1

From top to bottom: inferred magnetic flux density Φ = αB from the inversion of Map A (Sect. 2.1), Map B (Sect. 2.2), and Map C (Sect. 2.3). White areas correspond to regions where all three polarization profiles (Stokes Q, U, and V) are below the 4.5σ-level.

2.3. Map C

This map was recorded on January 17, 2007 between 07:05 UT and 09:58 UT. It originally consists of 1048 slits positions, of which 149 remain after temporally averaging. The center of slit was located at approximately the following coordinates on the solar surface: X =  −3.0′′ and Y = 697.1′′. This corresponds to a heliocentric position of μ = cosΘ ≈ 0.7 and to a latitude of Λ ≈ 40°. Here, the noise level is slightly higher than in Map A: σ = 3.2 × 10-4. We note that some consecutive slit positions in this map show very high noise in the Stokes profiles. Although we could not relate this effect to the South Atlantic Anomaly (increased flux of cosmic rays at certain orbits of the satellite) we have removed these slit positions from our analysis, which reduced the effective number of slit positions to 120.

Figure 1 displays the magnetic flux density Φ of maps A, B, and C as obtained through the inversion of the full Stokes vector with the VFISV (Very Fast Inversion of the Stokes Vector) inversion code (Borrero et al. 2010). For better visualization the maps in these figures are obtained from the inversion of the original data (i.e. not temporally averaged). This avoids pixelization in the X-axis of these plots. However, for the remainder of the paper, our discussions and figures are based only on the temporally averaged (67.1 s) data.

Although the previous figure only shows the total magnetic flux density, it is worth mentioning that the VFISV code also retrieves the three components of magnetic field vector B: B is the module of B, γ is the inclination of B with respect to the observer’s line-of-sight, and φ is the azimuth of B in the plane that is perpendicular to the observer’s line-of-sight. In addition, VFISV retrieves the magnetic filling factor α as well as the line-of-sight component of the velocity vector Vlos and a set of thermodynamic parameters T. We note that the magnetic flux density is defined as Φ = αB. For a more detailed overview on Milne-Eddington inversion codes, which include not only the magnetic field vector but also the thermodynamic and kinematic parameters relevant to the line formation, we refer the reader to del Toro Iniesta (2003), Borrero et al. (2010) and references therein.

3. Stokes profiles at different positions on the solar disk

The inversions carried out in the previous section could be employed to obtain histograms of the magnetic flux density Φ, module of the magnetic field vector B, and the inclination of the magnetic field vector with respect to the observer’s line-of-sight (γ) at different positions on the solar disk. However, in this paper we aimed to infer properties about the distribution of the magnetic field vector by directly studying the histograms of the Stokes profiles. Figure 2a presents distribution histograms of the maximum signals of the Stokes V(λ) (dashed lines) and Stokes Q(λ) and U(λ) (solid lines) normalized to the average quiet-Sun intensity over the entire map: Iqs. The colors indicate each of the different maps studied: red for map A (Sect. 2.1), green for map B (Sect. 2.2), and blue for map C (Sect. 2.3). Figure 2b displays the cumulative histogram of the pixels in each map that have a S/N equal to or higher than a given value. The colors and the line-styles are as in Fig. 2a. For instance: 31.6% of the pixels in map A posses signals in Q or U (solid-red line) that are above 4.5 times the noise level. To limit our analysis to the internetwork regions, we excluded from these figures the pixels in maps A, B, and C with a magnetic flux density Φ > 500 Mx cm-2.

thumbnail Fig. 2

Left panel: histogram of the number of pixels as a function of the maximum value in their polarization signals (normalized to the average continuum intensity). The vertical solid-black line represents the 3σ-level. Right panel: histogram of pixels with polarization signals above a certain S/N. This last histogram is cumulative, while the first one is not. In both panels the color lines indicate the same: red (map A at Λ = 0°), green (map B at Λ = 30°), blue (map C at Λ = 40°). The solid lines refer to the linear polarization profiles (Stokes Q and U) while the dashed lines correspond to the circular polarization profiles (Stokes V). For comparison we also display, in solid and dashed-black curves, the expected behavior due to noise.

We now focus on some of the features of the histograms in Fig. 2. A very striking one is the peak at around the 3σ-value for the linear polarization (solid lines) in Fig. 2a, and the fact that the amount of pixels with maximum linear polarization signals below and above this 3σ-level quickly drops. One might wonder how a peak can appear at around 3σ if the probability that photon noise (taken as a normal random distribution) will produce such a high value is only about 0.27%? The answer to this question is that photon noise has a probability of 0.27% to produce a signal at the 3σ-level at one particular wavelength, but Fig. 2 shows the maximum of the signal over all wavelengths. Indeed, it is possible to employ the binomial distribution to find the lower bound of the probability K(δ,N) that after N wavelengths, one of them will have a signal stronger than or equal to δ-times the noise level σ: K(δ,N)1[1p(δ)]N,\begin{equation} K(\delta,N) \geq 1-[1-p(\delta)]^N , \end{equation}(1)where p(δ) is the probability that at one single wavelength position, a normally distributed random variable (with a standard deviation σ) will yield a signal δ-times above the standard deviation is given by p(δ)=1erf(δ/2),\begin{equation} p(\delta) = 1-\textrm{erf}\left(\delta/\!\sqrt{2}\right) , \end{equation}(2)where erf denotes the so-called error function. Since the spectral line is sampled in N = 112 spectral positions for each Stokes parameter, the probability of finding a wavelength where the noise yields a signal at the 3σ is much higher than the 0.27% mentioned above. In particular, there is a K(N = 112,   δ = 3) ≥ 26.12% probability that the noise in the circular polarization (Stokes V) will yield a signal at the 3σ-level. Because the linear polarization consists of two Stokes parameters (Q and U), this probability is even higher: K(N = 224,   δ = 3) ≥ 45.42%. Given these high probabilities, it is not surprising that the histograms of the Stokes profiles peak close to the 3σ-level. This is certainly the case of the linear polarization (solid lines in Fig. 2a). It is noteworthy that, at the 3σ-location, Stokes V (circular polarization) only presents a local maximum (dashed lines in Fig. 2a) that actually disappears for maps B (dashed green) and C (dashed red). Indeed, the peak for the histogram of Stokes V appears to be located in the range of 2 − 3  ×  10-3 (in units of the quiet Sun continuum intensity), which is about 8 times above the noise level σ.

These differences between linear and circular polarization in Fig. 2a can be explained if we consider a distribution of B ∥  (component of the magnetic field vector that is parallel to the observer’s line-of-sight) that features a peak at a value such that the corresponding Stokes V profile would have a maximum at around 8σ. At the same time, the distribution of B (component of the magnetic field vector that is perpendicular to the observer’s line-of-sight) should feature a low probability of finding values of B that produce Stokes Q and U profiles above 3σ such that there is a peak at this level. We additionally investigated this by synthesizing Stokes profiles with different values of B ∥  and B and comparing them with the maximum signal of the resulting Stokes profiles. To obtain a good estimation we carried out this experiment with two different semi-empirical models that represent granules and intergranules (Borrero & Bellot Rubio 2002). Results are presented in Fig. 3. This figure shows that the distribution of B ∥  that is responsible for the Stokes V signal in Fig. 2a (dashed-lines) must posses a peak at around B ∥  ≈ 5 − 7 G, which is the value needed to produce a majority of Stokes V signals at the 8σ-level. In addition, Fig. 3 shows that the distribution of B must have very low probabilities for B ≥ 40 − 70 G, otherwise the peak in the histograms for Q and U in Fig. 2a (solid lines) would be significantly shifted above the 3σ-level.

thumbnail Fig. 3

Amount of circular polarization (normalized to the quiet-Sun continuum intensity) generated as a function of the component of the magnetic field vector that is parallel to the observer’s line-of-sight: B (dashed lines). Amount of linear polarization (normalized to the quiet-Sun continuum intensity) signals generated as a function of the component of the magnetic field vector that is perpendicular to the observer’s line-of-sight: B (solid lines). These curves were obtained by performing a synthesis of the Stokes profiles employing different values of B and B and using two different atmospheric semi-empirical models: red for intergranules and blue for granules. The horizontal black lines indicate the 3σ and 8σ levels assuming a noise of σ = 3 × 10-4 as in our observed maps (Sect. 2).

An important feature of Figs. 2a and b is that the distribution of the linear polarization signals (solid lines) is different at different positions on the solar surface. In addition, the distribution of the circular polarization signals (dashed lines) also changes slightly, but comparatively less than the linear polarization.

4. Results from theoretical distributions

In the previous section we have determined some general properties about the distribution of the magnetic field vector in the internetwork by considering some details from the observed histograms of the Stokes profiles. We now continue along these lines by investigating the sources of the differences in the distribution of the polarization signals at different positions on the solar disk. Following the notation introduced in Paper II, we refer to X as the set of physical parameters that affect the formation of the Stokes profiles: X=[T,Vlos,B],\begin{eqnarray} \ve{X} = [\ve{T}, V_{\rm los}, \ve{B}, \alpha] , \label{equation:x} \end{eqnarray}(3)where T refers to the thermodynamic parameters, Vlos to the line-of-sight-velocity, B to the magnetic field vector, and finally α refers to the magnetic filling factor (see also Sect. 2). Considering a Milne-Eddington atmosphere (see del Toro Iniesta 2003; Borrero et al. 2010) is equivalent to assuming that with the exception of source function, none of the thermodynamic, kinematic, and magnetic parameters change with optical depth τc in the photosphere. The source function, however, is considered to vary linearly with optical depth: S(τc) = S0 + τcS1, where S0 corresponds to the source function at the observer’s location, and S1 is the derivative of the source function with optical depth. In addition to S0 and S1, the other thermodynamic parameters in T are the Doppler width of the spectral line ΔλD, the damping parameter a, and quotient of the absorption coefficient in the continuum and in the line center η0.

We now refer to \hbox{$\mathcal{P}(\ve{X})\df \ve{X}$} as the probability of finding a pixel within the observed field-of-view where each of the physical parameters in X has values between Xi and Xi + dXi. Furthermore, we assume that the magnetic parameters are statistically independent of the thermodynamic and kinematic parameters, thereby allowing us to write 𝒫(X)dX=𝒫1(B)𝒫2(T,Vlos)dBdTdαdVlos.\begin{eqnarray} \mathcal{P}(\ve{X})\df \ve{X} = \mathcal{P}_1(\ve{B},\alpha) \mathcal{P}_2(\ve{T},V_{\rm los}) \df \ve{B} \df \ve{T} \df \alpha \df V_{\rm los} . \label{equation:pdftot} \end{eqnarray}(4)We now turn our attention to the probability distribution function of the magnetic parameters \hbox{$\mathcal{P}_1$}, which can be rewritten as 𝒫1(B)dBdα=𝒫1(Bi,Bj,Bk)dBidBjdBkdα,\begin{eqnarray} \begin{split} \mathcal{P}_1(\ve{B},\alpha)\df\ve{B} \df \alpha = & \mathcal{P}_1(B_i,B_j,B_k,\alpha) \df B_i \df B_j \df B_k \df \alpha , \end{split} \label{equation:pdfgeneric} \end{eqnarray}(5)which indicates the probability of finding a pixel whose magnetic field vector has the coordinates between the following values: Bi and Bi + dBi, Bj and Bj + dBj, Bk and Bk + dBk, and finally where the filling factor of the magnetic field has a value between α and α + dα. The reason for our choice of nomenclature in the three components of the magnetic field vector Bi, Bj and Bk is that the probability distribution function will generally not be expressed in spherical coordinates in the observer’s reference frame (the frame needed to solve the radiative transfer equation). In general, we therefore must perform a transformation of variables that will express Eq. (5) into the observer’s reference frame and into spherical coordinates: 𝒫1(Bi,Bj,Bk)dBidBjdBkdα=|J|𝒫1(B,γ,φ,α)×dBdγdφdα,\begin{eqnarray} \begin{split} \mathcal{P}_1(B_i,B_j,B_k,\alpha) \df B_i \df B_j \df B_k \df \alpha = \,& |J| \mathcal{P}_1(B,\gamma,\phi,\alpha) \\ & \times \df B \df\gamma \df\phi \df \alpha , \end{split} \label{equation:coortransform} \end{eqnarray}(6)where  | J |  is the determinant of the Jacobian matrix for the transformation between the two reference frames. This transformation might also introduce a dependence on the heliocentric angle Θ, thereby allowing us to evaluate the probability distribution function at different positions on the solar disk. The total probability must be equal to one: 01dα0dB0πdγ02π|J|𝒫1(B,γ,φ,α)dφ=1.\begin{equation} \int_0^1 \df \alpha \int_0^\infty \df B \int_0^\pi \df\gamma \int_0^{2\pi} |J| \mathcal{P}_1(B,\gamma,\phi,\alpha)\df\phi = 1 . \label{equation:normalization} \end{equation}(7)We then use this distribution function to evaluate the percentage of pixels (from the total) that posses a given magnetic field vector. In our simulations we employed a total of 2 × 106 pixels. To solve the radiative transfer equation, we need the probability distribution function of the thermodynamic and kinematic parameters \hbox{$\mathcal{P}_2$} (see Eq. (4)), in addition to the probability distribution function of the magnetic parameters \hbox{$\mathcal{P}_1$}. Hereafter we take, as \hbox{$\mathcal{P}_2(\ve{T},V_{\rm los}) {\textrm d}\ve{T} {\textrm d}V_{\rm los}$}, the distribution obtained from the results of map A (called map C in Paper I). We assume that this distribution does not depend with the position on the solar disk. Finally, in all our tests in this section we assume that the magnetic filling factor is unity: α = 1. We take this approach not to add a new degree of freedom that will make our subsequent analysis more cumbersome. Section 5 will address all the assumptions and simplifications of this section in more detail.

With this we now have all the necessary ingredients needed by the VFISV (Borrero et al. 2010) code to solve the radiative transfer equation in order to obtain theoretical Stokes profiles Q, U, and V. To these profiles we then add noise, assuming a normally distributed random variable (Leva 1992) with a standard deviation σ = 3 × 10-4 (in units of the quiet-Sun continuum intensity). This noise level is similar to that of the observed maps in Sect. 2. Once noise has been added, we select the peak values of the Stokes profiles and construct histograms like those derived from the observations in Figs. 2a, b. The theoretical and observed histograms are then compared with different theoretical distributions of the magnetic field in a attempt to match the observations.

The approach described in this section is indeed an inversion. However, it is not the same kind of inversion as those in Sect. 2 or in Paper I. First of all, the observables here are the histograms of the peak values in Stokes Q, U, and V (Figs. 2a, b), whereas in Paper I the observables were the full Stokes vector (I, Q, U, and V) including their wavelength dependence (not only the peak values) at each individual pixel. The model parameters also differ: while before the model parameters corresponded to the three components of the magnetic field vector for each individual pixel, here the model parameters correspond to a parametrized theoretical distribution function for the magnetic field vector that includes all pixels. As we show below, this parametrized theoretical distribution function possesses a very limited number of free parameters (1–3), which can be tuned to simultaneously fit all pixels in the field-of-view and at different heliocentric angles. This is a great advantage compared to traditional inversion of Stokes profiles (e.g. Borrero & Kobel 2011; Paper I), where there were about ten free parameters (those describing a Milne-Eddington atmosphere) for each pixel in the map. Another important difference is that, while regular inversion codes for the radiative transfer equation (Ruiz Cobo & del Toro Iniesta 1992; Borrero et al. 2010) are automatized, the procedure followed here is completely manual (i.e.: trial and error).

4.1. Isotropic distribution functions

thumbnail Fig. 4

Panel a) histogram of the number of pixels with a given maximum value in their polarization signals (normalized to the average continuum intensity in the quiet Sun). The vertical solid-black line corresponds to the 3σ-level. Panel b) histogram of pixels with polarization signals above a certain S/N. This last histogram is cumulative, but the first one is not. The solid and dashed lines refer to the linear and circular polarization profiles, respectively. These two panels can be readily compared to the observed histograms in Fig. 2. They were obtained employing an isotropic distribution function for the magnetic field vector (Eq. (13)) with three different mean values for the magnetic field strength B0: 20 (red), 30 (blue), and 40 (green) Gauss. The module of the magnetic field vector B and the inclination of this vector with respect to the observer’s line-of-sight γ arising from these isotropic distributions are presented in panels c) and d).

In this section we employ a theoretical distribution where the magnetic field vector is isotropic. We first define as isotropic a probability distribution function where the magnetic field vector has no preferred orientation. In the local reference frame on the solar surface this can be expressed as having a probability of finding a magnetic field vector that depends only on its module B=Bx2+By2+Bz2\hbox{$B=\sqrt{B_{x}^2+B_{y}^2+B_{z}^2}$}: 𝒫1(B)dB=Af(B)dBxdBydBz,\begin{equation} \mathcal{P}_1(\ve{B})\df \ve{B} = A f(B)\df B_{x} \df B_{y} \df B_{z} , \label{equation:isocarte} \end{equation}(8)where A is just a normalization constant. In this reference frame,  { ex,ey,ez } , the ez-axis is perpendicular to the solar surface. To take into account that the observer’s line-of-sight forms an angle Θ (heliocentric angle) with respect to the ez-axis, we perform a variable change into a new coordinate system {ex,ey,ez}\hbox{$\{\expr, \eypr, \ezpr \}$} that is rotated by an angle Θ around the ey-axis. Indeed, because the distribution is isotropic, it does not matter around which axis we consider the rotation. With this transformation the new ez\hbox{$\ezpr$}-axis is aligned with the observer’s line-of-sight. In this case, the relationship between the old coordinates of the magnetic field vector with the new ones is given by {\begin{eqnarray} \begin{cases} B_{x} = {B'}_{z} \sin\Theta + B{'}_{x} \cos\Theta , \\ B_{y} = B{'}_{y} , \\ B_{z} = B{'}_{z} \cos\Theta - B{'}_{x} \sin\Theta , \\ |J| = 1. \end{cases} \label{equation:rotation} \end{eqnarray}(9)Since this is a simple rotation, the determinant of the Jacobian matrix for the transformation is unity. In addition, the module of the magnetic field vector is the same in the old and new reference frames: B=Bx2+By2+Bz2=Bx2+By2+Bz2\hbox{$B=\sqrt{B_{x}^2+B_{y}^2+B_{z}^2} = \sqrt{B_{x}^{'2}+B_{y}^{'2}+B_{z}^{'2}}$}. We can therefore rewrite Eq. (8) in the observer’s reference frame as 𝒫1(B)dB=Af(B)dBxdBydBz,\begin{equation} \mathcal{P}_1(\ve{B})\df \ve{B} = A f(B)\df B{'}_{x} \df B{'}_{y} \df B{'}_{z} , \label{equation:isoobs} \end{equation}(10)which is functionally identical to Eq. (8) since it was defined to be isotropic and therefore independent of the viewing angle Θ. As previously mentioned, to solve the radiative transfer equation we need to express Eq. (10) in spherical coordinates. We therefore perform now an additional variable change where {\begin{eqnarray} \begin{cases} B{'}_{x} = B\sin\gamma\cos\phi , \\ B{'}_{y} = B\sin\gamma\sin\phi , \\ B{'}_{z} = B\cos\gamma , \\ |J| = B^2\sin\gamma . \end{cases} \label{equation:spherical} \end{eqnarray}(11)and thus, 𝒫1(B)dB=Af(B)B2sinγdBdγdφ.\begin{equation} \mathcal{P}_1(\ve{B})\df \ve{B} = A f(B) B^2 \sin\gamma \df B \df\gamma \df\phi . \label{equation:iso} \end{equation}(12)Since the distribution is isotropic, it must also be independent of Θ once it is expressed in the observer’s reference frame both in Cartesian coordinates or spherical ones (Eqs. (10) and (12), respectively). In addition, isotropism manifests itself as a sinγ-dependence, which comes from the determinant of the Jacobian matrix or, in order words, the volume-element expressed in spherical coordinates: dB = B2sinγdBdγdφ. For the distribution of the module of the magnetic field we consider an exponential function in the form f(B) ≈ exp( − B). After normalization, the resulting expression for the theoretical distribution function is 𝒫1(B,γ,φ)dBdγdφ=278πB2B03exp(3BB0)sinγdBdγdφ,\begin{equation} \mathcal{P}_1(B,\gamma,\phi) \df B \df\gamma \df\phi = \frac{27}{8\pi}\frac{B^2}{B_0^3} \exp\left(\frac{-3B}{B_0}\right) \sin\gamma \df B \df\gamma \df\phi , \label{equation:isotropic} \end{equation}(13)where B0 represents the mean value of the magnetic field vector module: B0=B=00π02πB·𝒫1(B,γ,φ)dBdγdφ.\begin{equation} B_0 = \langle B\rangle = \int_0^\infty \int_0^\pi \int_0^{2\pi} B \cdot \mathcal{P}_1(B,\gamma,\phi) \df B \df\gamma \df\phi . \end{equation}(14)Figures 4c–4d display the probability distribution functions for module B and inclination of the magnetic field vector with respect to the observer’s line-of-sight γ, corresponding to Eq. (13), and employing three different values of B0 = 20,30, and 40 G. The histograms for the resulting Stokes profiles after solving the radiative transfer equation with these distributions are given in Figs. 4a, b, which can be readily compared to Figs. 2a, b. This comparison shows that for B0 = 20 G, the linear polarization produced by an isotropic distribution (solid lines) is too low compared to the observed one. For B0 = 40 G, the amount of linear polarization is comparable to the observed one, but, in this case the isotropic distribution produces too much circular polarization (dashed lines). In addition, an isotropic distribution produces far too few Q and U profiles, about 50% of the total, with S/N > 3 (solid lines in Fig. 4b), while the observed distribution shows that about 70–80% of the profiles are above this level (solid lines in Fig. 2b). The misfit between the theoretical histograms and the observed ones becomes even clearer when we consider that all curves in Fig. 4 are independent of the position on the solar disk (as they should, because they correspond to an isotropic distribution; Eq. (13)), while the observed ones (Fig. 2) do change. Choosing a different f(B) function in Eq. (12) does not alter this because the histograms of the Stokes profiles would still be independent of the position on the solar disk. However, it is possible to obtain histograms that vary with the position on the solar disk by postulating that the mean magnetic field depends upon B0(Θ) or B0(Λ). In this case, the different curves in Fig. 4 will certainly change with the position on the solar disk, yet, in doing so we would implicitly introduce a dependence with Θ and/or Λ in the original distribution, given by Eq. (13), and thus it would no longer be isotropic.

4.2. Triple-Gaussian distribution functions

thumbnail Fig. 5

Same as Fig. 4 but for the triple-Gaussian theoretical distribution given by Eq. (19). In this case the mean value of the three components of the magnetic field vector are different: Bx0 = By0 = 22 G, and Bz0 = 5 G. Red lines correspond to Θ = 0° (disk center) while green and blue lines correspond to Θ = 45°. Solid and dashed lines indicate the circular (Stokes V) and linear (Stokes Q and U) polarization.

We now consider a distribution function in the local reference frame on the solar surface,  { ex,ey,ez } , in which each component of the magnetic field is statistically independent of the rest and shows a Gaussian-like dependence: 𝒫1(B)dB=dBxdBydBzπ3Bx0By0Bz0expBx2πBx02×expBy2πBy02expBz2πBz02·\begin{eqnarray} \begin{split} \mathcal{P}_1 (\ve{B})\df \ve{B} & = \frac{\df B_{x} \df B_{y} \df B_{z}}{\pi^3 B_{x0}B_{y0} B_{z0}} \exp{\left\{-\frac{B_{x}^2}{\pi B_{x0}^2}\right\}} \\ & \times \exp{\left\{-\frac{B_{y}^2}{\pi B_{y0}^2}\right\}} \exp{\left\{-\frac{B_{z}^2}{\pi B_{z0}^2}\right\}} \cdot \end{split} \label{equation:triplelocal} \end{eqnarray}(15)Again, with this definition, the ez-axis is perpendicular to the solar surface. The distribution is of course normalized to one: $𝒫1(Bx,By,Bz)dBxdBydBz=1.\begin{equation} \iiint_{-\infty}^{\infty} \mathcal{P}_1(B_{x},B_{y},B_{z}) \df B_{x} \df B_{y} \df B_{z} = 1 . \label{equation:norma} \end{equation}(16)In the previous equation, Bx0, By0, and Bz0 refer to the mean of the absolute value for each component of the magnetic field: Bx0=$+|Bx|𝒫1(Bx,By,Bz)dBxdBydBz.\begin{equation} B_{x0} = \iiint_{-\infty}^{+\infty} | B_{x} | \mathcal{P}_1 (B_{x},B_{y},B_{z})\df B_{x} \df B_{y} \df B_{z} . \end{equation}(17)Similar equations can be written for By0 and Bz0. As in the previous section, we perform a rotation of angle Θ around the ey-axis to transform Eq. (15) into the observer’s reference frame: {ex,ey,ez}\hbox{$\{\expr, \eypr, \ezpr \}$}. To guarantee that the result of this transformation does not depend on the choice of rotation axis in the XY-plane we must ensure that Bx0 = By0 such that Eq. (15) can be written in terms of Bx2+By2\hbox{$B_{x}^2+B_{y}^2$}. With this, we can now rewrite the theoretical distribution function as 𝒫1(B)dB=dBxdBydBzπ3Bx0By0Bz0expBy2πBy02×exp(BzsinΘ+BxcosΘ)2πBx02×exp(BzcosΘBxsinΘ)2πBz02·\begin{eqnarray} \begin{split} \mathcal{P}_1 (\ve{B})\df \ve{B} = \,& \frac{\df B{'}_{x} \df B{'}_{y} \df B{'}_{z}}{\pi^3 B_{x0}B_{y0} B_{z0}} \exp{\left\{-\frac{{B'}_{y}^{2}}{\pi B_{y0}^2}\right\}} \\ & \times \exp{\left\{-\frac{(B{'}_{z} \sin\Theta + B{'}_{x} \cos\Theta)^2}{\pi B_{x0}^2}\right\}} \\ & \times \exp{\left\{-\frac{(B{'}_{z} \cos\Theta - B{'}_{x} \sin\Theta)^2}{\pi B_{z0}^2}\right\}} \cdot \end{split} \label{equation:tripleobs} \end{eqnarray}(18)And finally, to solve the radiative transfer equation we need to express Eq. (18) in spherical coordinates (Eq. (11)). This transformation yields 𝒫1(B,γ,φ)dBdγdφ=B2sinγdBdγdφπ3Bx0By0Bz0×exp(Bsinγsinφ)2πBy02×exp(BcosγsinΘ+BsinγcosφcosΘ)2πBx02×exp(BcosγcosΘBsinγcosφsinΘ)2πBz02·\begin{eqnarray} \begin{split} \mathcal{P}_1 & (B,\gamma,\phi)\df B \df\gamma\df\phi =\frac{B^2 \sin\gamma\df B\df\gamma\df\phi}{\pi^3 B_{x0}B_{y0} B_{z0}} \\ & \times \exp{\left\{-\frac{(B\sin\gamma\sin\phi)^2}{\pi B_{y0}^2}\right\}} \\ & \times \exp{\left\{-\frac{(B\cos\gamma\sin\Theta+B\sin\gamma\cos\phi\cos\Theta)^2}{\pi B_{x0}^2}\right\}} \\ & \times \exp{\left\{-\frac{(B\cos\gamma\cos\Theta-B\sin\gamma\cos\phi\sin\Theta)^2}{\pi B_{z0}^2}\right\}} \cdot \end{split} \label{equation:triplesphe} \end{eqnarray}(19)With this theoretical distribution function of the magnetic field vector we can again (see Sect. 4.1) synthesize the Stokes profiles, add noise, and finally construct the theoretical histograms of the resulting Stokes profiles. Equations (18) and (19) now entail an explicit dependence on Θ, unlike the case of an isotropic distribution (Eq. (13)). Interestingly, if Bx0 = By0 = Bz0, the distribution becomes isotropic, because Eq. (15) can be re-written in terms of the module of the magnetic field vector B (see Eq. (8)). In the following we employ Bx0 = By0 = 22G;Bz0 = 5 G, and therefore the resulting distribution is non-isotropic. However, as indicated above, we impose Bx0 = By0 so that the distribution is symmetric in the XY-plane. The distributions for B and γ corresponding to the aforementioned values are presented in Figs. 5c and 5d for two different heliocentric angles: Θ = 0° (red curves), Θ = 45° (blue). The reason for employing Bz0 < Bx0,By0 is because in the isotropic case (see Fig. 4a) the value of B0 that produced a reasonable fit to the linear polarization at disk center yielded too much circular polarization and thus, to simultaneously fit both, we need to decrease the vertical component of the magnetic field vector while keeping the horizontal component at the appropriate level.

The first feature to notice is that the distribution of B (Fig. 5c) is independent of the heliocentric angle Θ. This was to be expected because the module of the magnetic field vector B is invariant with respect to rotations. The second feature to notice is that the distribution of γ shows a peak at γ = 90° at disk center (red curve). This is a consequence of having imposed Bz0 < Bx0,By0. This peak smoothes out as we move toward the poles (blue curves).

The theoretical histograms for the Stokes profiles resulting from the previous distribution (Eq. (19)) are displayed in Figs. 5a and 5b. As mentioned above, the values of Bx0, By0, and Bz0 were selected such that, for signals above 2 × 10-3, the theoretical histograms at disk center (Θ = 0°; red curves) are comparable to the observed ones (Fig. 2). However, at larger heliocentric angles (Θ = 45°; blue curves), the mismatch between theoretical and observed histograms is very clear, in particular for the circular polarization (dashed lines). Our theoretical distribution function (Eq. (19)) clearly predicts that the circular polarization should increase as Θ increases (blue dashed lines in Fig. 5). This happens as a consequence of having imposed Bx0 > Bz0, which means that away from disk center, the component of the magnetic field vector that is aligned with the observer’s line-of-sight, Bz = B, increases due to the contribution from Bx (see Eq. (9)). Because the slope of the V − B curve is so steep (see Fig. 3), a small increase in Bz = B translates into a large increase in the Stokes V signal. Therefore, the generated amount of circular polarization is much larger at larger heliocentric angles (dashed-blue curves; Θ = 45°) than at disk center (dashed-red curves; Θ = 0°). Likewise, since Bz0 < Bx0, the component of the magnetic field that is perpendicular to the observer’s line-of-sight B decreases as Θ increases. In this case, however, the linear polarization does not decrease as much as the circular polarization increases. This is due to the gentler slope of the Q − B curve (Fig. 3).

Interestingly, Figs. 2a, b show that the observed amount of circular polarization does not increase, as required by the change in the viewing angle, for larger heliocentric angles (i.e. as Θ increases) in accordance to the observed drop in linear polarization. We conclude therefore that the differences in the histograms of the observed Stokes profiles at different positions on the solar disk cannot be produced by the change in the viewing angle Θ. This problem is not unique to the theoretical distribution given by Eq. (19), but it will indeed affect any theoretical distribution that is prescribed in the local reference frame and features Bz0 < Bx0,By0. Prescribing a theoretical distribution in the local reference frame means that the underlying distribution is always the same regardless of the position on the solar disk. The only reason that this distribution changes is because the angle Θ, between the observer’s line of sight and the vector perpendicular to the solar surface, varies with the position on the solar disk.

4.3. Other distribution functions

To avoid the problem described in Sect. 4.2 we now prescribe a theoretical distribution of the magnetic field vector that, already in the local reference frame, depends on the latitude Λ. This will imply, unlike the previous distributions in Sects. 4.1 and 4.2, that the underlying distribution (i.e. in the local reference frame) is different at different positions on the solar disk. In particular, the Λ dependence means that the distribution function varies toward the poles, but not toward the limbs. The expression chosen in this section is 𝒫1(B)dB=dBxdBydBz192πβ5β3(1+ζsinΛ)5×[(BxcosΛ+BzsinΛ)2+By2]32[BzcosΛBxsinΛ]2×exp{|BzcosΛBxsinΛ|β}×exp|[(BxcosΛ+BzsinΛ)2+By2]12|(1+ζsinΛ)β·\begin{eqnarray} \begin{split} & \mathcal{P}_1(\ve{B}; \Lambda)\df \ve{B} = \frac{\df B_x \df B_y \df B_z}{192\pi\beta_\perp^5\beta_\parallel^3} (1+\zeta\sin\Lambda)^5 \\ & \times \left[(B_x\cos\Lambda+B_z\sin\Lambda)^2+B_y^2\right]^{\frac{3}{2}} \left[B_z\cos\Lambda-B_x\sin\Lambda\right]^2\\ & \times \exp\left\{-\frac{|B_z\cos\Lambda-B_x\sin\Lambda|}{\beta_\parallel}\right\} \\ & \times \exp\left\{-\frac{|[(B_x\cos\Lambda+B_z\sin\Lambda)^2+B_y^2]^{\frac{1}{2}}|(1+\zeta\sin\Lambda)}{\beta_\perp}\right\} \cdot \end{split} \label{equation:othercartesian} \end{eqnarray}(20)This theoretical distribution is normalized to one as it verifies Eq. (16). As in the two previous sections, we now perform a rotation of angle Θ around the ey-axis to express Eq. (20) in the observer’s reference frame, and then we transform into spherical coordinates in the observer’s reference frame. This is done by applying Eqs. (9) and (11). A close inspection (e.g. make for instance Λ = 0) of Eq. (20) shows that again the result does not depend on the choice of rotation axis in the XY-plane. Although trivial to obtain, the resulting expression for the distribution function in spherical coordinates in the observer’s reference frame is too long to be written here (we provide it for a simplified case below). It is noteworthy to mention, however, that even though Eq. (20) depends on Λ, once this equation is expressed in the observer’s reference frame, it depends both on Λ and Θ: \hbox{$\mathcal{P}_1(\ve{B}; \Lambda, \Theta)$}.

thumbnail Fig. 6

Same as Fig. 4 but for the theoretical distribution given by Eq. (20). The parameters employed are β ∥  = 3 G and β = 7 G. Red lines correspond to Λ = 0° (disk center), while green and blue lines correspond to Λ = 30° and Λ = 40°, respectively. Solid and dashed lines indicate the circular (Stokes V) and linear (Stokes Q and U) polarization.

The theoretical distributions functions of B and γ resulting from Eq. (20) are shown in Figs. 6c and d, respectively. Here we considered the following values: ζ = 0.3, β = 3, and β = 7. Unlike Sects. 4.1 and 4.2, β and β do not correspond to the mean values of B and B. For this reason we refer to them as β and β instead of B∥0 and B⊥0. In this new distribution, the module of the magnetic field vector B is different at different latitudes, with a slightly higher mean value at disk center (red line; Λ = 0°) than closer to the poles (green and blue lines; Λ = 30 − 40°). This is a direct consequence of prescribing a distribution of the magnetic field vector in the local reference frame that depends on Λ. Otherwise, as it occurred in Sects. 4.1 and 4.2, the distribution of B would be the same, since the module of the magnetic field vector is invariant with respect to rotations.

We now use the distribution given by Eq. (20) to solve the radiative transfer equation and produce theoretical histograms of the Stokes profiles. These are displayed in Figs. 6a and b. The shape of the curves for the circular (dashed) and linear (solid) polarization in Fig. 6a are very similar to the observed ones (Fig. 2a). Of particular interest is the fact that, in agreement with the observations, the histograms of the circular polarization are very similar in the three considered positions on the solar disk, and feature a peak at signals  ≈ 3 − 4 × 10-3. It is also important to mention that the linear polarization decreases from Λ = 0° toward Λ = 30−40°, which also agrees with the observations.

To understand how Eq. (20) can fit the observed histograms of the Stokes profiles at different latitudes, it is convenient to re-write this theoretical distribution function assuming that the observer’s line-of-sight is along the zero meridian. In this case the rotation angle Θ is equal to Λ, thereby simplifying the probability distribution function from Eq. (20) into (in spherical coordinates in the observer’s reference frame) 𝒫1(B,γ,φ)dBdγdφ=B7sin4γcos2γ192πβ3β5(1+ζsinΛ)5×exp{|Bsinγ|(1+ζsinΛ)β}×exp{|Bcosγβ}dBdγdφ.\begin{eqnarray} \begin{split} \mathcal{P}_1 (B,\gamma,\phi)\df B \df\gamma\df\phi = \,& \frac{B^7\sin^4\gamma\cos^2\gamma}{192 \pi \beta_\parallel^3 \beta_\perp^5} (1+\zeta\sin\Lambda)^5 \\ & \times \exp{\left\{-\frac{|B\sin\gamma|(1+\zeta\sin\Lambda)}{\beta_\perp}\right\}} \\ & \times \exp{\left\{-\frac{|B\cos\gamma}{\beta_\parallel}\right\}} \df B \df \gamma \df \phi. \end{split} \label{equation:othersphe} \end{eqnarray}(21)In this Eq. (21), the exponential term that refers to Bcosγ (i.e. component of the magnetic field vector that is parallel to the observer’s line-of-sight) does not depend on Λ and therefore the amount of circular polarization generated by this distribution will not change with latitude (see dashed lines in Figs. 6a and b). Moreover, the exponential term that contains Bsinγ (i.e. component of the magnetic field vector that is perpendicular to the observer’s line-of-sight) decreases as Λ increases and thus, this distribution function produces less linear polarization at higher latitudes (see solid lines in Figs. 6a and b).

Despite all these similarities between the theoretical and observed distribution of Stokes profiles there are still some significant differences, namely: a) there is a clear deficit of pixels with sufficient amount of linear polarization in the range of signals  ≈ 2 − 4 × 10-3 (Fig. 6a) or S/N > 5 (Fig. 6b) at all latitudes; and b) the theoretical histograms show a peak in the circular polarization at the 3σ-level at high latitudes (dashed-blue and dashed-green curves in Fig. 6a), while the observed histograms show no peak in the circular polarization at the 3σ-level at any disk position (dashed curves in Fig. 2a). Overall, however, the fit between the observed and the theoretical histograms of the Stokes profiles is clearly better for this distribution (Eq. (20)) than for the case of an isotropic distribution (Eq. (13)) or a triple Gaussian (Eq. (19)).

5. Assumptions and limitations

In the previous subsections we have focused on the effect that the probability distribution function of the magnetic field vector \hbox{$\mathcal{P}_1(\ve{B},\alpha)$} has on the observed histograms of the Stokes profiles at different positions on the solar disk. In our theoretical analysis a number of simplifying assumptions were made to keep the problem tractable. Although they have already been pointed out in Sect. 4, we summarize them here to briefly discuss their implications.

  • We have assumed that the thermodynamic and magnetic pa-rameters are statistically independent of each other. This al-lowed us to write the total probability distribution function inEq. (4) as the product of two distinct probability dis-tribution functions. However, as dictated by the Lorentz-forceterm in the momentum equation in magnetohydrodynamics,the magnetic field affects the thermodynamic structure of thesolar atmosphere. It is therefore clear that this assumption doesnot fully hold in the solar atmosphere. For instance, if we take thecommonly accepted picture of intergranular lanes harboringmore vertical and stronger magnetic fields than the granular cells,and we consider that intergranular cells have a smoother variationof the temperature with optical depth (see i.e. Fig. 3in Borrero & Bellot Rubio 2002), we could then have postulated acorrelation between the magnetic field vector B and the gradient of the source function with optical depth S1, which is contained in T (Eq. (3)). Indeed, the higher S1 is, the stronger will be the polarization profiles Q, U, and V (see Eq. (9.45) in del Toro Iniesta 2003). These correlations could potentially cause the observed histograms of the Stokes profiles (Fig. 2) to vary with the heliocentric angle, even if the underlying distribution of the magnetic field vector does not depend on Θ. Therefore it is important to investigate what an effect they have before conclusively proving that the distribution of the magnetic field vector is not isotropic (Sect. 4.1), or that the differences in the observed histograms of the Stokes profiles are not due to the viewing angle (Sect. 4.2). Unfortunately, the aforementioned correlations are not known for the solar internetwork simply because it is not clear how magnetic fields are distributed here. In the future we will explore this question by employing 3D numerical simulations of the solar atmosphere, because they provide correlations between B and T that are compatible with the MHD equations.

  • We have also assumed that the probability distribution function of the thermodynamic and kinematic parameters, \hbox{$\mathcal{P}_2(\ve{T},V_{\rm los})$}, does not depend on the position on the solar disk. Kinematic parameters (i.e. line-of-sight velocity Vlos) do not influence our study, since they have no effect on the amplitude of the Stokes profiles in Fig. 2. The same can be argued about other thermodynamic parameters in T, such as the source function at the observer’s S0 (affects only Stokes I), and damping parameter a (affects mostly the line width but not its amplitude). By far, the most important thermodynamic parameters affecting the amplitude of the Stokes profiles under the Milne-Eddington approximation are the gradient of the source function with optical depth S1 and the continuum-to-line-center absorption coefficient η0. In a 1D atmosphere both these parameters are known to decrease as Θ increases, because the line-of-sight samples a thinner vertical-portion of the atmosphere. However, since the dependence of the polarization profiles with S1 and η0 are identical (see Eqs. (8.14), (8.15) and (9.44) in del Toro Iniesta 2003), one would expect the same drop with increasing heliocentric angle angle in the amplitude of the circular polarization profiles (Stokes V) and linear polarization profiles (Stokes Q and U). However, Fig. 2 shows that the linear and circular polarization profiles (solid-color and dashed-color lines) behave differently, and therefore we can rule out the variations of S1 and/or η0 with Θ as being responsible for the observed histograms in the Stokes profiles. Of course, this would change in a 3D atmosphere, where the line-of-sight pierces through different inhomogeneous atmospheric layers, thereby opening the door for the possibility of S1 and or η0 to affect the linear and circular polarization profiles differently1.

  • Adopting a Milne-Eddington atmosphere also implies that we are assuming that the magnetic field vector B does not vary with optical depth τc in the photosphere. This can have important consequences, since at larger heliocentric angles the spectral line samples higher atmospheric layers than at disk center, where the probability distribution function of the magnetic field vector can be different. Employing the widely used 1D HOLMU model (Holweger & Müller 1974) we calculated that the continuum level τc = 1 rises by approximately 20 Km from disk center Θ = 0° (map A; Sect. 2.1) to Θ = 45° (maps B and C; Sects. 2.2 and 2.3). Since this vertical shift of the continuum level is rather small, we could argue that the histograms of the Stokes profiles in Fig. 2 are not affected by this effect. However, the value of 20 km should be considered only as a lower limit since a 1D model does not take into account the horizontal inhomogeneities present in the solar atmosphere. To properly account for this effect, more sophisticated 3D models should be employed.

  • Finally, we have considered α = 1 in our analysis (see Sect. 4). This is equivalent to considering that, at the resolution of the Hinode/SP instrument (0.32′′; Sect. 2), the magnetic structures are spatially resolved. This is, of course, highly unlikely, and therefore it would be important to drop this assumption in the future. Its importance can only be quantified with additional assumptions about the scale-distribution of the magnetic structures in the solar photosphere. This topic is, in itself, as controversial as the distribution of the magnetic field strength and inclination, which is the reason why we have refrained from addressing it here. Although employing 3D MHD simulations would certainly help to drop the α = 1 assumption, we are cautious about it since it is not clear whether these simulations are reliable at the smallest physical scales (Sánchez Almeida 2006).

6. Discussion and conclusions

The histograms of the observed Stokes profiles at different positions on the solar disk (Fig. 2) are clearly different from each other. One possible interpretation for this is that the distribution of the magnetic field vector in the solar internetwork is not isotropic. We explored this possibility in Sect. 4.1, where we employed an isotropic probability distribution of the magnetic field vector. This distribution yielded, as expected, the same distribution of Stokes profiles at all positions on the solar disk (Fig. 4). Martínez González et al. (2008) have also presented similar histograms but employing the Stokes profiles from the Fe I line pair at 1.56 μm (observed with the TIP2 instrument; Martínez Pillet et al. 1999). Their histograms (see their Fig. 2) showed no clear variation with the heliocentric angle, which lead them to conclude that the distribution of the magnetic field vector in the quiet Sun was isotropic. Interestingly, these authors also mentioned after a more detailed analysis that there could indeed be a dependence of the histograms with the heliocentric angle (as indeed we find here).

In addition to Martínez González et al. (2008), a number of works have also argued in favor of an isotropic distribution of magnetic fields in the internetwork. In particular, Asensio Ramos (2009) and Stenflo (2010), employing two different approaches, both concluded that for very weak magnetic fields (B → 0) the distribution becomes isotropic. With our present data we cannot argue against or in favor of this interpretation. The main reason for this is that, as discussed in Sect. 3, any distribution for the magnetic field vector where B has a peak below 40–70 will produce linear polarization profiles that are dominated by noise (3σ-level or S/N = 3). Therefore our current approach (described in Sect. 4) cannot be employed to discern the underlying distribution of the magnetic field vector from these profiles dominated by noise. However, it can be employed to establish that the number of pixels that would follow this hypothetically isotropic distribution cannot be much larger than 30% of the pixels in the internetwork, since this is the amount of pixels that show a peak at the 3σ-level in the polarization profiles (see Fig. 2a). For signals above  > 2 × 10-3, the histograms of the Stokes profiles deviate significantly from the ones predicted by an isotropic distribution, and thus we can establish that here the distribution of the magnetic field vector cannot be isotropic.

We can use a different argument to further clarify the previous point. Our theoretical distributions in Sect. 4 apply to all possible values of the module of the magnetic field vector. However, we could have employed distributions pieced together in the following form: 𝒫1(B)dB={\begin{eqnarray} \mathcal{P}_1(\ve{B})\df \ve{B} = \begin{cases} \mathcal{P}_a(\ve{B})\df\ve{B}, & \textrm{if}\; B<B^{*} \\ \mathcal{P}_b(\ve{B})\df \ve{B} , & \textrm{if}\; B>B^{*} , \end{cases} \label{equation:piecepdf} \end{eqnarray}(22)where \hbox{$\mathcal{P}_a$} could hypothetically correspond to an isotropic distribution for weak fields: B < B ∗ . This would explain the 3σ-peak in the linear polarization in Fig. 2a (dashed lines). In addition, \hbox{$\mathcal{P}_b$} could be a distribution, valid for larger fields B > B ∗ , that would fit the tails of the histogram. The distribution given by Eq. (22) does not need to be discontinuous because it could be prescribed such that \hbox{$\mathcal{P}_a(B^{*}) = \mathcal{P}_b(B^{*})$}.

In Sect. 4.2 we employed a triple Gaussian (one for each component of the magnetic field vector) distribution function and found that, under this assumption and at disk center, the best fit to the observed histograms of the Stokes profiles is produced by a distribution in which the mean value of the magnetic field vector component that is parallel to the solar surface is lower than the mean value of the magnetic field vector component that is perpendicular to the solar surface: Bz0<Bx02+By02\hbox{$B_{z0} < \sqrt{B_{x0}^2+B_{y0}^2}$}. This yields a distribution function where the magnetic field vector is highly inclined, in agreement with previous findings from Orozco et al. (2007a,b) and Lites et al. (2007, 2008). However, this distribution does not fit well the histograms of the Stokes profiles at other positions on the solar disk. In fact, in that section we found that it is not possible to fit the observed histograms for the Stokes profiles at different heliocentric angles employing a theoretical distribution function for the magnetic field vector prescribed in the local reference frame that only changes due to the viewing angle Θ. The reason for this is that, for an underlying distribution where the magnetic field vector is mostly horizontal at Θ = 0° (disk center), the amount of linear polarization slightly decreases when Θ increases, while the amount of circular polarization would significantly increase as Θ increases. However, the observed histograms of the Stokes profiles (Fig. 2) show that, although the amount of linear polarization decreases when Θ increases, the circular polarization does not particularly increase (see also discussion in Lites et al. 2008). This cannot be explained in terms of a simple rotation of the viewing angle Θ, and therefore we interpreted this fact, in Sect. 4.3, as proof that the underlying (i.e. in the local reference frame) distribution of the magnetic field vector must depend on the position on the solar disk.

Under the assumption that the distribution of the underlying magnetic field vector depends on the latitude Λ (see Sect. 4.3)), we were able to find a theoretical distribution of the magnetic field vector (Eq. (20)) that fits quite well the observed histograms of the Stokes profiles at different positions on the solar disk (Fig. 6). Among other properties, this distribution features a magnetic field whose mean value decreases toward the poles. We note here that this does not mean that this is the real distribution for the magnetic field vector present in the quiet Sun. One reason for this is that the fit is far from perfect (see discrepancies mentioned in Sect. 4.3), but most importantly, that we do not know whether this solution is unique because there can be other theoretical distributions that fit the observed Stokes profiles equally well, or even better. More work is indeed needed to confirm or rule out Eq. (20) as the real distribution of the magnetic field vector present in the Sun. In particular, a better fit to the observed histograms of the Stokes profiles is desirable. In addition, it is important to have maps at more latitudes to further constrain the possible distribution functions.

Moreover, it is important to bear in mind that the conclusions above are not necessarily the only possible interpretations, because postulating a probability distribution function of the thermodynamic and kinematic parameters, \hbox{$\mathcal{P}_2(\ve{T},V_{\rm los})$}, that varies with the heliocentric angle Θ, or postulating a correlation between the thermodynamic (T) and magnetic (B) parameters might also help explain the observed differences between the histograms of the Stokes profiles at different positions on the solar disk. Another effect that has not been accounted for is that magnetic field vector can vary with optical depth in the solar photosphere. Since the Stokes profiles sample increasingly higher atmospheric layers as the heliocentric angle increases, the distribution of the magnetic field vector can be different for different values of Θ, even if the probability distribution of the magnetic field vector is the same at all positions on the solar disk at a fixed geometrical depth. All these effects could be properly accounted for by means of 3D MHD simulations of the solar photosphere (Schüssler & Vögler 2008; Steiner et al. 2008, 2009; Danilovic et al. 2010).

In the future we expect to employ such simulations to either rule out or confirm our results in this paper. Consequently, our conclusions at this point should be regarded as preliminary only. Instead, the main purpose in this paper is to illustrate the methodology detailed in Sects. 3 and 4 to study the distribution of the magnetic field vector in the quiet Sun, by directly inverting the histograms of the Stokes profiles in entire maps instead of inverting the Stokes profiles at each spatial position in a given map. Our method has great potential to investigate several aspects of the photospheric magnetism in the solar internetwork. For instance, it can be used, as in Sect. 4.3, to confirm whether the mean value of the distribution of the magnetic field vector changes from disk center toward the poles (cf. Zwang 1987; Ito et al. 2010). This will have important consequences for theoretical models that explain the torsional oscillations in the butterfly diagram in terms of a geostrophic flow model (Spruit 2003), which requires a significant amount of magnetic flux at high latitudes at the beginning of the sunspot cycle. In addition, and although in this work we have restricted ourselves to variations in latitude (Λ), additional observations from disk center toward the solar limbs could be employed to investigate whether the properties of the magnetic field in the internetwork change also in longitude. This is already predicted by non-axisymmetric dynamo models (Moss 1991; Moss et al. 1999; Bigazzi & Ruzmaikin 2004; Charbonneau 2005) and can provide important clues about the strength of the differential rotation (Rüdiger & Elstner 1994; Zhang et al. 2003).


1

Outside the Milne-Eddington approximation, a similar argument could be made in terms of the stratification with optical depth of the temperature in the solar photosphere T(τc), instead of S1 and η0.

Acknowledgments

We would like to thank Luis Bellot Rubio, Mariam Martínez González and Oskar Steiner for fruitful discussions in the subject. Many thanks also to an anonymous referee, who pointed out an error in the probability distribution function of Sect. 4.2 in an early version of this manuscript. This work analyzes data from the Hinode spacecraft. Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Scientific operation of the Hinode mission is conducted by the Hinode science team organized at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner countries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), STFC (UK), NASA, ESA, and NSC (Norway). This work has also made use of the NASA ADS database.

References

  1. Asensio Ramos, A. 2009, ApJ, 701, 1032 [NASA ADS] [CrossRef] [Google Scholar]
  2. Asensio Ramos, A., Martínez González, M. J., López Ariste, A., et al. 2007, ApJ, 659, 829 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bellot Rubio, L. R., & Orozco Suárez, D. 2012, ApJ, 757, 19 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bigazzi, A., & Ruzmaikin, A. 2004, ApJ, 604, 944 [NASA ADS] [CrossRef] [Google Scholar]
  5. Borrero, J. M., & Bellot Rubio, L. R. 2002, A&A, 385, 1056 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Borrero, J. M., & Kobel, P. 2011, A&A, 527, A29 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Borrero, J. M., & Kobel, P. 2012, A&A, 547, A89 (Paper II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Borrero, J. M., Tomczyk, S., Kubo, M., et al. 2011, Sol. Phys., 273, 267 [Google Scholar]
  9. Charbonneau, P. 2005, Liv. Rev. Sol. Phys., 2 (cited on 2010) [Google Scholar]
  10. Danilovic, S., Schüssler, M., & Solanki, S. K. 2010, A&A, 513, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Domínguez Cerdeña, I., Sánchez Almeida, J., & Kneer, F. 2003, ApJ, 582, 55 [Google Scholar]
  12. Domínguez Cerdeña, I., Sánchez Almeida, J., & Kneer, F. 2006, ApJ, 636, 496 [NASA ADS] [CrossRef] [Google Scholar]
  13. Harvey, J. W., Bratson, D., Henney, C. J., & Keller, C. U. 2007, ApJ, 659, L177 [NASA ADS] [CrossRef] [Google Scholar]
  14. Henney, C. J., & Harvey, J. W. 2002, Sol. Phys., 207, 199 [NASA ADS] [CrossRef] [Google Scholar]
  15. Holweger, H., & Müller, E. A. 1974, Sol. Phys., 39, 19 [NASA ADS] [CrossRef] [Google Scholar]
  16. Ichimoto, K., Lites, B. W., Elmore, D., et al. 2008, Sol. Phys., 249, 233 [NASA ADS] [CrossRef] [Google Scholar]
  17. Ito, H., Tsuneta, S., Shiota, D., Tokumaru, M., & Fukiri, K. 2010, ApJ, 719, 131 [NASA ADS] [CrossRef] [Google Scholar]
  18. Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [NASA ADS] [CrossRef] [Google Scholar]
  19. Leva, J. L. 1992, ACM Transactions of Mathematical Software, 18, 454 [CrossRef] [Google Scholar]
  20. Lites, B. W., Socas-Navarro, H., Kubo, M., et al. 2007, PASJ, 59, 571 [Google Scholar]
  21. Lites, B. W., Kubo, M., Socas-Navarro, H., et al. 2008, ApJ, 672, 1237 [NASA ADS] [CrossRef] [Google Scholar]
  22. López Ariste, A., Martínez González, M. J., & Ramírez Vélez, J. C. 2007, A&A, 464, 351 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Martínez Pillet, V., Collados, M., Sánchez Almeida, J., et al. 1999, ASPC, 183, 264 [Google Scholar]
  24. Martínez González, M., Asensio Ramos, A., López Ariste, A., & Manso-Sainz, R. 2008, A&A, 479, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Moss, D. 1999, MNRAS, 306, 300 [NASA ADS] [CrossRef] [Google Scholar]
  26. Moss, D., Brandenburg, A., & Tuominen, I. 1991, A&A, 347, 576 [NASA ADS] [Google Scholar]
  27. Orozco Suárez, D., & Bellot Rubio, L. R. 2012, ApJ, 751, 1 [NASA ADS] [CrossRef] [Google Scholar]
  28. Orozco Suárez, D., Bellot Rubio, L. R., del Toro Iniesta, J. C., et al. 2007a, ApJ, 670, L61 [NASA ADS] [CrossRef] [Google Scholar]
  29. Orozco Suárez, D., Bellot Rubio, L. R., del Toro Iniesta, J. C., et al. 2007b, PASJ, 59, 837 [Google Scholar]
  30. Rüdiger, G., & Elstner, D. 1994, A&A, 281, 46 [NASA ADS] [Google Scholar]
  31. Ruiz Cobo, B., & del Toro Iniesta, J. C. 1992, ApJ, 398, 375 [NASA ADS] [CrossRef] [Google Scholar]
  32. Sánchez Almeida, J. 2005, A&A, 438, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Sánchez Almeida, J. 2006, A&A, 450, 1199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Sánchez Almeida, J., & Martínez González, M. J. 2011, in Proc. Solar Polarization Workshop 6, eds. J. R. Kuhn, D. M. Harrington, H. Lin, S. V. Berdyugina, J. Trujillo-Bueno, S. L. Keil, & T. Rimmele (San Francisco: ASP), 2011, 451 [Google Scholar]
  35. Schüssler, M., & Vögler, A. 2008, A&A, 481, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Steiner, O., Rezaei, R., Schaffenberger, W., & Wedemeyer-Böhm, S. 2008, ApJ, 680, L85 [NASA ADS] [CrossRef] [Google Scholar]
  37. Steiner, O., Rezaei, R., & Schlichenmaier, R. 2009, in Proceedings for the Second Hinode Science Meeting, eds. B. Lites, M. Cheung, T. Nagara, J. Mariska, & K. Reeves, ASP Conf. Ser., 415, 67 [Google Scholar]
  38. Shimizu, T., Nagata, S., Tsuneta, S., et al. 2008, Sol. Phys., 249, 221 [NASA ADS] [CrossRef] [Google Scholar]
  39. Spruit, H. 2003, Sol. Phys., 213, 1 [Google Scholar]
  40. Stenflo, J. O. 2010, A&A, 517, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Suematsu, Y., Tsuneta, S., Ichimoto, K., et al. 2008, Sol. Phys., 249, 197 [NASA ADS] [CrossRef] [Google Scholar]
  42. del Toro Iniesta, J. C. 2003, Introduction to Spectropolarimetry (Cambridge, UK: Cambridge University Press) [Google Scholar]
  43. Tsuneta, S., Suetmatsu, Y., Ichimoto, K., et al. 2008, Sol. Phys., 249, 167 [NASA ADS] [CrossRef] [Google Scholar]
  44. Zhang, K., Chan, K. H., Zou, J., Liao, X., & Schubert, G. 2003, ApJ, 596, 663 [NASA ADS] [CrossRef] [Google Scholar]
  45. Zwang, C. 1987, ARA&A, 25, 83 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

From top to bottom: inferred magnetic flux density Φ = αB from the inversion of Map A (Sect. 2.1), Map B (Sect. 2.2), and Map C (Sect. 2.3). White areas correspond to regions where all three polarization profiles (Stokes Q, U, and V) are below the 4.5σ-level.

In the text
thumbnail Fig. 2

Left panel: histogram of the number of pixels as a function of the maximum value in their polarization signals (normalized to the average continuum intensity). The vertical solid-black line represents the 3σ-level. Right panel: histogram of pixels with polarization signals above a certain S/N. This last histogram is cumulative, while the first one is not. In both panels the color lines indicate the same: red (map A at Λ = 0°), green (map B at Λ = 30°), blue (map C at Λ = 40°). The solid lines refer to the linear polarization profiles (Stokes Q and U) while the dashed lines correspond to the circular polarization profiles (Stokes V). For comparison we also display, in solid and dashed-black curves, the expected behavior due to noise.

In the text
thumbnail Fig. 3

Amount of circular polarization (normalized to the quiet-Sun continuum intensity) generated as a function of the component of the magnetic field vector that is parallel to the observer’s line-of-sight: B (dashed lines). Amount of linear polarization (normalized to the quiet-Sun continuum intensity) signals generated as a function of the component of the magnetic field vector that is perpendicular to the observer’s line-of-sight: B (solid lines). These curves were obtained by performing a synthesis of the Stokes profiles employing different values of B and B and using two different atmospheric semi-empirical models: red for intergranules and blue for granules. The horizontal black lines indicate the 3σ and 8σ levels assuming a noise of σ = 3 × 10-4 as in our observed maps (Sect. 2).

In the text
thumbnail Fig. 4

Panel a) histogram of the number of pixels with a given maximum value in their polarization signals (normalized to the average continuum intensity in the quiet Sun). The vertical solid-black line corresponds to the 3σ-level. Panel b) histogram of pixels with polarization signals above a certain S/N. This last histogram is cumulative, but the first one is not. The solid and dashed lines refer to the linear and circular polarization profiles, respectively. These two panels can be readily compared to the observed histograms in Fig. 2. They were obtained employing an isotropic distribution function for the magnetic field vector (Eq. (13)) with three different mean values for the magnetic field strength B0: 20 (red), 30 (blue), and 40 (green) Gauss. The module of the magnetic field vector B and the inclination of this vector with respect to the observer’s line-of-sight γ arising from these isotropic distributions are presented in panels c) and d).

In the text
thumbnail Fig. 5

Same as Fig. 4 but for the triple-Gaussian theoretical distribution given by Eq. (19). In this case the mean value of the three components of the magnetic field vector are different: Bx0 = By0 = 22 G, and Bz0 = 5 G. Red lines correspond to Θ = 0° (disk center) while green and blue lines correspond to Θ = 45°. Solid and dashed lines indicate the circular (Stokes V) and linear (Stokes Q and U) polarization.

In the text
thumbnail Fig. 6

Same as Fig. 4 but for the theoretical distribution given by Eq. (20). The parameters employed are β ∥  = 3 G and β = 7 G. Red lines correspond to Λ = 0° (disk center), while green and blue lines correspond to Λ = 30° and Λ = 40°, respectively. Solid and dashed lines indicate the circular (Stokes V) and linear (Stokes Q and U) polarization.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.