An adaptive gaussian quadrature for the Voigt function

TL;DR

This paper introduces an adaptive Gaussian quadrature method tailored for efficiently computing the Voigt function, crucial for advanced radiation transfer models involving non-Maxwellian particle velocity distributions.

Contribution

It presents a novel adaptive quadrature scheme specifically designed for the numerical evaluation of the Voigt function in complex radiation transfer calculations.

Findings

Validated the method with standard Voigt profile computations

Demonstrated improved efficiency for non-LTE radiation transfer models

Suitable for generalized redistribution in frequency functions

Abstract

We evaluate an adaptive gaussian quadrature integration scheme that will be suitable for the numerical evaluation of generalized redistribution in frequency functions. The latter are indispensable ingredients for "full non-LTE" radiation transfer computations i.e., assuming potential deviations of the velocity distribution of massive particles from the usual Maxwell-Boltzmann distribution. A first validation is made with computations of the usual Voigt profile.