EL PROBLEMA DE n CUERPOS

C. Sim6

¢1 .- Ecuaciones y definición del problema . Integrales primeras : Teo-

rema s de Bruns, Painlev é y Poincaré .

Desde el punto de vista matemático consideraremos el proble-

ma de n cuerpos como el de n masas no nulas, puntuales, sometidas
a la ley de la atracción universal de Newton . Si xi es la posición
del i-ésimo cuerpo respecto al c .d .m . y mi su masa, se tiene

x = Z mi i=1=n (1)
i i¢3 =1 Ix j -x i 13

como sistema de ecuaciones diferenciales que rige el movimiento de

los n cuerpos . (a)
Desde el punto de vista físico éste es un modelo idealizado .
Por un lado los cuerpos celestes no son puntuales . Por otro las
fuerzas en juego pueden ser de otro tipo distinto al l/r2 newto
niano (sin embargo muchos de los resultados de la teoría siguen
siendo válidos para otras fuerzas, por ejemplo del tipo l/ra
q e 11,31 , o exigiendo s61o comportamiento l/r2 para distancias
grandes) . Finalmente (1) implica una hipótesis conservativa . Las
pérdidas de energía por efecto de mareas, por ejemplo, no se tie
nen en cuenta . Volveremos a este punto más adelante (en¢5) .
Otro aspecto a tener en cuenta es que condiciones iniciales
despreciables matemáticamente, por ejemplo un conjunto de medida
nula en el espacio físico, pueden ser fundamentales físicamente.
Supongamos que el inicio del cosmos fuera debido a una explosión "
Se tendría entonces momento angular total nulo, lo que res-
tringe el movimiento a una variedad de medida de bebesgue nula .
A falta de una mejor comprensión de las leyes físicas que ri-
gen el movimiento de los astros tomaremos (1) como modelo . Este

(a) Se supondrá en lo que sigue que la constante de Newton vale 1 .

38
 es un sistema de 6 n ecuaciones de primer orden . Son
clásicas las in
tegrales primeras del centro de masas y de su velocidad, del
momen-
to angular C = E mi x i n x i y de la energia h = F, 1 mi z2
i _ 2 1
1
mi mj
_ E . En total se dispone de 10 integrales primeras .
1<i<j<n I x i -x j1

Además, debido a que la acción de S 1 en el espacio fásico transforma

toda solución de (1) en otra solución (eliminación . del nodo ) y a que
es posible eliminar el tiempo como variable independiente, se
tiene
una reducción a un sistema de 6 n -12 ecuaciones . As! el problema de
2 cuerpos queda resuelto y el de 3, reducido a 6ó orden . Tal
reducción
puede verse en Whittaker .
Es natural preguntarse por la existencia de más integrales prime
ras de (1) . Por el momento no se conoce ninguna otra (con carácter
global ; vease Losco para formas invariantes) y las experiencias num_é
ricas efectuadas parecen indicar que no existen más integrales pri-
meras (ver§5) .
Resultados clásicos en esta dirección son los célebres

Teorema (Bruns 1887) : Toda integral primera de (1) algebraica'en x, .x,

es función algebraica de las integrales clásicas .
Teorema (Painlevé 1896-1898) :Toda integral primera de (1) algebraica
en z es función algebraica de las clásicas .
Teorema (Poincaré 1889) : Sea F = F o + PF + 2 F2
1 ~l + . . . el hamiltoniano del
problema de n cuerpos desarrollado en potencias de un pequeño pará
metro u (por ejemplo, la masa de uno de los cuerpos) . Suponemos
F an_a
lítico en 4, F o = FO(p)(independiente de q), Fi , funciones de p,q,
i>0
2n- periódicas en q . Suponemos, además, que para cada cuerpo
p4gi son
variables canónicas asociadas a los elementos elípticos . Entonces el
sistema
qi =
= _bF/a g i , bF/api
Pi no tiene más integrales primeras

analíticas en P ,p,q y 2n -periódicas en q que las clásicas .

Este teorema no es excesivamente sorprendente ya
que se sabe(MoseB'

39
que para cierto conjunto denso de sistemas hamiltonianos no existen
1
integrales primeras de clase le distintas del hamiltoniano .

12 .- 2 cuerpos : Solución ; efemérides ; determinación de órbitas . Coli -

sione s : Regu larizaci ón LC y KS . Espacio de órbit as . Transferencia .

El problema más sencillo, único resuelto y el que sirve de base

a las técnicas de perturbaciones y a la mayoría de problemas de la
Astrodinámica es el de 2 cuerpos . Para estudiar el movimiento (Keple

riano) de uno de los cuerpos respecto al c .d .m común (u, homot .ética-

mente, respecto al otro cuerpo) se introducen 2 ángulos :I(inclinaci6n)
y Q(argumento del nodo) para fijar la posición del plano de la órbi
ta respecto a un sistema de referencia dado . La órbita en dicho pla-
no es una cónica de semieje a y ex
centricidad e, con un foco en el c .
d .m . . Se define el ángulo w (argu
mento del periastro) como el que
forman la dirección del nodo y la
del punto de la órbita más cercano
al c .d .ii1 . (periastro) .
Los elementos a,e,I,n , w,rr o , se
llaman Keplerianos . El sexto u 0 0

anomalía media en la época está de

finido módulo 2n y varia linealme-te con t . Tos demás son constantes .
La posición en la órbita se determina mediante el ángulo V, anomalía

2= 1+ ée
verdadera, formada por el periastr 1 sici6n actual del astro .
Este se relaciona mediante tg tg E/2 con la llamada

anomalia excéntrica E , ángulo auxiliar, que a su vez se relaciona

con el tiempo mediante M =M . + ta -1/3 = E + e sin E (mod 2n) .
Esta ecuación se conoce como ecuación de Kepler . Su resolución
impulsó en parte la contribución de Cauchy a la teoría de la variable
compleja . Una expresión de E-M como serie en sin KIT con coeficientes

40
funciones de e fue dada por Lagrange en 1769 . Los coeficientes re
sultan ser funciones de Bessel , cuyo estudio sistemático se abor
darla por primera vez 55 años más tarde .
El proceso de determinar la posición de un astro del que se co
nocen los elementos se llama cálculo de efemérides . Sin embargo,
para el cálculo efectivo de la posición deberemos tener en cuenta
la presencia de perturbaciones que alteran el problema de 2 cuer-
pos .
El problema inverso : conocida la posición del astro calcular
sus elementos, se conoce como determinación de órbitas . Debido a
la dificultad de conocer x y x en un instante dado, dicho cál
culo se convierte en un problema de valores en la frontera .
Es típica la determinación de la órbita de un asteroide o un
cometa mediante observaciones angulares, esto es : para cada instan
te de tiempo se conoce en qué dirección está respecto al observa
dor, pero no se saben ni su distancia ni su velocidad .
Los problemas motivados por la disponibilidad de un elevado
número de observaciones no compatibles debido a los errores de ob
servaci6n llevaron a Gauss a desarrollar en 1801 el método de mí-
nimos cuadrados .
Otro ejemplo típico de determinación de órbitas en el problema
de 2 cuerpos es el del rendez-vous o interceptación : los datos de
que se parte son t l , x 1 , t2 y x2 .
Notable interés presenta el estudio de las órbitas de colisión (es-

t-11-s, en que los dos cuerpos chocan) . En el espacio fásico (6 dimensional

llenan una variedad 4 dimensional, por lo que su medida es nula . Desta-
quemos que (1) deja de tener sentido si x 1 _x 2 : la trayectoria ha
llegado a una singularidad .
Es natural preguntarse si es posible encontrar un cambio de va
riables x { t que permita regularizar (1), esto es, transforma-
YO}
la en una ecuación sin singularidades, y en caso afirmativo inter-

41
pretar físicamente el resultado .
2
Para movimiento en R la respuesta fue dada en 1906 por Levi-
Civita . La transformación

1
Z 1 - Z21 Zl l 'x1
L (Z) Z =I Z2 Z 11 I = Ix 2' , ds = Ix1 dt
Z2I

que introduce un nueva variable espacial Z ( la raíz cuadrada de

x considerados x,Z e T ) y una nueva temporal, s, transforman (1)
en Z" + 2 Z = 0 (2), donde 1 = ds y h es energía total . Se obser
va de paso que (2) es la ecuación de un movimiento armónico . I
3
Se buscó en vano una matriz en R que generalizara la matriz
de Levi-Civita L(Z), o sea, cuya primera columna fuera Z y las
demás fueran ortogonales entre sí y a Z . Esta imposibilidad es de
bida al teorema de Poincaré de no existencia de campos vectoria-
les continuos en S 2 _
Sin embargo si existen dichos campos (que además engendran el
e .tangente) en S 3 . (Un célebre resultado establece que S l ,S 3 ,S7
son las únicas esferas con esa propiedad) .
En R4 se sustituye L(Z) por

Z 1 - Z 2 - Z3 Z4
1 _ Z4_ R4 -~ R4
Z2 Z Z3
L (Z) = La aplicación Z i L (Z)
3 4 1 2 Z
Z 2 Z Z
Z4_ Z3 Z2_ Z1

fue dada por Hopf en 1931 y resucitada por Kustaanheimo y Stie-

fel en 1965 que la aplicaron a la regularización . Se sumerge R 3
en R 4 haciendo x 4 = 0 . El cambio L (Z) Z =x, ds = Ixi dt convierte
(1) en (2) . La fibra de cada punto de R3 por la aplicación KS es
una circunferencia de centro el origen y radio YIXI Tanto en el ca
so LC como en el KS la regularización equivale físicamente a una co
lisi6n elástica : El movimiento Kepleriano se prolonga después de la
colisión conservando energía y cambiandó de signo la cantidad de mol
vimiento . -I
Otra pregunta que podemos formularnos acerca del problema de
dos cuerpos es la de qué estructura topol6gica puede darse al conjun-j

42
 P

to de las órbitas 9 . Este puede obtenerse del espacio fásico identi

ficando todos los puntos de una misma órbita o bien, equivalentemen
te, definiendo una distancia entre órbitas, ya sea a partir de las
integrales primeras o geométricamente, extendiendo adecuadamente la
métrica de Hausdorff a conjuntos no acotados .Para que se tenga un
espacio separado deben identificarse la órbita de llegada a una co-
lisión con la correspondiente de salida (Souriau) . Moser probó (1970)
que para h < 0 fijado, el conjunto de órbitas es homeomorfo a S 2 XS 2
(y que la subv . del espacio de las fases correspondiente es T 3
.S ) .
Generalizando este resultado se tiene
Teorema (Simó 1974) El conjunto de las órbitas del problema de 2
cuerpos es homeomorfo a S2 x S2 x R- fi a, -a,h)l a e S2, h? 0} . Este
resultado implica la no existencia de "elementos" sin singularida-
des en el problema de 2 cuerpos, pues una carta es insuficiente .
Otra consecuencia debida al teorema de Poincaré de no existencia de

