Comparison of inversion codes for polarized line formation in MHD simulations

I. Milne-Eddington codes

¹ Kiepenheuer Institut für Sonnenphysik, Schöneckstr. 6, 79104 Freiburg, Germany
e-mail: borrero@kis.uni-freiburg.de; rrezaei@kis.uni-freiburg.de
² High Altitude Observatory (NCAR), 3090 Center Green Dr., Boulder, CO 80301, USA
e-mail: lites@hao.ucar.edu; empel@hao.ucar.edu
³ Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, 37077 Göttingen, Germany
e-mail: lagg@mps.mpg.de

Received: 10 July 2014
Accepted: 29 August 2014

Abstract

Milne-Eddington (M-E) inversion codes for the radiative transfer equation are the most widely used tools to infer the magnetic field from observations of the polarization signals in photospheric and chromospheric spectral lines. Unfortunately, a comprehensive comparison between the different M-E codes available to the solar physics community is still missing, and so is a physical interpretation of their inferences. In this contribution we offer a comparison between three of those codes (VFISV, ASP/HAO, and HeLIx⁺). These codes are used to invert synthetic Stokes profiles that were previously obtained from realistic non-grey three-dimensional magnetohydrodynamical (3D MHD) simulations. The results of the inversion are compared with each other and with those from the MHD simulations. In the first case, the M-E codes retrieve values for the magnetic field strength, inclination and line-of-sight velocity that agree with each other within σ_B ≤ 35 (Gauss), σ_γ ≤ 1.2°, and σ_v ≤ 10 m s^-1, respectively. Additionally, M-E inversion codes agree with the numerical simulations, when compared at a fixed optical depth, within σ_B ≤ 130 (Gauss), σ_γ ≤ 5°, and σ_v ≤ 320 m s^-1. Finally, we show that employing generalized response functions to determine the height at which M-E codes measure physical parameters is more meaningful than comparing at a fixed geometrical height or optical depth. In this case the differences between M-E inferences and the 3D MHD simulations decrease to σ_B ≤ 90 (Gauss), σ_γ ≤ 3°, and σ_v ≤ 90 m s^-1.