Capítulo 10 — Conexión con Funciones Especiales

Funciones Gamma, Beta, Zeta de Riemann y Polinomios Ortogonales
Prof. Mariano Miguel Lanzi

Las funciones especiales constituyen el lenguaje natural de la física matemática. La función Gamma generaliza el factorial a argumentos complejos; la función Beta conecta integrales de Euler con productos de Gamma; la función Zeta de Riemann codifica la distribución de los números primos en sus propiedades analíticas. Todas ellas se definen y estudian de manera natural mediante integrales en el plano complejo, continuación analítica y el teorema del residuo. Este capítulo final explora estas conexiones profundas entre el análisis complejo y las funciones que aparecen en toda la física teórica.

10.1 Función Gamma $\Gamma(z)$

La función Gamma, introducida por Euler en 1729, extiende el factorial a todo el plano complejo. Su definición integral para $\Re(z) > 0$ es:

$$\Gamma(z) = \int_0^\infty t^{z-1} e^{-t}\,dt, \qquad \Re(z) > 0$$

Integrando por partes se obtiene la ecuación funcional fundamental:

$$\Gamma(z+1) = z\,\Gamma(z)$$

Como $\Gamma(1) = 1$, se sigue por inducción que $\Gamma(n+1) = n!$ para todo entero no negativo $n$. La ecuación funcional permite continuar analíticamente $\Gamma(z)$ a todo el plano complejo excepto los enteros no positivos $z = 0, -1, -2, \ldots$, donde presenta polos simples.

Desarrollo en Polos

En la vecindad de $z = -n$ (con $n = 0, 1, 2, \ldots$), el residuo de la función Gamma está dado por:

$$\operatorname{Res}(\Gamma, -n) = \frac{(-1)^n}{n!}$$

Esta fórmula se obtiene iterando la ecuación funcional $\Gamma(z) = \Gamma(z+n+1)/[z(z+1)\cdots(z+n)]$. La función Gamma es meromorfa en todo $\mathbb{C}$, sin ceros en ninguna parte del plano.

Representación de Weierstrass

La función Gamma admite una representación como producto infinito que exhibe explícitamente todos sus polos:

$$\frac{1}{\Gamma(z)} = z\,e^{\gamma z}\prod_{n=1}^\infty \left(1 + \frac{z}{n}\right) e^{-z/n}$$

donde $\gamma \approx 0.57721$ es la constante de Euler-Mascheroni. Esta expresión muestra que $1/\Gamma(z)$ es una función entera con ceros simples en $z = 0, -1, -2, \ldots$

10.2 Función Beta $B(p,q)$

La función Beta de Euler se define para $\Re(p) > 0$, $\Re(q) > 0$ mediante la integral:

$$B(p,q) = \int_0^1 t^{p-1} (1-t)^{q-1}\,dt$$

La conexión fundamental con la función Gamma es la identidad de Euler:

$$B(p,q) = \frac{\Gamma(p)\,\Gamma(q)}{\Gamma(p+q)}$$

Demostración vía Variable Compleja

La demostración más elegante utiliza el producto de dos integrales de Gamma y el cambio a coordenadas polares en el plano $uv$:

$$\Gamma(p)\Gamma(q) = \int_0^\infty \int_0^\infty u^{p-1}v^{q-1}e^{-(u+v)}\,du\,dv$$

Con la sustitución $u = r\cos^2\theta$, $v = r\sin^2\theta$, el jacobiano produce un factor $r\sin\theta\cos\theta$. La integración en $r$ da $\Gamma(p+q)$ y la integración angular produce la integral que define $B(p,q)/2$, estableciendo la identidad.

La función Beta es simétrica: $B(p,q) = B(q,p)$. La identidad de Euler permite además evaluar integrales que aparecen frecuentemente en mecánica estadística, teoría cuántica de campos y el modelo estándar de partículas.

10.3 Función Zeta de Riemann $\zeta(s)$

La función Zeta de Riemann, introducida por Bernhard Riemann en su célebre memoria de 1859, se define inicialmente para $\Re(s) > 1$ mediante la serie:

$$\zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s}, \qquad \Re(s) > 1$$

Euler descubrió la profunda conexión de esta serie con los números primos a través del producto de Euler:

$$\zeta(s) = \prod_{p\ \text{primo}} \frac{1}{1 - p^{-s}}, \qquad \Re(s) > 1$$

La continuación analítica de $\zeta(s)$ a todo el plano complejo (excepto un polo simple en $s = 1$ con residuo 1) se realiza mediante la fórmula de reflexión de Riemann:

$$\zeta(s) = 2^s \pi^{s-1} \sin\left(\frac{\pi s}{2}\right) \Gamma(1-s)\,\zeta(1-s)$$

