The Voigt and complex error function: Humlíček's rational approximation generalized

Abstract

Accurate yet efficient computation of the Voigt and complex error function is a challenge since decades in astrophysics and other areas of physics. Rational approximations have attracted considerable attention and are used in many codes, often in combination with other techniques. The 12-term code "cpf12" of Humlíček (1979) achieves an accuracy of five to six significant digits throughout the entire complex plane. Here we generalize this algorithm to a larger (even) number of terms. The $n=16$ approximation has a relative accuracy better than $10^{-5}$ for almost the entire complex plane except for very small imaginary values of the argument even without the correction term required for the cpf12 algorithm. With 20 terms the accuracy is better than $10^{-6}$. In addition to the accuracy assessment we discuss methods for optimization and propose a combination of the 16-term approximation with the asymptotic approximation of Humlíček (1982) for high efficiency.

9 pages, 5 figures

References in corpus (4)

• Molecfit: A general tool for telluric absorption correction. I. Method and application to ESO instruments
• Recommended isolated-line profile for representing high-resolution spectroscopic transitions (IUPAC Technical Report)
• Voigt Profile Fitting to Quasar Absorption Lines: An Analytic Approximation to the Voigt-Hjerting Function
• Algorithm 985: Simple, efficient, and relatively accurate approximation for the evaluation of the Faddeyeva function

Cited by in corpus (10)

• EXOPLINES: Molecular Absorption Cross-Section Database for Brown Dwarf and Giant Exoplanet Atmospheres
• PyExoCross: a Python program for generating spectra and cross-sections from molecular line lists
• Notes: An assessment of some closed-form expressions for the Voigt function III: Combinations of the Lorentz and Gauss functions
• Optab: Public code for generating gas opacity tables for radiation hydrodynamics simulations
• Computational Aspects of Speed-Dependent Voigt and Rautian Profiles
• SVEEEETIES: Singular vector expansion to estimate Earth-like exoplanet temperatures from infrared emission spectra
• A single-domain implementation of the Voigt/complex error function by vectorized interpolation
• Super-accuracy calculation for the width of a Voigt profile
• Efficient application of the Voigt functions in the Fourier transform
• Assessment of a Physics-based Retrieval of Exoplanet Atmospheric Temperatures from Infrared Emission Spectra