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:
- <math>
Li_s(z) \equiv \sum_{k=1}^\infty {z^k \over k^s}. <math>
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.
Contents |
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. zs=Exp(s Ln(z))).
Depending on the parameter s, the polylogarithm may be multi-valued. The principal branch of the polylogarithm is chosen to be that for which Lis(z) is real for z real, ≤ 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):
- <math>\textrm{Im}(Li_s(z)) = -{{\pi \mu^{s-1}}\over{\Gamma(s)}}.<math>
Going across the cut, if δ is an infinitesimally small positive real number, then:
- <math>\textrm{Im}(Li_s(z+i\delta)) = {{\pi \mu^{s-1}}\over{\Gamma(s)}}.<math>
The derivatives of the polylogarithm are:
- <math>z{\partial Li_s(z) \over \partial z} = Li_{s-1}(z)<math>
- <math>{\partial Li_s(e^\mu) \over \partial \mu} = Li_{s-1}(e^\mu).<math>
Particular values
See also "#Relationship to other functions" section below.
For integer values of s, we have the following explicit expressions:
- <math>Li_{1}(z) = -\textrm{ln}\left(1-z\right)<math>
- <math>Li_{0}(z) = {z \over 1-z}<math>
- <math>Li_{-1}(z) = {z \over (1-z)^2}<math>
- <math>Li_{-2}(z) = {z(1+z) \over (1-z)^3}<math>
- <math>Li_{-3}(z) = {z(1+4z+z^2) \over (1-z)^4}.<math>
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:
- <math>Li_{1}\left(1/2\right) = \textrm{ln}(2)<math>
- <math>Li_{2}(1/2) = {1 \over 12}[\pi^2-6(\ln 2)^2]<math>
- <math>Li_{3}(1/2) = {1 \over 24}[4(\ln 2)^3-2\pi^2\ln 2+21\,\zeta(3)]<math>
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:
- <math>
Li_{s+1}(z) \equiv {1 \over \Gamma(s+1)} \int_0^\infty {t^s \over e^t/z-1} dt. <math>
- 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:
- <math>
-Li_{s+1}(-z) \equiv {1 \over \Gamma(s+1)} \int_0^\infty {t^s \over e^t/z+1} dt. <math>
- 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.
- 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:
- <math>
Li_s(e^\mu)={{-\Gamma(1-s)}\over{2\pi i}}\oint_H {{(-t)^{s-1}}\over{e^{t-\mu}-1}}dt. <math>
- 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:
- <math>
Li_s(e^\mu)=-{{\Gamma(1-s)}\over{2\pi i}}\oint_H {{(-t)^{s-1}}\over{e^{t-\mu}}-1}dt + 2\pi i R <math>
- where R is the residue of the pole:
- <math>
R = {{\Gamma(1-s)(-\mu)^{s-1}}\over{2\pi}}. <math>
- The square relationship is easily seen from the defining equation (see also Clunie, Schrödinger):
- <math>
Li_s(-z) + Li_s(z) = 2^{1-s} ~ Li_s(z^2). <math>
- 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
- <math>Li_s(1) = \zeta(s)~~~~~~~~~~~~~(\textrm{Re}(s)>1). <math>
- The polylogarithm is related to Dirichlet eta function and
the Dirichlet beta function:
- <math>
- <math>
- The polylogarithm is equivalent to the Fermi-Dirac integral (GNU)
- <math>
- The polylogarithm is a special case of the Lerch Transcendent (Erdélyi § 1.11-14)
- <math>Li_s(z)=z~\Phi(z,s,1).<math>
- The polylogarithm is related to the Hurwitz zeta function by:
- <math>
- <math>\textrm{Re}(s)>1, \textrm{Im}(x)\ge 0, 0 \le \textrm{Re}(x) < 1<math>
- <math>\textrm{Re}(s)>1, \textrm{Im}(x)\le 0, 0 < \textrm{Re}(x) \le 1.<math>
- Using the relationship between the Hurwitz zeta function and the Bernoulli polynomials:
- <math>
- <math>
- <math>
- The polylogarithm with pure imaginary μ may be expressed in terms of Clausen functions Cis(θ) and Sis(θ) (Lewin, 1958 Ch VII § 1.4, Abramowitz & Stegun § 27.8)
- <math>
- The Inverse Tangent Integral Tis(z) (Lewin, 1958 Ch VII § 1.2) can be expressed in terms of polylogarithms:
- <math>
- The Legendre chi function χs(z) (Lewin, 1958 Ch VII § 1.1, Boersma) can be expressed in terms of polylogarithms:
- <math>
- The polylogarithm may be expressed as a series of Debye Functions Zn(z) (Abramowitz & Stegun § 27.1)
- <math>
- <math>
Series representations
- We may represent the polylogarithm as a power series about μ=0 as follows: (Robinson) Consider the Mellin transform:
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- <math>
- The defining equation may be extended to negative values of the parameter s using a Hankel contour integral (Wood, Gradshteyn & Ryzhik § 9.553):
- <math>
- <math>
- For negative integer s, the polylogarithm may be expressed as a series involving the Eulerian numbers
- <math>
- Another explicit formula for negative integer s is (Wood):
- <math>
Limiting behavior
The following limits hold for the polylogarithm (Wood):
- <math>
\lim_{|z|\rightarrow 0} Li_s(z) = \lim_{s \rightarrow \infty} Li_s(z) = z <math>
- <math>
\lim_{\mathrm{Re}(\mu) \rightarrow \infty} Li_s(e^\mu) = -{\mu^s \over \Gamma(s+1)} ~~~~~~(s\ne -1, -2,-3,\ldots) <math>
- <math>
\lim_{\mathrm{Re}(\mu) \rightarrow \infty} Li_{n}(e^\mu) = -(-1)^ne^{-\mu} ~~~~~~(n=1,2,3,\ldots) <math>
- <math>
\lim_{|\mu|\rightarrow 0} Li_s(e^\mu) = \Gamma(1-s)(-\mu)^s~~~~~~(s<1) <math>
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 <math>\rho=\left(\sqrt{5}-1\right)/2<math> as the reciprocal of the golden ratio. Then two simple examples of results from ladders include
- <math>Li_2(\rho^6)=4Li_2(\rho^3)+3Li_2(\rho^2)-6Li_2(\rho)+\frac{7\pi^2}{30}<math>
given by Coexeter in 1935, and
- <math>Li_2(\rho)=\frac{\pi^2}{10} - \log^2\rho<math>
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
- Template:Book reference
- 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/
- Template:Journal reference issue
- Template:Journal reference 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/
- Template:Book reference pp. 323-326
- Template:Journal reference
- Template:Book reference
- Template:Journal reference issue
- Template:Book reference
- GNU Scientific Library - Reference Manual http://www.gnu.org/software/gsl/manual/gsl-ref.html#SEC117
- Template:Book reference
- Template:Journal reference
- Template:Book reference
- Template:Book reference
- Template:Book reference
- Template:Journal reference
- Template:Journal reference
- Template:Book reference
- Template:Book reference (see § 1.2, pp 23-24)
- Template:Journal reference http://prola.aps.org/abstract/PR/v83/i3/p678_1
- Template:Book reference
- Template:Journal reference issue
- Template:Book reference edition
- 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