Fonction polylogarithme

La fonction polylogarithme (aussi connue sous le nom de fonction de Jonquière) est une fonction spéciale qui peut être définie pour tout s et |z| < 1 par :

$\operatorname {Li} _{s}(z)=\sum _{k=1}^{\infty }{z^{k} \over k^{s}}.$

Le paramètre s et l'argument z sont pris sur l'ensemble ℂ des nombres complexes. Les cas particuliers s = 2 et s = 3 sont appelés le polylogarithme d'ordre 2 ou dilogarithme et le polylogarithme d'ordre 3 ou trilogarithme respectivement. Le polylogarithme apparaît aussi dans la forme fermée de l'intégrale de la distribution de Fermi-Dirac et la distribution de Bose-Einstein et est quelquefois connue comme l'intégrale de Fermi-Dirac ou l'intégrale de Bose-Einstein.

Par prolongement analytique, on peut également donner un sens au polylogarithme pour |z| ≥ 1.

Différentes fonctions polylogarithmes dans le plan complexe :

Li-3(z)	Li-2(z)	Li-1(z)	Li0(z)	Li1(z)	Li2(z)	Li3(z)

Propriétés

Dans le cas important où le paramètre s est un nombre entier, il sera représenté par n (ou -n lorsqu'il est négatif). Il est souvent pratique de définir μ = ln(z) où ln est la branche principale du logarithme naturel, c’est-à-dire $-\pi <\Im (\mu )\leq \pi$. Ainsi, toute l'exponentiation sera supposée être à valeur unique. (e.g. z^s = e^{s ln(z)}).

Dépendant du paramètre s, le polylogarithme peut être à valeurs multiples. La branche principale du polylogarithme est celle pour laquelle Li_s(z) est réel pour z réel, 0 ≤ z ≤ 1 et est continu excepté sur l'axe des réels positifs, où une coupure est faite de z = 1 à l'infini telle que la coupure place l'axe réel sur le demi-plan le plus bas de z. En termes de μ, ceci s'élève à –π < arg(–μ) < π. Le fait que le polylogarithme puisse être discontinu en μ peut causer une certaine confusion.

Pour z réel et supérieur ou égal à 1, la partie imaginaire du polylogarithme est (Wood 1992) :

$\Im (\operatorname {Li} _{s}(z))=-{{\pi \mu ^{s-1}} \over {\Gamma (s)}}.$

En traversant la coupure :

$\lim _{\delta \rightarrow 0^{ }}\Im (\operatorname {Li} _{s}(z \mathrm {i} \delta ))={{\pi \mu ^{s-1}} \over {\Gamma (s)}}.$

Les dérivées du polylogarithme s'expriment également avec le polylogarithme :

$z{\partial \operatorname {Li} _{s}(z) \over \partial z}=\operatorname {Li} _{s-1}(z),\qquad {\partial \operatorname {Li} _{s}(\mathrm {e} ^{\mu }) \over \partial \mu }=\operatorname {Li} _{s-1}(\mathrm {e} ^{\mu }).$

Valeurs particulières

Voir aussi la section « Relation de parenté avec les autres fonctions » ci-dessous.

Pour les valeurs entières de s, on peut écrire les expressions explicites suivantes :

$\operatorname {Li} _{1}(z)=-\ln(1-z)$

$\operatorname {Li} _{0}(z)={z \over 1-z}$

$\operatorname {Li} _{-1}(z)={z \over (1-z)^{2}}$

$\operatorname {Li} _{-2}(z)={z(1 z) \over (1-z)^{3}}$

$\operatorname {Li} _{-3}(z)={z(1 4z z^{2}) \over (1-z)^{4}}.$

Le polylogarithme, pour toutes les valeurs entières négatives de s, peut être exprimé comme une fraction rationnelle en z (voir les représentations en série ci-dessus). Certaines expressions particulières pour les demi valeurs entières de l'argument sont :

$\operatorname {Li} _{1}\left({\frac {1}{2}}\right)=\ln(2)$

$\operatorname {Li} _{2}\left({\frac {1}{2}}\right)={\pi ^{2} \over 12}-{1 \over 2}\ln ^{2}2$

$\operatorname {Li} _{3}\left({\frac {1}{2}}\right)={1 \over 24}[4(\ln 2)^{3}-2\pi ^{2}\ln 2 21\,\zeta (3)]$

où $\zeta$ est la fonction zêta de Riemann. Aucune formule similaire de ce type n'est connue pour des ordres plus élevés (Lewin 1991, p. 2) ; par ailleurs, les seules valeurs connues de Li₂ exprimables à l'aide des fonctions élémentaires sont les huit valeurs suivantes :

$\operatorname {Li} _{2}(0)=0\ ;\ \operatorname {Li} _{2}(1)={\dfrac {\pi ^{2}}{6}}\ ;\ \operatorname {Li} _{2}(-1)={\dfrac {-\pi ^{2}}{12}}\ ;\ \operatorname {Li} _{2}(1/2)={\pi ^{2} \over 12}-{1 \over 2}\ln ^{2}2$

ainsi que

$\operatorname {Li} _{2}\left({3-{\sqrt {5}} \over 2}\right)={\dfrac {\pi ^{2}}{15}}-\ln ^{2}\left({1 {\sqrt {5}} \over 2}\right)\ ;\ \operatorname {Li} _{2}\left({-1 {\sqrt {5}} \over 2}\right)={\dfrac {\pi ^{2}}{10}}-\ln ^{2}\left({1 {\sqrt {5}} \over 2}\right)$ ;

$\operatorname {Li} _{2}\left({1-{\sqrt {5}} \over 2}\right)=-{\dfrac {\pi ^{2}}{15}} {1 \over 2}\ln ^{2}\left({1 {\sqrt {5}} \over 2}\right)\ ;\ \operatorname {Li} _{2}\left({-1-{\sqrt {5}} \over 2}\right)=-{\dfrac {\pi ^{2}}{10}} {1 \over 2}\ln ^{2}\left({1 {\sqrt {5}} \over 2}\right)$.

Expressions alternatives

• L'intégrale de la distribution de Bose-Einstein est exprimée en termes de polylogarithme :

$\operatorname {Li} _{s 1}(z)\equiv {1 \over \Gamma (s 1)}\int _{0}^{\infty }{t^{s} \over \mathrm {e} ^{t}/z-1}~\mathrm {d} t$

Celle-ci converge pour $\Re (s)>0$ et tous les z excepté pour les z réels et supérieurs ou égaux à 1. Le polylogarithme dans ce contexte est quelquefois connu comme l'intégrale de Bose ou de Bose-Einstein.

• L'intégrale de la distribution de Fermi-Dirac est aussi exprimée en termes de polylogarithme :

$-\operatorname {Li} _{s 1}(-z)\equiv {1 \over \Gamma (s 1)}\int _{0}^{\infty }{t^{s} \over \mathrm {e} ^{t}/z 1}~\mathrm {d} t.$

Celle-ci converge pour $\Re (s)>0$ et tous les z excepté pour les z réels et strictement inférieurs à -1. Le polylogarithme dans ce contexte est quelquefois connu comme l'intégrale de Fermi ou l'intégrale de Fermi-Dirac. (GNU)

• Le polylogarithme peut plutôt être généralement représenté par une intégrale sur un contour de Hankel (en) (Whittaker et Watson 1927).
  • Tant que le pôle t = μ de l'intégrande n'est pas relié à l'axe réel positif, et $s\neq 1,2,3\ldots$, on a :

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

    où H représente le contour de Hankel. L'intégrande possède une coupure le long de l'axe réel de zéro à l'infini, l'axe réel étant sur la moitié inférieure de la feuille ($\Im (t)\leq 0$).

  • Pour le cas où μ est réel et positif, nous pouvons simplement ajouter la contribution limitante du pôle :

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

    où R est le résidu du pôle :

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

• La relation carrée est facilement vue à partir de l'équation (voir aussi Clunie 1954 et Schrödinger 1952)

  $\operatorname {Li} _{s}(-z) \operatorname {Li} _{s}(z)=2^{1-s}~\operatorname {Li} _{s}(z^{2})$

La fonction de Kummer obéit à une formule de duplication très similaire.

Relation de parenté avec les autres fonctions

• Pour z = 1, le polylogarithme se réduit à la fonction zêta de Riemann

$\operatorname {Li} _{s}(1)=\zeta (s)\quad ({\textrm {Re}}(s)>1)$

• Le polylogarithme est relié à la fonction êta de Dirichlet et la fonction bêta de Dirichlet :

$\operatorname {Li} _{s}(-1)=-\eta (s)$

où η(s) est la fonction êta de Dirichlet.

Pour des arguments imaginaires purs, nous avons :

$\operatorname {Li} _{s}(\pm \mathrm {i} )=2^{-s}\eta (s)\pm \mathrm {i} \beta (s)$

où β(s) est la fonction bêta de Dirichlet.

• Le polylogarithme est équivalent à l'intégrale de Fermi-Dirac (GNU)

$F_{s}(\mu )=-\operatorname {Li} _{s 1}(-\mathrm {e} ^{\mu })\,$

• Le polylogarithme est un cas particulier de la fonction transcendante de Lerch (Erdélyi et al. 1981, § 1.11-14)

$\operatorname {Li} _{s}(z)=z~\Phi (z,s,1)$

• Le polylogarithme est relié à la fonction zêta de Hurwitz par :

$\operatorname {Li} _{s}(\mathrm {e} ^{2\pi \mathrm {i} x}) (-1)^{s}\operatorname {Li} _{s}(\mathrm {e} ^{-2\pi \mathrm {i} x})={(2\pi \mathrm {i} )^{s} \over \Gamma (s)}~\zeta \left(1-s,x\right)$

où Γ(s) est la fonction Gamma d'Euler. Ceci est valable pour ${\textrm {Re}}(s)>1,{\textrm {Im}}(x)\geq 0,0\leq {\textrm {Re}}(x)<1$

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

(l'équation équivalente d'Erdélyi et al. 1981, § 1.11-16 n'est pas correcte si on suppose que les branches principales du polylogarithme et le logarithme sont utilisés simultanément).

Cette équation fournit le prolongement analytique de la représentation en série du polylogarithme au-delà de son cercle de convergence |z| = 1.

• En utilisant la relation entre la fonction zêta de Hurwitz et les polynômes de Bernoulli :

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

qui reste valable pour tous les x réels et n entier positif, il peut être remarqué que :

$\operatorname {Li} _{n}(\mathrm {e} ^{2\pi \mathrm {i} x}) (-1)^{n}\operatorname {Li} _{n}(\mathrm {e} ^{-2\pi \mathrm {i} x})=-{(2\pi \mathrm {i} )^{n} \over n!}B_{n}(x)$

sous les mêmes contraintes sur s et x que ci-dessus. (L'équation correspondante d'Erdélyi et al. 1981, § 1.11-18 n'est pas correcte). Pour les valeurs entières négatives du paramètre, on a pour tous les z (Erdélyi et al. 1981, § 1.11-17) :

$\operatorname {Li} _{-n}(z) (-1)^{n}\operatorname {Li} _{-n}\left({\frac {1}{z}}\right)=0,\ n=1,2,3\ldots$

• Le polylogarithme avec un μ imaginaire pur peut être exprimé en termes de fonctions de Clausen Ci_s(θ) et Si_s(θ) (Lewin 1958, ch. VII, § 1.4 et Abramowitz & Stegun, § 27.8)

$\operatorname {Li} _{s}(\mathrm {e} ^{\pm \mathrm {i} \theta })=\operatorname {Ci} _{s}(\theta )\pm \mathrm {i} \operatorname {Si} _{s}(\theta )$

• L'arc tangente intégral Ti_s(z) (Lewin 1958, ch. VII, § 1.2) peut être exprimé en termes de polylogarithmes :

$\operatorname {Li} _{s}(\pm iy)=2^{-s}\operatorname {Li} _{s}(-y^{2})\pm \mathrm {i} \,\operatorname {Ti} _{s}(y)$

• La fonction chi de Legendre χ_s(z) (Lewin 1958, ch VII, § 1.1 et Boersma et Dempsey 1992) peut être exprimée en termes de polylogarithmes :

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

• Le polylogarithme peut être exprimé comme une série de fonctions de Debye Z_n(z) (Abramowitz & Stegun, § 27.1)

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

Une expression remarquablement similaire relie la fonction de Debye au polylogarithme :

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

Représentations en séries

On peut représenter le polylogarithme comme une série de puissances pour μ = 0 comme suit (Robinson 1951). Considérons la transformation de Mellin :

$M_{s}(r)=\int _{0}^{\infty }{\textrm {Li}}_{s}(f\mathrm {e} ^{-u})u^{r-1}~\mathrm {d} u={1 \over \Gamma (s)}\int _{0}^{\infty }\int _{0}^{\infty }{t^{s-1}u^{r-1} \over \mathrm {e} ^{t u}/f-1}~\mathrm {d} t~\mathrm {d} u$

Le changement de variables t = ab, u = a(1-b) permet à l'intégrale d'être séparée :

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

pour f = 1 on a, à travers la transformation inverse de Mellin :

$\operatorname {Li} _{s}(\mathrm {e} ^{-u})={1 \over 2\pi \mathrm {i} }\int _{c-\mathrm {i} \infty }^{c \mathrm {i} \infty }\Gamma (r)\zeta (s r)u^{-r}~\mathrm {d} r$

où c est une constante à droite des pôles de l'intégrande.

Le chemin d'intégration peut être converti en un contour fermé, et les pôles de l'intégrande sont ceux de Γ(r) à r = 0 , –1, –2, … et de ζ(s r) à r = 1–s. Sommer les résidus donne, pour $|\mu |<2\pi$ et $s\neq 1,2,3,\ldots$

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

Si le paramètre s est un entier positif n, ainsi que le k = n–1^e terme, la fonction gamma devient infinie, bien que leur somme ne l'est pas. Pour un entier k > 0, nous avons :

$\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}-\ln(-\mu )\right)$

et pour k = 0 :

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

Ainsi, pour s = n où n est un entier positif et $|\mu |<2\pi$, nous avons :

$\operatorname {Li} _{n}(\mathrm {e} ^{\mu })={\mu ^{n-1} \over (n-1)!}\left(H_{n-1}-\ln(-\mu )\right) \sum _{k=0,k\neq n-1}^{\infty }{\zeta (n-k) \over k!}~\mu ^{k},\ (n=2,3,4,\ldots )$

$\operatorname {Li} _{1}(\mathrm {e} ^{\mu })=-\ln(-\mu ) \sum _{k=1}^{\infty }{\zeta (1-k) \over k!}~\mu ^{k}~~~~~~~~~~(n=1)$

où H_n-1 est un nombre harmonique :

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

Le problème des termes contient maintenant -ln(-μ) qui, lorsqu'ils sont multipliés par μ^k, tendront vers zéro quand μ tend vers zéro, excepté pour k = 0. Ceci reflète le fait qu'il existe une vraie singularité logarithmique en Li_s(z) en s = 1 et z = 1, puisque :

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

En utilisant la parenté entre la fonction zêta de Riemann et les nombres de Bernoulli B_k

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

on obtient pour les valeurs entières négatives de s et $|\mu |<2\pi$ :

$\operatorname {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 )$

puisque, excepté pour B₁, tous les nombres de Bernoulli impairs sont égaux à zéro. Nous obtenons le terme n = 0 en utilisant $\zeta (0)=B_{1}=-{\frac {1}{2}}$. Encore, l'équation équivalent d'Erdélyi et al. 1981, § 1.11-15 n'est pas correcte si nous supposons que les branches principales du polylogarithme et le logarithme sont utilisées simultanément, puisque $\ln \left({\frac {1}{z}}\right)$ n'est pas uniformément égal à –ln(z).

L'équation définie peut être étendue aux valeurs négatives du paramètre s en utilisant une intégrale sur un contour de Hankel (en) (Wood 1992 et Gradshteyn et Ryzhik 1980) :

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

où H est le contour de Hankel qui peut être modifié pour qu'il entoure les pôles de l'intégrande, à t – μ = 2kiπ, l'intégrale peut être évaluée comme la somme des résidus :

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

Ceci restera valable pour ${\textrm {Re}}(s)<0$ et tous les z excepté pour z = 1.

Pour les entiers négatifs s, le polylogarithme peut être exprimé comme une série impliquent les nombres eulériens

$\operatorname {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 )$

où $\left\langle {n \atop i}\right\rangle$ sont les nombres eulériens.

Une autre formule explicite pour les entiers négatifs s est (Wood 1992) :

$\operatorname {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 )$

où S(n,k) sont les nombres de Stirling de deuxième espèce.

Comportement aux limites

Les limites suivantes restent valables pour le polylogarithme (Wood 1992) :

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

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

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

$\lim _{|\mu |\rightarrow 0}\operatorname {Li} _{s}(\mathrm {e} ^{\mu })=\Gamma (1-s)(-\mu )^{s},\ (s<1)$

Échelles de polylogarithmes

Leonard Lewin a découvert une remarquable généralisation d'un grand nombre de relations classiques sur les polylogarithmes pour des valeurs particulières. Celles-ci sont maintenant appelées les échelles de polylogarithmes. Définissons $\rho ={\frac {({\sqrt {5}}-1)}{2}}$ comme l'inverse du nombre d'or. Alors, deux exemples simples des résultats issus des échelles incluent

$\operatorname {Li} _{2}(\rho ^{6})=4\operatorname {Li} _{2}(\rho ^{3}) 3\operatorname {Li} _{2}(\rho ^{2})-6\operatorname {Li} _{2}(\rho ) {\frac {7\pi ^{2}}{30}}$

donné par Coxeter en 1935, et

$\operatorname {Li} _{2}(\rho )={\frac {\pi ^{2}}{10}}-\log ^{2}\rho$

donné par Landen. Les échelles de polylogarithmes apparaissent naturellement et profondément en K-théorie.

Notes et références

• (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne)
• (en) J. Boersma et J. P. Dempsey, « On the evaluation of Legendre's chi-function », Math. Comp., vol. 59, n^o 199,‎ 1992, p. 157-163
• (en) J. Clunie, « On Bose-Einstein functions », Proc. Phys. Soc. Section A, vol. 67,‎ 1954, p. 632-636
• (en) A. Erdélyi (en), W. Magnus, F. Oberhettinger (de) et F. G. Tricomi, Higher Transcendental Functions, vol. 1, New York, Krieger, 1981
• (en) I. S. Gradshteyn et I. M. Ryzhik, Tables of Integrals, Series, and Products, New York, Academic Press, 1980, § 9.553
• (en) L. Lewin, Dilogarithms and Associated Functions, London, Macdonald, 1958
• (en) Leonard Lewin, Structural Properties of Polylogarithms, 1991 [détail de l’édition]
• (en) J. E. Robinson, « Note on the Bose-Einstein integral functions », Physical Review, 2^e série, vol. 83,‎ 1951, p. 678-679 (lire en ligne)
• (en) E. Schrödinger, Statistical Thermodynamics, Cambridge, 1952
• (en) E. T. Whittaker et G. N. Watson, A Course of Modern Analysis, Cambridge, 1927, § 12.22 et 13.13
• (en) David C. Wood, Technical Report 15-92, Canterbury, UK, University of Kent Computing Laboratory, juin 1992 (lire en ligne)

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Polylogarithm » (voir la liste des auteurs).

Bibliographie

• (en) David Bailey, Peter Borwein et Simon Plouffe, « On the Rapid Computation of Various Polylogarithmic Constants », Math. Comp., vol. 66, n^o 218,‎ 1997, p. 903–913 (lire en ligne)
• (en) David H. Bailey et David J. Broadhurst, « A Seventeenth-Order Polylogarithm Ladder », 1999, « math.CA/9906134 », texte en accès libre, sur arXiv.
• (en) Jonathan M. Borwein, David M. Bradley, David J. Broadhurst et Petr Lisonek, « Special Values of Multidimensional Polylogarithms », Trans. Amer. Math. Soc., vol. 353,‎ 2001, p. 907-941 (lire en ligne), « math.CA/9910045 », texte en accès libre, sur arXiv.
• (en) B. C. Berndt, Ramanujan's Notebooks, Part IV, Springer, 1994, p. 323-326
• (en) B. Fornberg et K. S. Kölbig, « Complex zeros of the Jonquiére or polylogarithm function », Math. Comp., vol. 29, n^o 130,‎ 1975, p. 582-599
• (en) GNU Scientific Library, Reference Manual
• (en) Eugen Jahnke et Fritz Emde (de), Tables of Functions with Formulae and Curves, Dover, 1945
• (en) K. S. Kölbig, J. A. Mignaco et E. Remiddi, « On Nielsen's generalized polylogarithms and their numerical calculation », BIT, vol. 10,‎ 1970, p. 38-74
• (en) K. S. Kölbig, « Nielsen's Generalized Polylogarithms », SIAM J. Math. Anal., vol. 17,‎ 1986, p. 1232-1258
• (en) L. Lewin, Polylogarithms and Associated Functions, North-Holland, 1981
• (en) J. McDougall et E. C. Stoner (en), « The computation of Fermi-Dirac functions », Philos. Trans. Royal Soc., Series A, vol. 237,‎ 1939, p. 67-104 (lire en ligne)
• (en) B. Markman, « The Riemann Zeta Function », BIT, vol. 5,‎ 1965, p. 138-141
• (de) N. Nielsen (de), « Der Eulersche Dilogarithmus und seine Verallgemeinerungen », Nova Acta Leopoldina, vol. 90,‎ 1909, p. 123-211
• (en) A. P. Prudnikov, O. I. Marichev (de) et Yu. A. Brychkov, The Generalized Zeta Function, Bernoulli Polynomials, Euler Polynomials, and Polylogarithms, Vol. 3: More Special Functions, Gordon and Breach, 1990, p. 23-24, § 1.2
• (en) C. Truesdell, « On a function which occurs in the theory of the structure of polymers », Ann. of Math., 2^e série, vol. 46, n^o 1,‎ 1945, p. 144-1457
• (en) Don Zagier, « Appendix A : Special Values and Functional Equations of Polylogarithms », dans L. Lewin (éd.), Structural Properties of Polylogarithms, AMS, 1991
• (en) Don Zagier, « The Dilogarithm Function », dans P. Cartier, B. Julia, P. Moussa, P. Vanhove (éds.), Frontiers in Number Theory, Physics and Geometry II, Springer-Verlag, 2006 (lire en ligne)