Veronique Bommier

The Télescope Héliographique pour l’Etude du Magnétisme et des Instabilités Solaires (THEMIS) has been operating in the Canary Islands since 1998. A total of 187 publications are listed in the THEMIS database. The telescope was upgraded in 2019 with adaptive optics and was fully operational in 2024. When operated in polarimetric mode, the telescope is calibration-free and has a high polarimetric sensitivity, which enables important results to be obtained. We will summarize a few of these results, obtained mainly during coordinate campaigns with the multi-spacecraft, outlined as follows: the horizontal magnetic field in prominences, the existence of flux rope in flare regions, and the magnetic field interchange reconnection between jets and filaments.

Context. All magnetic field vector measurements lead to ambiguous results. Aims. We demonstrate that carrying out observations in two different lines belonging to the same multiplet, but with different absorption coefficients (and thus formed at two different depths, e.g., Fe I 6302.5 Å and 6301.5 Å) allow for the azimuth ambiguity remaining from the Zeeman signal interpretation to be resolved. Methods. The measurement based on an interpretation of the Zeeman effect relates to the magnetic field, H, and not the divergence-free magnetic induction, B. We analyzed how the anisotropy of the photosphere, which is strongly stratified due to gravity and density at the star surface, affects the estimation of div H. We also considered how the ambiguity resolution ought to be performed in these conditions. Results. We obtained two ambiguity-resolved field vector maps at two different, but nonetheless close altitudes. This allowed us to derive the current density full vector via curl H = J. The process revealed the horizontal component of the current density, which is generally found to be markedly larger than the better known vertical one. We observed some systematical trends and present some examples in this paper, such as circular currents wrapping spots clockwise about a positive polarity spot and counterclockwise about a negative polarity spot. We also detail the strong horizontal current components crossing active region neutral lines. Conclusions. Finally, we remark that the Na I D1 and D2 lines form another such line pair. We propose them as a way to access the low chromosphere where they are formed. However, a better spatial resolution for both observations and analysis would likely be necessary in such a medium.

In the present paper we consider the full nonlocal thermodynamic equilibrium (non-LTE) radiation transfer problem. This formalism allows us to account for deviation from equilibrium distribution of both the radiation field and the massive particles. In the present study, two-level atoms with broadened upper level represent the massive particles. In the absence of velocity-changing collisions, we demonstrate the analytic equivalence of the full non-LTE source function with the corresponding standard non-LTE partial frequency redistribution (PFR) model. We present an iterative method based on operator splitting techniques that can be used to numerically solve the problem at hand. We benchmark it against the standard non-LTE transfer problem for a two-level atom with PFR. We illustrate the deviation of the velocity distribution function of excited atoms from the equilibrium distribution. We also discuss the dependence of the emission profile and the velocity distribution function on elastic collisions and velocity-changing collisions.

Whitepaper #333 in the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. Main topics: space weather applications; basic research. Additional topics: space weather research/operations/research loop; space weather research/applications/operations pipeline; […]

The problem of bias, meaning over- or under-estimation, of the component perpendicular to the line-of-sight [ $B_{\perp }$ B ⊥ ] in vector magnetic-field maps is discussed. Previous works on this topic have illustrated that the problem exists; here we perform novel investigations to quantify the bias, fully understand its source(s), and provide mitigation strategies. First, we develop quantitative metrics to measure the $B_{\perp }$ B ⊥ bias and quantify the effect in both local (physical) and native image-plane components. Second, we test and evaluate different options available to inversions and different data sources, to systematically characterize the impacts of these choices, including explicitly accounting for the magnetic fill fraction [ $f\!\!f$ f f ]. Third, we deploy a simple model to test how noise and different models of the bias may manifest. From these three investigations we find that while the bias is dominantly present in under-resolved structures, it is also present in strong-field, pixel-filling structures. Noise in the spectropolarimetric data can exacerbate the problem, but it is not the primary cause of the bias. We show that fitting $f\!\!f$ f f explicitly provides significant mitigation, but that other considerations such as the choice of $\chi ^{2}$ χ 2 -weights and optimization algorithms can impact the results as well. Finally, we demonstrate a straightforward “quick fix” that can be applied post facto but prior to solving the $180^{\circ}$ 180 ∘ ambiguity in $B_{\perp }$ B ⊥ , and which may be useful when global-scale structures are, e.g., used for model boundary input. The conclusions of this work support the deployment of inversion codes that explicitly fit $f\!\!f$ f f or, as with the new neural-net, that are trained on data that did so.

