Función polilogarítmica

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda

El polilogaritmo (también conocido como función de Jonquière) es una función especial \operatorname{Li}_{s}(z) definida por la siguiente serie:


\operatorname{Li}_s(z) = \sum_{k=1}^\infty {z^k \over k^s}.

Esta no es, en general, una función elemental, aunque esté relacionada con la función logarítmica. La definición dada arriba es válida para todo número complejo s y z tal que \vert z\vert < 1. Para obtener el polilogaritmo en el resto del plano complejo, hay que extender la definición mediante una continuación analítica.

Varias funciones polilogarítmicas en el plano complejo
Complex polylogminus3.jpg
Complex polylogminus2.jpg
Complex polylogminus1.jpg
Complex polylog0.jpg
Complex polylog1.jpg
Complex polylog2.jpg
Complex polylog3.jpg

\operatorname{Li}_{-3}(z)

\operatorname{Li}_{-2}(z)

\operatorname{Li}_{-1}(z)

\operatorname{Li}_{0}(z)

\operatorname{Li}_{1}(z)

\operatorname{Li}_{2}(z)

\operatorname{Li}_{3}(z)

El caso especial s = 1 nos da la relación de estas funciones con el logaritmo (\operatorname{Li}_{1}(z) = -\operatorname{ln}(1 - z)) mientras que los casos especiales s = 2 y s = 3 se denominan dilogaritmo (o función de Spence) y trilogaritmo respectivamente. El nombre de la función proviene del hecho de que podría der definida como integrales iteradas de la misma función:


\operatorname{Li}_{s+1}(z) = \int_0^z \frac {\operatorname{Li}_s(t)}{t}dt

así, el dilogaritmo es una integral del logaritmo, el trilogaritmo del dilogaritmo y así continuamente. Para valores enteros negativos de s, el polilogaritmo es una función racional.

El polilogaritmo también aparece en la forma cerrada de la integral de la distribución de Fermi-Dirac y de la distribución de Bose-Einstein, denominándose a veces como la integral de Fermi-Dirac o la integral de Bose-Einstein. El polilogaritmo no debe confundirse con las funciones polilogarítmicas ni con la función logaritmo integral, la cual tiene una notación similar.

Propiedades[editar]

En el caso en el que el parámetro s sea un entero, este estará representado por n (o -n cuando sea negativo). Suele ser conveniente definir \mu = \operatorname{ln}(z) donde \operatorname{ln}(z) es la rama principal del logaritmo complejo \operatorname{Ln}(z), de tal manera que -\pi < \operatorname{Im}(\mu)<\pi. Además, toda exponencialización se considerará univaluada: z^{s} = e^{s\operatorname{ln}(z)}.

Dependiendo del parámetro s, el polilogaritmo puede ser multivaluado. La rama principal del polilogaritmo se escoge que sea aquella para la que \operatorname{Li}_{s}(z) sea real para z\in [0,1] y sea continua excepto en el eje positivo real, donde hay un corte en el intervalo z\in [1,\infty) tal que el corte coloca a los puntos del eje real en el semiplano inferior de z. En término de \mu, la rama principal está definida para aquellos valores de \mu tales que \operatorname{arg}(\mu)\in (-\pi,\pi]. El hecho de que el polilogaritmo sea discontinuo en \mu puede dar lugar a confusión.

Para z real y z ≥ 1, la parte imaginaria del polilogaritmo es:

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

Si se atraviesa el corte, esto es, tomando un parámetro infinitesimal δ positivo y real, entonces se tiene:

\textrm{Im}(\operatorname{Li}_s(z+i\delta)) = {{\pi \mu^{s-1}}\over{\Gamma(s)}}.

Las derivadas del polilogaritmo son:

z{\partial \operatorname{Li}_s(z) \over \partial z} = \operatorname{Li}_{s-1}(z)
{\partial \operatorname{Li}_s(e^\mu) \over \partial \mu} = \operatorname{Li}_{s-1}(e^\mu).

