Polylogarithm

The polylogarithm (also known as Jonquiére's function) is a special function and may be defined for all s and |z|<1 by:

$Li_{s}(z)\equiv \sum _{k=1}^{\infty }{z^{k} \over k^{s}}$

Both the parameter s and the argument z are taken to be complex numbers. The special cases s=2 and s=3 are called the dilogarithm and trilogarithm respectively. The polylogarithm also arises in the closed form of the integral of the Fermi-Dirac distribution and the Bose-Einstein distribution and is sometimes known as the Fermi-Dirac integral or the Bose-Einstein integral. Polylogarithms should not be confused with polylogarithmic functions nor with the offset logarithmic integral which has a similar notation.

The polylogarithm is actually defined over a larger range of z than the above definition allows by the process of analytic continuation.

Properties

In the important case where the parameter s is an integer, it will be represented by n (or -n when negative). It is often convenient to define μ=Ln(z) where Ln is the principal branch of the natural logarithm so that -π < Im(μ) ≤ π. Also, all exponentiation will be assumed to be single valued. (e.g. z^s=Exp(s Ln(z))).

Depending on the parameter s, the polylogarithm may be muli-valued. The principal branch of the polylogarithm is chosen to be that for which Li_s(z) is real for z real, 0 ≤ z ≤ 1 and is continuous except on the positive real axis, where a cut is made from z=1 to ∞ such that the cut puts the real axis on the lower half plane of z. In terms of μ, this amounts to -π < arg(-μ) ≤ π. The fact that the polylogarithm may be discontinuous in μ can cause some confusion.

For z real and z ≥ 1 the imaginary part of the polylogarithm is (Wood):

${\textrm {Im}}(Li_{s}(z))=-{{\pi \mu ^{s-1}} \over {\Gamma (s)}}$

Going across the cut, if δ is an infinitesimally small positive real number, then:

${\textrm {Im}}(Li_{s}(z+i\delta ))={{\pi \mu ^{s-1}} \over {\Gamma (s)}}$

The derivatives of the polylogarithm are:

$z{\partial Li_{s}(z) \over \partial z}=Li_{s-1}(z)$

${\partial Li_{s}(e^{\mu }) \over \partial \mu }=Li_{s-1}(e^{\mu })$

Particular values

See also "#Relationship to other functions" section below.

For integer values of s, we have the following explicit expressions:

$Li_{1}(z)=-{\textrm {ln}}\left(1-z\right)$

$Li_{0}(z)={z \over 1-z}$

$Li_{-1}(z)={z \over (1-z)^{2}}$

$Li_{-2}(z)={z(1+z) \over (1-z)^{3}}$

$Li_{-3}(z)={z(1+4z+z^{2}) \over (1-z)^{4}}$

The polylogarithm for all negative integer values of s can be expressed as a ratio of polynomials in z (See series representations below). Some particular expressions for half-integer values of the argument are:

$Li_{1}\left(1/2\right)={\textrm {ln}}(2)$

$Li_{2}(1/2)={1 \over 12}[\pi ^{2}-6(\ln 2)^{2}]$

$Li_{3}(1/2)={1 \over 24}[4(\ln 2)^{3}-2\pi ^{2}\ln 2+21\,\zeta (3)]$

where ζ is the Riemann zeta function. No similar formulas of this type are known for higher orders (Lewin, 1991 p2).

Alternate Expressions

• The Integral of the Bose-Einstein distribution is expressed in terms of a polylogarithm:

  $Li_{s+1}(z)\equiv {1 \over \Gamma (s+1)}\int _{0}^{\infty }{t^{s} \over e^{t}/z-1}dt$

  This converges for Re(s) > 0 and all z except for z real and ≥ 1. The polylogarithm in this context is sometimes referred to as a Bose integral or a Bose-Einstein integral.
• The integral of the Fermi-Dirac distribution is also expressed in terms of a polylogarithm:

  $-Li_{s+1}(-z)\equiv {1 \over \Gamma (s+1)}\int _{0}^{\infty }{t^{s} \over e^{t}/z+1}dt.$

  This converges for Re(s)>0 and all z except for z real and <(-1). The polylogarithm in this context is sometimes referred to as a Fermi integral or a Fermi-Dirac integral. (GNU)
• The polylogarithm may be rather generally represented by a Hankel contour integral (Whittaker & Watson § 12.22, § 13.13). As long as the t=μ pole of the integrand does not lie on the non-negative real axis, and s ≠ 1,2,3..., we have:

  $Li_{s}(e^{\mu })={{-\Gamma (1-s)} \over {2\pi i}}\oint _{H}{{(-t)^{s-1}} \over {e^{t-\mu }-1}}dt$

  Where H represents the Hankel contour. The integrand has a cut along the real axis from zero to infinity, with the real axis being on the lower half of the sheet (Im(t) ≤ 0). For the case where μ is real and non-negative, we can simply add the limiting contribution of the pole:

  $Li_{s}(e^{\mu })=-{{\Gamma (1-s)} \over {2\pi i}}\oint _{H}{{(-t)^{s-1}} \over {e^{t-\mu }}-1}dt+2\pi iR$

  where R is the residue of the pole:

  $R={{\Gamma (1-s)(-\mu )^{s-1}} \over {2\pi }}$

• The square relationship is easily seen from the defining equation: (See also Clunie, Schrödinger):

  $Li_{s}(-z)+Li_{s}(z)=2^{1-s}~Li_{s}(z^{2})$

  Note that Kummer's function obeys a very similar duplication formula.

Relationship to other functions

• For z=1 the polylogarithm reduces to the Riemann zeta function

  $Li_{s}(1)=\zeta (s)~~~~~~~~~~~~~({\textrm {Re}}(s)>1)$

• The polylogarithm is related to Dirichlet eta function and the Dirichlet beta function:

  $Li_{s}(-1)=\eta \left(s\right)$

  where η(s) is the Dirichlet eta function. For pure imaginary arguments, we have:

  $Li_{s}(\pm i)=2^{-s}\eta (s)\pm i\beta (s)$

  where β(s) is the Dirichlet beta function.
• The polylogarithm is equivalent to the Fermi-Dirac integral (GNU)

  $F_{s}(\mu )=-Li_{s+1}(-e^{\mu })\,$

• The polylogarithm is a special case of the Lerch Transcendent (Erdélyi § 1.11-14)

  $Li_{s}(z)=z~\Phi (z,s,1)$

• The polylogarithm is related to the Hurwitz zeta function by:

  $Li_{s}(e^{2\pi ix})+(-1)^{s}Li_{s}(e^{-2\pi ix})={(2\pi i)^{s} \over \Gamma (s)}~\zeta \left(1-s,x\right)$

  where Γ(s) is the gamma function. This holds for

  ${\textrm {Re}}(s)>1,{\textrm {Im}}(x)\geq 0,0\leq {\textrm {Re}}(x)<1$

  and also for

  ${\textrm {Re}}(s)>1,{\textrm {Im}}(x)\leq 0,0<{\textrm {Re}}(x)\leq 1$

  (Note that Erdélyi's equivalent Equation § 1.11-16 is not correct if we assume that the principal branches of the polylogarithm and the logarithm are used simultaneously.) This equation furnishes the analytical continuation of the series representation of the polylogarithm beyond its circle of convergence |z|=1.
• Using the relationship between the Hurwitz zeta function and the Bernoulli polynomials:

  $\zeta (-n,x)=-{B_{n+1}(x) \over n+1}$

  which holds for all x and n=0,1,2,3,... it can be seen that:

  $Li_{n}(e^{2\pi ix})+(-1)^{n}Li_{n}(e^{-2\pi ix})=-{(2\pi i)^{n} \over n!}B_{n}\left({x}\right)$

  under the same constraints on s and x as above. (Note that the corresponding equation Erdélyi § 1.11-18 is not correct) For negative integer values of the parameter, we have for all z (Erdélyi § 1.11-17):

  $Li_{-n}(z)+(-1)^{n}Li_{-n}\left(1/z\right)=0~~~~~n=1,2,3\ldots$

• The polylogarithm with pure imaginary μ may be expressed in terms of Clausen functions Ci_s(θ) and Si_s(θ) (Lewin, 1958 Ch VII § 1.4, Abramowitz & Stegun § 27.8)

  $Li_{s}(e^{\pm i\theta })=Ci_{s}(\theta )\pm iSi_{s}(\theta )$

• The Inverse Tangent Integral Ti_s(z) (Lewin, 1958 Ch VII § 1.2) can be expressed in terms of polylogarithms:

  $Li_{s}(\pm iy)=2^{-s}Li_{s}(-y^{2})\pm i\,Ti_{s}(y)$

• The Legendre chi function χ_s(z) (Lewin, 1958 Ch VII § 1.1, Boersma) can be expressed in terms of polylogarithms:

  $\chi _{s}(z)={1 \over 2}~[Li_{s}(z)-Li_{s}(-z)]$

• The polylogarithm may be expressed as a series of Debye Functions Z_n(z) (Abramowitz & Stegun § 27.1)

  $Li_{n}(e^{\mu })=\sum _{k=0}^{n-1}Z_{n-k}(-\mu ){\mu ^{k} \over k!}~~~~~~(n=1,2,3,\ldots )$

  A remarkably similar expression relates the Debye function to the polylogarithm:

  $Z_{n}(\mu )=\sum _{k=0}^{n-1}Li_{n-k}(e^{-\mu }){\mu ^{k} \over k!}~~~~~~(n=1,2,3,\ldots )$

Series representations

• We may represent the polylogarithm as a power series about μ=0 as follows: (Robinson) Consider the Mellin transform:

  $M_{s}(r)=\int _{0}^{\infty }{\textrm {Li}}_{s}(fe^{-u})u^{r-1}\,du={1 \over \Gamma (s)}\int _{0}^{\infty }\int _{0}^{\infty }{t^{s-1}u^{r-1} \over e^{t+u}/f-1}~dt~du$

  The change of variables t=ab, u=a(1-b) allows the integrals to be separated:

  $M_{s}(r)={1 \over \Gamma (s)}\int _{0}^{1}b^{r-1}(1-b)^{s-1}db\int _{0}^{\infty }{a^{s+r-1} \over e^{a}/f-1}da=\Gamma (r){\textrm {Li}}_{s+r}(f)$

  for f=1 we have, through the inverse Mellin transform:

  $Li_{s}(e^{-u})={1 \over 2\pi i}\int _{c-i\infty }^{c+i\infty }\Gamma (r)\zeta (s+r)u^{-r}dr$

  where c is a constant to the right of the poles of the integrand. The path of integration may be converted into a closed contour, and the poles of the integrand are those of Γ(r) at r=0,-1,-2,..., and of ζ(s+r) at r=1-s. Summing the residues gives, for $|\mu |<2\pi$ and s ≠ 1,2,3,...

  $Li_{s}(e^{\mu })=\Gamma (1-s)(-\mu )^{s-1}+\sum _{k=0}^{\infty }{\zeta (s-k) \over k!}~\mu ^{k}$

  If the parameter s is a positive integer, n, both the k=n-1 term and the gamma function become infinite, although their sum does not. For integer k>0 we have:

  $\lim _{s\rightarrow k+1}\left[{\zeta (s-k)\mu ^{k} \over k!}+\Gamma (1-s)(-\mu )^{s-1}\right]={\mu ^{k} \over k!}\left(\sum _{m=1}^{k}{1 \over m}-{\textrm {Ln}}(-\mu )\right)$

  and for k=0:

  $\lim _{s\rightarrow 1}\left[\zeta (s)+\Gamma (1-s)(-\mu )^{s-1}\right]=-{\textrm {Ln}}(-\mu )$

  So, for s=n where n is a positive integer and $|\mu |<2\pi$ we have the following:

  $Li_{n}(e^{\mu })={\mu ^{n-1} \over (n-1)!}\left(H_{n}-{\textrm {Ln}}(-\mu )\right)+$

  $\sum _{k=0,k\neq n-1}^{\infty }{\zeta (n-k) \over k!}~\mu ^{k}~~~~~~~~~~~~~~~~~~~~~~(n=2,3,4,\ldots )$

  $Li_{1}(e^{\mu })=-{\textrm {Ln}}(-\mu )+\sum _{k=1}^{\infty }{\zeta (1-k) \over k!}~\mu ^{k}~~~~~~~~~~(n=1)$

  Where H_n is a harmonic number:

  $H_{n}\equiv \sum _{k=1}^{n}{1 \over k}$

  The problem terms now contain -Ln(-μ) which, when multiplied by μ^k will tend to zero as μ tends to zero, except for k=0. This reflects the fact that there is a true logarithmic singularity in Li_s(z) at s=1 and z=1 Since:

  $\lim _{\mu \rightarrow 0}\Gamma (1-s)(-\mu )^{s-1}=0~~~~~({\textrm {Re}}(s)>1)$

  Using the relationship between the Riemann zeta function and the Bernoulli numbers B_k:

  $\zeta (-n)=(-1)^{n}{B_{n+1} \over n+1}~~~~~~~~~~~(n=0,1,2,3,\ldots )$

  we obtain for negative integer values of s and $|\mu |<2\pi$ :

  $Li_{-n}(z)={n! \over (-\mu )^{n+1}}-\sum _{k=0}^{\infty }{B_{k+n+1} \over k!~(k+n+1)}~\mu ^{k}~~~~~~~~~~~(n=1,2,3,\ldots )$

  since, except for B₁, all odd Bernoulli numbers are zero. We obtain the n=0 term using ζ(0)=B₁=-1/2. Note again that Erdélyi's equivalent Equation § 1.11-15 is not correct if we assume that the principal branches of the polylogarithm and the logarithm are used simultaneously, since Ln(1/z) is not uniformly equal to -Ln(z).
• The defining equation may be extended to negative values of the parameter s using a Hankel contour integral (Wood, Gradshteyn & Ryzhik § 9.553):

  $Li_{s}(e^{\mu })=-{\Gamma (1-p) \over 2\pi i}\oint _{H}{(-t)^{s-1} \over e^{t-\mu }-1}dt$

  where H is the Hankel contour which can be modified so that it encloses the poles of the integrand, at $t-\mu =2k\pi i$ and the integral can be evaluated as the sum of the residues:

  $Li_{s}(e^{\mu })=\Gamma (1-s)\sum _{k=-\infty }^{\infty }(2k\pi i-\mu )^{s-1}$

  This will hold for ${\textrm {Re}}(s)<0$ and all z except z=1.
• For negative integer s, the polylogarithm may be expressed as a series involving the Eulerian numbers

  $Li_{-n}(z)={1 \over (1-z)^{n+1}}\sum _{i=0}^{n-1}\left\langle {n \atop i}\right\rangle z^{n-i}~~~~~~~~~~~~~(n=1,2,3,\ldots )$

  where $\left\langle {n \atop i}\right\rangle$ are Eulerian numbers:
• Another explicit formula for negative integer s is (Wood):

  $Li_{-n}(z)=\sum _{k=1}^{n+1}{(-1)^{n+k+1}(k-1)!S(n+1,k) \over (1-z)^{k}}~~~~~~~~~~(n=1,2,3,\ldots )$

  where S(n,k) are Stirling numbers of the second kind.

Limiting behavior

The following limits hold for the polylogarithm (Wood):

$\lim _{|z|\rightarrow 0}Li_{s}(z)=\lim _{s\rightarrow \infty }Li_{s}(z)=z$

$\lim _{\mathrm {Re} (\mu )\rightarrow \infty }Li_{s}(e^{\mu })=-{\mu ^{s} \over \Gamma (s+1)}~~~~~~(s\neq -1,-2,-3,\ldots )$

$\lim _{\mathrm {Re} (\mu )\rightarrow \infty }Li_{n}(e^{\mu })=-(-1)^{n}e^{-\mu }~~~~~~(n=1,2,3,\ldots )$

$\lim _{|\mu |\rightarrow 0}Li_{s}(e^{\mu })=\Gamma (1-s)(-\mu )^{s}~~~~~~(s<1)$

Polylogarithm ladders

Leonard Lewin discovered a remarkable and broad generalization of a number of classical relationships on the polylogarithm for special values. These are now called polylogarithm ladders. Define $\rho =\left({\sqrt {5}}-1\right)/2$ as the reciprocal of the golden ratio. Then two simple examples of results from ladders include

$Li_{2}(\rho ^{6})=4Li_{2}(\rho ^{3})+3Li_{2}(\rho ^{2})-6Li_{2}(\rho )+{\frac {7\pi ^{2}}{30}}$

given by Coexeter in 1935, and

$Li_{2}(\rho )={\frac {\pi ^{2}}{10}}-\log ^{2}\rho$

given by Landen. Polylogarithm ladders occur naturally and deeply in K-theory.

History

Don Zagier remarked that "The dilogarithm is the only mathematical function with a sense of humor."

References

• Abramowitz, M. and Stegun, I. A. (ed.), Handbook of Mathematical Functions, National Bureau of Standards, 1964; reprinted Dover Publications, 1965.
• Bailey, D.; Borwein, P.; and Plouffe, S. "On the Rapid Computation of Various Polylogarithmic Constants." http://www.cecm.sfu.ca/~pborwein/PAPERS/P123.ps.
• Bailey, D. H. and Broadhurst, D. J. "A Seventeenth-Order Polylogarithm Ladder." 20 Jun 1999. http://xxx.lanl.gov/abs/math.CA/9906134/
• Boersma, J. and Dempsey, J. P. "On the evaluation of Legendre's chi-function", Mathematics of Computation, 59, 199, pp. 157-163, 1992.
• Borwein, J. M.; Bradley, D. M.; Broadhurst, D. J.; and Losinek, P. "Special Values of Multidimensional Polylogarithms." CECM-98:106, 14 May. http://www.cecm.sfu.ca/preprints/1998pp.html\#98:106
• Borwein, J. M.; Bradley, D. M.; Broadhurst, D. J.; and Losinek, P. "Special Values of Multidimensional Polylogarithms." 8 Oct 1999. http://xxx.lanl.gov/abs/math.CA/9910045/
• Berndt, B. C. Ramanujan's Notebooks, Part IV. New York: Springer-Verlag, pp. 323-326, 1994.
• Clunie, J., "On Bose-Einstein functions", Proceedings of the Physical Society, Section A, 67, pp. 632-636, 1954.
• Erdélyi, A.; Magnus, W.; Oberhettinger, F.; and Tricomi, F. G. Higher Transcendental Functions, Vol. 1. New York: Krieger, 1981.
• Fornberg, B. and Kölbig, K. S., "Complex zeros of the Jonquiére or polylogarithm function", Mathematics of Computation, 29, 130, pp. 582-599, 1975.
• Gradshteyn, I.S. and Ryzhik, I.M., Tables of Integrals, Series, and Products, Academic Press, New York 1980.
• GNU Scientific Library - Reference Manual http://www.gnu.org/software/gsl/manual/gsl-ref.html#SEC117
• Jahnke, E. and Emde, F., Tables of Functions with Formulae and Curves, Dover, 1945.
• Kölbig, K. S., Mignaco, J. A., and Remiddi, E., "On Nielsen's generalized polylogarithms and their numerical calculation", BIT, 10, pp. 38-74, 1970.
• Lewin, L. Dilogarithms and Associated Functions. London: Macdonald, 1958.
• Lewin, L. Polylogarithms and Associated Functions. New York: North-Holland, 1981.
• Lewin, Leonard. (Ed.). Structural Properties of Polylogarithms. Providence, RI: Amer. Math. Soc., 1991. ISBN 0-8218-4532-2
• McDougall, J. and Stoner, E. C. "The computation of Fermi-Dirac functions", Philosophical Transactions of the Royal Society, Series A, 237, pp. 67-104, 1939.
• Markman, B., "The Riemann Zeta Function", BIT,5, pp. 138-141, 1965.
• Nielsen, N. Der Euler'sche Dilogarithms. Leipzig, Germany: Halle, 1909.
• Prudnikov, A. P.; Marichev, O. I.; and Brychkov, Yu. A. "The Generalized Zeta Function, Bernoulli Polynomials, Euler Polynomials, and Polylogarithms ." §1.2 in Integrals and Series, Vol. 3: More Special Functions. Newark, NJ: Gordon and Breach, pp. 23-24, 1990.
• Robinson, J. E. "Note on the Bose-Einstein integral functions", Physical Review, Series 2, 83, pp. 678-679, 1951. http://prola.aps.org/abstract/PR/v83/i3/p678_1
• Schrödinger, E., Statistical Thermodynamics, Cambridge, 1952.
• Truesdell, C. "On a function which occurs in the theory of the structure of polymers", Annals of Mathematics, Series 2, 46, No 1, pp. 144-1457, 1945.
• Whittaker, E. T., and Watson, G. N., A Course of Modern Analysis, Cambridge, 1927.
• Wood, David C., Technical Report 15-92, University of Kent computing Laboratory, University of Kent, Canterbury, UK June, 1992. http://www.cs.kent.ac.uk/pubs/1992/110/
• Zagier, D. "Special Values and Functional Equations of Polylogarithms." Appendix A in Structural Properties of Polylogarithms (Ed. L. Lewin). Providence, RI: Amer. Math. Soc., 1991.