campos continuos en S 2 es que no existe ninguna aplicación continua

9->-R 3 que asigne a cada órbita un punto de la misma . Esto es, ningún
elemento del tipo Mp puede ser asignado de .forma continua a una ór-
bita .
Finalmente un problema que no se halla en la naturaleza pero
aparece en astrodinámica es el de transferencia optimal de órbitas .
Se trata de pasar de una órbita Kepleriana a otra de manera que el
consumo de energía en dicho paso sea mínimo . El problema está resue_1
to en el caso de que ambas órbitas sean coplanarias, pero no en el
caso general . Además existen casos en los que la transferencia 6pti
ma desde el punto de vista energético requiere tiempo infinito (las
llamadas transferencias por el infinito, a través de 2 órbitas in-
termedias) .

3 .- Tres cuerpos : Pequeño T . de Sundman Teorema de Sundman Solu-

ciones conocidas : soluciones periódicas Posiciones de equilibrio
relativo . Interés en el cálc ulo de variedades invariantes .
 El paso de 2 a- 3 cuerpos representa la aparición de multitud
de diferencias . Podemos decir que en el problema general(masas
cualesquiera) el conocimiento de soluciones del problema de 2
cuerpos no representa ninguna aportación .
Un primer teorema que se aplica a un n4 cualquiera de cuer
pos es llamado t .del colapso total o pequeño t.de Sundman : Condi
ci6n necesaria para que los n cuerpos colisionen es que el momen
to angular total sea nulo . Por tanto si c9£0 y n = 3 s61o puede
haber colisiones dobles . Sundman probó que la presencia del ter
cer cuerpo no es obstáculo para que el movimiento sea regulariza
ble, y de hecho, si no se anula c, di6 la solución al problema
de 3 cuerpos .

Teorema (Sundman(1913) Si c y¿0, las coordenadas cartesianas y

el tiempo son funciones analíticas de w que convergen si Iwl < 1,
donde w = th s= ~ t (U+ 1) dt siendo U la energía
b4 to

potencial y b un número positivo función de las masas y de las

condiciones iniciales para t =t ., En particular el desarrollo es
válido para toda t real .
A pesar de que este teorema resuelve el problema de 3 cuerpos
excepto en un caso de medida nula (c =0), no es práctico dada la
lenta convergencia de las series y el hecho de que no suministra
información del comportamiento cualitativo .
La obtención de soluciones para intervalos de tiempo relati-
vamente pequeños no ofrece dificultad (sea desde el punto de vista
analítico o numérico .Véa3e sin embargofi0) .Sin embargo si queremos pre
decir algo para tiempo no acotado tendremos dificultades .
Estas desaparecen si la órbita es periódica . En este caso el
conocimiento de la misma para un tiempo finito permite conocerla
para todo tiempo .
Póincaré conjeturó que dada una órbita x(t) de condiciones

44
iniciales X o , Xo para t = 0 VE > 0 , d T > 0 existen condicio-
nes X*, Xó tales que IXO- Xó~ < 6 IX o - X* I <g dan una ór
bita periódica x * (t) de manera que dt e [O,T] Ix(t)-x*(t)~ < t .
Esta conjetura, a la que había que añadir h < 0 es la que ha impul
sado la búsqueda de órbitas periódicas . Actualmente no se conoce
todavía la certeza o falsedad de dicha conjetura .
En el caso de que 2 de las masas fueran pequeñas (caso pla-
netario) Poincaré dividió las órbitas en 3 especies :
1) I=0, er0
2) I =0, e grande
3) I grande, e-0

La búsqueda efectiva de soluciones periódicas del problema de 3