The problem of bias, meaning over- or underestimation, of the component perpendicular to the line-of-sight, Bperp, in vector magnetic field maps is discussed. Previous works on this topic have illustrated that the problem exists; here we perform novel investigations to quantify the bias, fully understand its source(s), and provide mitigation strategies. First, we develop quantitative metrics to measure the Bperp bias and quantify the effect in both local (physical) and native image-plane components. Second we test and evaluate different inversion options and data sources, to systematically characterize the impacts of choices, including explicitly accounting for the magnetic fill fraction ff. Third we deploy a simple model to test how noise and different models of the bias may manifest. From these three investigations we find that while the bias is dominantly present in under-resolved structures, it is also present in strong-field pixel-filling structures. Noise in the magnetograms can exacerbate the problem, but it is not the primary cause. We show that fitting ff explicitly provides significant mitigation, but that other considerations such as choice of chi^2 weights and optimization algorithms can impact the results as well. Finally, we demonstrate a straightforward "quick fix" that can be applied post-facto but prior to solving the 180deg ambiguity in Bperp, and which may be useful when global-scale structures are, e.g., used for model boundary input. The conclusions of this work support the deployment of inversion codes that explicitly fit ff or, as with the new SyntHIA neural-net, that are trained on data that did so.

Aims. We present 24 synoptic maps of solar filaments, in which the average unambiguous magnetic field vectors of 296 prominences were determined with Pic-du-Midi observations between 1974 and 1982. This was the ascending phase of cycle 21. Methods. The magnetic field was determined by interpreting the Hanle effect, which is observed in the He I D₃ line. Previous results for the prominence field polarity and prominence chirality were applied to solve the fundamental ambiguity. The measurements were averaged in each prominence for accuracy reasons. Results. The result is twofold. First, alternating field directions can be observed from one neutral line to the next. Second, a general field alignment is found along a solar north-south field that is distorted by the differential rotation effect.

Context. This paper presents a numerical application of a self-consistent theory of partial redistribution in nonlocal thermodynamical equilibrium conditions, developed in previous papers of the series. Aims. The code was described in IV of this series. However, in that previous paper, the numerical results were unrealistic. The present paper presents an approximation able to restore the reliability of the outgoing polarization profiles. Methods. The convergence of the results is also proved. It is demonstrated that the step increment decreases like 1/N^α, with α > 1. Results. Thanks to these additions, the results series behaves like a Riemann series, which is absolutely convergent. However, convergence is not fully reached in line wings within the allocated computing time. Development of efficient acceleration methods would be desirable for future work. Conclusions. Agreement between the computed and observed linear polarization profiles remains qualitative only. The discrepancy is assigned to the plane parallel atmosphere model, which is insufficient to describe the chromosphere, where these lines are formed. As all the integrals are numerical in the code, it could probably be adapted to more realistic and higher dimensional model atmospheres. However, this is time consuming for lines with a hyperfine structure, as in the Na I D lines. The net linear polarization observed in Na I D₁ with the Zürich Imaging Polarimeter ZIMPOL mounted on the McMath-Pierce telescope at Kitt Peak is not confirmed by the present calculations and could be an artefact of instrumental polarization. The presence of instrumental polarization could be confirmed by the higher linear polarization degree observed by this instrument in the Na I D₂ line center with respect to the present calculation result where the magnetic field is not accounted for. At this precise point, the Hanle effect acts as a depolarizing effect in the second solar spectrum. The observed linear polarization excess is found to be of the same order of magnitude in both line centers, namely 0.1%, which is also comparable to the instrumental polarization compensation level of this experiment.

This paper presents a numerical application of a self-consistent theory of partial redistribution in non-LTE conditions, developed in previous papers of the series. The code was described in a previous paper of this series. However, in that previous paper (number IV of the series), the numerical results were unrealistic. The present paper presents an approximation, which was able to restore the reliability of the outgoing polarization profiles. The convergence of the results is also proved. It is demonstrated that the step increment decreases like 1/N^a, with a > 1, so that the results series behaves like a Riemann series, which is absolutely convergent. However, agreement between the computed and observed linear polarization profiles remains qualitative only. The discrepancy is assigned to the plane parallel atmosphere model, which is insufficient to describe the chromosphere, where these lines are formed. As all the integrals are numerical in the code, it could probably be adapted to more realistic and higher dimensioned model atmospheres. However, it is time consuming for lines having an hyperfine structure as the Na I D lines are. The net linear polarization observed in Na I D1 with the polarimeter ZIMPOL mounted on the McMath-Pierce telescope at Kitt Peak is not confirmed by the present calculations and could be an artefact of instrumental polarization. The presence of instrumental polarization could be confirmed by the higher linear polarization degree observed by this instrument in Na I D2 line center, with respect to the present calculation result, where the magnetic field is not accounted for, when the Hanle effect acts as a depolarizing effect in the Second Solar Spectrum. The observed linear polarization excess is found of the same order of magnitude 0.1% in both line centers, which is also comparable to the instrumental polarization compensation level of this experiment.