Distintos valores[editar]

Polylogarithm plot negative.svg

Para valores enteros de s, se tienen las siguientes relaciones explícitas:

\operatorname{Li}_{1}(z)  = -\log\left(1-z\right)
\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}.
\operatorname{Li}_{-4}(z) = {z(1+z)(1+10z+z^2) \over (1-z)^5}.

Para todos los valores negativos de s, se puede expresar el polilogaritmo como un cociente de polinomios en z, siendo por tanto funciones racionales. Algunos valores de polilogaritmos para argumentos semienteros son:

\operatorname{Li}_{1}\left(1/2\right) = \log 2
\operatorname{Li}_{2}(1/2) = {1 \over 12}[\pi^2-6(\log 2)^2]
\operatorname{Li}_{3}(1/2) = {1 \over 24}[4(\log 2)^3-2\pi^2(\log 2)+21\,\zeta(3)]

donde ζ es la función zeta de Riemann. No se conocen fórmulas similares para mayores órdenes.

Expresiones alternativas[editar]


\operatorname{Li}_{s+1}(z) = {1 \over \Gamma(s+1)}
\int_0^\infty {t^s \over e^t/z-1} dt.
Esta integral converge para \operatorname{Re}(s)>0 y para todo z exceptuando los z reales y ≥ 1. En este contexto se suele llamar al polilogaritmo integral de Bose o integral de Bose-Einstein.

-\operatorname{Li}_{s+1}(-z) = {1 \over \Gamma(s+1)}
\int_0^\infty {t^s \over e^t/z+1} dt.
Esta integral converge para \operatorname{Re}(s)>0 y para todo z exceptuando los z reales menores que -1. En este contexto se suele llamar al polilogaritmo integral de Fermi o integral de Fermi-Dirac.
  • El polilogaritmo puede ser también representado por una integral de contorno. Mientras que el polo t = μ del integrando no caiga en el eje real no negativo y s\neq \,1,\,2,\,3,\,\dots, se tiene:

\operatorname{Li}_s(e^\mu)={{-\Gamma(1-s)}\over{2\pi i}}\oint_H {{(-t)^{s-1}}\over{e^{t-\mu}-1}}dt.
donde H representa un contorno de Hanker. El integrando tiene un corte en el eje real desde el cero hasta el infinito, encontrándose el eje real en el semiplano inferior (es decir, siendo alcanzable el eje real desde el semiplano inferior de forma continua y desde el semiplano superior de forma discontinua). Para el caso en el que μ es real y no negativo, podemos simplemente añadir la contribución del polo:

