Fermi-Dirac integrals appear frequently in semiconductor problems, so an understanding of their properties is essential. The purpose of these notes is to collect in one place, some basic information about Fermi-Dirac integrals and their properties.
We also present Matlab functions (in a zipped file) that calculate Fermi-Dirac integrals (the "script F" defined by Blakemore (1982)) in three different ways.
The function, "FD_int_approx.m", evaluates Fermi-Dirac integrals using analytic approximations developed by Bednarczyk et al. (1978) and Aymerich-Humet et al. (1981, 1983).
The folder, "FD_int_Pulfrey", includes a set of Matlab files that calculate Fermi-Dirac integrals using the approximations proposed by Halen and Pulfrey (1985, 1986). The main function is "FDjx.m". This function gives a better accuracy, but the simulation time is a little longer.
Fermi-Dirac integrals can be evaluated accurately by numerical integration. The function, "FD_int_num.m", calculates Fermi-Dirac integrals using the composite trapezoidal rule. This approach provides very high accuracy, but the CPU time is considerably longer.
J. S. Blakemore, Solid-St. Electron, 25, 1067 (1982)
X. Aymerich-Humet, F. Serra-Mestres, and J. Millan, Solid-St. Electron, 24, 981 (1981)
X. Aymerich-Humet, F. Serra-Mestres, and J. Millan, J. Appl. Phys., 54, 2850 (1983)
P. V. Halen and D. L. Pulfrey, J. Appl. Phys., 57, 5271 (1985)
P. Van Halen and D. L. Pulfrey, J. Appl. Phys., 59, 2264 (1986)
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipies: The art of scientific computing, 3rd Ed., Cambridge University Press, 2007.