Voigt profile fitting to quasar absorption lines: an analytic approximation to the Voigt–Hjerting function

The Voigt–Hjerting function is fundamental in order to correctly model the profiles of absorption lines imprinted on the spectra of bright background sources by intervening absorbing systems. In this work, we present a simple analytic approximation to this function in the context of absorption-line profiles of intergalactic H  i absorbers. Using basic calculus tools, we derive an analytic expression for the Voigt–Hjerting function that contains only fourth-order polynomial and Gaussian functions. In connection with the absorption coefficient of intergalactic neutral hydrogen, this approximation is suitable for modelling Voigt profiles with an accuracy of 10⁻⁴ or better for an arbitrary wavelength baseline, for column densities up to   N _H I= 10 ²² cm⁻²  , and for damping parameters   a ≲ 10⁻⁴  , that is, the entire range of parameters characteristic to all Lyman transitions arising in a variety of H  i absorbing systems such as Lyman α (Lyα) forest clouds, Lyman limit systems and damped Lyα systems. We hence present an approximation to the Voigt–Hjerting function that is both accurate and flexible to implement in various types of programming languages and machines, and with which Voigt profiles can be calculated in a reliable and very simple manner.