cuerpos, se aborda a partir del problema restringido (ver§4), pe
ro una exploración completa dista de llevarse a cabo .
Un tipo particular de soluciones periódicas es aquél en que
se pasa de Xi(0) a Xi(t)
mediante una homotecia H . función de t,
t
idéntica para i =l`. n,(en particular H t puede ser un giro), Euler
(1765) obtuvo los valores de Xi(0) en el caso en que Xl ,X 2 ,X3
estuvieran alineadas, obteniendo 3 soluciones y Lagrange (1772)
determinó las otras 2 soluciones restantes (simétricas) en las
que los 3 cuerpos forman triángulo equilátero . (Naturalmente hay
soluciones de tipo homotético no periódicas, llevando a colapso
o a desintegración) .
El conjunto de condiciones iniciales que da órbitas homoté-
ticas periódicas es cerrado por homotecias . cada clase se llama
una posición de equilibrio relativo (p .e .r .) . Se demuestra que
en una p .e .r . los n cuerpos deben estar en un plano . El interés
del conocimiento de los p .e .r . en el caso general se ha visto au-
mentado por el hecho de que los p .e .r . coinciden con los puntos
mi m*
criticos del potencial -E restringido a la variedad
l$i<j<n I x i. - x .

45
l2 E mi l xil 2 = 1 .
i-1
Consideremos ahora la aplicación (E,J) que a cada punto del
espacio fásico del que se excluyen las colisiones le hace corres
ponder energia y momento angular . Dados h E R, c e R3 el conjun-
-1
to Ihc =(E,J) (h,c) es'invariante por el flujo .. El conocimien-
to de Ih c puede dar información sobre las órbitas con esos va-
lores de energía y momento .
Smale introdujo el concepto de conjunto de bifurcación : va-
lores de h,c donde Ih. cambia de estructura . Gracias a resulta-
dos de Easton y Cabral, se sabe que el conjunto de bifurcación
contiene a los puntos críticos .
En el caso plano del problema de n cuerpos los puntoscriticos
son :

U {(-V (c z).c) E R2 1 cyÉ 0, zEa .

2 E mi 1 Z i~ 2 -1~
a E P " e .r .

y el conjunto bifurcación es la unión de los puntos críticos

con los ejes h,c .
En el caso espacial los puntos críticos de (E,J) se obtie-
nen por unión del eje h con los puntos (h,c) tales que
(fi, Ic 1) es punto critico del caso plano .
El estudio del conjunto de bifurcación en el caso espacial
debe llevarse a cabo directamente, pero en todo caso requiere
conocer el de puntos críticos, y este el de p.e .r .
Esto ha sido llevado a cabo por C . Simó y J . Llibre . Se pue

2
de calcular Ih c en términos de antiimágenes de lc donde

1 (z) = -V(Z)2/ siendo E mi IZil 2 = l, V el poten-

C 2y c (z)

cial y Yc (z) la energía cinética de los cuerpos considerados

como sólido rígido con momento angular c . Sin entrar en detalles
digamos que en el caso de masas distintas se tienen 11 estruc-
turas de Ihc para el caso plano y 17 en el caso espacial .

46
 Destaquemos que mientras en el caso plano se tienen difeomorfis
mos,en el caso espacial se reducen a homeomorfismos . De hecho
Ihc
no son variedades topológicas en el caso espacial : La dimensión
lo-
cal varía .

¢4 .- El problema restr inqido : Ecuaciones . Curvas de velocidad nula .

Soluciones periódicas Aplicación de Poincaré E stabilidad . Puntos
fijos .
Un caso sencillo del problema de 3 cuerpos es aquel en que una
de las masas es infinitesimal de manera que su influencia sobre las
otras dos es despreciable . El movimiento de éstas es, pues, el de
un problema de dos cuerpos .Supondremos que éste es circular y que
el 3ó se mueve en el plano de los 2 primeros . Hablamos de movimien-
to restringido plano y circular . Si usamos ejes (sinódicos ) girato
rios tales que el eje X contenga siempre los 2 primarios se obtie-
ne

,y+ 2x =n y n =
2 L
(1-p)r2+ ur2] +

+ 1- ~L + AL-
r1 r2
que son las ecuaciones que rigen el movimien-
to del cuerpo infinitesimal . Tienen una inte-
gral primera (Jacobi 1836) : 0 -
2
Pijadas unas condiciones iniciales conocemos
(x 2+ y 2 ) = C .

C . Haciendo ik = Y = 0 se tienen las llamadas

curvas de velocidad nula que delimitan-regiones (en coordenadas si-
n6dicas), que no puede ocupar el Ser cuerpo . En la figura se muestra
la forma de esas curvas para valores
distintos de C . El movimiento sólo es
posible en el interior de las mismas
Como para la luna se tiene un caso aná
logo a C l , ésta no puede escapar de la
proximidad de la tierra para caer hacia
el sol ( en el modelo de 3 cuerpos restringido)
 Los puntos x,y en los que n es estacionario corresponden a po
siciones de equilibrio en el sistema sinódico y de equilibrio re-
lativo en el sidéreo . Se encuentran las ya conocidas posiciones de
Euler y Lagrangé . El problema restringido de 3 cuerpos es suficien
temente sencillo como para que s=haya abordado de manera exhaustiva
la búsqueda de órbitas periódicas . El primer estudio sistemático se
inició en 1900 por,Strómgren y sus colaboradores con el llamado pro
blema de Copenhage Iv caso en quely = 2 . Dada una órbita periódica
podemos buscar condiciones iniciales próximas a las de dicha órbita
que nos proporcionen otra órbita periódica . Por ejemplo, si bus-
camos órbitas periódicas simétricas, las c .i . pueden ser x o ,y o =0 1
xo =0 1 YO , o equivalentemente x o ,c o . Si la órbita periódica puede
prolongarse(curva en el espacio, x,c, por ejemplo), hablamos de
familia de órbitas periódicas . (Para sistemas hamiltonianos si
es periódica y 6 es una aplicación de Poincaré sobre una variedad
2n - 2' dimensional E transversal a y contenida en la variedad de
energía constante a la que pertenece y, entonces esta órbita se ha
lla en una familia 1 -paramétrica de órbitas periódicas si 1 no
esde
valor. propio ~ .~~i °- matriz de ,Tiví'avurviTiia,
D8 (m) d Siendo üi

(Véase Robinson) .
En el caso de Copenhagen se hallaron 13 familias de órbitas
periódicas, algunas de ellas con varias subclases . Posteriormente
(Hénon, Bartlett) se han halladol0 familias más de órbitas perió-
dicas simétricas, se han encontrado familias no simétricas y se
han estudiado los casos ~ X1/2 .
Una herramienta básica para la obtención de órbitas periódicas
es la aplicación de Poincaré . Sobre H = cte ., cortemos el flujo trars
versalmente por una vaiedad E.Se define la aplicación de Poincaré9
(quizá sobre un subconjunto de E) como la aplicación que a un punto
m de E le hace corresponder la próxima intersección A(m) de E con
el flujo pasando por m . Es claro que g(m) = m es la condición de ór
bita periódica (simple), o, en general 8 (m) =m .

48
 La aplicación de Poincaré es simpléctica (conserva la forma
dp A dq de E) y en particular preserva medida . Los teoremas de pun
to fijo pueden servir para garantizar existencia de órbitas periódi
cas (pero téngase en cuenta que por preservar volumen no pueden
aplicarse teoremas de contracción) .

Una primera indicación de la estabilidad (orbital) de la ór-

bita periódica la suministra el análisis del espectro de De .Si
existen valores propios reales el punto es hiperbólico en lo que
hace referencia a la parte lineal (y también a la nolineal) . En ca
so contrario es elíptico (linealmente) y la estabilidad frente a
perturbaciones no lineales puede obtenerse en ciertos casos de la
Teoría KAM (ver §5) .

¢5 .- Totalidad de soluciones en el problema restringido : Integrabi -

lidad . Teorema de Liouville-Arnol'd . Teorema KAM Movimientos gua-
siperi6dicos . Aplicaciones : Teorema de Leontovich Sistema solar .
Difusión de Arnol'd .
Hénon ha continuado la exploración numérica en el caso del pro
blema de Copenhagen mediante aplicación de Poincaré respecto al plá
no x,x, tratando todo tipo de órbitas (y no sólo las periódicas).Ios
resultados son muy instructivos .
Para valores C grandes (> 4 .5) las suces mas imágenes de x,x
por 8 parecen llenar una curva cerrada de un modo denso correspon-
diendo a trayectorias quasi periódicas .
Sin embargo para número de rotación racional (nó de rotación= frac-
ción de vuelta girada alrededor
de un punto fijo en cada aplica
ci6n de 6, en promedio), apare-
cen islas que envuelven puntos
elípticos y que quedan -separadas
por puntos hiperbólicos . A medi

49
da que nos acercamos al punto fijo de 9 (correspondiente a la
órbita periódica) la anchura de las islas ( y de las zonas hi- .
perb6licas) disminuye siendo no apreciable .

Rara valores menores de C aparecen regiones en las que las sucesivas

imágenes de (x,x)por 6 parecen aleatc ..as y, en cierto modo lle
nan todo el posible espacio de variación de (x,~c) formando lo
que se llama'Un mar erg6dico en el que permanecen islas integra-
bles" .
Si un sistema hamiltoniano posee n integrales primeras in-
dependientes F i en involuci6n, el teorema de Liouville-Arnol'd
asegura que si Fi= cte . i =l- n define una variedad compacta,
ésta es un toro n-dimensional (en .el caso no compacto será (Si )k x

xR n-k en el que el flujo viene dado por qi = wi = cte .

Al cortar dicho toro por E aparece esta variedad foliada por
curvas invariantes . Esta es la característica de sistemas inte-
grables . Por este motivo la presunta aparición de esas curvas ha
motivado la búsqueda de una 21 integral (la lá es la de Jacobi) .
En problemas más generales, esto es, con mayor nó de grados de
libertad se habla de la búsqueda de la 3a integral (siendo las
2 primeras las de energía y momento) . Es claro que esta búsglle
da es formal ya que la aparición de puntos elípticos e hiperb6
licos de 6k entre curvas invariantes de ® es general . En .este
aspecto una aproximación mediante la forma normal de Birkhoff
es útil .
En el extremo opuesto de los sistemas integrables están los
erg6dicos .La aplicación 6 es erg6dica si para toda f eL l (E), la
media de f sobre E (supuesto E compacta) es igual a la media "tem
poral" sobre la trayectoria por 6 . Esto es : b xe E

1 N k 1 f C .
lim Ef (6 x ) =
N , c, Ñ
k =0 l. (E) E

50
 Se tiene de inmediato que si A es erg6dica todo subconjun
to propio de E invariante por 8 tiene medida nula .

En la práctica parece que todo sistema hamiltoniano tiene

parte de "componente integrable" y parte de "componente erg66
dica" . Sin embargo no está claro como abordar el problema .
Si un sistema hamiltoniano H está cerca de un sistema in
tegrable H o : H = H o +e H1 , podemos preguntarnos si retiene al-
go de las características de integrabilidad .
La respuesta la da el teorema de Kolmogoroff-Arnol'd-Moser

(1954-1962-1963) probado por Arnol'd por el caso analítico,

le 333
Moser para el caso y luego por Rüsmann para le s (se

sabe que para 1?4 es falso), que asegura (tanto para sis-
temas hamiltonianos como para aplicaciones globalmente ca-
n6nicas) que si las perturbaciones son suficientemente peque
ñas y se cumple una condición de no degeneración de H
o (rela
tiva a su jacobiano y hessiano ), para la mayor_ parte de los
toros invariantes estos sufren una deformación pero no desapá
recen . Esta " mayor parte" está formada por aquellos toros pá
ra los que Hke Zn , dk o e Z I(k,w) -k o i > Cw /Ik1 v donde

Cw es una constante que depende de w , V > n +2 es un entero

adecuado y ¡K¡ = E lK i l . Observemos que en dimensión 1 hay
que eliminar como posibles valores de w los racionales y los
nºs de Lioville .

Sobre los toros invariantes el flujo es el de un movimiento

quasi periódico : qi = w i cuyas órbitas son densas . La prime-

ra referenlia a las funciones quasi-peri6dicas puede encontrar

se en la teoría clásica del movimiento de los astros a base de

51
eoiciclos desarrollada ya por los ?strónomos babil6nicos .
Una aplicación del teorema KAM fue dada en 1962 por Leonto
vich : En el problema plano circular restringido de 3 cuerpos,
si los puntos triangulares de equilibrio relativo son estables
por la parte lineal (lo que ocurre si ~ < 2 (1-T23/27)) lo

son también (con la posible exclusión un nó finito de valores

de las masas) al considerar la parte nolineal .
Arnol'd aplicó su teorema al sistema solar con el siguien
te resultado : si las masas, excentricidades e inclinacio- :
nes de los planetas son suficientemente pequeños, para la ma
yor parte de condiciones iniciales, en el sentido de medida de
Lebesgue, el sistema solar es estable .
Destaquemos sin embargo que M . Hénon ha obtenido cotas de
lo que significa en este caso "suficientemente pequeño", ha-
llando que dichos valores deberían ser 0(10- 50 ) como=mucho .

Sin embargo podemos citar aquí la curiosa experiencia de

P . Nacozy : Si se integra numéricamente el movimiento alrede
dor del sol de Júpiter y Saturno con sus masas actuales (no
despreciables) durante un periodo largo de tiempo, no hay
evolución apreciable .
Si se incrementan las masas y se repite la experiencia se
halla idéntico resultado hasta que las masas son 29 veces
las actuales . Entonces Saturno es expulsado rápidamente (en
unos pocos centenares de años) . Parece pues que las cotas
para las que es válida la teoría de KAM pueden ser notable
mente aumentadas .

52
 Finalmente destaquemos que H = cte. da una variedad 2n -1 di-

mensional y los toros invariantes de un sistema próximo a integra

ble son n-dimensionales . Si n = 2 dichos toros ."separan" la varie
dad de energía constante . Si un movimiento se inicia en un toro
invariante queda siempre sobre él y si entre 2 toros no puede
escapar del recinto que delimitan .
Esto no es así si n ? 3 en cuyo caso los toros "no separan"
y el sistema puede evolucionar destruyéndose la estabilidad . Un
ejemplo fue dado por Arnol'd mediante las llamadas cadenas de
transición de toros bigotudos . Se tiene así la llamada inestabi-
lidad topol6gica o difusión de Arnol'd .
Según Moser podemos imaginar el espacio fásico como una es-
ponja estable . Los minúsculos agujeros corresponden a las zonas
situadas fuera de los toros invariantes y la inestabilidad corres
ponde a la existencia de finísimos tubos comunicando los agujeros .
Hasta la actualidad no se sabe si la difusl6n de Arnol'd se
presenta en la mayoría de los sistemas hamiltonianos o si ocurre
lo contrario . En todo caso es muy difícil de determinar numéri-
camente .

§6 .- Colisión triple : Importancia física . Teorema de Siegel . Re-

qularizaci6n .
Otra característica esencialmente distinta del problema de
3 cuerpos respecto al de 2 es la posibilidad de colisiones tri-
ples . Se sabe (Siegel 1941) que las órbitas que conducen a coli
sión triple forman una variedad 13-dimensional en el espacio de
fases 18-dimensional . Sin embargo se cree que el conjunto de ór-
bitas que pasan cerca de la colisión tt#Je tiene considerable importan-

cia en la evolución del problema de tres cuerpos . De hecho pare

ce comprobado que gran parte de las expulsiones de estrellas de
una galaxia se debe a aproximación de una estrella a una binaria

53
fuerte (es decir, formada por estrellas muy cercanas, con elevada
energía potencial) . Se está así en una situación próxima a colisión
triple .
Un resultado clásico asegura que cuando K cuerpos se acercan a
un punto en el que colisionan, sus posiciones relativas a dicho pun
to tienden a una configuración homotética . Así, en el caso K =3 los
cuerpos tienden a formar triángulo equilatero (independiente de sus
masas) o a estar alineados . Un caso parecido se tiene cuando los
cuerpos se alejan unos de otros en una expansión parabólica o hipar
b6lica .
Siegel probó que con la excepción de un conjunto discreto de
valores de las masas las coordenadas de los tres cuerpos que coli-
sionan (referidos al c .d .m .) se expresan analíticamente (para t > 0)
-~
en función de las variables t1/2 , t 2 , t-~ 3 en el caso equilateral
y t2/3 , t-X2 en el colineal . Exceptuando un conjunto discreto de
son
masas,X2 ,X y irracionales y por tanto la singularidad para t = 0
es esencial (a diferencia de la colisión binaria, en la que el pa-
t2/3) .
rámetro de uniformizaci6n es
As¡ pues, no parece que sea posible una regularización analíti-
ca de la colisión triple en el caso general . De hecho McGehee, Mar
chal y Waldvogel han dado ejemplos de casos en los que tal regula-
rización es imposible.
Una regularización más débil que la analítica es la geométrica .
Sea {Xi,}{i } una sucesión de condiciones iniciales que tiende hacia
una c .i .(x,x) que lleva a colisión triple . Si V t e R* fijo, X(t,xi)
depende continuamente de Xi . podemos definir la prolongación de
X(t) en el caso de colisión triple como lim X(t,xi ) . Sin embargo

hay ejemplos en los que V s > 0 la e-bola de (x,x) contiene puntos

cuyas trayectorias se alejan despues del paso por las proximidades
de colisión triple y uno de los cuerpos puede salir de esta situa-
ción con energía arbitrariamente grande . Sin embargo éste es otro
campo en el que hay mucho camino por recorrer .

54
 7 .- Evolución final : Clasificación de tipos de movimientos .
Chazy (1922) clasificó los tipos de movimientos del problema
de tres cuerpos cuando tico en 7 tipos :

1) H : Movimiento hiperbólico, en que las distancias mutuas

Ixi - x j 1 son de la forma 0(t) .

2) HPi : Movimiento hiperbólico-parabólico . El cuerpo i se

aleja hiperbólicamente de los otros 2, y éstos se hallan

en un movimiento parabólico (esto es la distancia es

O(t 2/3 ))-

2/3
3) P : Movimiento triparab6lico : (X _ X j 1 = 0(t ),
i
4) HE i : Movimiento hiperbólico-elíptico . Análogo a HP i pero
hallándose los 2 cuerpos restantes en movimiento elípti-
co .
5) PEi : Parabólico elíptico .
6) L : Movimientos estables a la Lagange : sup {IX X .I} l}<
.
i_
ij, t>0
7) OS : Movimientos oscilatorios : lim sup {/X .- X .I} = ~,
t~oo i,j 1 J

lim sup {ix .- X .I} < . Del tipo 7) Chazy no disponía

t-, oo i, j 1 J

de ningún ejemplo . El primero fue suministrado por Sit-

nikov(1959) . (ver § 8) .
Si h > 0 é61o son posibles 1),2),4) . Dado que i puede ser 1,
2,3, hay un total de 7 posibles evoluciones finales para t - + m .
Las mismas se tienen para t-+-oo . Puede pasarse de cualquiera de
ellas para t-4 -co a cualquier otra t-> +oo, según ha sido probado
por Birkhoff, Sitnikov y Alexeyev, y el conjunto de condiciones
que permite el paso de una a otra tiene medida positiva (excep-
to en el tipo HPi, cuya medida total es nula) . Es decir, puede
darse captura parcial (H -+ HEi+ 0~) y su reciproco (desintegra-
ci6n completa), así como intercambio (HE . -> HE . ) . Si h 9 0 la
1-00 7+ 00

situación es mucho más complicada . Dejando de lado el caso PE i

55
(de medida nula) consideramos el HE i , L y OS (cuya medida no se
conoce aún) .

Se sabe que es posible el paso de cualquiera de estas situa-

ciones para t-> -oo a otra t~ + oo , . pero que las transiciones
HE . HE . , HE HE . L y L ocurren para conjun-
1 .} J.- J + m, - °° + °°
tos de c .i . de medida positiva mientras que los pasos
HE,1- -> L + , HE . --> OS+ o0 tienen lugar para conjuntos de medi-
oo ao 1 -ao

da nula . Existe multitud de criterios para reconocer que se es-

tá en un caso HE, cuando se procede ala integración numérica de
1
las ecuaciones del movimiento .
HE i son abiertos conexos en el espacio fásico, pero cada uno
de ellos tiene ramificaciones complicadas en el espacio de fases .
Birkhoff decía que HE i eran como tres corrientes de agua provi-
niendo de -oo que se ramificaban en un conjunto numerable de arro
yos que-posteriormente se juntaban para (a excepción de una can-
tidad de agua de medida nula que puede quedar en L o OS) pasar
a formar las 3 corrientes HE .
1 + oo'

Szebehely ha conjeturado que si la relación entre 2 masas pe-

queñas es suficientemente cercana a 1 y si c2 h < 0 y suficiente-

mente cercano a cero, casi todas las soluciones son del tipo HE i
tanto para t-s + oo como para t -> - oo .

8 .- Movimiento Quasi-al eatorio : Problema de Sitnikov .

Consideremos dos masas m, =m 2 = 1/2 en movimiento elíptico en

el plano x,y y una tercera masa m, n: 0 en el eje Z, pasando por
el c .d .m . de n~ y m2 . Sea 2n el periodo de m, y m2 . La ecuación
-3/2
que rige
es el movimiento de m 3 z = -z ( Z 2 + r2 (t)) , donde

r(t) es la distancia de n~ al origen . Si la órbita es circular

el movimiento de m3 es periódico .

56
 Sean tk , kE Z, ceros de z . Suponemos tk < tk + 1, e introducimos

Sk = tk 2 l - tk . A cada solución z(t) podemos asociar la Buce

si6n (infinita por ambos lados) { Sk} .

Teorema .(Sitnikov) Dada una excentricidad e> p suficientemente pe-

queña existe un entero m = m(e) tal que cualquier sucesión (Sk} con
Sk k m corresponde a una solución de la ecuación diferencial ante-
rior .
El teorema puede extenderse a valores de m3 no nulos y de e
cualesquiera (con excepción de un conjunto discreto) .
Si, por ejemplo Sk-soo para k ,,, , Sk acotado para k ., -oo ten
dremos una transición de movimiento tipo Lagrange a oscilatorio .
Si Sk vale infinito para un valor de k finito, se tendrá un escape .
Aparentemente el movimiento, siendo determinista, tiene un compor-
tamiento aleatorio ya que la sucesión ( SO puede ser cualquiera
(con Sk > m) .

Básicamente la demostración del teorema usa el shift de Ber-

noulli : Sea A un alfabeto finito o numerable . Considera el espacio
X de las sucesiones doblemente infinitaslsk }, Sk EA . En dicho espa
cio se introducé una topología d:fmiendo como base de entornos de

S _ { Skti a . Uj={SI Sk = Sk , ¡k¡ < j } . Si P es una medida en

A¡ S E P(Sk) = 1, puede dotarse a X de la medida producto de la forma

keA
habitual . La aplicación a : X , X I (a (s)) k = S k- 1 se llama el shift
de Bernoulli . Trivialmente preserva medida, es ergódica y además es
un homeomorfismo . Definiendo la entropia como E - P(Sk) 1 n P(Sk)
SkCA
se tiene que dos shifts son isomorfos si y sólo si tienen la mis-
ma entropia .
En el caso de que A sea finito Q equivale a la llamada transfor

maci6n del .panadero .

La idea esencial es incluir el shift como subsistema de un homeo
morfismo ~ = Y , Y , esto es :
 á
3 T tal que X X conmuta, siendo
T1 1T
T un homeomorfismo entre X y su
y -r Y
imagen .

La teoría de los sistemas dinámicos quasi-aleatorios (=_ siste-

mas con entropía topol6gica positiva) ha sido desarrollada por
Alexeyev .
Se sabe (Markus, Meyer) que con carácter genérico (una pro-
piedad es genérica si es cierta para una intersección numera-
ble de abiertos densos), los sistemas hamiltonianos no son
ni integrables ni erg6dicos (en cada variedad H = cte) . Podemos
conjeturar sin embargo que un comportamiento del tipo del pro-
blema de Sitnikov (relacionados con la difusión de Arnol'd) si
es genérico .

9 .- Particularidades del problema de cuatro cuerpos : P e .r . y

evolución final .
El problema de cuatro cuerpos ofrece un nuevo salto en cuan
to a dificultad respecto al de 3 cuerpos . Sin embargo parece
que no haya ya variación cualitativa de las situaciones que pue
den darse al aumentar el número de cuerpos (aunque no hay resul-
tados en este sentido) .
En cuanto a p .e .r . se conocían clásicamente algunas solucio
nes . Si las masas son iguales las soluciones eran del tipo cua-
drado, triángulo equilátero con una masa en el .centro o colinea-
les . El autor , en colaboración con J . Llibre, ha encontrado
recientemente otra familia de soluciones en las que tres de las
masas se hallan en los vértices de un triángulo isósceles y la
cuarta en el interior . En total se conocen 50 p .e .r . y se tenTan
fundados motivos para creer que eranlasúnicas . Véanse sin embar-.
go resultados de Palmore . Smale conjeturó que la energía potencial
1
sobre E m i ( zi )2 =1 es funci6n de Morse . En este caso K dado
2
que dicho conjunto es homeoformo a Pn-Z ( T ) -D (D = colisiones) , la
conjetura implica que el ná de p .e .r . es finito para
todo n, y que para n = 4 la suma al-

58
ternada del n4 de p .e .r . que tienen indice por e impar es 2(valor
dela caracteristica de Euler-Poincaré de P2 (T) -D) . Esto se cumple
en nuestro caso ya que hay 12 p .e .r . alineadas y 8 de tipo equilá-
tero con indice 2, 24 de tipo isósceles con indice 3 y 6 de tipo
cuadrado con indice 4 .
A diferencia del caso n=3 en qie el número de las p . e . r . es inde-
pendiente de las masas, sabemos que si una de las masas es suficien
temente pequeña hay 38 p .e .r . si las otras 3 son suficientemente
parecidas y 34 en caso contrario . Para los valores en que ocurre
la transición (que forman una variedad de codim . 1 en el espacio de
las masas) el potencial restringido deja de ser f . de Morse . Un
análisis completo está en estudio .
Respecto a evolución final aparecen dos tipos nuevos de compor-
tamiento :
a) La expansión super-hiperbólica . S61o - se presenta si n ? 4 .
Si d(t) es el diámetro del sistema, se tiene d(t)/t ~ . . Como
mínimo hay 2 distancias mutuas ¡Xi -X j 1 tales que lim IX i - X j i/t =+ a

lim IXi - X j l /t = o .

b) La expansión infinita en tiempo finito (Mather-McGehee) . S61o

se presenta si n>
_ 4 . Para ello imaginemos un cuerpo C de masa peque
ña oscilando entre una bina-
ria AB y el cuarto cuerpo D A9
que se alejan de forma apro-
ximadamente lineal del c .d .m .
Las colisiones con D son de
tipo elástico . . Con la binaria AB son de tipo explosivo, extrayen-
do C tanta energía como se quiera, que posteriormente comunicará a
D y colapsando la binaria hacia energías potenciales mayores .
Tanto en el caso a) como en b) la demostración se lleva a ca-
bo por un proceso análogo al del T . de Sitnikov .
Un notable esfuerzo se dedica últimamente a la obtención de

59
órbitas periódicas del problema general de 4_ . cuerpos (e incluso pa
ra mayor número de cuerpos) por parte de Hadjidemetricu y su escue
la . Partiendo de un problema con 2 cuerpos no nulos y los restantes
nulos talesque .las masas nulas describen órbitas periódicas del pro
blema restringido con periodos conmensurables, se usa el método (b la con
tinuaci6n analítica prolongando la familia de órbitas periódicas a
masas no nulas .

§10 .- valores mayores de n . Soluciones particulares . Teorema de Moul -

ton. Experiencias numéricas . Caso n < 1000 . RecTular ización. Caso
ñk 10004

Para valores de n > 4 se sabe que el problema de n cuerpos pue-

de descomponerse en subsistemas o cgmulos. Si cuando t -+ oo esos sub

sistemas se alejan unos de otros en movimiento del tipo parabólico

o hiperbólico los centros de masas de cada subsistema tienden
hacia una solución homotética . Poco se sabe del movimiento interno
en cada cúmulo (excepto si está formado por 2 o 3 cuerpos o si su
movimiento es del tipo l) .
Para n > 4 y masas cualesquiera se dispone de algunas experien
cias en la búsqueda de órbitas periódicas (ver¢9) . Si todas las ma
sas excepto una son pequeñas (caso planetario) es posible la apli-
cación del teorema KAM (ver§5) .
Se desconoce el nó y forma de los p.e .r . (ver¢4) . Sin embargo
si se conocen para el caso alineado . En efecto, Moulton probó que
en el problema de n cuerpos existen para todos los valores de las
masas n!1/2 p .e .r . alineadas . El teorema se demuestra por inducción
respecto al nó de masas y por prolongación respecto al valor de la
última masa introducida .
Por lo que antecede sólo queda el recurso de las experiencias
numéricas . Este ha sido utilizado in extenso y ha contribuido en
gran parte al refinamiento de los métodos de resolución de e .d .o .

60
 Son corrientes los programas que integran las ecuaciones del
movimiento para varios cientos de cuerpos, sea directamente o por
métodos de Montecarlo . En todo caso debe pensarse que la integra
ci6n directa requiere pasos variables para cada cuerpo . En caso
contrario cada paso de integración significaría la evaluación de

(Z) distancias . En este problema jueganun papel preponderante

las binarias . Estas exigen un paso de integración muy pequeño y en

las proximidades de una colisión pueden tenerse errores notables .
Por ello es imprescindible efectuar un paso a variables regulari-
zadas para tratar este problema .
Una técnica debida a Heggie permite la regularización de todas
las colisiones binarias . Para ello se consideran = xi- x j como
yij
variables independientes y se efectua la regularización KS para
cada una de ellas . Sin embargo el nó de grados de libertad crece
de 3n a 2n(n-1), por lo que deja de ser práctico si n-> 4 . Indique
mos que en este caso el Hamiltoniano es polinomial .Para valores de
n mayores debe recurrirse a procedimientos aproximados mas groseros .
Una técnica adecuada puede ser la de estimación del campo creado
por las partículas alejadas y el estudio correcto sólo para los
más próximos . Otro método consiste en la aproximación de los tér-
minos per cantidades previamente almacenades en memoria . Esta
r12
1)
técnica ha sido usada por Hohl y le permite tratar cientos de mi
les de partículas . De esta manera puede conseguirse simular el com
portamiento de galaxias . Sin embargo el movimiento es inestable y
el escape de estrellas es continuo .

¢11 .- Técnicas de perturbaciones . Teorías planetarias . Movimientos

de la luna . Pequeños denominadores . Resotiancias .

En los problemas prácticos de cálculo de efemérides de astros

61
del sistema solar se rec=e a la integración numérica y a las teorías

de perturbaciones . En todo caso si la órbita real se aleja lentamen

te de una órbita Kepleriana, puede tomarse ésta como órbita de re-
ferencia y.dar sólo las desviaciones con respecto a la misma . En
ciertos casos es más .interesante tomar como órbita de referencia
otras soluciones, sean las del problema de los dos centros (como
en la teoría de Vinti de un satélite artificial en el campo terres
tre) o las órbitas intermedias de Garfinkel o Aksnes para

el mismo problema, o bien la órbita de variación de Hill para el

movimiento de la luna, etc .
Otra técnica consiste en facilitar en cada momento la órbita
Kepleriana ( o no) osculatriz . Este se define como la órbita keple

riana cuya posición y velocidad coinciden en cada momento con la

.
órbita real . Por supuesto que los elementos a,e,I,ui,n, M o - ta-1/3
serán funciones del tiempo . Sea c el vector formado por dichos
elementos . Se tendrá á = e Cp (a " t), donde el factor 6 indica que
seria ceros¡ no hubiera la perturbación, y que ésta se supone
pequeña . Partiendo de un valor inicial a se obtiene a, median
te a~ =a 0 + e tp ( a0, T ) dT . El métodm(le Picard deja de ser prác
J 0
tico ya que puede cp4e l , t) ser expresión di f f cil , . Se aplica un desa--
c
llo de Taylor limitado cp(a 1 ,t) = cp(a 0 ,t) + Dcp (a,,, t) E ~ tp(Q o ,T)dT
" n
dg .
con lo que a 2 =a 1 + e2 Dcp( ao . T) T cp(ao .Q) .
Jc0 J0
La iteración da lugar al método clásico de teoría de perturba-
ciones iniciado por Lagrange, Laplace y Poisson .
Los métodos de series de Lie, promedios, von Zeipel, etc . con-
sisten en técnicas iterativas para eliminar la perturbación del

hamiltoniano inicial H = Ho + CH, + e8H2 + . . . .con objeto de convertir

lo en un hamiltoniano independiente de las variables angulares .
Las sucesivas iteraciones y su convergencia se ven dificultadas

62
por la apari ción . d e pequeños denominadores . Por ejemplo, en el mé
todo clásico se tiene en cada paso que la función a integrar res-
pecto al tiempo es del tipo E At (p) exp i(k,q) y si (k,4) es
ke23

pequeño dicho término aparece multiplicado por un factor grande .

Este hecho implica que la mayor parte de las series utilizadas en
mecánica celeste tienen caracter asintótico y no son válidas uni-
formemente .
En las teorías planetarias se suele usar como pequeño parámetro
la masa de los cuerpos perturbadores . El estudio iniciado por
Laplace de la variación de e y w para los distintos planetas debi
do a sus perturbaciones (al menos en una aproximación de primer or
den) condujo a introducir la ecuación o polinomio característico .
La teoría de la luna es un ejemplo de la dificultad que encie-
rran los cálculos en mecánica Celeste .Por una parte ni su masa ni
la perturbación solar son pequeñas, con lo que los desarrollos que
se obtienen convergen lentamente . Por otra parte la precisión de
las medidas actuales (del orden de decenas de Cm . en la distancia
mediante uso de reflectores laser) obliga a utilizar teorías de o_r
den muy elevado . A título de ejemplo digamos que una teoría acep-
table de la luna que diera correctamente su posición con error de
unos pocos metros por siglo . requeriría tener un cuenta la influen-
cia de todos los planetas(excepto Plutón) y trabajar con desarro-
llos de Fourier conteniendo millones de términos .
Un tal trabajo requiere el uso de ordenadores en el manejo de
expresiones analíticas y en procesos analíticos de derivación, inte
graci6n y multiplicación de series .
El ajuste de las desviaciones de la posición de la luna respec
to a las que predice la teoría newtoniana mediante precisas y muy
numerosas observaciones de ocultaciones de estrellas por la misma
ha permitido constatar la disminución de la cons-

63
tante de la gravitáci6n y en breve se dispondrá de valores de la
velocidad de variación .
Finalmente indiquemos que la apárici6n de pequeños denominado-
rp `- está relacionada con el fen6menr -e resonancia . Hay que desta-
car las numerosas resonancias que se producen en el sistema solar,
tanto entre planetas como con satélites, asteroides o con anillos
de S£Lturno . Los asteroides evitan ciertos valores del movimiento
medio mientras que muestran preferencia por otros . Se está lejos
de tener una explicación satisfactoria de estos hechos, pero en to
do caso la resonancia parece tener un papel importante en la evo-
luci6n del sistema solar hacia su estado actual .

