What polarized hydrogen light can reveal about the Sun’s magnetic atmosphere

A magnetic clue with several possible causes

How can astronomers measure magnetic fields in the Sun’s chromosphere, the structured atmospheric layer above its visible surface? One clue comes from Hα, a spectral line produced by hydrogen. Scattering gives this light a small degree of linear polarization: a preference for particular orientations of its electric field.

Magnetic fields can modify that polarization through the Hanle effect. But motion and uneven illumination also change the signal. An asymmetric polarization profile—a different shape on either side of its peak—therefore cannot automatically be read as evidence for a particular magnetic structure.

Jaume Jaume Bestard, Jiří Štěpán and Javier Trujillo Bueno investigate these overlapping influences in an arXiv preprint posted on November 9, 2026. Its submission record reports acceptance for publication in Astronomy & Astrophysics. The work is a computational study compared with existing solar observations, not a new observing campaign.

Following light through a three-dimensional atmosphere

The researchers used a solar-atmosphere simulation produced with Bifrost, a code that models plasma and magnetic fields together. It includes nonequilibrium hydrogen ionization, allowing the ionization state to depart from an instantaneous equilibrium value as conditions change.

They then used PORTA to calculate how radiation interacts with hydrogen and travels through the atmosphere. Their non-local thermodynamic equilibrium treatment means the atomic excitation cannot be determined from local temperature alone: radiation arriving from elsewhere matters too.

The calculation couples hydrogen transitions responsible for Hα and two ultraviolet lines, Lyα and Lyβ. It includes atomic polarization, magnetic depolarization, horizontal radiation transport and Doppler shifts from plasma motion. The radiation calculation assumes statistical equilibrium for atomic excitation while retaining the simulation’s nonequilibrium ionization state.

For comparison with measurements from the Zurich Imaging Polarimeter, ZIMPOL, the team blurred and sampled the synthetic signals to represent observational limitations, including spatial and spectral resolution and polarimetric accuracy. This matters because fine structure visible in a simulation may disappear when averaged by a telescope.

Broad agreement, with visible shortcomings

The modeled polarization follows the observed overall trend from the center of the solar disk toward its edge, or limb, with broadly comparable amplitudes. The calculations also produce substantial variation across the simulated atmosphere and asymmetric profiles resembling features seen in observations.

Nevertheless, the synthetic intensity profiles are narrower and deeper than observed. The authors associate this discrepancy partly with missing small-scale motions and other limitations of the atmospheric model. Even after observational degradation, synthetic polarization structures generally have larger amplitudes and sharper spatial variations than the measurements.

One important mismatch is time: the calculation uses a single atmospheric snapshot, whereas observations average over several minutes. A changing atmosphere can smooth polarization signals during an exposure. The modeled magnetic environment also differs from that sampled by the coronal-hole observations used for the center-to-limb comparison.

Profile shapes track magnetic differences statistically

The team grouped synthetic polarization profiles into symmetric single peaks, asymmetric single peaks and multiple peaks. The quantity analyzed, Q/I, expresses one component of linear polarization relative to total intensity.

The multiple-peaked group preferentially samples stronger magnetic fields. At the reference optical depth used to characterize line formation, its mean field is approximately 22.5 gauss, compared with about 18 gauss for symmetric single peaks and 16 gauss for asymmetric single peaks.

The amplitude distributions differ too. Both single-peaked groups have distributions peaking around 0.25% polarization, while the multiple-peaked group peaks near 0.05%. These are typical locations within distributions of simulated signals, not universal thresholds for interpreting solar measurements.

The association of stronger fields with weaker polarization is consistent with Hanle depolarization. However, the classification excludes profiles with substantial negative polarization, and the fraction labeled asymmetric depends strongly on the chosen threshold.

A diagnostic possibility, not a magnetic shortcut

The study establishes correlations within this model, not a unique conversion from profile shape to magnetic strength. Nearly symmetric profiles still occur at fields capable of producing a Hanle effect.

Crucially, the researchers did not repeat the calculation with the magnetic field removed while preserving the same temperature, density and velocities. They therefore cannot cleanly separate magnetic contributions from atmospheric motion and three-dimensional structure. Hα remains a promising probe, but its interpretation requires those influences to be considered together.