A highly efficient Voigt program for line profile computation

Mofreh R. Zaghloul, Jacques Le Bourlot

Research output: Contribution to journalArticlepeer-review

Abstract

Evaluation of the Voigt function, a convolution of a Lorentzian and a Gaussian profile, is essential in various fields such as spectroscopy, atmospheric science, and astrophysics. Efficient computation of the function is crucial, especially in applications where the function may be called for an enormous number of times. In this paper, we present a highly efficient novel algorithm and its Fortran90 implementation for the practical evaluation of the Voigt function with accuracy in the order of 10−6. The algorithm uses improved fits based on Chebyshev subinterval polynomial approximation for functions in two variables. The algorithm significantly outperforms widely-used competitive algorithms in the literature, in terms of computational speed, making it highly suitable for real-time applications and large-scale data processing tasks. The substantial improvement in efficiency positions the present algorithm and computer code as a valuable tool in relevant scientific domains. The algorithm has been adopted and implemented in the Meudon PDR code at Paris Observatory and is recommended for similar applications and simulation packages.

Original languageEnglish
Article number109234
JournalJournal of Quantitative Spectroscopy and Radiative Transfer
Volume330
DOIs
Publication statusPublished - Jan 2025

Keywords

  • Bivariate Chebyshev polynomial approximation
  • Fortran90 implementation
  • Numerical methods
  • Superior efficiency
  • Voigt spectral line profile

ASJC Scopus subject areas

  • Radiation
  • Atomic and Molecular Physics, and Optics
  • Spectroscopy

Fingerprint

Dive into the research topics of 'A highly efficient Voigt program for line profile computation'. Together they form a unique fingerprint.

Cite this