CONCLUSION : Wintner decia que toda generación de matemáticos ataca

el problema de 3 cuerpos ( y por extensión el de n) de una manera
peculiar y con herramientas propias . Hoy en día podemos afirmar que
métodos utilizados abarcan de la topología al estudio probabilís-
tico,de la teoría de puntos críticos a las aproximaciones diofánti
cas y del análisis no lineal a las técnicas numéricas . Además la
precisión alcanzada por la técnica en las mediciones ha permiti-
do que las medidas sean mucho más precisas que las teorías,obli-
gando a rehacerlas y a estudiar las bases físicas en que se apoyan .
Concluyamos con una anécdota . V . Szebehely autor del mejor
texto sobre el problema restringido y de numerosísimos artículos
y persona sumamente afable, ha felicitado las Navidades del 75 a
sus amigos y conocidos que trabajan en Mecánica Celeste diciendo :
Deseo fervientemente que usted "resuelva" en 1976 el problema de
3 cuerpos .

64
BIBLIOGRAFIA

Abraham,R ., Marsden,J .E . : Foundations of Mechanics , Benjamin,1967

Alexeyev,V .M . : "Quasirandom dynamical systems . I,II,III", Math .
USSR Sbornik,5,(1968),73-128 ; 6,(1968),505-560 ; 7,(1969)
1-43 .
Alexeyev,V .M . :" Sur 1'allure finale du mouvement dans le probléme
des trois corps", Actes Congrés Inter . Math ., Nice,1970,
Tome 2, 893-907, Gauthier-Villars .
Arnol'd,V .I . : "Small denominators and problems of stability of mo-
tion in Classical and Celestial Mechanics", Russ . Math .
Surveys, 18,(1963),85-192 .
Arnol'd,V .I . : "Instability of Dynamical Systems with Several De-
grees of Freedom", Soviet Math . Dokl ., 5,(1964),581-585 .
Arnol'd,V .I ., Avez,A . : Ergodic Problems of Classical Mechanics .
Benjamin, 1968 .
Bartlett,J .H . : "The Restricted Problem of Three Bodies, I,II", Mat .
Fys . Skr . Vid . Selsk . , 2, n4 7,(1964) ; 3, nó 1,(1965) .
Birkhoff,G . : Dynamical Systems , Am . Math . Soc ., 1927 .
Cabral, H .E . : "On the Integral Manifolds of the N-Body Problem",
Inventiones Math ., 20, (1973) 50-72 .
Chazy,J . : "Sur l'allure finale du mouvement dans le probléme des
3 corps quand le temps croit indéfiniment", Annales de
1'Ecole Normale Supérieure, 32 série, 39,(1922) 29-130 .
Easton,R . :"The Topology of the Regularized Integral Surfaces of
the 3-Body Problem",J . of Diff . Equa- 12, (1972) 361-
384 .
Easton,R . : "Some Topology of n-Body Problems", J . of Diff . Equa .,
19, (1975) 258-269 .
Escobal, P .R . : Methods of Orbit Determination , J .wiley, 1965 .

