Computing the Fermi−Dirac Functions by Exponentially Convergent Quadratures
- Authors: Kalitkin N.N.1, Kolganov S.A.2
-
Affiliations:
- Keldysh Institute of Applied Mathematics
- National Research University of Electronic Technology
- Issue: Vol 10, No 4 (2018)
- Pages: 472-482
- Section: Article
- URL: https://journals.rcsi.science/2070-0482/article/view/202493
- DOI: https://doi.org/10.1134/S2070048218040063
- ID: 202493
Cite item
Abstract
Highly accurate specialized quadrature formulas are constructed for directly computing the Fermi−Dirac functions of the half-integer index. It is shown that the dependence of the error on the number of nodes is not power-law but exponential. The properties of such formulas are investigated. It is demonstrated that the factor of the convergence exponent is determined by the distance to the nearest pole of the integrand. This provides a very fast convergence of the quadratures. Simple approximations of the Fermi−Dirac functions of the integer and half-integer indices with an accuracy better than 1% are constructed; these approximations are convenient for physical estimates. In passing, asymptotic representations for Bernoulli numbers are found.
About the authors
N. N. Kalitkin
Keldysh Institute of Applied Mathematics
Author for correspondence.
Email: kalitkin@imamod.ru
Russian Federation, Moscow, 125047
S. A. Kolganov
National Research University of Electronic Technology
Email: kalitkin@imamod.ru
Russian Federation, Zelenograd, Moscow, 124498