Hipótesis de Riemann

La ecuación funcional revela que $\zeta(s)$ tiene ceros triviales en los enteros pares negativos $s = -2, -4, -6, \ldots$ (provienen del factor $\sin(\pi s/2)$). Todos los demás ceros deben estar en la banda crítica $0 \leq \Re(s) \leq 1$. La Hipótesis de Riemann —uno de los siete problemas del milenio del Clay Institute— conjetura que todos los ceros no triviales se encuentran sobre la línea crítica $\Re(s) = 1/2$. Su demostración tendría consecuencias revolucionarias en la teoría de números y en la comprensión de la distribución de los números primos.

Continuación Analítica vía Integral de Contorno

Una de las maneras de continuar $\zeta(s)$ al plano completo utiliza la integral de contorno de Hankel. Partiendo de la representación integral para $\Re(s) > 1$:

$$\zeta(s)\,\Gamma(s) = \int_0^\infty \frac{t^{s-1}}{e^t - 1}\,dt$$

se deforma el contorno a un contorno de Hankel alrededor del semieje real positivo. Esta deformación, justificada por el teorema de Cauchy, produce automáticamente la continuación a todo $s \neq 1$.

10.4 Conexión con Series de Fourier

Los coeficientes de Fourier de una función periódica admiten una interpretación natural como integrales de contorno en el plano complejo. Para una función $f(z)$ analítica en un anillo que contiene la circunferencia unidad, su serie de Laurent centrada en el origen:

$$f(z) = \sum_{n=-\infty}^\infty c_n z^n, \qquad c_n = \frac{1}{2\pi i}\oint_{|z|=1} \frac{f(z)}{z^{n+1}}\,dz$$

evaluada sobre $z = e^{i\theta}$ produce inmediatamente la serie de Fourier de $f(e^{i\theta})$:

$$f(e^{i\theta}) = \sum_{n=-\infty}^\infty c_n e^{in\theta}$$

Relación entre Coeficientes de Laurent y Fourier

Los coeficientes de Laurent $c_n$ de $f(z)$ en el anillo que contiene $|z| = 1$ son exactamente los coeficientes de Fourier de la función $g(\theta) = f(e^{i\theta})$:

$$c_n = \frac{1}{2\pi}\int_0^{2\pi} f(e^{i\theta})\,e^{-in\theta}\,d\theta$$

Esta identidad es fundamental: permite usar el teorema del residuo para calcular sumas de series de Fourier que serían intratables por integración real directa.

Recíprocamente, la convergencia de la serie de Laurent en un anillo se traduce en la convergencia uniforme de la serie de Fourier bajo condiciones de analiticidad, un resultado que conecta el análisis armónico con la teoría de funciones analíticas.

10.5 Fórmula de Reflexión y de Duplicación de Gamma

La fórmula de reflexión de Euler relaciona los valores de la función Gamma en puntos simétricos respecto de $z = 1/2$:

$$\Gamma(z)\,\Gamma(1-z) = \frac{\pi}{\sin(\pi z)}, \qquad z \notin \mathbb{Z}$$

Una demostración elegante utiliza la representación de Weierstrass y la factorización del seno como producto infinito. La fórmula revela inmediatamente que:

$$\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi}$$

La fórmula de duplicación de Legendre comprime el producto de dos evaluaciones de Gamma en una sola:

$$\Gamma(2z) = \frac{2^{2z-1}}{\sqrt{\pi}}\,\Gamma(z)\,\Gamma\left(z + \frac{1}{2}\right)$$

Esta fórmula, junto con la de reflexión, forma la base del cálculo de la función Gamma en argumentos semienteros, esenciales en volúmenes de esferas $n$-dimensionales y en la normalización de armónicos esféricos.

Generalización: Fórmula de Multiplicación de Gauss

Para cualquier entero $m \geq 2$ se tiene la fórmula de multiplicación:

$$\Gamma(mz) = \frac{m^{mz-1/2}}{(2\pi)^{(m-1)/2}} \prod_{k=0}^{m-1} \Gamma\left(z + \frac{k}{m}\right)$$

El caso $m = 2$ recupera la fórmula de duplicación de Legendre. La demostración utiliza la representación límite de la función Gamma: $\Gamma(z) = \lim_{n\to\infty} \frac{n!\,n^z}{z(z+1)\cdots(z+n)}$.

10.6 Relación con Polinomios Ortogonales

Los polinomios ortogonales clásicos —Hermite, Laguerre, Legendre, Chebyshev— admiten funciones generatrices que se expresan como funciones analíticas en el plano complejo. La expansión en serie de Laurent de estas funciones generatrices produce los polinomios como coeficientes, y el teorema del residuo proporciona representaciones integrales.