65
Gauss,K .F . : Theoria Motus , Dover (reprint) 1963
Hadjidemetriou, J .D . : "Families of Periodic Orbits in the General
N-Body Problem", en Long-Time Prediction in Dynamics ;
223-240, Ed . V. Szebehely, B.D . Tapley, Reidel, 1976 .
Hagihara, Y. : Celestial Mechanics , Vol . 1, MIT Press, 1970 ; Vol .
2, MIT Press, 1972 ; Vol. 3, Japan Society for the Pro
motion of the Science, 1974 ; Vol . 4, Japan Soc ., 1975 .
Heggie,D . : "A Global Regularisation of the Gravitational'N-Body
Problem", Celestial Mechanics , 10, (1974) 217-242 .
Hénon, M. : "Exploration numérique du probléme restreint . III,IV",
Bulletin Astronomique, série 3, 1, fasc . 1, 57,(1966),
fasc . 2, 49 (1966)
Hoh1,F . : "N-Body Simulations of disks", en Dynamics of Stellar
Systems, A . Hayli (editor), IAU Symposium, nó 69,
p . 349, Reidel, 1975 .
Kolmogoroff, A .N . : "General Theory of Dynamical Systems and Cla-
ssical Mechanics",Proceed . 1954, Inter . Cong . Math (en
ruso) . Traducido al inglés en Abraham, apéndice D .
Leontovich, A .M . : "On the Stability of the Lagrange Periodic Solu-
tions for the Reduced Problem of the Three Bodies",
Soviet Math . Doklady, 3, (1965) 425-430 .
Losco, L . : "Sur un invariant,intégral du probléme des n-corps
conséquence de 1'homogéneité du potentiel", en The
Stability of the Solar System and of Small Stellar Sys -Í
tems . Y. Kozai (Editor), IAU Symposium N4 62, 249-255,
Reidel, 1974 .
McGehee, R . : "Triple Collision in the Collinear Three-Body Pro- _
blem", Inventiones Math ., 27, (1974) 191-227 .
McGehee, R ., Mather, J .N . : "Solutions of the collinear four body
problem which become unbounded in finite time", Batte-
lle Recontres, 1974 Proceedings .
Marchal, Ch . : "Qualitative Methods and Results in Celestial Mecha-

