La fórmula de sumación de Poisson y parientes cercanosmatmat/ebook2017/V2017n02-ebook.pdfEsto es...
Transcript of La fórmula de sumación de Poisson y parientes cercanosmatmat/ebook2017/V2017n02-ebook.pdfEsto es...
MAT 2MATerials MATemàticsVersió per a e-book deltreball no. 2 del volum 2017www.mat.uab.cat/matmat
La fórmula de sumación de Poissony parientes cercanos
Fernando Chamizo Lorente
Dulcinea Raboso Paniagua
1. ¿Quién teme a las series que no con-vergen?
¿Recuerdas cómo te metían mie-
do en Cálculo I con ese teorema de
Riemann que afirmaba que las se-
ries condicionalmente convergentes
se podían reordenar para que la su-
ma fuera cualquier número? ¿Y qué
decir de las no convergentes? ¿Te acuerdas del vér-
tigo que te produjo 1 = 1− (1−1)− (1−1)−· · · =(1−1)+(1−1)+ · · · = 0? ¿No dejaste de confiar en
el profe de Física I cuando te dijo que la delta de
Dirac δ era una función que valía cero en todos los
puntos e infinito en el cero pero que integraba uno?
Pues, agárrate, la fórmula que resume la sumación
de Poisson es∞∑
n=−∞δ(x− n) =
∞∑n=−∞
e−2πinx. (1.1)
Los lectores que hayan superado el susto1· mere-
cen una explicación, primero intuitiva y después
con una prueba matemática de esas que nos dejan
satisfechos.
Si metemos el hacha de cortar infinitos en elsegundo miembro de (1.1), queda una humilde se-rie geométrica que podemos sumar y hasta dibujarporque es real como la vida misma
DN (x) =
N∑n=−N
e−2πinx =e−2πi (N+1) x − e2πiN x
e−2πi x − 1
=sin(π (2N + 1)x
)sin(π x)
= · · · -2 -1 1 2 · · ·(1.2)
El bulto grande que aparece en los enteros es
de altura 2N + 1 y de anchura media más o me-
nos 1/2N , entonces (1.1) luce bien visualmente:
-2 -1 1 2
hacia ∞
hacia ∞
hacia ∞
hacia ∞
hacia ∞
= lımN→∞
-2 -1 1 2
2N+1
1/2N
Y si integramos DN (x) contra una función que de-caiga bien, el resultado se aproxima por la sumade los valores de f en los enteros que es donde es-tán los bultos. Si confiamos en que N →∞ no noshace ninguna jugarreta, tenemos la forma habitualde la fórmula de sumación de Poisson
∞∑n=−∞
f(n) =
∞∑n=−∞
f(n)
con f(ξ) =
∫ ∞−∞
f(x) e−2πi ξ x dx . (1.3)
Esto es equivalente a (1.1) o lo que le da sentido.
La f es la famosa transformada de Fourier de f .Vale, la deducción de (1.3) como broma está
bien, pero ¿cuál es la demostración sin chapuzas?
Muy fácil: se define la función 1-periódica F (x) =∞∑
k=−∞f(x+ k), que bajo condiciones adecuadas de
regularidad debe ser igual a su serie de Fourier, portanto
F (x) =
∞∑n=−∞
(∫ 1
0
F (t) e−2πin t dt
)e2πinx
=
∞∑n=−∞
(∫ ∞−∞
f(t) e−2πin t dt
)e2πinx (1.4)
y basta tomar x = 0.
2. Identidades a tutiplén
Una vez que tenemos el martillo, busquemos
los clavos. La función más sencilla que sabemos
integrar es la exponencial. Pongámosle un menos
en el argumento para que decaiga, valor absolu-
to para que no se desmande en −∞ y una pizca
positiva de 2πα para decorar. Esto es, tomemos
f(x) = e−2π α |x|. El lado izquierdo de (1.3) es una
serie geométrica, bueno. . . dos, que sabemos su-
mar, lo que unido al cálculo de las integrales del
lado derecho da lugar a
eπα + e−πα
eπα − e−πα=α
π
∞∑n=−∞
1
n2 + α2. (2.1)
Una bella identidad que generaliza el bien cono-
cido π2/6 =∞∑n=1
n−2 del maestro sumador Euler.
¿Lo generaliza? Sí, se deja como ejercicio deducirlo
de (2.1) haciendo la pizca 2π α infinitesimal.
Nuestra prueba de (2.1) requiere α > 0 pero
Poisson nos ha dado más de lo esperado: por ex-
tensión analítica se cumple para todo α ∈ C con
α i 6∈ Z. Para α imaginario puro, (2.1) es el conoci-
do desarrollo de la cotangente, también del maes-
tro, y su desarrollo de Laurent permite calcular∞∑n=1
n−k cuando k es par.
Los valores absolutos tienen mala fama en el
cálculo infinitesimal por eso de que los pobrecillos
no son derivables. Teniendo en cuenta que los cua-
drados son positivos a más no poder, cambiamos
de variable y consideramos f(x) = e−2π αx2 . Ju-
gando con variaciones2· de∫ ∞−∞
e−x2dx =
√π ,
la fórmula (1.3) produce la identidad válida para
<(α) > 0
∞∑n=−∞
e−2π αn2=
1√2α
∞∑n=−∞
e−π n22α . (2.2)
Tomando α = 18π la segunda suma tiene términos
minúsculos para n 6= 0, eso da pie a que te burles de
tus amigos más crédulos con la fórmula falsa π =(12
15∑n=−15
e−n2/4
)2
cuyo error, de tamaño menor
que 10−16, es indetectable con una calculadora.
Tomemos ahora el intervalo I = [−t, t] con t >0 no entero. Si ft vale uno en I y 0 fuera, (1.3)
aplicado a ft nos da
2btc+ 1 = 2 t+∑n6=0
sen(2π n t)
π n(2.3)
con btc la parte entera de t. Pero ¡un momento! si t
es entero está claro que esto es falso porque la serie
se anula ¿Por qué 1.999 es tan diferente de 2 a efec-
tos de esta fórmula? Una manera de explicarlo es
que si t no es entero entonces se puede buscar una
g muy regular que aproxime muy bien a ft y tal
que las sumas que aparecen en (1.3) sean muy pa-
recidas para ft y g. La regularidad permite aplicar
el análisis de Fourier sin reparos en la demostra-
ción de (1.3). Por otro lado, si t es entero, la suma
del primer miembro es muy sensible a tomar fun-
ciones parecidas, por ejemplo, para ft+ε y ft−ε con
ε = 10−mucho la diferencia es de dos unidades. Otra
manera de verlo es notar que (2.3) no es más que el
desarrollo en serie de Fourier de la onda triangular
2btc − 2 t + 1 que, obedeciendo a la teoría, debe
converger al punto medio en las discontinuidades
de salto. La falta de regularidad está reflejada en la
ausencia de convergencia absoluta en la serie (2.3)
permitiendo que t = 1.999 y t = 2 den resultados
dispares. En principio es posible (pero no saludable
por las cuestiones de regularidad) considerar las se-
ries de Fourier como aplicaciones de la fórmula de
sumación de Poisson porque a fin de cuentas son
sumas sobre los enteros.
Una vez que hemos rozado el agridulce sabor de
la falta de convergencia absoluta, nos meteremos en
el barrizal de la absoluta no convergencia.
La integral∫ ∞−∞|x|−s cos(2π x) dx tiene senti-
do para cualquier 0 < s < 1 y se puede evaluarcomo 2 (2π)s−1 sin(π s/2) Γ(1 − s) donde Γ es la
función especial Γ(u) =
∫ ∞0
xu−1 e−x dx introdu-
cida por Euler (¡otra vez!) [Art64] para generalizar
el factorial: Γ(n) = (n − 1)! cuando n ∈ Z+. To-mando f(x) = |x|−s en (1.3) y cerrando los ojos seobtiene
∞∑n=−∞
|n|−s = 2 (2π)s−1 sin(π s
2
)Γ(1−s)
∞∑n=−∞
|n|s−1 .
(2.4)
Esto tiene muy mal aspecto, si una suma existe la
otra no y n = 0 da lugar a infinitos pero. . . había
una cosa llamada función ζ(s) que extendía analíti-
camente∞∑n=1
n−s y servía para estudiar los primos.
Si nos olvidamos de los n = 0 y de la convergencia,
lo anterior se reescribe como
ζ(s) = 2 (2π)s−1 sin(π s
2
)Γ(1− s) ζ(1− s) , (2.5)
que es la ecuación funcional correcta en forma asi-
métrica. Para el que esté ansioso por justificar este
misterio, está la bibliografía de (la vieja) guardia
[Gui41] [Ing90, III.4]. Un ingrediente principal es
que para s > 0, s 6= 1, se cumple
ζ(s) = lımN→∞
( N∑n=1
n−s − N1−s
1− s
), (2.6)
lo cual es obvio para s > 1 y permite cancelar
infinitos en el caso 0 < s < 1.
3. Más dimensiones y menos dimensio-nes
Extender (1.3) a más dimensiones es pura ru-
tina: se suma en ~n ∈ Zd y en la definición de f se
escribe ~ξ · ~x. Una vez que hemos subido la dimen-
sión la podemos bajar. No es broma del todo. Si
f tiene simetrías especiales, es posible agrupar va-
lores iguales dando lugar a una fórmula en menos
dimensiones aparentemente muy distinta.Por ejemplo, si d = 2 y f es radial, f(~n) =
g(n21 +n2
2) y g(n) aparecerá con multiplicidad r(n),el número de representaciones como suma de dos
cuadrados. Trabajando los detalles [IK04, §4.4], setiene una impresionante fórmula con función deBessel y todo:
∞∑n=1
r(n) g(n) = π
∫ ∞0
g(x) dx
+ π
∞∑n=1
r(n)
∫ ∞0
g(x) J0(2π√nx ) dx (3.1)
que es fundamental para estudiar el problema del
círculo de Gauss [LWL24].
Por si eso sabe a poco, generalicemos Zd con-
siderando un retículo en Rd. Esto es el conjunto
de combinaciones lineales enteras de d vectores li-
nealmente independientes o simplemente Λ = AZd
con A una matriz cuadrada no singular y se llama
retículo dual a Λ∗ =(A−1
)t Zd.
Z2 Λ = AZ2
Con la fórmula de cambio de variable tenemos
fácilmente la generalización d-dimensional en re-
tículos ∑~n∈Λ
f(~n) = |Λ|−1∑~n∈Λ∗
f(~n) (3.2)
donde |Λ| := det(A). Por ejemplo, si elegimos
f(~x) = e−2π α ‖~x‖2 , (2.2) se generaliza a∑~n∈Λ
e−2π α ‖~n‖2 = |Λ|−1 (2α)−d/2∑~n∈Λ∗
e−π ‖~n‖2
2α .
(3.3)
Vamos a ver una bella y rápida aplicación
[Sko02] de (3.3) al “kissing number problem” que
pregunta cuál es el máximo número de esferas que
pueden besar simultáneamente a otra, todas del
mismo tamaño. Aquí “besar” significa “tocar”. Sí,
está claro que besar no es lo mismo que tocar. . .
mejor no seguir por ahí.
En dimensiones 1 y 2 es sencillo calcular el nú-
mero de osculación o “kissing number” κ(d). Gráfi-
camente:
κ(1)=2 κ(2)=6
La dificultad aumenta en dimensión 3, y probar
κ(3) = 12 no es nada fácil. El problema es famoso
por el desacuerdo entre los matemáticos I. Newton
(sí, también físico) y D. Gregory, polémica3· que
no fue resuelta con rigor hasta mediados del siglo
XX.
Situando siete esferas como en la figura ante-
rior para d = 2, podemos colocar exactamente tres
esferas más en la parte superior que estén en con-
tacto con la esfera central. Si hacemos lo mismo
en la parte inferior, tenemos doce esferas tocando
la central que están simultáneamente en contacto
con otras. Esta manera de acomodar las esferas es
conocida como la configuración de Kepler.
Config. de Kepler Config. icosaedro
También podemos colocar las esferas besuquea-
doras en los vértices de un icosaedro regular pero, a
diferencia de lo que ocurre con la configuración de
Kepler, ahora las esferas se pueden mover libremen-
te sin llegar a tocar al resto. De hecho, el espacio
“extra” entre cualesquiera dos de las doce esferas
exteriores no es para nada despreciable y no está
claro si cabe una decimotercera apretujándolas un
poco.
Aquí nos restringimos al “lattice kissing num-
ber” `(d), que es lo mismo pero suponiendo que
los centros de las esferas están en un retículo. Uno
tendería a pensar que κ(d) = `(d) por eso de la
elegancia y el orden en matemáticas, de hecho los
seis valores conocidos4· de κ(d) coinciden con `(d)
pero se sabe indirectamente que κ(9) 6= `(9). Por
otro lado, está claro que κ(d) ≥ `(d).
Llamemos F (α) al lado izquierdo de (3.3) mul-
tiplicado por αd/2. Si hemos logrado poner ` esferas
de radio 1/2 con centros en Λ besando una esfera en
el origen, separando de la suma ‖~n‖ = 0, 1, mayor,
F (α) = αd/2(1 + ` e−2π α +
∑‖~n‖>1
e−2π α ‖~n‖2) .(3.4)
Lo crucial es que (3.3) asegura que F (α) es crecien-te. Si α > d
4π ,
0 ≤ F ′(α) =d
2αd/2−1 + αd/2−1 ` e−2π α
(d2− 2π α
)+ cosas negativas . (3.5)
Eligiendo de manera óptima α = d+24π , se concluye
`(d) ≤ d
2ed/2+1 . (3.6)
Esto no está nada mal si se compara con los re-
sultados conocidos cuando d → ∞. Por cierto, un
problema relacionado es el empaquetamiento de es-
feras en dimensiones altas y las mejores cotas su-
periores se obtienen con procedimientos parecidos
[PZ04] [CE03].
4. Sumas parciales con giros inespera-dos
En teoría analítica de números a menudo apa-
rece el problema de conseguir acotaciones no tri-
viales de sumas trigonométricas. Por ejemplo, al
escribir una función como superposición de ondas,
g(x) e2πi f(x) con g la amplitud y f la fase, si que-
remos estudiar el promedio de la función en los
entreros a, a+ 1, a+ 2, . . . b, tras sumar por partes
(separando la amplitud no oscilatoria), nos enfren-
tamos a sumas del tipo
b∑n=a
e2πi f(n) . (4.1)
El método de van der Corput [GK91] se basa
en que las sumas trigonométricas se pueden apro-
ximar por integrales y que estas a su vez pueden
ser estimadas. La acotación más básica permite ob-
tener acotaciones no triviales en rangos en los que
la derivada segunda de la fase es moderadamen-
te pequeña y la sumación de Poisson es un arma
fundamental para transformar cualquier suma su-
ficientemente regular en una suma de integrales os-
cilatorias.
Si en (1.3) consideramos e2πi f(x) multiplicada
por la función característica de [a, b] se llegaría a
b∑n=a
e2πi f(n) “=”∞∑
n=−∞
∫ b
ae2πi (f(x)−nx) dx .
(4.2)
Las comillas se deben a que se requiere cierta regu-
laridad y la discontinuidad no gusta mucho, siem-
pre puede llevar a complicaciones técnicas. Aunque
tampoco estamos tan lejos de la “verdad” si tene-
mos cierto control sobre alguna derivada de f en
[a, b]. Por ejemplo, si f ′ es monótona y |f ′| ≤ 1/2
tendríamosb∑
n=a
e2πi f(n) =
∫ b
ae2πi f(x) dx+ error , (4.3)
donde error es menor que cierta constante absoluta,
o si f satisface f ′′ ≥ λ > 0, entonces∣∣∣∣∣b∑
n=a
e2πi f(n)
∣∣∣∣∣ ≤ C (f ′(b)− f ′(a) + 1)λ−1/2 (4.4)
para cierta constante C.
Un ejemplo curioso es el que ocu-
rre al dibujar los números complejos
zN =
N∑n=1
eπi√n , (4.5)
cuyo resultado revela una estructura de espiral que
es inesperada si se desconoce el método de van der
Corput.
Escribiendo S(N ;α) =N∑n=1
e2πiα√n , con α >
0, la derivada de α√x decrece, de hecho para α
razonablemente pequeño después de algunos tér-
minos esta derivada es menor que uno, y según las
ideas del método de van der Corput podemos es-
perar una buena aproximación mediante∫ N
1e2πiα
√x dx =
√N
π αe2πiα
√N −πi/2 + error
(4.6)
donde, como antes, el error es menor que una cons-
tante. Esto sugiere que al tener en cuenta úni-
camente el término principal, la sucesión finita
{S(n;α)}Nn=1 se debería parecer a una espiral de
Arquímedes
t
2π2 α2(sin t,− cos t) con t ∈ [1, 2π α
√N ] .
La sorpresa es que en los dibujos muy rara vez se
observa la estructura de espiral. Por ejemplo, para
N = 1000 se obtienen las siguientes figuras para
los valores indicados de α,
α = 1 α = 65/64 α = 1.3
Esto es en cierto sentido una ilusión óptica,
puesto que si conectamos puntos consecutivos en
la sucesión con líneas, las espirales de Arquímedes
reaparecen y, como predice (4.6), un valor mayor
de α da una menor distancia entre giros sucesivos.
Los patrones que se observan llevan a consideracio-
nes aritméticas y se estudian con detalle para α2
entero en [CR15b].
5. Bienvenido al mundo modular
La variable compleja es muy rígida, en cuan-
to pones muchas condiciones te quedas sin funcio-
nes, como las compañías teatrales mediocres. Por
ejemplo, si uno busca las funciones F holomorfas
en C (enteras) que tengan dos periodos, digamos
F (z) = F (z + 1) y F (z) = F (z + i), sólo hay una
muy aburrida una vez que se ha especificado F (0),
la función constante F (z) = F (0).
Cuando uno considera funciones holomorfas en
el semiplano superior H = {z : =(z) > 0} e im-
pone condiciones relacionadas con ciertos grupos
fuchsianos, la teoría adquiere una riqueza y pro-
fundidad increíbles. Es la teoría de las formas mo-
dulares y si no has oído mencionarla, no eres de
este mundo (matemático).
Esencialmente una función modular5· es una
función holomorfa en H invariante por las trans-
formaciones z 7→ z + 1, z 7→ −1/z y se dice forma
modular si además es holomorfa en el infinito. La
traslación y la inversión (negativa) generan el gru-
po SL2(Z) pero las necesidades aritméticas llevan
a considerar distintos subgrupos, entre los que des-
tacan los grupos de congruencias Γ0(N).
En este contexto, se puede probar por ejemplo
que sólo hay una función holomorfa en H que cum-
ple
F (z) = F (z + 1), F (z) =−1
4 z2F(− 1
4 z
), (5.1)
una vez que especificamos el valor de F (i∞) :=
lım=(z)→+∞
F (z), supuesto finito al igual que
lım=(z)→0+
=(z)2017 F (z), donde puedes cambiar 2017
por cualquier N > 2, siempre da cero.
Vamos a ver cómo deducir de ello un resultado
clásico pero complicado de teoría de números gra-
cias a (1.3). Definamos F (z) como el lado izquierdo
de (2.2) elevado a la cuarta potencia con α = −i z.
Esto es,
F (z) :=( ∞∑n=−∞
e2πin2 z)4
=
∞∑n=0
r4(n) e2πin z,
(5.2)
donde r4(n) es el número de representaciones como
suma de cuatro cuadrados. Gracias a (2.2), la fun-
ción F satisface (5.1) y proclamamos, según lo an-
terior, que es la única con F (i∞) = 1. El alumno
genial y desconfiado de la última fila, un tal Ei-
senstein sin escalinata y con “s” y “e”, podría dar
un ejemplo desconcertante
G(z) =1
π2
∑m
∑n
(m,n)6=(0,0)
( 4
(m+ 4n z)2− 1
(m+ n z)2
).
(5.3)
A pesar de lo aparatoso de la fórmula es muy senci-
llo (sí, de verdad) comprobar que G verifica (5.1) si
no te distraes preocupándote con la convergencia.
Por otro lado, hallar G(i∞) se reduce a considerar
la contribución de n = 0 que es π−2∑m 6=0
3m−2 = 1.
Hay que admitir entonces la identidad delirante
F (z) = G(z) por la unicidad de las funciones que
cumplen (5.1).
Todavía hay más, con (1.3) aplicado a f(n) =
(n + w)−2 y un poco de pericia con las integrales
(de Cauchy, si puede ser), se tiene
∞∑m=−∞
1
(m+ w)2= −4π2
∞∑m=1
m e2πimw, (5.4)
que es válida para w ∈ H tal como está y exten-sible a −w ∈ H cambiando e2πimw por e−2πimw.Es posible también obtener esta fórmula como con-secuencia indirecta de (2.1). Sustituyendo en G yhaciendo los cálculos, separando enes positivas ne-gativas y nulas, la identidad delirante se convierteen
∞∑n=0
r4(n) e2πin z
= 1 + 8
∞∑n=1
∞∑m=1
(m e2πimnz − 4m e2πi 4mnz
). (5.5)
Si ahora comparamos los coeficientes de e2πin z en
ambos miembros, se tiene el golpe final, la fórmula
apabullante:
r4(n) = 8∑
4-d, d|n
d . (5.6)
Una identidad limpia, breve y profunda. Por ejem-
plo, sin escribir ni una cuaterna de cuadrados, se
deduce
r4(2017) = 8 (1 + 2017) = 16144 . (5.7)
¿Y qué ocurre para otro número de cuadrados?
¿Hay fórmulas tan limpias? Sólo en algunos casos y
eso está relacionado con las propiedades de las for-
mas modulares [Gro85] [Har20] [Ran77] [Iwa97]. El
lector abrumado encontrará en [Ven70] una prueba
elemental de (5.6) y de otras fórmulas relacionadas.
Las formas modulares han estado tradicional-
mente ligadas a la variable compleja. Algunas de
ellas aparecen naturalmente al aplicar (1.3) para
probar la ecuación funcional de ciertos objetos re-
levantes en aritmética, como la función ζ. En este
contexto, H. Maass se percató de que al tratar cuer-
pos algebraicos reales era natural introducir funcio-
nes con simetrías similares pero no holomorfas6·.
Si perdemos la variable compleja ¿qué se salva?
Una teoría endiabladamente complicada pero útil
y bella.
El semiplano H fue introducido por Poincaré
como modelo de geometría hiperbólica (curvatu-
ra −1) dotado de la métrica ds2 = y−2(dx2 +dy2
).
Esto significa que arriba las cosas están más cer-
ca de lo que parecen y abajo lo están menos. Las
formas modulares no holomorfas invariantes por la
acción de un grupo Γ, se pueden entender como
funciones en la superficie (de Riemann) obtenida
al identificar los bordes de una región llamada do-
minio fundamental. Por ejemplo, para SL2(Z) se
tiene una esfera en la que un punto se lleva a infi-
nito.
−→ −→
De hecho siempre se puede elegir como dominio
fundamental un polígono hiperbólico cuyos lados
son arcos de semicircunferencias o rectas verticales
(los dos tipos de geodésicas).
Una diferencia notable entre los dos primeros domi-
nios, correspondientes a SL2(Z) y Γ0(4), y el terce-
ro es que este último corresponde a un grupo com-
pacto que no tiene cúspides (vértices en la recta
real o en el infinito). Aunque el caso compacto es
menos común en teoría de números desde el pun-
to de vista de sus posibles aplicaciones aritméticas,
tiene una importante ventaja en relación a un te-
ma básico en el análisis armónico, la expresión de
una función como superposición de “tonos puros”:
autofunciones de un operador.
Un ejemplo familiar son las autofunciones e2πinx
del operador f 7→ f ′′, efectivamente las series de
Fourier en R/Z de toda la vida. En H hay un
operador natural que es invariante por las trans-
formaciones del grupo, con autovalores 0 = λ0 <
λ1 ≤ λ2 ≤ . . . correspondientes7· a autofunciones
normalizadas {uj(z)}∞j=0 (no holomorfas) llamadas
formas de Maass.
La primera autofunción es la función constante
u0(z) = |Γ\H|−1/2. Las imágenes representan las
primeras formas de Maass no triviales para Γ =
SL2(Z).
u1 u2 u3
Aunque para obtener el desarrollo espectral de
cualquier función automorfa, se deben considerar
otras funciones que, a diferencia de las formas de
Maass, no son de cuadrado integrable pero que si-
guen siendo autofunciones. Hablamos de las series
de Eisenstein (el alumno listo de última fila reapa-
rece) de las que se dice que pertenecen al espectro
continuo y que están relacionadas con las cúspi-
des de Γ, sus puntos del infinito. De este modo, al
sustituir R por H, obtenemos
f(z) =∑
aj uj(z)︸ ︷︷ ︸espectro discreto
+
(contribución del
espectro continuo
).
(5.8)
¿Y la ventaja de trabajar con grupos compactos?
pues en ese caso la parte continua (con aspecto
realmente feo) desaparece por lo que nos quedamos
únicamente con la parte discreta del desarrollo.
Aún más, si bien el mundo holomorfo nos ofrece
identidades falsas (o aproximadas) ligadas al desa-
rrollo de Fourier de formas modulares clásicas, el
no holomorfo no iba a ser menos. Por ejemplo,
∞∑n=1
r(n) r(3n+ 2)√n e−(log(n)/4)2 (5.9)
está muy cerca de 72 e9√π . De hecho, su error
relativo a pesar de que no es cero es menor que
3 · 10−7 y, por increíble que parezca, la precisión
del error en esta y otras aproximaciones similares
depende del tamaño de ciertos autovalores.
Para no dejar de sorprendernos, también tene-
mos una fórmula de sumación de Poisson en es-
te tipo de superficies con inmediatas consecuencias
aritméticas: la fórmula de Kuznetsov [Kuz80] que
da una interpretación en términos de autovalores
de la distribución del inverso multiplicativo cuan-
do variamos el módulo (publicidad poco sublimi-
nal: en [CR15a] dimos una prueba de la fórmula
de Kuznetsov que hasta un niño8· podría seguir).
Los últimos cuarenta años han mostrado que lo po-
co que sabemos de autovalores y autofunciones es
suficiente para decir cosas que no sabíamos de la
distribución de los inversos y la teoría espectral en
este contexto se ha desarrollado de la mano de la
teoría de números [Rab15].
Así podemos ver cómo distintas áreas de las
matemáticas conviven en estrecha relación, aunque
por supuesto lo anterior solo muestra la punta de
un iceberg, e icebergs hay muchos (como veremos
en la siguiente sección).
6. Análisis, aritmética y geometría
Comencemos respetando el orden alfabético
por el análisis. La ecuación de ondas en la recta
real es
utt = uxx con u = u(x, t) (6.1)
donde los subíndices indican derivadas. ¿De dón-
de sale esta ecuación? ¿Por qué es una ecuación
“de ondas”? La manera más rápida de responder
usando solo conceptos físicos básicos y trampas, a
partes iguales, es que viene de que la aceleración es
derivar dos veces con respecto al tiempo y la fuer-
za por unidad de longitud que opone una cuerda
elástica a ser deformada es proporcional a la cur-
vatura que queramos imponerle con una constante
de proporcionalidad dada por la tensión, y de aquí
se saca utt = uxx, escribiendo a = F/m con uni-
dades llamativas. ¿Y las trampas? Si a tu profe de
geometría diferencial le dices que la curvatura de
la gráfica de y = f(x) es f ′′(x), te mirará con ca-
ra de −273.15◦C (cero absoluto) pero si la gráfica
es muy plana (f ′ ≈ 0) entonces ambas cantidades
son muy similares, además tampoco la tensión es
muy constante sin esta hipótesis. Para los más exi-
gentes, un siglo después, el electromagnetismo de
J.C. Maxwell ofreció un contexto físico e importan-
te (¿qué harías sin tu móvil, tu WiFi, tu TV. . . ?)
donde la ecuación de ondas es exacta [Max65]. En
cualquier caso, piensa en (6.1) como un modelo pa-
ra cosas que oscilan como una cuerda. La ecuación
(6.1) en realidad tiene una solución general fácil9·
dada por la fórmula u(x, t) = f(x− t) + g(x+ t) y
por tanto toda onda solución se escribe como suma
de una que viaja a la derecha y otra a la izquierda,
ambas con velocidad uno.
El gran mantra del análisis armónico es que to-
da onda es superposición de ondas estacionarias.
Estas ondas especiales son, como su nombre indi-
ca, las que de algún modo no se mueven: para cada
punto fijado solo vemos oscilaciones arriba y abajo
en el tiempo de las de toda la vida, con un seno
o un coseno. Si optamos por las segundas, mate-
máticamente son las soluciones de (6.1) de la for-
ma u(x, t) = cos(ω t)φ(x) y podemos limitarnos
a ω ≥ 0 por la simetría. Si no imponemos nin-
guna condición de contorno, (6.1) nos lleva a que
φ(x) debe ser una combinación lineal de sen(ω x)
y cos(ω x). El lector avispado notará la parado-
ja que plantea que las ondas estacionarias no se
muevan y la velocidad uno, antes mencionada, de
las soluciones de (6.1). Pensémoslo con un ejemplo:
la onda estacionaria u(x, t) = cos(2 t) cos(2x) no
se mueve, pero se descompone en la semisuma de
cos(2 (x − t)
)y cos
(2 (x + t)
)que se mueven con
velocidad uno hacia la derecha e izquierda, respec-
tivamente. Mientras se inventa el papel animado,
aquí van unos pocos fotogramas superpuestos:
1 2 3 4 5 6
-1
-0.5
0.5
1
1 2 3 4 5 6
-1
-0.5
0.5
1
cos(2 t) cos(2 x) para t = 0, 1/2, 1, 3/2 cos(2 (x− t)
)para t = 0, 1/2, 1, 3/2
De alguna forma la interferencia entre las ondas
que viajan a derecha y a izquierda es tan perfecta
que el resultado es la inmovilidad.
Pensemos ahora que nuestras oscilaciones tie-
nen lugar en un anillo muy fino representado por
la circunferencia unidad, lo que equivale a poner
condiciones 1-periódicas en la x en (6.1) identifi-
cando 2π x con el ángulo. Para que la oscilación
no se rompa, se necesita que ω sea un múltiplo en-
tero de 2π, la longitud de la circunferencia.
ω = 56π ω = 59.4π
Debe conceder el lector que ω ∈ 2πN es una apari-
ción modesta pero honesta de la aritmética. Antes
de dar paso a la geometría, consideremos el siguien-
te monstruo no convergente
u(x, t) = 1 + 2∞∑n=1
cos(2π n t) cos(2π n (x− x0)
).
(6.2)
Formalmente resuelve (6.1) y según (1.1), con x 7→x0 − x y tomando partes reales, para t = 0 es una
delta de Dirac en la circunferencia sobre el punto de
ángulo 2π x0. Como las soluciones de la ecuación
de ondas tienen velocidad uno, esta delta de Dirac
reaparecerá en el mismo lugar para cada tiempo
entero, la redondez de la circunferencia hace indi-
ferente que las ondas viajen a derecha o izquier-
da. Si tomamos la traza, esto es, si integramos en
0 ≤ x = x0 ≤ 1 (sí, ya sabemos que bastaría to-
mar x = x0), se concluye que∞∑
n=−∞cos(2π n t) es
una delta de Dirac en cada entero. Hemos hecho
un negocio redondo, nunca mejor dicho: usando la
ecuación de ondas hemos deducido (1.1) de (1.1).
Es un buen momento para pasar a la geometría.
En general, la ecuación de ondas en una varie-
dad riemanniana M (un objeto geométrico en que
sepamos medir ángulos y distancias) es:
utt = ∆u con u = u(x, t) (6.3)
donde ∆u, denota el laplaciano de u, paraM = Rn
es la suma de las derivadas segundas respecto a las
coordenadas de x. El nombre no deriva de un santo-
ral rural sino del insigne P. Laplace que lo usó con
especial maestría en relación con la gravitación. Si
queremos saber por ejemplo cómo oscila una mem-
brana, debemos usar (6.3). Aquí nos centraremos
en el caso en que M es compacta, para fijar ideas
uno puede pensar en el caso de una superficie ce-
rrada como una esfera o un toro.
Las ondas estacionarias u(x, t) = cos(ω t)φ(~x)
para ser solución de (6.3) deben satisfacer −∆φ =
ω2 φ, esto es, que ω2 sea un autovalor de −∆ con
autofunción φ. Se sabe que tales autovalores for-
man un conjunto discreto 0 = λ0 < λ1 < λ2 < . . . .
El primo n dimensional del monstruo (6.2) es10·
u(x, t) =
∞∑k=0
cos(t√λk)∑m
φkm(x)φkm(x0) .
(6.4)donde φkm(x) son las autofunciones correspondien-tes al autovalor λk ortonormalizadas, esto significa∫Mφkm φk′m′ = 0 excepto si (k,m) = (k′,m′), en
cuyo caso da uno. ¿De dónde sale el parentesco en-
tre (6.4) y (6.2)? Cuando ∆ es la derivada segundade Cálculo I, ya habíamos visto que las φkm so-lo podían ser sen(2π k t) y cos(2π k t), que corres-ponden a
√λk = 2π k, si las ortonormalizamos,
adquieren unos factores de fantasía:
φk1(x) =√
2 sen(2π k t)
y φk2(x) =√
2 cos(2π k t) , (6.5)
con la excepción φ02(x) = 1 .
Haciendo las cuentas, con esta elección, (6.2) es
(6.4) para la circunferencia.
La similitud va más allá, el mantra toda onda
es superposición de ondas estacionarias (llamado
completitud espectral entre los expertos) asegura
que (6.4) para t = 0, es una delta de Dirac en el
punto x0.
Antes de dar el toque de gracia, dejaremos des-
cansar al lector con un ejemplo. ¿Qué ocurre si ti-
ramos un fuerte pellizco a una pelota elástica? En
el primer instante el estiramiento rebotará violen-
tamente y después se propagarán las oscilaciones
a lo largo de la esfera. Una delta de Dirac es un
pellizco superbestia y entonces suena bien que si
nos quedamos solo con los primeros términos de
(6.4) tendremos un pellizco solo bestia11·. En el
caso de la esfera se conoce que λk = k (k+ 1) y las
φkm son unas funciones raritas pero importantes
llamadas armónicos esféricos. Enseñándole a nues-
tro ordenador cómo calcularlas y limitando la suma
a k ≤ 20 se tienen las siguiente figuras, en las que
se ha escalado el tamaño del pellizco por motivos
estéticos:
t = 0 t = 0.1257 t = 0.2513 t = 0.3770
El rebote es, como preveíamos, muy rápido. Des-
pués de que las ondas se propagan hasta el polo
sur y causan un chichón allí, vuelven. Veamos eso
últimos instantes en que el pellizco vuelve a repro-
ducirse:
t = 4.7752 t = 5.2779 t = 6.0319 t = 6.2832
¡Visto y no visto! El pellizco asoma para t = 6.0319
con fuerza menor que al principio pero el rebote es
tan rápido que no es fácil capturar numéricamen-
te el momento de máxima altura. Más adelante le
tendremos que pedir algunas disculpas a nuestro
ordenador.Volvamos a (6.4) en total generalidad pero con-
servando en mente las figuras anteriores. Si inte-
gramos como antes en x = x0 (tomamos la traza),surge una conjetura natural y alucinante:
∞∑k=0
cos(t√λk)
?=∑
δ(t− tiempos de retorno de las ondas) . (6.6)
En el caso de la circunferencia, esto cuadra y en elde la esfera el valor t = 6.0319 para el que el pe-llizco reaparece con altura notable (aunque menorque la original) difiere solo en un 4 % de 2π que eslo que se tarda en volver al polo norte caminandoa velocidad uno por un meridiano, y ese 4 % suenaa error en la aproximación de pellizco por delta.¿Por qué por un meridiano? Porque es localmenteel camino más corto, el que se sigue por inercia sinhacer fuerza. Esto del camino más corto es lo quese llama geodésica en geometría. Así la conjetura
se vuelve más alucinante todavía:
∞∑k=0
cos(t√λk)
?=∑
δ(t− longitudes de geodésicas cerradas) (6.7)
donde se entiende que las geodésicas se pueden re-correr varias veces o en sentido contrario, lo quese indica con números negativos para conservar laparidad. ¡Increíble, una fórmula de sumación dePoisson que relaciona autovalores (objetos del aná-lisis) con geodésicas (objetos de la geometría)! De-masiado increíble, tanto que la conjetura anterior,así escrita, es falsa. Sin embargo, en 1956 antes deque nadie considerase estos razonamientos A. Sel-berg consiguió una fórmula de este tipo para cier-tas superficies de curvatura constante menos uno(no compactas en general). La fórmula, cuya prue-ba requiere la teoría de formas de Maass de la quehemos hablado, es complicadísima y ocupa variaslíneas en el original [Sel56]. Si acumulamos unos
renglones en puntos suspensivos, su forma es:∞∑k=1
f(√
λk − 1/4)
=1
2
∑n
∞∑k=1
`n f(k `n)
sinh(k `n/2)+ · · ·
(6.8)
donde `n son las longitudes de las geodésicas cerra-
das recorridas solo una vez. Aunque parezca menti-
ra, la fórmula tiene unas consecuencias aritméticas
nada triviales, como la similitud formal entre el
conjunto{
e`n}y el de los primos de toda la vida
pero con una función ζ que satisface, casi trivial-
mente, la hipótesis de Riemann. Y por si eso fuera
poco, en uno de los casos más simples las longitu-
des están relacionadas con las llamadas soluciones
fundamentales de la ecuación de Pell y sus multipli-
cidades, con números de clases de cuerpos cuadrá-
ticos reales; dos objetos realmente elusivos sobre
los que la fórmula de la traza da una información
sorprendente [Sar82].
Comparando la complejidad de (6.8) con la sim-
plicidad de (6.7) está claro que la conjetura fue de-
masiado optimista. De hecho ya en el caso de la
esfera encontramos dificultades: los√λk se pare-
cen muchísimo a k + 1/2 y entonces para t muy
cercano a 2π, el lado derecho de (6.7) debería pa-
recerse más a −δ(t− 2π) que a δ(t− 2π), hay un
signo que no cuadra. Esto se refleja en el numérico
y la figura correspondiente a t = 6.2832 ≈ 2π ocul-
ta un pellizco muy similar al de t = 0 pero hacia
dentro ¡el pico para t = 6.0319 no era tan impor-
tante como lo que no se ve después! Si vamos más
allá de la esfera, en general una geometría compli-
cada de las geodésicas causa interferencias mutuas
difíciles de cuantificar.
¿Se puede salvar algo? Los trabajos más re-
levantes a este respecto por orden histórico son
[CdV72] (hay también otros del mismo autor)
[Cha74] y [DG75]. En resumen, lo que se puede afir-
mar con bastante generalidad es que∞∑k=0
cos(t√λk)
no es singular fuera de los múltiplos enteros de las
longitudes de las geodésicas cerradas. Esto implica
que si f tiene soporte que no contiene estas longi-
tudes prohibidas,∞∑k=0
∫f(t) cos
(t√λk)no puede
depender de valores sueltos de f . Es algo así como
una versión un poco debilitada de que cuando en la
fórmula de Poisson clásica tomamos f con sopor-
te disjunto con Z, entonces∑f(n) es nulo. Bajo
hipótesis geométricas sobre las geodésicas cerradas
se puede ser mucho más concreto y llegar a algo
que sin ser una fórmula exacta como la original de
Poisson o (6.8), se acerca a ella.
7. Abel, Meyer, Barcelona
El Universo, la persona amada, las fórmulas bo-
nitas. . . siempre tendemos a pensar que todas las
cosas importantes son únicas, excepcionales ¿si no
por qué nos iban a gustar?
Bajo esta filosofía, consideremos dos ristras cre-
cientes de números reales {xn}∞n=−∞ y {yn}∞n=−∞
que no se amontonen, en el sentido de que haya una
cota inferior positiva para xn+1 − xn y yn+1 − yn.Ahora imponemos que para funciones buenas se
cumpla el análogo de (1.3)∞∑
n=−∞f(xn) =
∞∑n=−∞
f(yn) . (7.1)
En esta situación necesariamente {xn} = {yn} =
Z, no hay una fórmula de sumación de Poisson exó-
tica. Este resultado se prueba en [Cór88] generali-
zado a Rd, esto es, (3.2) es única.Veamos la demostración completa en dimen-
sión 1. Si escogemos f(x) = e−2π αx2 se sigue (2.2)con xn e yn, lo que implica 0 ∈ {xn} toman-do α → ∞, digamos x0 = 0. Consideremos µ =
xj0+1−xj0 la distancia mínima12· entre elementos
de {xn} y la función “tienda de campaña” φ(x) =
max(1− |x| , 0
). Su transformada de Fourier cum-
ple φ(0) = 1 y φ ≥ 0 (solo se anula en un conjuntodiscreto). Eligiendo f(x) = e2πi t x φ
((x − xj)/µ
)en (7.1), el primer miembro únicamente tiene elsumando no nulo f(xj) = e2πi t xj y simplificandoun poco se concluye
1 = µ
∞∑n=−∞
e−2πi xj yn φ(µ (yn − t)
)para todo t ∈ R. (7.2)
Recordando φ(0) = 1, φ ≥ 0, la única posibilidadpara que el caso j = 0, que da 1 = µ
∑nφ(µ (yn −
t)), sea coherente con j 6= 0 es que e−2πixj yn = 1.
En particular µ yn = (xj0+1 − xj0) yn es entero.Además t = yn en (7.2) implica µ ≤ 1. Integrando(7.2) sobre −1/2 ≤ t ≤ 1/2 se sigue
1 =
∞∑n=−∞
∫In
φ con In =[µ yn −
µ
2, µ yn +
µ
2
).
(7.3)
Los intervalos In son disjuntos pero por otra parte∫ ∞−∞
φ = φ(0) = 1 fuerza a que⋃In = R y enton-
ces necesariamente µ = 1 y {µ yn} = Z. Apelandoa (1.3), {yn} = Z implica {xn} = Z.
El caso d-dimensional es técnicamente más
complicado pero sigue el mismo argumento. En
[Cór89], A. Córdoba afirma que esto da “confir-
mación a un importante paradigma cristalográfico:
cuando se ven picos de Bragg en el espectro de Fou-
rier, se debe tener una estructura periódica, esto es,
un cristal” y aborda la situación más complicada
en que se admiten coeficientes de cierto tipo en la
versión d-dimensional de (7.1). Bajo ciertas condi-
ciones, se deduce que la única posibilidad es sumar
las fórmulas (3.2) para varios retículos. Varios au-
tores han contribuido a este problema pero su so-
lución completa en dimensió 1, debida a N. Lev y
A. Olevskii [LO15], es muy reciente. En pocas pa-
labras, siempre que se respete la condición de que
haya una distancia mínima positiva, las únicas po-
sibilidades para una variante (7.1) con coeficientes
multiplicando a f(xn) y f(yn), son combinaciones
lineales finitas de (1.3) escalando o trasladando los
enteros. En términos cristalográficos, los picos ais-
lados siempre representan superposiciones finitas
de estructuras periódicas.
Para Rd, d > 1, en [LO15] se requiere cierta
positividad de los coeficientes, por lo cual toda-
vía queda algún hueco para la investigación pero
aparte de eso, una vez que con mucha generalidad
se conoce que (1.3) es única ¿qué hacer ahora? El
multiverso, los amantes, los contraejemplos. . . to-
do ello también tiene su público, igual que atreverse
con la no unicidad de (1.3). ¿Por qué no jugar a las
dos cosas? Lev y Olevskii demuestran en [LO16]
(ver también [Kol16]) que hay dos conjuntos dis-
cretos {xn}, {yn} que no contienen ninguna copia
escalada de los enteros tales que (7.1) se cumple
con ciertos coeficientes. ¿Dónde está el truco? Un
conjunto discreto es el que solo consta de puntos
aislados y en contra de lo que parece intuitivo eso
no implica que sus elementos guarden una distan-
cia de seguridad. Por ejemplo, las sumas parciales
de la serie armónica [Gas17] forman un conjunto
discreto pero no hay una distancia mínima positi-
va.
Seguramente no hay que informar mucho a nin-
gún lector sobre el premio Abel. Después de la Me-
dalla Fields, es el premio más prestigioso en ma-
temáticas, otorgado anualmente por el gobierno de
Noruega desde 2003 a un reconocido matemático.
Para el que no se deje llevar por las veleidades so-
ciales del prestigio y el reconocimiento, también es-
tá una dotación económica que da miedo escribir.
Unas semanas antes de escribir este artículo se ha
anunciado que el galardonado en 2017 es Y. Meyer,
un matemático bien conocido particularmente por
su pionera contribución a las wavelets.
Poco antes de emocionarse por recibir el pre-
mio, estaba muy emocionado por un resultado que
había obtenido en relación con el problema que nos
ocupa. No es broma, leemos en [Mey16]: “En es-
te artículo se resuelve un problema importante de
análisis armónico: ¿Es única la fórmula de suma-
ción de Poisson o pertenece a una clase más am-
plia? Lo segundo es cierto. El método que se usa
para probar esta afirmación es sorprendente. Nues-
tras nuevas fórmulas de Poisson estaban ocultas
en un antiguo y casi olvidado artículo de 1959 de
A.P. Guinand. El papel desempeñado por la teoría
de números en este tema es fascinante”. En otro lu-
gar, escribe “Estos patrones eran totalmente ines-
perados y podrían abrir nuevas puertas en cristalo-
grafía”. Quizá es pertinente recordar la temprana
contribución de Meyer a los cuasicristales, estruc-
turas con gran sabor aritmético que tienen una pe-
riodicidad solo aproximada.La base de esta “nueva” fórmula radica en el
análogo de (3.1) cuando uno generaliza el problemadel círculo de Gauss a esferas. En ese caso, en vezde funciones de Bessel, se vuelven a obtener senosy cosenos y con ello una relación con la transfor-mada de Fourier unidimensional [IK04, Cor. 4.8].Aplicando una traslación, Meyer deduce la fórmu-la
∑~n∈Z3
e2πi~β·~n
‖~n+ ~α‖f(‖~n+ ~α‖
)= −i e−2πi ~α·
~β∑~n∈Z3
e−2πi ~α·~n∥∥~n+ ~β∥∥ f(∥∥~n+ ~β
∥∥)con ~α, ~β 6∈ Z3 (7.4)
válida para f impar, si eso molesta, basta cambiar
las f(x) por f(x) − f(−x), que es la forma que
aparece en [Mey16]. A pesar del sumatorio triple,
es una fórmula unidimensional, los {xn} y {yn} an-teriores son los valores tomados por ‖~n+ ~α‖ y por∥∥~n+ ~β
∥∥, que son números reales. Lo que Meyer re-
salta en su resultado principal es que, eligiendo ~α
y ~β, se consigue que cualquier subconjunto finito
de {xn} y {yn} sea linealmente independiente so-
bre los racionales. Esto significa que {xn} y {yn}son radicalmente distintos de los enteros o de una
unión de copias escaladas de ellos.
Muy bien, con esto tenemos a Meyer y a
Abel ¿y Barcelona? En enero de 2017, el Institut
d’Estudis Catalans albergó el evento “Abel in Bar-
celona” bajo los auspicios de la Societat Catalana
de Matemàtiques. Es allí donde tuvo lugar la úl-
tima reunión del comité cuya tarea fue seleccionar
el galardonado con el premio Abel 2017. Un gran
premio desde una bella ciudad.
8. Historias e historietas
Siméon Denis Poisson (1781–1840) es ese señorde la imagen de la primera página. Algunas de suscontribuciones que aparecen en los grados de físicao matemáticas son:
Pr(θ) =1− r2
1− 2 r cos θ + r2, −∆F = 4π ρ ,
{f, g} =∂f
∂qi∂g
∂pi− ∂g
∂qi∂f
∂pi, P (X = k) =
λk e−λ
k!.
Fue un investigador y profesor infatigable que vivió
los turbulentos años de la restauración monárquica
tras la revolución francesa. En la vida de Galois
desempeña un papel, posiblemente injusto, del lado
de los malos de la película por no haber apreciado
su famosa memoria y recomendar a la Academia de
Ciencias que no fuera admitida13·. El nombre de
Poisson es uno de los 72 de ilustres franceses que
aparece en la torre Eiffel. No es para menos, sus
contribuciones a la física y a las matemáticas son
ciertamente notables. Su producción es inmensa, se
estima que escribió entre 300 y 400 trabajos.
En su trabajo Sur les intégrales définies et sur
la sommation des séries, en la página 451 del vo-
lumen de 1823 del Journal de l’École Royale Poly-
technique, encontramos
En caligrafía moderna LATEXiana,
1
2F (0) +
∞∑n=1
F (2 l n) =1
2A0 +
∞∑n=1
An
con An =1
l
∫ ∞0
cos(nπ z
l
)F (z) dz.
(8.1)
Para deducir (1.3) tomamos l = 1/2 y f la exten-
sión par de F , f(x) = F (|x|). El primer miembro
de (8.1) es
1
2f(0) +
∞∑n=1
f(n) =1
2
∞∑n=−∞
f(n). (8.2)
Mientras que el segundo miembro es
∫ ∞0
f(x) dx+ 2
∞∑n=1
∫ ∞0
cos(2π nx)f(x) dx
=
∞∑n=−∞
∫ ∞0
cos(2π nx)f(x) dx
=
∞∑n=−∞
1
2f(n) . (8.3)
Con esto queda probado (1.3) a partir de la fór-
mula original (8.1) para funciones pares. ¿Y para
las impares? Muy fácil: 0 = 0. ¿Y para el resto?
Todas las funciones son suma de una par y otra
impar. Poisson no da referencias a trabajos ante-
riores, por lo cual es natural suponer que es en este
trabajo donde obtuvo su fórmula originalmente.
¿Cómo deduce (8.1) Poisson? Esencialmente
como dedujimos (1.3) de (1.4) pero de una ma-
nera más enrevesada, considera la serie de cosenos
en lugar de la de Fourier, teniendo cuidado con las
discontinuidades que induce en los extremos. Ya
es tarde para que mande una fe de erratas pero
de hecho hay un pequeño descuido en su expresión
de partida justamente por la contribución de esos
términos.
Para terminar, un reto para el lector. Unas pá-
ginas antes de llegar a su fórmula, Poisson da sen-
tido a algunas series divergentes curiosas:
¿Qué estaba usando para llegar a tan estrafalarias
identidades? Pista: empieza por “núcleo de”.
Agradecimientos. Este artículo tiene como baseun coloquio impartido el 14 de diciembre de 2016 en laUniversitat Autònoma de Barcelona por el primer au-tor quien aprovecha la oportunidad para agradecer aManuel Castellet, Xavier Xarles y al resto de los miem-bros del Departament de Matemàtiques de la Univer-sitat Autònoma de Barcelona atenderle tan inmereci-damente bien durante su visita. Agradecemos tambiéna Carlos Pastor sus finas observaciones relativas a unaversión preliminar de este trabajo.
Fernando Chamizo está parcialmente financiado
por el proyecto MTM2014-56350-P y por el Programa
Centro de Excelencia Severo Ochoa SEV-2015-0554 del
MINECO.
Referencias
[Art64] E. Artin. The Gamma Function. Transla-
ted by Michael Butler. Athena Series: Se-
lected Topics in Mathematics. Holt, Ri-
nehart and Winston, New York-Toronto-
London, 1964.
[CdV72] Y. Colin de Verdière. Spectre du lapla-
cien et longueurs des géodésiques pério-
diques. C. R. Acad. Sci. Paris Sér. A-B,
275:A805–A808, 1972.
[CE03] H. Cohn and N. Elkies. New upper
bounds on sphere packings. I. Ann. of
Math. (2), 157(2):689–714, 2003.
[Cha74] J. Chazarain. Formule de Poisson
pour les variétés riemanniennes. Invent.
Math., 24:65–82, 1974.
[Cór88] A. Córdoba. La formule sommatoire de
Poisson. C. R. Acad. Sci. Paris Sér. I
Math., 306(8):373–376, 1988.
[Cór89] A. Córdoba. Dirac combs. Lett. Math.
Phys., 17(3):191–196, 1989.
[Cor00] F. Corbalán. Galois. Revolución y mate-
máticas. Nivola libros y ediciones, S.L.,
2000.
[CR15a] F. Chamizo and D. Raboso. On the
Kuznetsov formula. J. Funct. Anal.,
268(4):869–886, 2015.
[CR15b] F. Chamizo and D. Raboso. Van der Cor-
put method and optical illusions. Indag.
Math. (N.S.), 26(5):723–735, 2015.
[DG75] J. J. Duistermaat and V. W. Guillemin.
The spectrum of positive elliptic opera-
tors and periodic bicharacteristics. In-
vent. Math., 29(1):39–79, 1975.
[Gas17] A. Gasull. Sumes harmóniques. Mat.
Mat., 12017(1), 2017.
[GK91] S. W. Graham and G. Kolesnik. van der
Corput’s method of exponential sums, vo-
lume 126 of London Mathematical Society
Lecture Note Series. Cambridge Univer-
sity Press, Cambridge, 1991.
[Gro85] E. Grosswald. Representations of inte-
gers as sums of squares. Springer-Verlag,
New York, 1985.
[Gui41] A. P. Guinand. On Poisson’s summation
formula. Ann. of Math. (2), 42:591–603,
1941.
[Har20] G. H. Hardy. On the representation of
a number as the sum of any number of
squares, and in particular of five. Trans.
Amer. Math. Soc., 21(3):255–284, 1920.
[IK04] H. Iwaniec and E. Kowalski. Analytic
number theory, volume 53 of American
Mathematical Society Colloquium Publi-
cations. American Mathematical Society,
Providence, RI, 2004.
[Ing90] A. E. Ingham. The distribution of prime
numbers. Cambridge Mathematical Li-
brary. Cambridge University Press, Cam-
bridge, 1990. Reprint of the 1932 original,
With a foreword by R. C. Vaughan.
[Iwa97] H. Iwaniec. Topics in classical automorp-
hic forms, volume 17 of Graduate Studies
in Mathematics. American Mathematical
Society, Providence, RI, 1997.
[Kol16] M. N. Kolountzakis. Fourier pairs of dis-
crete support with little structure. J.
Fourier Anal. Appl., 22(1):1–5, 2016.
[Kuz80] N. V. Kuznecov. The Petersson conjectu-
re for cusp forms of weight zero and the
Linnik conjecture. Sums of Kloosterman
sums. Mat. Sb. (N.S.), 111(153)(3):334–
383, 479, 1980.
[LO15] N. Lev and A. Olevskii. Quasicrystals
and Poisson’s summation formula. In-
vent. Math., 200(2):585–606, 2015.
[LO16] N. Lev and A. Olevskii. Quasicrystals
with discrete support and spectrum. Rev.
Mat. Iberoam., 32(4):1341–1352, 2016.
[LWL24] J. E. Littlewood, A. Walfisz, and E. Lan-
dau. The lattice points of a circle. Proc.
Roy. Soc. London. Ser. A., 106(739):478–
488, 1924.
[Max65] J. C. Maxwell. A dynamical theory of
the electromagnetic field. Philosophical
Transactions of the Royal Society of Lon-
don, 155:459–513, 1865.
[Mey16] Y. F. Meyer. Measures with locally finite
support and spectrum. Proc. Natl. Acad.
Sci. USA, 113(12):3152–3158, 2016.
[PZ04] F. Pfender and G. M. Ziegler. Kis-
sing numbers, sphere packings, and some
unexpected proofs. Notices Amer. Math.
Soc., 51(8):873–883, 2004.
[Rab15] D. Raboso. When the modular world
becomes non-holomorphic. In Trends in
number theory, volume 649 of Contemp.
Math., pages 221–244. Amer. Math. Soc.,
Providence, RI, 2015.
[Ran77] R. A. Rankin. Modular forms and fun-
ctions. Cambridge University Press,
Cambridge-New York-Melbourne, 1977.
[Sar82] P. Sarnak. Class numbers of indefinite bi-
nary quadratic forms. J. Number Theory,
15(2):229–247, 1982.
[Sel56] A. Selberg. Harmonic analysis and dis-
continuous groups in weakly symmetric
Riemannian spaces with applications to
Dirichlet series. J. Indian Math. Soc.
(N.S.), 20:47–87, 1956.
[Sko02] N.-P. Skoruppa. Quick asymptotic upper
bounds for lattice kissing numbers. Mat-
hematika, 49(1-2):51–57 (2004), 2002.
[Spi65] M. Spivak. Calculus on manifolds. A mo-
dern approach to classical theorems of ad-
vanced calculus. W. A. Benjamin, Inc.,
New York-Amsterdam, 1965.
[Ven70] B. A. Venkov. Elementary number
theory. Translated from the Russian
and edited by Helen Alderson. Wolters-
Noordhoff Publishing, Groningen, 1970.
[Wei79] A. Weil. Scientific works. Collected pa-
pers. Vol. III (1964–1978). Springer-
Verlag, New York-Heidelberg, 1979.
Fernando Chamizo Lorente
Departamento de Matemáticas e IC-
MAT
Universidad Autónoma de Madrid
Dulcinea Raboso Paniagua
Doctora por la
Universidad Autónoma de Madrid
Publicat el 9 de juny de 2017
Notes
1Si esto te parecen niñerías y te gustan las emo-
ciones fuertes, abre un libro de teoría cuántica de
campos, allí verás pocas cantidades finitas distintas
de los números de página (vale, es una exageración,
pero si eres matemático nunca habrás visto tantas
integrales divergentes juntas: la energía del vacío
es mínima y ya vale infinito). ´
2Según se recoge en [Spi65, 3-41], Lord Kelvin
afirmó que esta evaluación de la integral es tan evi-
dente para un matemático como que dos y dos son
cuatro. No te preocupes, aparte de despoblar las
facultades de matemáticas lo habría hecho con las
escuelas de ingenieros aeronáuticos por su descon-
fianza en la navegación aérea. ´
3Mientras que para el primero el límite era 12,
Gregory pensaba que se podría añadir una esfera
más. ´
4Son κ(1) = 2, κ(2) = 6, κ(3) = 12, κ(4) = 24,
κ(8) = 240, κ(24) = 196560, como dijo J.H. Con-
way del último, es que hay un montón de espa-
cio allí arriba. Fue Gauss el primero que probó
`(3) = 12, dando un retículo sencillo que sustituye
al icosaedro antes mencionado. ´
5Las funciones modulares de peso k respon-
den a la misma definición pero exigiendo f(z) =
z−k f(−1/z). ´
6En palabras de A. Weil [Wei79, p.463], fue ne-
cesario Maass para sacarnos del gueto de las fun-
ciones holomorfas. ´
7Las formas modulares holomorfas de alguna
manera se corresponden con λ = 0 ya que son
funciones armónicas por las ecuaciones de Cauchy-
Riemann. ´
8(con un máster en matemáticas) ´
9Las complicaciones comienzan cuando uno po-
ne condiciones de contorno en intervalos acotados
o cuando se consideran dimensiones mayores donde
en vez de dos direcciones, derecha e izquierda, hay
infinitas. ´
10En la literatura especializada no se le llama
ni primo ni monstruo porque no sería serio, sino
solución fundamental de la ecuación de ondas. ´
11Si quieres algo intermedio pero mucho más vio-
lento busca en la red simulaciones de colisiones de
un asteroide con la Tierra. Con un poco de imagi-
nación se dan un aire a las figuras de más adelante,
¿verdad? ´
12Si solo hay ínfimo, el argumento también fun-
ciona (ejercicio fácil). ´
13Provisionalmente: “Se puede pues esperar a que
el autor haya publicado su trabajo completo para
formarse una opinión definitiva” [Cor00]. ´