\operatorname{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
donde R es el residuo del polo:

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

\operatorname{Li}_s(-z) + \operatorname{Li}_s(z) = 2^{1-s} ~ \operatorname{Li}_s(z^2).
Nótese que la función de Kummer obedece una fórmula de duplicación muy parecida. Este es un caso especial de la fórmula de multiplicación, para cualquier entero p:

\sum_{m=0}^{p-1}\operatorname{Li}_s(ze^{2\pi i m/p}) = p^{1-s}\,\operatorname{Li}_s(z^p)
La cual puede ser porbada usando la definición mediante series de la función polilogarítmica y la ortogonalidad de los términos exponenciales (véase transformada discreta de Fourier).

Relación con otras funciones[editar]


\operatorname{Li}_s(1) = \zeta(s)~~~~~~~~~~~~~(\textrm{Re}(s)>1).

\operatorname{Li}_s(-1) = -\eta\left(s\right)
Cuando se aplica sobre argumentos imaginarios puros se tiene:

\operatorname{Li}_s(\pm i) = 2^{-s}~\eta(s)\pm i~\beta(s)\,
  • El polilogaritmo es equivalente a la integral de Fermi-Dirac

F_s(\mu)=-\operatorname{Li}_{s+1}(-e^\mu).\,
\operatorname{Li}_s(z)=z~\Phi(z,s,1).

\operatorname{Li}_s(e^{2\pi i x})+(-1)^s~\operatorname{Li}_s(e^{-2\pi i x})={(2\pi i)^s \over \Gamma(s)}~\zeta\left
(1-s,x\right)
donde \Gamma(s) es la función gamma. Esta última igualdad es válida para
\textrm{Re}(s)>1, ~\textrm{Im}(x)\ge 0, ~0 \le \textrm{Re}(x) < 1
y para
\textrm{Re}(s)>1, ~\textrm{Im}(x)\le 0, ~0 <   \textrm{Re}(x) \le 1.

Esta ecuación proporciona la continuación analítica de la representación mediante series del polilogaritmo más allá del círculo de convergencia \vert z\vert = 1.

Por otra parte, para todo s\in\mathbb{C} y para todo ~z\not\in~]0;1[, la fórmula de inversión es:


\operatorname{Li}_s(z)+(-1)^s~\operatorname{Li}_s(1/z) = {(2\pi i)^s \over \Gamma(s)}~\zeta\left
(1\!-\!s,~\frac{1}{2}+{\ln(-z)\over2i\pi}\right),
mientras que para todo s\in\mathbb{C} y para todo z~\not\in~]1;+\infty[

\operatorname{Li}_s(z)+(-1)^s~\operatorname{Li}_s(1/z) = {(2\pi i)^s \over \Gamma(s)}~\zeta\left
(1\!-\!s,~\frac{1}{2}-{\ln(-1/z)\over2i\pi}\right),

Más abajo se tiene esta fórmula para el caso en el que s sea un entero, en el cual se simplifica bastante.



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

que se cumple para todo x y n = \,0\,1,\,2,\,3,\dots se obtiene la siguiente relación:


\operatorname{Li}_{n}(e^{2\pi i x})+ (-1)^n~\operatorname{Li}_{n}(e^{-2\pi i x}) 
= -{(2 \pi i)^n\over n!} B_n\left({x}\right)

Con las mismas condiciones sobre x que arriba, para valores enteros negativos del parámetro s, se tiene que, para todo z\in\mathbb{C}:


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

Más generalmente, n=0,\pm1,\pm2,\pm3,\cdots:


\operatorname{Li}_{n}(z)+ (-1)^n~\operatorname{Li}_{n}\left(1/z\right) = -\frac{(2i\pi)^n}{n!}\,
B_n\left(\frac{1}{2}+{\ln(-z)\over 2i\pi}\right)  \qquad z~\not\in~]0;1[,

\operatorname{Li}_{n}(z)+ (-1)^n~\operatorname{Li}_{n}\left(1/z\right) = -\frac{(2i\pi)^n}{n!}\,
B_n\left(\frac{1}{2}-{\ln(-1/z)\over 2i\pi}\right)  \qquad z~\not\in~]1;+\infty[.
  • El polilogaritmo sobre valores de μ imaginarios puros puede expresarse en término de las funciones de Chausen \operatorname{Ci}_{s}(\theta) y \operatorname{Si}_{s}(\theta).

\operatorname{Li}_s(e^{\pm i \theta}) = Ci_s(\theta) \pm i~Si_s(\theta).

\operatorname{Li}_s(\pm iy) = 2^{-s}~\operatorname{Li}_s(-y^2) \pm i~\operatorname{Ti}_s(y).

\chi_s(z) = {1 \over 2}~[\operatorname{Li}_s(z)-\operatorname{Li}_s(-z)].

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

Otra expresión bastante similar relaciona la función de Debye con el polilogaritmo:


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


Representaciones en forma de serie[editar]

  • Se puede representar al polilogaritmo como una serie de potencias en torno a \mu = 0 como sigue: Comenzemos con la transformada de Mellin siguiente:

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.

El cambio de variables t = ab, u = a(1 - b) transforma a la integral doble en una integral separable:


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).