66
 nics", en Long-Time Predictions in Dynami cs, Ed . V .
Szebehly, B .D . Tapley, 181-208, Reidel, 1976 .
Markus,L .,Meyer,K .R . : "Generic Hamiltonian Dynamical Systems are
neither Integrable nor Ergodic", Memoirs of the A .M .S .,
144, 1974 .
Moser,J . : "Nonexistence of Integrals for Canonical Systems of
Differential Equations" Comm . on Pure and Appl .Math .,
8, (1955) 409-436 .
Moser, J . : "On invariant curves of area-preserving mappings of
an annulus", Nach . Akad . Wiss . GSttingen . Math . Phys .
Klass ., 1, (1962), 1-20 .
Moser,J . : "Regularization of Kepler's Problem.and the Averaging
Method on a Manifold", Comm . on Pure & Appl . Math .,
23, (1970) 609-636 .
Moser,J . : Stable and Random Motion in Dynamical Systems , Prince-
ton Univ . Press, 1973 .
Moulton, F .R . : "The Straight line solutions of the problem of
N-bodies", Annals of Math ., 12, (1910), 1-17 .
Palmore,J . : "Classifying Relative Equilibria,I,II", Bull . A .M .S .,
79,(1973) 904-908 ; 81, (1975), 489-491 .
Palmore,J . : "Classifying Relative Equilibria,III", Letters in
Math . Phys ., 1, (1975) 71-73 .
Poincaré,H . : Les méthosies nouvelles de la Mécanique Céleste .
Gauthier-Villars, 1892, 1893, 199 .
Robinson,C . : "Lectures on Hamiltonian Systems", Monografías de
Matemática, Guanabara, 1971 .
Rüsmann, H . : "Uber invariante Kurven differenzierbarer Abbildun-
gen eines Kreisringes% Narichten Akéd . Wiss . G2ittin-
gen Math . Phys . Klass ., 1970 .
Saari, D .G . : "Expanding gravitational systems", T .A .M .S ., 156,
(1971) 1-22 .
Siegel, C .L ., Moser,J .K . : Lectures on Celestial Mechanics , Sprin-
ger-Verlag, 1971 .
Simó, C . : "La variedad de órbitas Keplerianas y la teoría gene- I
ral de perturbaciones" Tesis, Universidad de Barcelo-
na, 1974 .
Simó, C . : "Aspects Topologiques en Mécanique Céleste", Table ron-
de de Mec . Cel ., Institut Henri Poincaré, Paris, 21,
mai 1975 .
Simó, C . : "Posiciones de equilibrio relativo del problema de
3 + 1 cuerpos y su evolución", comunicación en la
Reunión Anual de Matemáticos Españoles, Málaga, 1976 .
Sitnikov, K .A . : "The existence of oscillatory motions in 3-body
problem", Soviet Phys . Doklady, 5, (1961), 647-650 .
Smale,S . : "Topology and Mechanics", Inventiones Math ., 10, (1970)1
305-331 ; 11, (1970) 45-64 .
Souriau,J .M . : Structure des Svstémes dynamictues , Dunod, 1970
Sternberg,S . : Celestial Mechanics , Benjamin, 1969
Stiefel, E .L ., Scheifele, G. : Linear and Regular Celestial Me~
chanics , Springer-Verlag, 1971
Szebehely, V . : Theory of Orbits . The restricted problem of 3 bo-
dies , Academic-Press, 1967 .
Waldvogel, J . : "Triple Collision" en Long-Time Predictions in Dy -
namics, 241-258, Ed . V . Szebehely, B.D . .Tapley, Rei-
del, 1976 .
Whittaker, E .T . : A Treatise on the Analvtical Dynamics of parti -
des and rigid bodies , Cambridge Univ . Press, 1973 .
Wintner, A . : The Analvtical Foundations of Celestial Mechanics ,
Princeton Univ . Press, 1941 .

68