Polinomios de Hermite $H_n(x)$

La función generatriz de los polinomios de Hermite es una función entera en la variable compleja $t$:

$$e^{2xt - t^2} = \sum_{n=0}^\infty H_n(x)\,\frac{t^n}{n!}$$

Por la fórmula integral de Cauchy para derivadas (o equivalentemente, por la fórmula de Laurent):

$$H_n(x) = \frac{n!}{2\pi i} \oint_C \frac{e^{2xz - z^2}}{z^{n+1}}\,dz$$

donde $C$ es cualquier contorno cerrado simple alrededor del origen. Esta representación es la base para deducir relaciones de recurrencia, ortogonalidad y comportamiento asintótico mediante el método del punto silla.

Polinomios de Laguerre $L_n^{(\alpha)}(x)$

La función generatriz de los polinomios de Laguerre generalizados:

$$\frac{e^{-xt/(1-t)}}{(1-t)^{\alpha+1}} = \sum_{n=0}^\infty L_n^{(\alpha)}(x)\,t^n$$

define una función analítica en $|t| < 1$ con un punto de ramificación en $t = 1$. La representación integral de contorno resultante:

$$L_n^{(\alpha)}(x) = \frac{1}{2\pi i} \oint_C \frac{e^{-xz/(1-z)}}{(1-z)^{\alpha+1}\,z^{n+1}}\,dz$$

permite estudiar el comportamiento de estos polinomios en el límite $n \to \infty$, relevante para el problema de Kepler en mecánica cuántica (átomo de hidrógeno).

Método del Punto Silla (Steepest Descent)

Para obtener el comportamiento asintótico de $H_n(x)$ cuando $n \to \infty$, se escribe la representación integral como $\oint e^{n\,\Phi(z)} dz$ con fase $\Phi(z) = 2x z/n - z^2/n - \log z$. El método del punto silla (o steepest descent) deforma el contorno para que pase por el punto silla donde $\Phi'(z_0) = 0$, siguiendo la dirección de descenso más pronunciado de $\Re(\Phi)$. La contribución dominante proviene de la vecindad de $z_0$, y el teorema de Cauchy garantiza que la deformación no altera el valor de la integral. Este método, de origen puramente complejo, es una de las herramientas más poderosas de la física matemática para el análisis asintótico.

Fórmula de Rodríguez en el Plano Complejo

La fórmula de Rodríguez para polinomios ortogonales se deriva de manera unificada usando la representación integral de Cauchy para la $n$-ésima derivada de una función analítica. Para los polinomios de Legendre $P_n(x)$:

$$P_n(x) = \frac{1}{2^n n!}\frac{d^n}{dx^n}(x^2-1)^n = \frac{1}{2\pi i}\oint_C \frac{(z^2-1)^n}{2^n (z-x)^{n+1}}\,dz$$

Esta es la fórmula de Schläfli, que proporciona una representación integral de contorno válida para $x$ complejo. Es la base de la fórmula de Laplace para $P_n(\cos\theta)$ y del análisis asintótico de los armónicos esféricos.

Ortogonalidad y Teorema del Residuo

Las relaciones de ortogonalidad $\int_a^b p_n(x)p_m(x)w(x)dx = 0$ para $n \neq m$ pueden demostrarse sistemáticamente insertando las representaciones de contorno de $p_n$ y $p_m$, intercambiando integrales, y evaluando la integral de contorno resultante mediante el teorema del residuo. Esta técnica, conocida como método de las funciones generatrices complejas, unifica el tratamiento de familias de polinomios aparentemente dispares.

Cuestionario

1. ¿Cuál es la definición integral de la función Gamma $\Gamma(z)$ para $\Re(z) > 0$?

2. ¿Cuál es la ecuación funcional fundamental de la función Gamma?

3. ¿Dónde se encuentran los polos simples de $\Gamma(z)$?

4. ¿Cuál es la relación entre la función Beta $B(p,q)$ y la función Gamma?

5. ¿Para qué valores de $s$ converge la serie $\zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s}$ originalmente definida por Riemann?

6. ¿Qué expresa el producto de Euler para $\zeta(s)$?

7. ¿Cuál es el valor de $\Gamma(1/2)$?

8. ¿Cuál es la fórmula de reflexión de Euler para la función Gamma?

9. ¿Qué representación integral proporciona la fórmula de Schläfli para los polinomios de Legendre $P_n(x)$?

10. ¿Cómo se conectan los coeficientes de Laurent de una función analítica con sus coeficientes de Fourier sobre la circunferencia unidad?