Para f = 1 tenemos, usando la transformada inversa de Mellin:


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

donde c es una constante situada a la derecha de los polos del integrando. El camino de integración puede deformarse en un contorno cerrado, y los polos del integrando son los de la función gamma \Gamma(r) en r = \,0\,-1,\,-2,\dots y el polo de la función zeta de Riemann \zeta(s + r) situado en r = 1 - s. Sumando estos residuos, para \vert\mu\vert<2\pi y s\neq \,1,\,2,\,3,\dots


\operatorname{Li}_s(e^\mu) =
\Gamma(1-s)(-\mu)^{s-1} +
\sum_{k=0}^\infty {\zeta(s-k) \over k!}~\mu^k.

Si el parámetro s es un entero positivo n, tanto el término k = n - 1 como la función gamma se hacen infinito, aunque su suma no. Para enteros k > 0 se tiene:


\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)

y para k = 0:


\lim_{s\rightarrow 1}\left[
\zeta(s)+\Gamma(1-s)(-\mu)^{s-1}\right]
= -\ln(-\mu).

Así, para s = n, donde n es un entero positivo y \vert\mu\vert<2\pi, se tiene lo siguiente:


\operatorname{Li}_{n}(e^\mu) =
{\mu^{n-1} \over (n-1)!}\left(H_{n-1}-\ln(-\mu)\right) +

\sum_{k=0,k\ne n-1}^\infty {\zeta(n-k) \over k!}~\mu^k, 
~~~~~~~~~~~~~~~~~~~~~~n=2,3,4,\ldots

\operatorname{Li}_{1}(e^\mu) =-\ln(-\mu)+
\sum_{k=1}^\infty {\zeta(1-k) \over k!}~\mu^k, 
~~~~~~~~~~(n=1)

donde H_{n} es un número armónico:


H_n = \sum_{k=1}^n{1\over k}.

El término problemático contiene a -\operatorname{ln}(-\mu) el cual, cuando se multiplica por \mu^{k} tiende a cero cuando μ tiende a cero, excepto en el caso en el que k = 0. Esto refleja el hecho de que el polilogaritmo \operatorname{Li}_{s}(z) contiene una singularidad logarítmica cuando s = 1 y z = 1, entonces:


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

Usando la relación existente entre la función zeta de Riemann y los números de Bernoulli B_{k}:


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

se obtiene, para órdenes s negativos enteros del polilogaritmo y \vert\mu\vert<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

Como, exceptuando B_{1}, todos los números de Bernoulli con índice impar son cero, se obtiene el término n = 0 usando \zeta(0) = B_{1} = -1/2.


  • Como se ha dicho anteriormente, el polilogaritmo puede ser extendido a valores negativos del parámetro s usando una integral de contorno:

\operatorname{Li}_s(e^\mu)=-{\Gamma(1-s) \over 2\pi i}\oint_H{(-t)^{s-1} \over e^{t-\mu}-1}dt

donde H es el contorno de Hankel, s\neq \,1,\,2,\,3,\dots, y el polo t = \mu del integrando no cae en el eje no negativo real. El contorno de Hanker puede deformarse para que contenga los polos del integrando en t-\mu=2k\pi i y la integral puede ser evaluada como la suma de los residuos de los polos que contenga con la fórmula integral de Cauchy.


\operatorname{Li}_s(e^\mu)=\Gamma(1-s)\sum_{k=-\infty}^\infty (2k\pi i-\mu)^{s-1}.

Esto se cumplirá para \operatorname{Re}(s)<0 y para todo μ excepto cuando e^{\mu} = 1. Sumando estas series se obtiene:


\operatorname{Li}_s(e^\mu)=-\sum_{k=0}^\infty \frac{1}{k!} 
\left[1-\frac{2}{2^{s-k}}\right]\zeta(s-k) (\mu-\pi i)^k

