Inversion codes for the Radiative Transfer equation
At the Leibniz Institut für Sonnenpysik we have a strong tradition in the development and maintenance of inversion codes for the radiative transfer equation that are widely used to analyse spectropolarimetric data (i.e. Stokes vector) in order to infer the physical properties of the solar atmosphere (temperature, magnetic field, plasma velocity, etc). We are currently maintaining the VFISV inversion code as well as actively developing the FIRTEZ-DZ inversion code.
Very Fast Inversion for the Stokes Vector (VFISV)
This inversion code was originally developed between 2004-2010 in order to analyse data from the Helioseismic and Magnetic Imager (HMI) instrument on-board the Solar Dynamics observatory (SDO). HMI is a Michelson-type instrument with Lyot elements to perform the wavelength tunning. At a later stage, the code was adapted to analyse data from spectrograph-type (i.e. spectropolarimeters) instruments as well as for Fabry-Perot-type instruments. The code is also being employed for quick-look analysis of GRIS data at KIS Science Data Centre.
VFISV focuses on speed and therefore it is a Milne-Eddington (M-E) inversion code. M-E codes retrieve the average kinetic and magnetic properties of the solar atmosphere over the region in which the analysed spectral lines are formed. More information about VFISV can be found here. The code is currenly being maintained by Dr. J. M. Borrero. This inversion code is free and open-software licensed under the GPL2 licence. There are currently two different versions available:
- Version for Fabry-Perot instruments: https://gitlab.leibniz-kis.de/borrero/vfisv_fpi
- Version for Spectrograph instruments: https://gitlab.leibniz-kis.de/borrero/vfisv_spec
- Manual: https://gitlab.leibniz-kis.de/borrero/vfisv_fpi/-/blob/master/manual_vfisv.pdf
If you use VFISV in your work/research, we kindly ask that you cite the original paper: Borrero et al. (2011), Sol.Phys, 273, 267
Forward and Inverse solver for the Radiative Transfer in geometrical height "Z".
The FIRTEZ-DZ inversion code is capable of retrieving, not the average kinematic and magnetic properties like M-E inversion codes do, but rather the full vertical stratification of these physical properties. FIRTEZ-DZ is the only inversion code capable of achieving this in the geometrical scale z instead of the optical depth scale. This is done by combining the radiative transfer equation for polarized light with the magnetohydrostatic momentum equation in ideal magnetohydrodynamics (MHD). The code is under heavily development and currently maintained by Dr. Adur Pastor Yabar and Dr. J. M. Borrero. Original funding for its development was provided by the DFG.
FIRTEZ-DZ is publicly available as free and open software under the GPL2 licence. It can be downloaded/cloned here:
- Inversion code source: https://github.com/apy-github/firtez-dz
- Manual: https://github.com/apy-github/firtez-dz/blob/master/manual/code_usermanual.pdf
- MHS module source code: https://gitlab.leibniz-kis.de/borrero/firtez-mhs
Note: The MHS module is not yet implemented inside of FIRTEZ-DZ so they need to be used separately and iteratively. Please contact us if you want to do this. In the future we hope to have the MHS module integrated inside FIRTEZ.
If you use FIRTEZ-DZ in your work/research, we kindly ask that you cite the original paper: Pastor Yabar et al. (2019), A&A, 629, 24. If you use the MHS module as well, please cite Borrero et al. (2021), A&A, 674, 190.
Instrumental Polarisation Analysis Tool
Polarimetric data acquired from astronomical instruments require calibration processes to compute the physical parameters accurately. A key aspect of such calibration is the estimation of instrumental polarisation (IP). In Stokes polarimetry, instrumental polarization, that is the polarisation cross-talks introduced by the observing optical system, can be represented by a 4×4 Mueller matrix. The usual approach to compute IP is to device telescope specific models by means of analytical approach or ray tracing software. PyAstroPol is a targeted Python library developed to simplify the modelling process and enable smoother integration of models into the data processing pipelines.
The tool is developed and maintained by Dr. Hemanth Pruthvi. This work was part of the HELLRIDE development (Project: Jets in the Solar Atmosphere). The source code is hosted at Github and the corresponding article has been published in the JOSS.