Nótese que esta suma es más compacta escrita en términos de la función eta de Dirichlet.



\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

donde \left\langle{n\atop i}\right\rangle son los números Eulerianos.


  • Otra fórmula explícita para enteros negativos s es:

\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)

donde S(n,k) son los números de Stirling de segunda especie.


Límites del polilogaritmo[editar]

Los siguientes límites se cumplen para el polilogaritmo:


\lim_{|z|\rightarrow 0} \operatorname{Li}_s(z) = 0

\lim_{s \rightarrow \infty}
\operatorname{Li}_s(z) = z

\lim_{\mathrm{Re}(\mu) \rightarrow \infty} \operatorname{Li}_s(e^\mu) = -{\mu^s \over \Gamma(s+1)}
~~~~~~(s\ne -1, -2,-3,\ldots)

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

\lim_{|\mu|\rightarrow 0} \operatorname{Li}_s(e^\mu) =  \Gamma(1-s)(-\mu)^{s-1}~~~~~~(s<1)

Dilogaritmo[editar]

El dilogaritmo es el polilogaritmo con s=2. Una representación integral del dilogaritmo es:


\operatorname{Li}_2 (z) = -\int_0^z{\ln (1-t) \over t} dt.

La identidad de Abel para el dilogaritmo está dada por

\ln(1-x)\ln(1-y)= 
\mbox{Li}_2 \left( \frac{x}{1-y} \right)
+\mbox{Li}_2 \left( \frac{y}{1-x} \right)
-\mbox{Li}_2 \left(x \right)
-\mbox{Li}_2 \left(y \right)
-\mbox{Li}_2 \left( \frac{xy}{(1-x)(1-y)} \right)
.

Es inmediato que esta igualdad se cumple en los casos x = 0 y y=0, y es fácilmente verificable por diferenciación ∂/∂x ∂/∂y. Para y = 1-x la identidad se reduce a la fórmula de reflexión de Euler:

\mbox{Li}_2 \left(x \right)+\mbox{Li}_2 \left(1-x \right)= 
\frac{\pi^2}{6}-\ln(x)\ln(1-x)

donde se ha usado \operatorname{Li}_{2}(1) = \frac{\pi^{2}}{6}.

Haciendo el cambio de variables u = x/(1-y), v = y/(1-x), la identidad de Abel se transforma así:

\mbox{Li}_2(u)+\mbox{Li}_2(v)-\mbox{Li}_2(uv)=
\mbox{Li}_2 \left( \frac{u-uv}{1-uv} \right)
+\mbox{Li}_2 \left( \frac{v-uv}{1-uv} \right)
+\ln \left( \frac{1-u}{1-uv} \right)\ln\left( \frac{1-v}{1-uv} \right)
,

llamándose a esta la Identidad Pentagonal.

Nota histórica:Don Zagier remarcó que "El dilogaritmo es la única función matemática con sentido del humor."

Escaleras de polilogaritmos[editar]

Leonard Lewin descubrió una generalización de varias relaciones clásicas para valores especiales del polilogaritmo, estas son las escaleras de polilogaritmoseee. Sea \rho=\left(\sqrt{5}-1\right)/2 el recíproco del número áureo. Entonces, dos ejemplos de estas relaciones son los siguientes:

\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}
\operatorname{Li}_2(\rho)=\frac{\pi^2}{10} - \log^2\rho

Estas relaciones aparecen de forma natural en teoría K y Geometría algebraica.

Monodromía[editar]

El polilogaritmo tiene dos puntos de ramificación, uno en z=1 y otro en z=0. El segundo punto (el situado en z=0) no es visible en la hoja principal del polilogaritmo. Este sólo se hace visible cuando el polilogaritmo es extendido mediante continuación analítica a sus otras hojas. El grupo de monodromía para el polilogaritmo consiste en las clases de homotopía de los caminos que rodean los puntos de ramificación. Denotando a estos puntos como m_0 y m_1, el grupo de monodromía tiene la siguiente representación de grupo:

\langle m_0, m_1\vert w=m_0m_1m^{-1}_0m^{-1}_1,\, wm_1=m_1w\rangle

Para el caso especial de dilogaritmo, se tiene además que wm_0=m_0w, y el grupo de monodromías es el grupo de Heisenberg (identificando m_0,m_1 y w con x,y,z).


Referencias[editar]

  • Erdélyi, A., Magnus, W., Oberhettinger, F., and Tricomi, F. G. (1981). Higher Transcendental Functions, Vol. 1. New York: Krieger. 
  • Fornberg, B.; Kölbig, K. S. (1975). «Complex zeros of the Jonquiére or polylogarithm function». Mathematics of Computation 29 (130):  pp. 582–599. doi:10.2307/2005579. 
  • Gradshteyn, I.S. and Ryzhik, I.M. (1980). Tables of Integrals, Series, and Products. Academic Press, New York. ISBN 0-12-294760-6. 
  • GNU Scientific Library. «Reference Manual».
  • Jahnke, E. and Emde, F. (1945). Tables of Functions with Formulae and Curves. Dover. 
  • Kölbig, K. S.; Mignaco, J. A. , and Remiddi, E. (1970). «On Nielsen's generalized polylogarithms and their numerical calculation». BIT 10:  pp. 38–74. doi:10.1007/BF01940890. 
  • Lewin, L. (1958). Dilogarithms and Associated Functions. Macdonald - London. 
  • Lewin, L. (1981). Polylogarithms and Associated Functions. North-Holland-New York. ISBN 0-444-00550-1. 
  • Lewin, Leonard (Ed.) (1991). Structural Properties of Polylogarithms. Amer. Math. Soc. - Providence, RI. ISBN 0-8218-1634-9. 
  • McDougall, J.; Stoner, E. C. (1939). «The computation of Fermi-Dirac functions». Philosophical Transactions of the Royal Society, Series A 237:  pp. 67–104. doi:10.1098/rsta.1938.0004. 
  • Markman, B. (1965). «The Riemann Zeta Function». BIT 5:  pp. 138–141. 
  • Nielsen, N. (1909). Der Euler'sche Dilogarithms. Halle - Leipzig, Germany. 
  • Prudnikov, A. P.; Marichev, O. I.; and Brychkov, Yu. A. (1990). Integrals and Series, Vol. 3 (More Special Functions): The Generalized Zeta Function, Bernoulli Polynomials, Euler Polynomials, and Polylogarithms. Gordon and Breach - Newark, NJ.  (see § 1.2, pp 23-24)
  • Robinson, J. E. (1951). «Note on the Bose-Einstein integral functions». Physical Review, Series 2 83:  pp. 678–679. http://prola.aps.org/abstract/PR/v83/i3/p678_1. 
  • Rogers, L. J. (1907). On Function Sum Theorems connected with the series \sum_{n=1}^\infty\frac{x^n}{n^2}. 4.  pp. 169-189. 
  • Schrödinger, E. (1952). Statistical Thermodynamics. Cambridge. 
  • Truesdell, C. (1945). «On a function which occurs in the theory of the structure of polymers». Annals of Mathematics, Series 2 46 (1):  pp. 144–1457. 
  • Vepstas, Linas (February 2007). «An efficient algorithm for accelerating the convergence of oscillatory series, useful for computing the polylogarithm and Hurwitz zeta functions». ArXiv math.CA (math.CA/0702243). http://arxiv.org/abs/math.CA/0702243. 
  • Whittaker, E. T., and Watson, G. N. (1962). A Course of Modern Analysis (Fourth edition edición). Cambridge University Press. 
  • Wood, David C. (June, 1992). «Technical Report 15-92». University of Kent computing Laboratory, University of Kent, Canterbury, UK.