Floquet

4.
Estabilidad
4.1. Introducci on
Una vez obtenidas las soluciones tipo breather en una red de osciladores, cabe
preguntarse si son estables, es decir, si persisten bajo perturbaciones.
Ya en su prueba de existencia, MacKay y Aubry [132] apuntaron que era muy
probable que los breathers con un sitio excitado (1-site breathers o 1sB para abre-
viar) fueran linealmente estables, es decir, que una perturbaci on aplicada sobre
ellos crezca, como mucho, linealmente en el tiempo. Estos resultados fueron pos-
teriormente demostrados por Marn y Aubry [19, 142] y extendidos para sistemas
innitos por MacKay y Sepulchre [133].
Para el an alisis de estabilidad de las soluciones peri odicas se utiliza la teora
de Floquet y las propiedades de la ecuaci on de Hill. En este captulo se recopi-
lan algunos resultados que son importantes para el desarrollo de esta tesis. Para
mayor claridad, las demostraciones de la mayora de los resultados se encuentran
en el apendice A de este trabajo. Para una explicaci on m as completa de ciertos
resultados, el lector puede referirse a la tesis doctoral de Jose Luis Marn en la que
se encuentra una magnca recopilaci on de las bases de la teora de estabilidad
en redes anarm onicas. Adem as, en este captulo y en los Apendices A y B expon-
dremos algunos resultados, que, al menos, no han podido ser encontrados por el
autor. Estos se pueden resumir en:
Demostraci on de la hermiticidad y simplectitud del operador de Newton.
Estructura de los autovectores del operador de Floquet cuando no existe
degeneracion.
Estructura de los autovectores del operador de Floquet en el caso degenerado.
Tambien se exponen resultados recientemente publicados por el autor, tales
como:
Relaci on entre los autovalores del operador de Floquet y los del jacobiano.
Demostraci on de la estabilidad de multibreathers para acoplos peque nos
bas andose en la teora de perturbaciones degenerada.
28 4. Estabilidad
Toda la teora que se expone en esta tesis se reere a estabilidad lineal. Existe,
adem as, un estudio de estabilidad en el sentido de Nekhoroshev (es decir, estabi-
lidad exponencial) realizado por Bambusi [27], seg un el cual, para acoplamientos
debiles, los breathers linealmente estables conservan la energa localizada durante
tiempos exponencialmente largos.
4.2. Ecuaciones din amicas perturbadas. Operador de
Newton
Recordemos las ecuaciones dinamicas de una red hamiltoniana de osciladores
acoplados:
F
n
(u) u
n
+ V
(u
n
) + C
i
W
(u
n
u
n+i
) = 0, (4.1)
donde u
n
se considera una funcion d-dimensional y n Z
d
. Si se introduce una
perturbaci on en la red, es decir, u u + , esta cumplir a la ecuaci on:
n
(t) + V
(u
n
(t))
n
(t) C
i
W
(u
n+i
(t) u
n
(t))(
n+i
(t)
n
(t)) +
+C
i
W
(u
n
(t) u
ni
(t))(
n
(t)
ni
(t)) = 0. (4.2)
Esta ecuacion describe ademas los modos lineales del sistema cerca de la so-
lucion u = {u
n
(t)} de las ecuaciones din amicas de la red (4.1). Si restringimos la
interacci on al caso de potencial de acoplo lineal y de pr oximos vecinos, nos queda
una expresi on m as simplicada:
n
(t) + V
(u
n
(t))
n
(t) + C(2
n
(t)
n+1
(t)
n1
(t)) = 0. (4.3)
La ecuaci on (4.2) se puede expresar de manera m as compacta:
N(u(t)) = 0, (4.4)
donde {
n
(t)} y
N(u(t))
d
2
dt
2
+ V
(u(t)) + CJ(u(t)), (4.5)

con [V
(u)]
nm
=
nm
V (u
n
), siendo J la matriz de acoplo. El operador N recibe
el nombre de operador de Newton. Si C
2
, el estudio del espectro del operador
de Newton nos ofrece informaci on sobre la estabilidad. Si, adem as, E
2
s
(
b
), el
operador de Newton es el jacobiano (o derivada de Frechet) correspondiente a las
ecuaciones din amicas del sistema (4.1), es decir,
u
F = N(u(t)).
La ecuaci on de autovalores del operador de Newton es, l ogicamente,
4.2. Ecuaciones din amicas perturbadas. Operador de Newton 29
N(u(t)) = E, (4.6)
por lo que la ecuaci on de las perturbaciones (4.2) se corresponde con el caso par-
ticular E = 0. Sin embargo, el estudio completo del espectro de N ofrece informa-
cion adicional de gran interes. En el se basa la teora de bandas de Aubry [19] y el
teorema de estabilidad de multibreathers [12].
El operador de Newton posee una serie de propiedades, que se demuestran
Apendice A. Estas propiedades tienen consecuencias importantes sobre sus auto-
valores y autofunciones, y, en consecuencia, sobre la estabilidad de las soluciones.
Pasemos a enumerarlas:
1. N conmuta con el operador de inversi on temporal. Esto implica que
las autofunciones de N pueden ser escogidas como simetricas o antisimetricas
en el tiempo, es decir,
(t) = (t). (4.7)
2. N conmuta con el operador de traslaci on de un periodo. Seg un el
teorema de Bloch, las autofunciones de N se pueden expresar de la forma
denominada funci on de Bloch:
(t) = e
jt/T
v(t), (4.8)
donde v(t) es una funci on peri odica de periodo T = 2/
b
con C.
3. N es real. Esto implica que dos autofunciones conjugadas tienen sus auto-
valores conjugados, es decir,
N
= E
. (4.9)
4. N es hermtico. Esta propiedad se basa en el hecho de que N es simplecti-
co, es decir, el producto simplectico de dos autofunciones cualesquiera se
mantiene constante en el tiempo. Si llamamos:
(t) =
_
(t)
(t)
_
, donde (t)

(t), (4.10)
el producto simplectico se dene como la suma de las areas de los subes-
pacios denidos por cada par de variables conjugadas, o, dicho en terminos
matem aticos:
[
1
(t),
2
(t)] =
n
(
1n
(t)
2n
(t)
2n
(t)
1n
(t)). (4.11)
Esta propiedad implica que los autovalores de N son reales.
30 4. Estabilidad
4.3. An alisis de Floquet
4.3.1. El operador de Floquet. Autovalores y autofunciones
Cuando la soluci on del sistema u(t) es periodica, el an alisis de estabilidad
requiere la integraci on a lo largo de un solo periodo. Este modo de estudiar la
estabilidad se denomina an alisis de Floquet [97]. Para ello, se considera el siguiente
mapa:
_
(T)
(T)
_
= F
o
_
(0)
(0)
_
, (4.12)
o, en distinta notaci on:
(T) = F
o
(0). (4.13)
Toda la informaci on sobre la estabilidad se encuentra en el operador F
o
, tam-
bien llamado operador de Floquet. A su representaci on se le llama matriz mono-
droma o, simplemente, monodroma
1
.
Si el sistema tiene N partculas, la monodroma tendr a 2N autovectores. Estos
seran los correspondientes al operador de Newton con E = 0 y t = 0 o t = . Por
la propiedad 2 del operador de Newton (teorema de Bloch) se tiene, trivialmente,
que:
(T) = e
j
(0), (4.14)
es decir, que los autovalores de la monodroma son de la forma = e
j
, ( C).
A se le conoce con el nombre de multiplicador de Floquet, mientras que j recibe
el nombre de exponente de Floquet, siendo el argumento de Floquet.
Como se demuestra en el apendice A, el operador de Floquet es real y sim-
plectico. Esto implica que si es autovalor, tambien lo seran
, 1/ y 1/
. Pero
este hecho tiene a un otra consecuencia importante: como, seg un la teora de Flo-
quet, una orbita u(t) es estable si ninguno de los multiplicadores de Floquet tiene
m odulo mayor que uno, una condici on necesaria y suciente para la estabilidad
lineal de nuestro sistema es que todos los multiplicadores de Floquet esten sobre
la circunferencia unidad. En otras palabras, los argumentos de Floquet deben
ser reales para que el sistema sea linealmente estable.
En la gura 4.1 se ilustra este hecho con un ejemplo de breather estable y otro
de breather inestable.
4.3.2. Signatura de Krein
Un concepto muy util para el estudio de la estabilidad lineal es el de signatura
de Krein, la cual se dene para cada argumento de Floquet de la siguiente forma
1
Del griego oo: uno, oo: ciclo.
4.3. An alisis de Floquet 31
(a) (b)
1 0.5 0 0.5 1
1
0.5
0
0.5
1
Re()
I
m
(
)
1 0 1
1.5
1
0.5
0
0.5
1
1.5
Re()
I
m
(
)
Fig. 4.1: Multiplicadores de Floquet para un breather en un potencial de Morse
con frecuencia
b
= 0,8 y acoplamiento (a) C = 0,05 y (b) C = 0,14.
El primer caso corresponde a un breather estable (todos los autovalores
estan en la circunferencia), mientras que el segundo corresponde a un
breather inestable (existen dos autovalores con m odulo distinto de uno).
Las cruces (+) indican signatura de Krein positiva, las aspas () sig-
natura de Krein negativa y los crculos () signatura de Krein nula (ver
Seccion 4.3.2).
32 4. Estabilidad
[15]:
() = sgn([Re((t)), Im((t))]) = sgn
_
j
i
(
i
(t)
i
(t)
i
(t)
i
(t))
_
, (4.15)
donde (t) es la autofunci on correspondiente a . Como el operador de Newton es
simplectico, la signatura de Krein se conserva y est a bien denida. Si = 0 o ,
las autofunciones son reales (pues el multiplicador correspondiente tambien lo es)
y la signatura de Krein es cero. Denimos la signatura de Krein de un conjunto
de puntos, interior a una curva cerrada, como la suma de la signatura de Krein de
los puntos de dentro de la misma, es decir, K =
(). Este valor se conserva

siempre que ning un punto cruce la curva contorno.
Cuando los par ametros del sistema varan, la unica manera de que un sistema
se desestabilice es que colisione y se salga del crculo unidad un cuadruplete de
autovalores (si = 0 y = ) o una pareja (si = 0 o = ). En ese caso, antes
de una inestabilidad deben coincidir los autovalores implicados. Sin embargo, no
siempre que coincidan dos autovalores se produce una inestabilidad. Para saber
cuando se produce una colisi on se introduce el criterio de Krein [15], seg un el cual
solo es posible una inestabilidad si la signatura de Krein de los autovalores que
coinciden son distintas. Este criterio es crucial para obtener que los breathers con
un s olo sitio excitado (1SB) son estables, siempre que la constante de acoplo sea
los sucientemente peque na, tal como veremos en secciones posteriores.
4.3.3. Bifurcaciones
Cuando un breather se perturba, la soluci on inicial, u(t), que era de periodo T =
2/
b
, se transforma en la nueva soluci on u(t) = u(t)+(t), que no tiene porque ser
periodica, ya que la perturbaci on (t) no lo es. Recordemos que por el teorema de
Bloch (4.8), la perturbaci on puede expresarse como (t) = exp(jt/T)v(t), donde
v(t) tiene la misma periodicidad del breather. As, la soluci on perturbada se puede
expresar de la forma:
u(t) = u(t) + e
j
b
t/2
v(t), (4.16)
es decir, la solucion perturbada oscilar a con dos frecuencias: una de ellas ser a la
del breather,
b
, mientras que la segunda ser a
b
= (/2)
b
mod
b
.
Cuando se produce una colisi on de los autovalores, la soluci on original u(t)
puede volverse inestable y, en un entorno innitamente cercano a la bifurcaci on, la
soluci on estable emergente ser a igual a la perturbada. Por ello, existen tres tipos
distintos de bifurcaciones debidas a cambios de estabilidad:
1. Bifurcaci on arm onica. Se produce cuando una pareja de autovalores co-
lisiona en = 0. En ese caso,
b
= 0, por lo que u(t) = u(t) + v(t) y las
soluciones emergentes tienen la misma frecuencia.
4.4. Estabilidad lineal de los breathers 33
2. Bifurcaci on subarm onica. Se produce cuando una pareja de autovalores
colisiona en = . En ese caso,
b
= (2k + 1)
b
/2 con k N, por lo que
u(t) oscilar a con periodo 2T. En otras palabras, se produce una duplicaci on
de periodo.
3. Bifurcaci on oscilatoria o de Hopf. Se produce cuando un cuadruplete de
autovalores se sale del crculo, es decir, dos parejas conjugadas de autovalores
se salen en = 0 o = . A este tipo de colisiones tambien se les conoce
como Krein crunches (o colisiones de Krein). Si el angulo en el que se
produce la colisi on es distinto de p/q con p, q Z, la frecuencia de la
solucion emergente
b
es inconmensurable con la frecuencia de la soluci on
original
b
, por lo que la soluciones emergentes ser an cuasi-periodicas. En
otras palabras, se produce una bifurcaci on oscilatoria.
Las soluciones cuasi-peri odicas en sistemas hamiltonianos no se pueden cal-
cular con los metodos desarrollados en el captulo 3. Es m as, el teorema de
existencia las prohbe, pues siempre existen tres n umeros p, q, r Z tales que
p
b
+ q
b
= r
ph
, siendo
ph
las frecuencias de los fonones. Sin embargo,
este tipo de soluciones s se encuentran en el contexto de la ecuacion DNLS
[111, 112] y del modelo de Holstein adiab atico [19]. Recientemente, Bambu-
si y Vella [28] han demostrado la existencia de breathers cuasi-peri odicos,
bas andose en el concepto de lmite anti-continuo. Esta prueba de existencia
es valida siempre que el sistema tenga una integral de movimiento distinta
del hamiltoniano.
En sistemas disipativos se pueden obtener facilmente soluciones cuasiperio-
dicas [145], ya que las resonancias con los fonones no producen la desaparici on
del breather.
Un esquema de como se visualizan estas bifurcaciones en el crculo unidad se
encuentra en la gura 4.2.
4.4. Estabilidad lineal de los breathers
Una vez vistas las herramientas para el an alisis de la estabilidad lineal, se puede
proceder al estudio de la estabilidad de las soluciones tipo breather. Para ello, se
parte de la forma del espectro del operador de Floquet en el lmite anti-continuo
y se observa c omo vara cuando se introduce el acoplamiento entre las partculas.
4.4.1. An alisis de estabilidad lineal en el lmite anti-continuo
Supongamos una red nita de N osciladores en la que se induce un breather
tal que en el lmite anti-continuo posea un solo sitio excitado. En esta situaci on,
los osciladores se desacoplan de forma que el excitado se describe mediante una
ecuacion de Hill:
34 4. Estabilidad
(a) (b)
(c) (d)
(e) (f)
Fig. 4.2: Esquemas para bifurcaci on arm onica, (a) y (b), subarm onica, (c) y (d), y
oscilatoria, (e) y (f). Las guras de la izquierda muestran los autovalores
antes de colisionar y las de la derecha los muestra despues de la colisi on,
es decir, en una situacion de inestabilidad.
o
(t) + V
(u
o
(t))
o
(t) = 0, (4.17)
mientras que los osciladores en reposo se describen por:
n
(t) +
2
o
n
(t) = 0, (4.18)
donde
o
es la frecuencia de las oscilaciones de peque na amplitud.
La ecuaci on (4.17) tiene como soluci on
o
(t) = u(t) y
o
(t) = u
u/
b
,
como se puede comprobar derivando las ecuaciones din amicas respecto al tiempo
y respecto a la frecuencia. Estas soluciones se denominan modo de fase y modo
de crecimiento, respectivamente. El primer modo aparece debido a que el sistema
es invariante bajo traslaci on temporal. El segundo, en cambio, debe su nombre a
que una perturbaci on en su direccion hace saltar la orbita del oscilador a otra de
distinta frecuencia. Como se puede ver en el apendice A, el modo de crecimiento
tiene la forma u
= (t) + t u(t), donde (t) es una funcion peri odica. Por tanto,
el modo de crecimiento diverge de manera lineal. Por este motivo, se dice que el
modo de crecimiento es un modo marginal, pues una perturbaci on en su direccion
solo diverge linealmente, lo cual es estable en sentido fsico (pero no en sentido
matem atico), ya que en el sistema completo da lugar a otra soluci on pr oxima
estable.
Como se demuestra en el apendice A, ambos modos se corresponden con un
bloque de Jordan en el operador de Floquet con un multiplicador doblemente
degenerado cuyo valor es +1. Este bloque de Jordan no es diagonalizable, por lo
que el modo de crecimiento no es autovector y no se puede obtener mediante las
rutinas numericas usuales de diagonalizaci on de matrices.
En el caso de la ecuacion (4.18), tenemos el autovector
n
(t) = e
j
o
T
n
(0),
que estar a N 1 veces degenerado. El argumento de Floquet correspondiente
sera =
o
T mod 2. Estos modos
n
(t) se corresponden con las oscilaciones
lineales.
En el caso de un multibreather con p sitios excitados, tendremos que el modo de
fase y de crecimiento estar an p veces degenerados, mientras que los modos lineales
lo estar an N p veces.
Teniendo en cuenta estos resultados, se puede calcular la signatura de Krein
de los autovalores. Los autovalores en +1 tienen signatura de Krein nula, mientras
que los autovalores en =
o
T tienen signatura = +1 si (0, ) y = 1
si (, 2) (ver Apendice A). Recordemos que la signatura de Krein es nula en
= 0 y = . Esto implica que, para que un breather sea estable, debe cumplirse
la siguiente condici on de resonancia:
n
2
b
=
o
, n = 0, 1, 2, 3, ... (4.19)
Si n es par, se reobtienen las condiciones de no resonancia impuestas por el
teorema de MacKayAubry (resonancias arm onicas). Si n es impar, las resonancias
se llaman subarm onicas, las cuales son equivalentes a las resonancias parametricas
36 4. Estabilidad
de la ecuaci on de Hill [137]. En resumen, siempre que se cumpla la condici on de
no resonancia (4.19), la soluci on pr oxima al lmite anti-continuo es estable.
4.4.2. Efecto del acoplamiento en la estabilidad
1site breathers
Cuando las partculas del sistema se acoplan, el autovalor doble en +1 se man-
tiene, pues u(t) sigue siendo soluci on del problema de autovalores. El resto de los
modos seguir an correspondiendo a oscilaciones de peque na amplitud, siempre que
el acoplo sea bajo. Por ello, los autovalores situados en = 0 para C = 0 se expan-
dir an formando una banda: la banda de fonones. En el caso de un acoplamiento
de pr oximos vecinos, la banda viene dada por la relaci on de dispersi on:
ph
=
_
2
o
+ 4C sin
2
(q/2), < q < . (4.20)
En sistemas innitos estas bandas forman un conjunto denso [133]. La signatura
de Krein de la banda del semiplano superior es distinta de la del semiplano inferior.
Por ello, en principio, s olo hay dos modos de que el sistema se desestabilice
2
:
Las bandas se solapan en +1. Esto hara que el breather dejara de existir ya
que no se cumplira la condici on de no resonancia del teorema de MacKay
Aubry (la frecuencia del breather resonara con la de los fonones).
Las bandas se solapan en -1. En este caso, se tendra una resonancia sub-
arm onica, la cual dara lugar a una bifurcaci on de duplicaci on de periodo,
tal como vimos en la secci on anterior.
Sin embargo, pueden ocurrir dos fen omenos que modican este panorama. El
primero se debe al hecho de que la condici on de no resonancia subarm onica es
necesaria pero no suciente para la desestabilizaci on del breather. Esto ocurre,
por ejemplo, cuando se considera una red con un potencial de Morse onsite y un
acoplamiento arm onico y atractivo.
El otro fen omeno, el cual es de importancia capital para la posible movilidad de
los breathers, consiste en que, tal como muestran Baesens et al [24], existen ciertos
tipos de redes en las cuales, a partir de un valor umbral del acoplamiento, aparecen
modos exponencialmente localizados, los cuales est an aislados del conjunto denso
de valores que forma la banda de fonones. No obstante, pueden existir tambien
modos localizados dentro de la banda de fonones.
El origen fsico de estos modos localizados radica en el hecho de que al aumentar
el acoplamiento del sistema, el n umero de partculas que se excitan va aumentando,
lo cual produce una disminucion de los grados de libertad correspondientes a las
oscilaciones de peque na amplitud.
2
Este hecho es independiente de que la constante de acoplo sea positiva o negativa.
Pues bien, estos modos localizados pueden colisionar en = 0 produciendo
bifurcaciones arm onicas. Estas bifurcaciones no destruyen al breather pues la con-
dicion de no resonancia del teorema de existencia se sigue cumpliendo. Tambien
pueden hacerlo para otros valores de dando lugar al escenario de bifurcaciones
expuesto con anterioridad.
El valor en el que se localizan los modos depende del modelo y, por el momento,
no es posible predecirlo. Adem as, existen casos en los que estos modos localizados
ni siquiera aparecen.
En la gura 4.3 se pueden encontrar ejemplos de estos efectos.
Evoluci on de los autovalores
A partir de la relaci on de dispersion de los fonones se puede conocer de manera
cualitativa la evoluci on con el acoplo de los multiplicadores de Floquet, teniendo
en cuenta la frecuencia del breather y el signo de la constante de acoplamiento.
La banda de fonones se expande entre los valores
o
y
1
, que vienen dados por:
o
=
2
o
b
(mod 2);
1
=
2
o
1 +
4C
2
o
(mod 2). (4.21)
Como se puede observar,
1
>
o
si C > 0 y
1
<
o
si C < 0. Supongamos un
potencial duro y un acoplamiento positivo. En ese caso, habr a dos comportamientos
diferentes dependiendo de la frecuencia. As, si
b
> 2
o
,
o
< , lo cual implica
que el borde inferior de la banda est a a la derecha del borde superior. Si
o
<
b
< 2
b
,
o
> , y el borde inferior estar a a la izquierda. En consecuencia,
cuando se localiza un modo, tendr a un patr on de vibraci on similar al del breather.
En el caso de un potencial duro, la localizaci on provendr a entonces del borde
superior de la banda. As, en el primer caso, los modos localizados colisionaran en
1, mientras que en el segundo, lo haran en +1.
En el caso de un potencial blando (con acoplamiento positivo), el comporta-
miento es similar. Si
b
/
o
(2/2k +1, 1/k) con k = 1, 2, 3, ..., entonces
o
< . Si
b
/
o
(1/k +1, 2/2k +1), entonces
o
> . Ahora, los modos localizados saldr an
del borde inferior de la banda, por lo que, en el primer caso, las colisiones tienen
lugar en +1, mientras que en el segundo, las colisiones ocurren en 1.
En el caso de acoplamiento negativo, la posici on relativa de
1
y
o
se invierte
con respecto al caso C > 0. Esto implica que la expansi on de la banda sea contraria
al caso anterior. Adem as, los autovalores que se localizan lo hacen por el borde
superior de la banda si el potencial es blando y por el borde inferior, si es duro.
Esto se debe a que los patrones de vibraci on se intercambian cuando vara el signo
del acoplamiento.
En las guras 4.4 a la 4.7 se ilustran los diferentes casos con diversos ejemplos.
38 4. Estabilidad
(a) (b)
1 0.5 0 0.5 1
1
0.5
0
0.5
1
Re()
I
m
(
)
1 0.5 0 0.5 1
1
0.5
0
0.5
1
Re()
I
m
(
)
(c) (d)
1 0.5 0 0.5 1
1
0.5
0
0.5
1
Re()
I
m
(
)
1 0.5 0 0.5 1
1
0.5
0
0.5
1
Re()
I
m
(
)
Fig. 4.3: Multiplicadores de Floquet en distintos casos. (a) Breather con frecuen-
cia
b
= 0,8 en el lmite anti-continuo. (b) Breather con la misma fre-
cuencia, potencial de sine-Gordon y acoplamiento C = 0,10. Se puede
observar la expansi on de la banda de fonones hacia 1. (c) Breather
con los mismos par ametros que en el caso anterior pero con potencial
de Morse. Se puede observar, adem as de la expansi on de la banda de
fonones, la existencia de modos localizados que se mueven hacia +1. (d)
Breather con potencial de Morse, frecuencia
b
= 0,6 y acoplamiento
C = 0,04. Se observa ahora una expansi on de la banda hacia +1 y la
existencia de modos localizados que se mueven hacia 1.
Multibreathers
Recordemos que en el caso de los multibreathers en el lmite anti-continuo,
existen p parejas de autovalores en +1, mientras que hay Np parejas conjugadas
en = 0. Cuando el sistema se acopla, el efecto sobre los autovalores correspon-
dientes a oscilaciones de peque na amplitud es similar al caso de los 1sB, de forma
que la banda de fonones evoluciona de manera similar.
Sin embargo, las p parejas de autovalores en +1 tienen un comportamiento
distinto. De entrada, sabemos que hay una pareja (la correspondiente a los modos
de fase y crecimiento) que siempre permanecer a en ese punto. La inc ognita es saber
que ocurre con las p 1 parejas restantes, ya que cuando el sistema se acopla, se
rompe la degeneraci on. Obviamente, se pueden presentar dos casos:
Que los autovalores se salgan del crculo, volviendose el sistema inestable.
Que los autovalores se queden en el crculo, y se comporten de manera pare-
cida a los modos localizados que mostramos en el caso de los 1sB. Adem as,
los modos localizados pueden colisionar con los modos extendidos de la ban-
da de fonones (o con otros modos localizados que hayan partido de esta) y
producir bifurcaciones oscilatorias debidas a colisiones de Krein (es decir, en
= 0 o = ).
Este comportamiento se puede predecir mediante el Teorema de Estabilidad de
Multibreathers (TEM), que se encuentra explicado en la Secci on 4.6 y demostrado
en el Apendice B. Este teorema puede aplicarse facilmente siempre que el poten-
cial de acoplamiento sea arm onico y los sitios excitados en el multibreather sean
adyacentes.
A continuaci on, se exponen los resultados obtenidos mediante la aplicaci on del
TEM para 2site breathers (2sB), es decir, breathers con dos sitios adyacentes
excitados en el lmite anti-continuo. Los resultados coinciden con los predichos por
un teorema de Aubry
3
. As, existen varios casos cuando se introduce el acopla-
miento en la red, dependiendo de que el potencial sea blando o duro, que el 2sB
este en fase ( = {. . . , 1, 1, . . .}) o en anti-fase ( = {. . . , 1, 1, . . .}), o que el
acoplamiento sea atractivo (C > 0) o repulsivo (C < 0).
Los resultados para acoplamiento atractivo se resumen en la siguiente tabla (el
caso de acoplamiento repulsivo es el opuesto):
Codigo Potencial blando Potencial duro
1 1 Inestable Estable
1 -1 Estable Inestable
Este comportamiento se verica en las guras 4.84.11.
3
Este teorema no esta publicado pero s esta enunciado y demostrado para acoplamientos
armonicos en la tesis de Jose Luis Marn [137]
40 4. Estabilidad
El TEM tambien se puede aplicar para predecir el comportamiento de 3-site
breathers, es decir, breathers con tres sitios adyacentes excitados en el lmite anti-
continuo. La siguiente tabla muestra el comportamiento para C > 0:
Codigo Potencial blando Potencial duro
1 1 1 Inestable Estable
1 -1 1 Estable Inestable
1 1 -1 Inestable Inestable
Se puede establecer, en general, que si en un n-site breather (es decir, un mul-
tibreather con n sitios adyacentes excitados), todos los valores del c odigo distintos
de cero son iguales a +1, el comportamiento es identico al del 2sB con codigo
= {. . . , 1, 1, . . .}, mientras que si existe una alternancia entre 1 y 1, el com-
portamiento ser a equivalente al 2sB con = {. . . , 1, 1, . . .}. En cualquier otro
caso, el multibreather es inestable para cualquier tipo de potencial y signo de la
constante de acoplo.
(a) (b)
0.05 0.1 0.15 0.2 0.25 0.3 0.35
0
0.5
1
1.5
2
2.5
3
C
0.05 0.1 0.15 0.2 0.25 0.3 0.35
0.8
0.9
1
1.1
1.2
1.3
|
|
C
(c) (d)
0 0.02 0.04 0.06 0.08 0.1
0
0.5
1
1.5
2
2.5
3
C
0 0.02 0.04 0.06 0.08 0.1
0.5
1
1.5
2
2.5
3
|
|
C
Fig. 4.4: Evoluci on con el acoplo de los argumentos de Floquet (izquierda) y del
m odulo de los autovalores (derecha) para un breather con potencial de
Morse y acoplamiento positivo. Las guras superiores corresponden a
una frecuencia
b
= 0,8 y las inferiores a
b
= 0,6. Como se puede com-
probar, el breather con frecuencia
b
= 0,8 es estable hasta C 0,13,
momento en el que un modo localizado experimenta una bifurcaci on
arm onica. Sin embargo, para C 0,24, el breather recupera su estabili-
dad. Adem as, se observa a partir de C 0,17 la formaci on de burbujas
de inestabilidad va bifurcaciones oscilatorias y subarm onicas. Estas
ultimas tienen su origen en el tama no nito del sistema. En el caso
b
= 0,6, dos modos se localizan y evolucionan hacia = . Se produce
una peque na colisi on en C 0,053 y el sistema sigue siendo estable
hasta que dos modos localizados colisionan en = 0 para C 0,07.
42 4. Estabilidad
(a) (b)
0.08 0.06 0.04 0.02 0
0
0.5
1
1.5
2
2.5
3
C
0.08 0.06 0.04 0.02 0
0.98
0.985
0.99
0.995
1
1.005
1.01
1.015
1.02
|
|
C
(c) (d)
0.15 0.1 0.05 0
0
0.5
1
1.5
2
2.5
3
C
0.15 0.1 0.05 0
0.6
0.8
1
1.2
1.4
1.6
|
|
C
m odulo de los autovalores (derecha) para un breather con potencial de
Morse y acoplamiento negativo. Las guras superiores corresponden a
una frecuencia
b
b
= 0,6. El breather con
frecuencia
b
= 0,8 es siempre estable hasta C 0,08 y la banda
de fonones se expande hacia = 0. En el caso
b
= 0,6, el breather es
estable hasta que las bandas se solapan en C 0,055. En ambos casos,
solo se observan modos localizados cuando el acoplamiento es alto y ya
esta cerca de producirse la resonancia con los fonones.
(a) (b)
0.05 0.1 0.15 0.2 0.25 0.3
0
0.5
1
1.5
2
2.5
3
C
0.05 0.1 0.15 0.2 0.25 0.3
0.98
0.985
0.99
0.995
1
1.005
1.01
1.015
1.02
|
|
C
(c) (d)
0 0.2 0.4 0.6 0.8 1
0
0.5
1
1.5
2
2.5
3
C
0 0.2 0.4 0.6 0.8 1
0.98
0.985
0.99
0.995
1
1.005
1.01
1.015
1.02
|
|
C
Fig. 4.6: Evoluci on con el acoplo de los argumentos de Floquet (izquierda) y
del modulo de los autovalores (derecha) para un breather con potencial
4
duro y acoplamiento positivo. Las guras superiores corresponden a
una frecuencia
b
b
= 2,5. Las dos familias de
breathers son siempre estables. Se puede observar que en el caso
b
= 1,5
la banda de fonones evoluciona hacia = 0, mientras que para
b
= 2,5
lo hace hacia = .
44 4. Estabilidad
(a) (b)
0.05 0.1 0.15 0.2 0.25 0.3 0.35
0
0.5
1
1.5
2
2.5
3
C
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
0
0.5
1
1.5
2
2.5
3
C
(c)
0 0.5 1 1.5
0
0.5
1
1.5
2
2.5
3
C
Fig. 4.7: Evoluci on con el acoplo de los argumentos de Floquet para un breather
(a) con potencial
3
y
b
= 0,8, (b) con potencial de sine-Gordon y
b
= 0,8 y (c) con potencial
4
duro onsite y de acoplo y
b
= 3. En
el primer caso, se observa la existencia de dos modos localizados que se
salen de la banda de fonones pero que posteriormente vuelven a entrar
en ella. En el segundo caso, no existen modos localizados. En el tercer
caso, existe un modo localizado que bifurca en = 0 produciendo una
inestabilidad.
(a) (b)
0.05 0.1 0.15 0.2 0.25 0.3 0.35
0
0.5
1
1.5
2
2.5
3
C
0.05 0.1 0.15 0.2 0.25 0.3 0.35
0.6
0.8
1
1.2
1.4
1.6
1.8
2
2.2
|
|
C
(c) (d)
0 0.02 0.04 0.06 0.08 0.1
0
0.5
1
1.5
2
2.5
3
C
0 0.02 0.04 0.06 0.08 0.1
0.8
0.9
1
1.1
1.2
1.3
|
|
C
del modulo de los autovalores (derecha) para dos 2sB con potencial de
Morse y acoplamiento positivo. Las guras superiores corresponden a
un 2sB en fase con frecuencia
b
= 0,8, mientras que las inferiores lo
hacen con un 2sB en anti-fase de frecuencia
b
= 0,75. Como se puede
comprobar, el 2sB en fase es inestable hasta C 0,13, momento en el
que entra en el crculo unidad un autovalor que estaba antes fuera de
el. Sin embargo, para C 0,24, el 2sB recupera su estabilidad. Como
vemos, este comportamiento es el opuesto al caso del 1sB con las mismas
caractersticas (Figura 4.4a). Este hecho, denominado intercambio de
estabilidad es muy importante a la hora de producir un breather m ovil.
El 2sB en anti-fase es estable hasta C 0,0885, momento en el que
experimenta una bifurcaci on oscilatoria. Los dos modos que colisionan
tienen su origen en un modo que se localiza de la banda de fonones y en
uno de los modos que estaban en = 0 en el lmite anti-continuo.
46 4. Estabilidad
(a) (b)
0.07 0.06 0.05 0.04 0.03 0.02 0.01 0
0
0.5
1
1.5
2
2.5
3
C
0.07 0.06 0.05 0.04 0.03 0.02 0.01 0
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
|
|
C
(c) (d)
0.1 0.08 0.06 0.04 0.02 0
0
0.5
1
1.5
2
2.5
3
C
0.1 0.08 0.06 0.04 0.02 0
0.5
1
1.5
2
2.5
3
3.5
|
|
C
m odulo de los autovalores (derecha) para dos 2sB con potencial de Morse
y acoplamiento negativo. Las guras superiores corresponden a un 2sB
en fase con frecuencia
b
= 0,8, mientras que las inferiores lo hacen con
un 2sB en anti-fase de frecuencia
b
= 0,75. El 2sB en fase es estable
hasta que un modo localizado, cuyo autovalor estaba inicialmente en
1, colisiona con la banda de fonones en C 0,028, produciendose
una bifurcaci on oscilatoria, recuperando su estabilidad en C 0,054.
El 2sB en anti-fase, en cambio, es siempre inestable.
(a) (b)
0 0.05 0.1 0.15 0.2
0
0.5
1
1.5
2
2.5
3
C
0 0.05 0.1 0.15 0.2
0.92
0.94
0.96
0.98
1
1.02
1.04
1.06
1.08
|
|
C
(c) (d)
1 0.8 0.6 0.4 0.2 0
0
0.5
1
1.5
2
2.5
3
C
1 0.8 0.6 0.4 0.2 0
0
200
400
600
800
1000
1200
|
|
C
del modulo de los autovalores (derecha) para dos 2sB con potencial
4
duro, acoplamiento positivo y frecuencia
b
= 1,5. Las guras superio-
res corresponden a un 2sB en fase, mientras que las inferiores lo hacen
con un 2sB en anti-fase. Como se puede comprobar, el 2sB en fase es
estable hasta que un modo localizado colisiona con la banda de fonones
en C 0,14, produciendose una bifurcaci on oscilatoria. En cambio, el
2sB en anti-fase es inestable para todo valor de C.
48 4. Estabilidad
(a) (b)
1 0.8 0.6 0.4 0.2 0
0
0.5
1
1.5
2
2.5
3
C
0.4 0.3 0.2 0.1 0
0
1
2
3
4
5
6
7
|
|
C
(c) (d)
1 0.8 0.6 0.4 0.2 0
0
0.5
1
1.5
2
2.5
3
C
0.4 0.3 0.2 0.1 0
0
1
2
3
4
5
6
7
|
|
C
m odulo de los autovalores (derecha) para dos 2sB con potencial
4
du-
ro, acoplamiento negativo y frecuencia
b
= 1,5. Las guras superiores
corresponden a un 2sB en fase, mientras que las inferiores lo hacen
con un 2sB en anti-fase. Como se puede comprobar, el 2sB en fase es
siempre inestable. En cambio, el 2sB en anti-fase es estable hasta que
para C 0,26 comienza una cascada de bifurcaciones arm onicas.
4.5. Teora de bandas de Aubry 49
4.5. Teora de bandas de Aubry
4.5.1. Introducci on
Hemos visto que el an alisis de Floquet nos permite conocer la estabilidad de un
breather e incluso predecir algunos fen omenos. Sin embargo, existen otros hechos
de los que no puede dar cuenta. Adem as, hay ciertas situaciones en los que no se
puede saber con exactitud si un breather es estable o no porque pueden aparecer
errores numericos cuando hay varios autovalores con valores similares. Por ello,
Aubry [19] introdujo un an alisis m as completo de la estabilidad, bas andose en el
hecho de que la ecuaci on din amica de la perturbaci on (4.2) es equivalente a una
red de osciladores parametricos. Estos osciladores estan descritos por ecuaciones
de Hill, las cuales son equivalentes a la ecuaci on de Schr odinger para un electr on
en un potencial peri odico, salvo que la periodicidad en la ecuaci on de Hill es
en el tiempo. Es bien sabido que las soluciones para la ecuaci on de Schr odinger
en el caso anteriormente mencionado tienen una estructura de bandas, la cual
aparecer a tambien en la ecuacion de Hill.
Como vimos en la secci on 4.2, la ecuaci on de autovalores del operador de New-
ton es:
N(u(t)) = E. (4.22)
El teorema de Bloch implica que los autovalores forman un conjunto de bandas
E
(), las cuales no se solapan. Adem as, E
() es una funci on continua, no cons-

tante y de periodo 2. Debido a que las autofunciones del operador de Newton se
pueden elegir como funciones de Bloch (ecuaci on (4.8)), el par ametro , tambien
llamado par ametro de fase, cumplir a la ecuaci on:
(, T) = e
j
(E)
(, 0), (4.23)
El par ametro se elige de forma que este en la primera zona de Brioullin, es
decir, (, ]. As, al conjunto de autovalores E
() tal que pertenezca a la

primera zona de Brioullin, se le denomina banda .
Las bandas estar an asociadas a soluciones estables si
R. Ademas, como
a una autofunci on (t) le corresponde otra
(t) con el mismo autovalor E y ar-

gumento de Floquet
(E), las bandas ser an simetricas respecto a
= 0 (ver
Apendice A).
Los valores (E) se obtienen diagonalizando la matriz F
E
, que se obtiene inte-
grando las ecuaciones (4.6):
_
(T)
(T)
_
= F
E
_
(0)
(0)
_
. (4.24)
La monodroma se corresponder a, l ogicamente, con F
o
. Esto quiere decir que
los argumentos de Floquet ser an
(0). Es decir, para que una soluci on sea estable,

50 4. Estabilidad
deben existir 2N bandas (contando la degeneraci on) que corten o sean tangentes
al eje E = 0.
Este analisis de bandas permite mejorar el criterio de la signatura de Krein para
prever si cuando dos autovalores coinciden en el crculo unidad se salen o no del
crculo unidad. Se sabe, por lo visto en secciones anteriores, que el criterio es una
condici on necesaria pero no suciente para que se produzca una colisi on. Aubry
[19] demuestra que la signatura de Krein es el opuesto del signo de la pendiente
de la banda que corresponde a cada autovalor en E = 0, es decir,
(
) = sgn
_
dE
()
d
_
. (4.25)
Como se muestra en [19], la condici on necesaria y suciente para que se produz-
ca la colisi on es que los autovalores pertenezcan a la misma banda, ya que en ese
caso las pendientes en los puntos de bifurcacion tendr an distinto signo. Este hecho
se ve ilustrado en la gura 4.12, en la cual se muestra de manera esquem atica la
forma de las bandas en los distintos escenarios de bifurcaciones.
La teora de bandas de Aubry se ha utilizado para entender c omo el tama no
nito de la red afecta a la estabilidad de las soluciones [141] y es la base para el
estudio de la estabilidad de multibreathers.
4.5.2. Estructura de bandas de un breather
Supongamos una red en la que se induce un multibreather con p sitios excitados.
En el lmite anti-continuo, la ecuaci on de autovalores del operador de Newton se
desacopla de forma que quedan p ecuaciones de Hill:
o
(t) + V
(u(t))
o
(t) = E
o
(t), (4.26)
que se corresponden con los osciladores excitados. Por otro lado, los osciladores en
reposo se describen mediante las N p ecuaciones:
n
(t) +
2
o
n
(t) = E
n
(t). (4.27)
Estas tienen por solucion:
n
(t) = e
j
2
o
E
n
(0), (4.28)
por lo que, teniendo en cuenta (4.23), el par ametro de fase se dene de la forma:
(E) =
2
b
_
2
o
E, (4.29)
de forma que este restringido a la primera zona de Brioullin. As, las bandas que
describen los osciladores en reposo vienen dadas por:
(a) (b)
E

(c) (d)
E
o
E
o
Fig. 4.12: Comportamiento de las bandas en diferentes tipos de bifurcaciones.
La evoluci on del sistema es tal que cuando cambian los par ametros,
el valor de E disminuye. (a) corresponde a una bifurcaci on arm onica
(se pierde la interseccion en = 0); (b) lo hace con una bifurcaci on
subarm onica (se pierde la intersecci on en = ); (c) y (d) muestran
dos casos en los que dos autovalores de la monodroma coinciden en
o
.
El primero se corresponde con una colisi on de Krein pues los modos
que coinciden pertenecen a la misma banda. El segundo representa un
cruce de los autovalores: no hay bifurcaci on pues los modos pertenecen
a distintas bandas.
52 4. Estabilidad
E
() =
2
o

2
b
_

2
_
2
. (4.30)
En el caso de los osciladores excitados, las bandas se obtienen integrando
numericamente las ecuaciones (4.26). Sin necesidad de realizar la integraci on se
puede deducir que existe una banda tangente al eje E = 0 debido a que u(t) es so-
lucion de (4.26), de forma que el punto de tangencia es (E, ) = (0, 0) (ver Apendice
A). Esta banda tiene curvatura positiva si el potencial es blando, y negativa si es
duro [50, 121].
La estructura de bandas est a acotada para E > 0 y no lo est a para E < 0.
Este hecho es valido tanto para osciladores en reposo como excitados.
Cuando el sistema se acopla, las bandas pierden su degeneraci on. De todas las
bandas correspondientes a los osciladores excitados, habr a una que permanezca
tangente al eje E = 0 debido a la existencia de un modo de fase global.
La gura 4.13 muestra, con un ejemplo concreto, la estructura de bandas tanto
en el lmite anti-continuo como para un valor determinado de la constante de
acoplo.
Puede ocurrir que cuando vare el acoplo, varas bandas se mezclen. Esto
puede tener lugar entre bandas correspondientes a fonones o entre una banda
correspondiente a un modo localizado y las bandas de los fonones. En ambos
casos, se pueden abrir huecos que originan burbujas de inestabilidad [137, 141].
4.5.3. Relaci on entre las bandas de Aubry y los autovalores del
jacobiano
Como vimos en la secci on 4.2, el operador de Newton coincide con el jacobiano
si las autofunciones son peri odicas reversibles y en el tiempo, es decir, si
E
2
s
(
b
). Sin embargo, las teoras de Floquet y de Aubry exigen simplemente que
tenga segunda derivada continua, es decir, C
2
. Esto indica que existen ciertas
condiciones en las que las autofunciones del operador de Floquet F
E
coinciden con
las autofunciones del jacobiano, ya que C
2
E
2
s
(
b
).
Cuando se dan estas condiciones? Pues bien, los autovectores del operador de
Newton, seg un el teorema de Bloch, pueden ser expresados por la ecuaci on (4.8),
es decir:
(t) = e
j
(E)t/T
v(t), (4.31)
donde v(t) es una funcion peri odica de periodo T y reversible en el tiempo. Por
tanto, si
(E) = 0, (t) sera peri odica y reversible en el tiempo, es decir,

E
2
s
(
b
). En otras palabras, los autovalores del jacobiano coinciden con los cortes
de las bandas con el eje = 0, siempre que estas correspondan a modos peri odicos
y reversibles en el tiempo. Debido a este ultimo hecho, el modo de fase y el de
crecimiento no aparecen en el espectro del jacobiano (ver Apendice A).
(a) (b)
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
1
1.5
2
E
E
o
E
1
E
2
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
1
E
(c) (d)
3 2 1 0 1 2 3
6
4
2
0
2
4
E
E
o
E
1
E
2
3 2 1 0 1 2 3
6
4
2
0
2
4
E
Fig. 4.13: Diagramas de bandas en el lmite anti-continuo (izquierda) y para un

acoplo C = 0,05 (derecha) correspondientes a un breather con potencial
de Morse y frecuencia
b
= 0,8 (arriba) y potencial
4
duro y frecuencia
b
= 1,5 (abajo). Los valores E
o
, E
1
y E
2
indican los cortes con el eje E
que se corresponden con autovalores del jacobiano (por claridad, s olo
se han especicado en el lmite anti-continuo). Observese que el corte
con el eje E para E = 0 no aparece en el jacobiano, tal como se predice
en el Apendice A).
54 4. Estabilidad
Este resultado ha sido expuesto por primera vez en [55]. Recordemos que el
jacobiano se puede calcular mediante las derivadas de las ecuaciones din amicas
respecto a los coecientes de Fourier de la soluci on, F
k
n
/z
k
n
(ecuacion (3.13)).
Veamos que tambien estos resultados coinciden con los esperados seg un la teora
estandar de bifurcaciones.
Para que ocurra una bifurcaci on se debe anular un autovalor del jacobiano.
Esto esta de acuerdo con la teora de bandas, pues seg un esta, se debe perder una
interseccion con el eje E = 0 para que haya un cambio de estabilidad. Seg un la
teora de bandas, las bifurcaciones oscilatorias debidas a colisiones de Krein y las
bifurcaciones subarm onicas se deben a perdidas de interseccion con el eje E = 0
para valores de = 0, mientras que las bifurcaciones arm onicas se deben a perdidas
de interseccion en = 0. Esto implica que, mediante el estudio del espectro del
jacobiano, s olo pueden detectarse las bifurcaciones arm onicas. Este hecho tambien
se conrma por los resultados de las teoras de bifurcaciones. En efecto, cuando
un autovalor del jacobiano se anula, es porque se deja de cumplir el teorema de la
funcion implcita, es decir, el sistema deja de tener solucion unica en el espacio de
trabajo. En el caso de una bifurcaci on arm onica, las soluciones emergentes siguen
perteneciendo a E
2
s
(
b
). Sin embargo, cuando la bifurcaci on es subarm onica, las
nuevas soluciones pertenecen a E
2
s
((2k + 1)
b
/2) E
2
s
(
b
), k N y cuando la
bifurcacion es oscilatoria, las soluciones dejan de ser peri odicas.
4.6. Teorema de Estabilidad de Multibreathers
4.6.1. Formulaci on del Teorema
Como se ha comentado anteriormente, existe un teorema, en virtud del cual, se
puede predecir la estabilidad o inestabilidad de un multibreather cuando se acopla
el sistema. Este teorema, el cual conjuga la teora de bandas de Aubry con la
teora de perturbaciones degeneradas, ha sido formulado recientemente por JFR
Archilla y se encuentra en la referencia [12], aplic andose a potenciales asimetricos
exclusivamente en la presente tesis. En esta secci on nos limitaremos a exponer el
teorema y varias aplicaciones del mismo, de modo que la demostraci on se encuentra
en el Apendice B.
El sistema considerado es una red homogenea de tipo KleinGordon con aco-
plamiento arm onico, aunque el teorema puede extenderse f acilmente a otro tipo
de acoplamientos. As, se pueden expresar las ecuaciones din amicas de la forma:
u
n
+ V
(u
n
) + C
N
m=1
J
nm
u
n
= 0 n = 1, . . . , N, (4.32)
donde C es la constante de acoplo y J es una matriz de acoplo, la cual es constante
e incluye las condiciones de contorno.
Recordemos que para un multibreather con p sitios excitados en el lmite anti-
continuo, existen p bandas degeneradas tangentes al eje E = 0 en el punto (0, 0).
4.6. Teorema de Estabilidad de Multibreathers 55
Pues bien, el Teorema de Estabilidad de Multibreathers permite conocer los auto-
valores E
m
del operador de Newton (4.6) cuando se rompe la degeneraci on debido
a la introducci on del acoplamiento. Estos autovalores vendr an dados por:
E
m
= C
m
, m = 0, . . . , p 1, (4.33)
siendo
m
los autovalores de la matriz de la perturbaci on Q, la cual se construye
siguiendo los siguientes pasos:
1. Se calcula el coeciente de simetra , el cual viene dado por:
=
k1
(1)
k
k
2
z
2
k
k1
k
2
z
2
k
> 0, (4.34)
donde z
k
son los coecientes de Fourier de la solucion de la ecuaci on din amica
de un oscilador aislado sometido al potencial onsite V . N otese que este
coeciente es igual a la unidad si V es simetrico. En el caso del potencial de
Morse, =
b
, tal como se demuestra en el Apendice B.
2. Se construye la matriz

J
de orden N cuyas componentes son iguales a:
nm
=
1
2
J
nm
[(
n
m
+ 1) + (
n
m
1)], (4.35)
donde
n
es la enesima componente del codigo del multibreather. Es decir,
nm
=
_
_
_
J
nm
si
n
=
m
J
nm
si
n
=
m
. (4.36)
3. Se considera la matriz

J, de orden p, que sera igual a

J
, salvo que se han

eliminado las las y columnas correspondientes a osciladores en reposo en el
lmite anti-continuo.
4. La matriz de la perturbaci on Q sera vendr a dada por:
Q
nm
=
_
_
_
J
nm
si n = m
m=n
J
nm
si n = m
. (4.37)
Debido a que la curvatura de las bandas consideradas es positiva si el potencial
es blando y negativa si es duro, para que el multibreather sea estable se necesita
que los autovalores E
m
sean negativos (positivos) si el potencial es blando (duro).
Un hecho a destacar es que siempre existe un autovalor E
m
nulo debido a la
existencia del modo de fase global. Sin embargo, si existe m as de un autovalor nulo,
la teora de perturbaciones de primer orden no se puede aplicar, estando indenida
la predicci on del TEM.
56 4. Estabilidad
4.6.2. Aplicaci on del Teorema
En este apartado, generalizamos para potenciales asimetricos la aplicaci on del
TEM para psite breathers expuesta en [12]. N otese la sencillez de la demostraci on
de las propiedades de estabilidad de 2site breathers, comparada con la que aparece
en [137].
2site breathers
Tal como se vio en la secci on 4.4.2, existen dos casos posibles. Estos vendr an
dados por el c odigo = {0, . . . , 0, 1,
1
, 0, . . . , 0}, donde
1
= 1 si los sitios ex-
citados del 2site breather vibran en fase y
1
= 1 si lo hacen en oposici on de
fase.
La matriz de la perturbaci on Q vendr a dada por:
Q =
1
2
[(1 +
1
) (1
1
)]
_
1 1
1 1
_
. (4.38)
Los autovalores de esta matriz ser an
0
= 0 y
1
=
_
_
_
2 si
1
= 1
2 si
1
= 1
. (4.39)
Esto implica que el signo del autovalor E no nulo ser a opuesto al de la constante
de acoplo C si el 2site breather est a en fase, mientras que ambos signos ser an
iguales si el 2site breather est a en anti-fase. Estos resultados estan de acuerdo
con los obtenidos en la Secci on 4.4.2.
3site breathers
En este caso, consideraremos los multibreathers dados por el codigo =
{0, . . . , 0,
1
, 1,
3
, 0, . . . , 0}. Habr a entonces tres casos diferentes:
1
=
3
= 1 y
3
= 1.
Si la interacci on es de primeros vecinos, la matriz Q es de la forma:
Q =
_
_
1

1
0
1

1
+
3

3
0
3

3
_
_
, (4.40)
donde:
n
=
1
2
[(1 +
n
) (1
n
)]. (4.41)
As, la ecuaci on caracterstica de la matriz Q es:
3
2(
1
+
3
)
2
+ 3
1
3
= 0. (4.42)
Como cabra esperar, un autovalor (
0
) es cero. Sin necesidad de resolver la
ecuacion caracterstica se puede saber cu ando los dos autovalores restantes tienen
distinto signo, es decir,
1
2
< 0. Este hecho implicara que el multibreather sera
inestable, independientemente de la forma del potencial y del signo de la constante
de acoplo, pues una banda subira y otra bajara. El valor de este producto es:
2
= 3
1
3
=
_
_
3 si
1
=
2
= 1
3
2
si
1
=
2
= 1
3 si
1
2
= 1
. (4.43)
Como podemos apreciar, siempre que
1
y
3
tengan distinto signo, el mul-
tibreather ser a inestable. En caso de que
1
y
3
fueran iguales, la estabilidad
dependera del tipo de potencial y del valor de la constante de acoplo.
El valor de los autovalores
1
y
2
se encuentra en la siguiente tabla:
Caso
1

2
1
=
3
= 1 1 3
1
=
3
= 1 3
3
= 1 1
_
2
+ + 1 1 +
_
2
+ + 1
Estos resultados implican los expuestos en la Seccion 4.4.2. Las guras 4.14 y
4.15 muestran, para 3site breathers, los autovalores de la matriz Q con potenciales
4
blando y Morse, respectivamente. En esas guras puede observarse tambien el
rango de validez de la teora de perturbaciones.
Multibreathers con un numero arbitrario de sitios
Supongamos un multibreather con p sitios adyacentes excitados (es decir, un p
site breather) y con interacci on de primeros vecinos. El TEM permite determinar
facilmente la estabilidad de aquellas soluciones en las que los sitios excitados vibran
en fase o en zigzag.
Si todos los osciladores est an en fase, la matriz Q = Q
f
, de orden p, sera de la
forma:
Q
f
=
_
_
1 1 0 . . . 0 0
1 2 1 0 . . . 0
. . . . . . . . . . . . . . . . . .
0 0 . . . 1 2 1
0 0 . . . 0 1 1
_
_
. (4.44)
58 4. Estabilidad
(a) (b)
0.02 0.04 0.06 0.08 0.1
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
E
C
0.02 0.04 0.06 0.08 0.1
0.25
0.2
0.15
0.1
0.05
0
E
C
(c)
0.02 0.04 0.06 0.08 0.1 0.12
0.2
0.1
0
0.1
0.2
0.3
E
C
Fig. 4.14: Autovalores de la matriz de perturbaci on para distintos 3site brea-
thers con potencial
4
. (a) corresponde al c odigo {1, 1, 1}, (b) a
{1, 1, 1} y (c) a {1, 1, 1}. Los puntos representan los valores cal-
culados mediante la teora de bandas, mientras que las lneas corres-
ponden a las predicciones del TEM.
(a) (b)
0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
0
0.05
0.1
0.15
0.2
E
C
0.01 0.02 0.03 0.04 0.05 0.06 0.07
0.2
0.15
0.1
0.05
0
E
C
(c)
0.01 0.02 0.03 0.04 0.05 0.06 0.07
0.2
0.15
0.1
0.05
0
0.05
0.1
0.15
0.2
E
C
Fig. 4.15: Autovalores de la matriz de perturbaci on para distintos 3site brea-
thers con potencial de Morse. (a) corresponde al c odigo {1, 1, 1}, (b) a
{1, 1, 1} y (c) a {1, 1, 1}. Los puntos representan los valores cal-
culados mediante la teora de bandas, mientras que las lneas corres-
ponden a las predicciones del TEM.
60 4. Estabilidad
La diagonalizaci on de esta matriz no es trivial a simple vista. Sin embargo, si
se escribe la ecuacion caracterstica Q x = x en funcion de las componentes de
los autovectores x = {x
n
},
_
_
x
1
x
2
= x
1
2x
n
x
n+1
x
n1
= x
n
, n = 1, p
x
p
x
p1
= x
p
, (4.45)
se pueden calcular los autovalores partiendo del hecho que que la ecuaci on central
es similar a la ecuaci on de dispersi on de los fonones en una red FPU. As, teniendo
en cuenta (2.8) se llega a que los autovalores son de la forma:
f
m
= 4 sen
2
m
2p
m = 0, . . . , p 1, (4.46)
es decir, son positivos, coincidiendo el resultado con el obtenido para 2site y 3site
breathers con c odigos = {. . . , 1, 1, . . . } y = {. . . , 1, 1, 1, . . . }, respectivamente.
En el caso de que los sitios excitados vibraran en zigzag, la matriz de la per-
turbaci on sera Q
z
= Q
f
, por lo que
z
m
=
f
m
, es decir:
z
m
= 4 sen
2
m
2p
m = 0, . . . , p 1. (4.47)
En consecuencia, los autovalores seran negativos. En consecuencia, se reobtie-
nen los resultados para los 2site y 3site breathers con c odigos = {. . . , 1, 1, . . . }
y = {. . . , 1, 1, 1, . . . }, respectivamente.
En el caso de que no se siguiera ninguno de los dos patrones anteriores, se
puede ver numericamente que algunos autovalores son positivos y otros negativos,
lo cual implica que el multibreather es inestable con independencia de la forma del
potencial onsite y del signo de la constante de acoplo.
En caso de que existiera en el c odigo un cero entre regiones con valores no nulos,
la matriz Q se dividira en bloques. Si quedara en el c odigo alg un valor distinto de
cero entre dos ceros, no sera v alida la teora de perturbaciones de primer orden y
el TEM no podra predecir la estabilidad o inestabilidad.
4.7. Estabilidad en sistemas disipativos
En esta seccion se exponen brevemente ciertos aspectos relativos a la estabilidad
de sistemas disipativos. La evoluci on de los autovalores del operador de Floquet
con el acoplamiento es bastante m as complicada que en el caso de los sistemas
hamiltonianos y su estudio est a fuera del alcance de esta tesis. Lo que se exponen
en este apartado son resultados b asicos que son de utilidad en el desarrollo posterior
de la tesis, los cuales est an basados en el trabajo de Marn et al [145]. Supongamos
4.7. Estabilidad en sistemas disipativos 61
una red de osciladores forzada y amortiguada, tal como la describimos en el captulo
3 (ecuaci on (3.23)):
u
n
+ u
n
+ V
(u
n
) + C(2u
n
u
n+1
u
n1
) = F sin
b
t. (4.48)
La ecuaci on de la perturbaci on ser a, por tanto:
n
+

n
+ V
(u
n
)
n
+ C(2
n
n+1
n1
) = 0. (4.49)
Como podemos observar, esta ecuaci on es independiente del valor de la fuerza
externa F. Sin embargo, los resultados sobre la estabilidad variar an dependiendo
de que el sistema este o no forzado.
La ecuaci on (4.49) se puede transformar en otra equivalente a la ecuaci on de una
perturbaci on en sistemas hamiltonianos (4.2) si hacemos un cambio de variables
de la forma [192]:
n
(t) = e
t/2
n
(t). (4.50)
Con este cambio, la ecuaci on (4.49) se transforma en:

n
+ (V
(u
n
)
2
/4)
n
+ C(2
n
n+1
n1
) = 0. (4.51)
Esta ecuacion tiene las propiedades vistas para sistemas hamiltonianos. Sin
embargo, las autofunciones del operador de Floquet de las ecuaciones originales
tendr an otras propiedades pues hay que deshacer el cambio (4.50). As, los auto-
valores deber an multiplicarse por el factor e
T/2
.
En consecuencia, cuando el acoplamiento es peque no, los multiplicadores de
Floquet, en lugar de estar en la circunferencia unidad, se situar an en otra concentri-
ca de radio r < 1. Este radio es, obviamente, r = e
T/2
. Ademas, si es autovalor,
tambien lo seran r/,
y r/
.
El hecho de que un autovalor abandone el crculo de radio r no implica ines-
tabilidad. El sistema ser a inestable cuando alguno de los autovalores se salga del
crculo unidad.
En el caso de un sistema forzado, no existen los modos de fase y crecimiento.
Esto se puede observar derivando la ecuaci on din amica (3.23). Fsicamente, este
hecho se debe a que tanto la frecuencia del breather como el origen de tiempos est an
jados por la fuerza externa. En el caso de que el sistema no estuviera forzado, los
modos de fase y crecimiento seguiran existiendo.
Respecto al diagrama de bandas, la ecuaci on es equivalente al caso hamiltonia-
no. No obstante, si el sistema est a forzado, la banda tangente al eje E = 0 deja de
existir, pues el modo de fase u(t) deja de ser soluci on de la ecuaci on (4.49).
En el lmite anti-continuo, pueden calcularse las bandas de los osciladores en
reposo, que vendr an dadas por:
62 4. Estabilidad
(a) (b)
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
1
E
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
E
Fig. 4.16: Diagramas de bandas en el lmite anti-continuo (a) y para un aco-

plo C = 0,05 (b) para un breather disipativo con un potencial de
sine-Gordon. La constante de amortiguamiento es = 0,02, la fuerza
externa tiene amplitud F = 0,02 y la frecuencia es
b
= 0,2. Como
se puede observar, no existe banda tangente al eje E = 0. Se nalar que
los cortes con el eje = 0 aparecen en el jacobiano ya el espacio de
trabajo no esta restringido a soluciones reversibles en el tiempo.
E
() =
2
o

2
4

2
b
_

2
_
2
. (4.52)
La gura 4.16 muestra un ejemplo de la estructura de bandas de un breather
con disipaci on y fuerza externa.
5. Estabilidad de dark breathers
5.1. Introducci on
Hasta ahora, hemos considerados breathers en los que la localizaci on consiste
en unos pocas partculas vibrando con gran amplitud mientras que las restantes
lo hacen con poca amplitud. Puede darse tambien el caso contrario, es decir, la
mayora de los osciladores de la red vibra con gran amplitud, excepto unos po-
cos que lo hacen con peque na amplitud. A este tipo de multibreather lo podemos
denominar dark breather o breather oscuro [8]. Este nombre tiene su origen en la
analoga con los solitones oscuros que aparecen en sistemas continuos [195, 196]
y discretos gobernados por la ecuaci on DNLS [110, 118, 120]. El estudio de dark
breathers en redes KleinGordon se ha abordado dos maneras. La primera, llevada
a cabo por Morgante et al [150, 151], los considera como desestabilizaciones de on-
das estacionarias no lineales. El segundo punto de vista [8] estudia las propiedades
de estabilidad de los dark breathers a partir del an alisis de bandas (Captulo 4)
aplicado a multibreathers.
En este captulo vamos a centrarnos en el an alisis de la estabilidad de dark
breathers bas andonos en la teora de bandas de Aubry y en el Teorema de Estabi-
lidad de Multibreathers. Los sistemas que se consideran son tanto hamiltonianos
como disipativos.
Por ultimo, quisieramos se nalar que a los breathers en los que la mayora de
las partculas vibran con peque na amplitud reciben en este captulo el nombre de
bright breathers (breathers brillantes) en analoga con la terminologa empleada en
solitones.
5.2. Estabilidad de dark breathers en sistemas
hamiltonianos
5.2.1. Presentaci on de las ecuaciones
En esta seccion se estudia numericamente la estabilidad de los dark breathers en
sistemas hamiltonianos con acoplamiento arm onico y atractivo, cuyas ecuaciones
din amicas vienen dadas por:
u
n
+ V
(u
n
) + C(2u
n
u
n+1
u
n1
) = 0, (5.1)
64 5. Estabilidad de dark breathers
con C > 0 y condiciones de contorno peri odicas. Tambien se consideran siste-
mas con acoplamiento repulsivo. Para darle un signicado fsico a este tipo de
acoplamiento, vamos a suponer que la interacci on es del tipo que aparece en el
estudio del modelo con acoplamiento dipolar del ADN de los Captulos 11 y 12,
pero restringida al primer vecino. As, la ecuaciones din amicas vienen dadas por:
u
n
+ V
(u
n
) + C(u
n+1
+ u
n1
) = 0, (5.2)
con C > 0, igualmente. El perl de varios dark breathers se encuentra en la gura
5.1.
Las ecuaciones linealizadas para el estudio de la estabilidad, en el caso atractivo,
son:
n
+ V
(u
n
)
n
+ C(2
n
n+1
n1
) = 0, (5.3)
mientras en el caso repulsivo son:
n
+ V
(u
n
)
n
+ C(
n+1
+
n1
) = 0. (5.4)
Se nalemos que las soluciones que vamos a estudiar en esta secci on vendr an
dadas por el c odigo {. . . 111101111 . . .}, es decir, 1site dark breathers en fase. Al
conjunto de osciladores excitados se le conoce tambien con el nombre de fondo.
5.2.2. An alisis de estabilidad
An alogamente al caso de un bright breather, el operador de Floquet tiene p
parejas de autovalores correspondientes a los osciladores excitados (fondo) cuyo
valor es +1, mientras que N p corresponden a los osciladores en reposo. La dife-
rencia con el caso de un dark breather es que en este la mayora de los osciladores
estan excitados. Como se nalamos anteriormente, vamos a considerar 1site dark
breathers, por lo que p = N 1.
Este hecho se traslada al an alisis de bandas de forma que existen N 1 dege-
neradas y tangentes al eje E = 0 en el punto (, E) = (0, 0) correspondientes al
fondo, mientras que existe otra banda que corta al eje y se corresponde al oscilador
en reposo. Es decir, el esquema es el mismo que el de las guras 4.13a y 4.13c.
La curvatura de las bandas del fondo es positiva si el potencial es blando, siendo
negativa si el potencial es duro.
Cuando se introduce el acoplo, se rompe la degeneraci on de las bandas corres-
pondientes al fondo, de forma que N 2 pierden la tangencia con el eje, y otra, la
que se corresponde con el modo de fase, permanece ja. Como se observa numeri-
camente, la perdida de tangencia de las bandas se produce porque estas suben
cuando el acoplo es atractivo o bajan cuando este es repulsivo. As, cuando se
introduce el acoplo, se pueden dar los siguientes casos:
Si el potencial es blando, el dark breather es estable para acoplo repulsivo e
inestable para acoplo atractivo.
5.2. Estabilidad de dark breathers en sistemas hamiltonianos 65
(a) (b)
8 6 4 2 0 2 4 6 8
0
0.2
0.4
0.6
0.8
1
A
m
p
l
i
t
u
d
Sitio
8 6 4 2 0 2 4 6 8
0.2
0
0.2
0.4
0.6
0.8
1
1.2
A
m
p
l
i
t
u
d
Sitio
(c) (d)
8 6 4 2 0 2 4 6 8
0.2
0
0.2
0.4
0.6
0.8
1
1.2
A
m
p
l
i
t
u
d
Sitio
8 6 4 2 0 2 4 6 8
0
0.2
0.4
0.6
0.8
1
A
m
p
l
i
t
u
d
Sitio
Fig. 5.1: Perl de dark breathers. (a) y (b) corresponden a breathers con potencial
de Morse onsite y acoplo atractivo y repulsivo, respectivamente (
b
=
0,8 y C = 0,01). (c) y (d) son an alogas a las guras anteriores, salvo que
el potencial onsite es
4
duro (
b
= 1,3 y C = 0,01).
(a)
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
1
1.5
2
E
(b) (c)
3 2 1 0 1 2 3
0
0.05
0.1
0.15
0.2
E
3 2 1 0 1 2 3
0.1
0.05
0
0.05
0.1
0.15
E
Fig. 5.2: Diagramas de bandas en el lmite anti-continuo (arriba) y para un aco-

plo C = 0,01 (abajo) para un dark breather con potencial de Morse y
frecuencia
b
= 0,8. (b) corresponde a un potencial de interacci on atrac-
tivo, mientras que (c) lo hace con un potencial de interacci on repulsivo.
Si el potencial es duro, el dark breather es estable para acoplo atractivo e
inestable para acoplo repulsivo.
Otro fen omeno que aparece cuando se introduce el acoplo es la hibridaci on de
las bandas correspondientes al fondo y al oscilador en reposo, con la consiguiente
aparici on de huecos. Esto hace que, al aumentar el acoplo, se produzcan colisiones
de Krein.
Ademas de estas bifurcaciones, existen otras de tipo subarm onico cuando el
potencial onsite no es simetrico. Esto se debe a que la banda naliza en = , por
lo que, cuando estas bajan lo suciente debido al acoplo, perder an su interseccion
en ese punto, y aparecen bifurcaciones. En el caso de potenciales simetricos, las
soluciones tambien tienen periodo T/2, por lo que las bandas no nalizan en ,
sino en 2. Por esta raz on no aparecen este tipo de bifurcaciones en el caso de
potenciales onsite simetricos.
Estos fenomenos se ilustran en las guras 5.2 y 5.3.
5.2. Estabilidad de dark breathers en sistemas hamiltonianos 67
(a)
3 2 1 0 1 2 3
6
4
2
0
2
4
E
(b) (c)
3 2 1 0 1 2 3
0.08
0.06
0.04
0.02
0
0.02
0.04
0.06
E
3 2 1 0 1 2 3
0.5
0.4
0.3
0.2
0.1
0
E
Fig. 5.3: Diagramas de bandas en el lmite anti-continuo (arriba), y para un aco-

plo C = 0,01 (abajo), para un dark breather con potencial
4
duro
y frecuencia
b
= 1,3. (b) corresponde a un potencial de interacci on
atractivo, mientras que (c) lo hace con un potencial repulsivo.
5.3. Aplicaci on del Teorema de Estabilidad de
Multibreathers
En este apartado veremos c omo los resultados anteriores pueden predecirse me-
diante la aplicaci on del Teorema de Estabilidad de Multibreathers expuesto en la
Seccion 4.6. Para ello, aplicaremos el TEM a phonobreathers, cuyo conocimien-
to permite establecer trivialmente los resultados para dark breathers. Adem as,
predeciremos la estabilidad de otros tipos de dark breathers los cuales no fueron
considerados en el apartado anterior.
5.3.1. Phonobreathers
Condiciones de contorno de extremos libres
Recordemos que un phonobreather es un multibreather con todos sus sitios
excitados. Tal como se vio en la Secci on 4.6, para que un nsite breather sea estable
es necesario que sus sitios excitados vibren en fase o en zigzag. Un phonobreather
es un caso particular de nsite breather en el que n = N, siendo N el n umero
de partculas de la red. Esto es, si un phonobreather no deriva del fon on con
q = 0 (vibraci on en fase) o q = (vibraci on en zigzag), ser a siempre inestable,
con independencia del signo de la constante de acoplo y de la forma funcional del
potencial onsite.
La matriz de la perturbaci on en el caso de un phonobreather con condiciones de
contorno de extremos libres es la misma que un Nsite breather. As, los autovalores
son positivos para q = 0 y negativos para q = . Las consecuencias de este hecho
se resumen en la siguiente tabla:
Interaccion atractiva Interaccion repulsiva
q Potencial blando Potencial duro Potencial blando Potencial duro
0 Inestable Estable Estable Inestable
Estable Inestable Inestable Estable
Existe un total acuerdo entre estos resultados y los obtenidos como consecuen-
cia del estudio de la inestabilidad modulacional [60] y de la ecuaci on DNLS [151].
Si las condiciones de contorno son de extremos jos, la matriz de la perturbaci on
es la misma que si los extremos estuvieran libres, de forma que los resultados son
identicos.
Condiciones de contorno peri odicas. Inestabilidades de paridad
Si las condiciones de contorno son peri odicas, la matriz Q es de la forma:
5.3. Aplicaci on del Teorema de Estabilidad de Multibreathers 69
Q
f
=
_
_
2 1 0 . . . 0 1
1 2 1 0 . . . 0
. . . . . . . . . . . . . . . . . .
0 0 . . . 1 2 1
1 0 . . . 0 1 2
_
_
, (5.5)
siempre que el phonobreather vibre en fase. En este caso, la ecuaci on de autovalores
se puede escribir como:
_
_
2x
1
x
2
x
N
= x
1
2x
n
x
n+1
x
n1
= x
N
, n = 1, N
2x
N
x
N1
x
1
= x
N
. (5.6)
Esta expresi on es identica a la ecuaci on de dispersi on para los fonones en una
cadena FPU con condiciones de contorno peri odicas, cuya soluci on es:
f
m
= 4 sen
2
m
N
, m = 0, . . . , N 1. (5.7)
En este caso, todos los autovalores son positivos, tal como ocurre con el pho-
nobreather con condiciones de contorno de extremos libres.
En el caso de que el phonobreather vibre en zigzag, aparecen nuevos fen omenos.
La matriz Q tiene la forma:
Q
z
=
_
_
( 1) 1 0 . . . 0
1 2 1 0 . . . 0
. . . . . . . . . . . . . . . . . .
0 0 . . . 1 2 1
0 . . . 0 1 ( 1)
_
_
, (5.8)
con =
1
2
[(
1
N
+ 1) + (
1
N
1)].
Si el n umero de partculas (N) es par, el codigo del phonobreather es de la forma
= {1, 1, 1, ..., 1, 1, 1}, por lo que = 1 y Q
z
= Q
f
, obteniendose que
z
m
=
f
m
.
Sin embargo, si N es impar, = {1, 1, 1, ..., 1, 1, 1}, por lo que = 1/
y Q
z
= Q
o
. Se encuentra que hay un autovalor positivo
+
(cuyo valor es
4/( + 2) si N es lo sucientemente grande, tal como se demuestra a continua-
cion), mientras que el resto son negativos. Esto implica que el phonobreather es
siempre inestable si el n umero de partculas del sistema es impar, produciendose
el fenomeno que denominamos inestabilidad de paridad. El origen fsico de estas
inestabilidades radica en el hecho de que el patr on de vibraci on del sistema se
rompe al unir las partculas de los extremos.
Para calcular el valor del autovalor
+
responsable de la inestabilidad de pari-
dad, renumeramos las partculas del sistema de forma que la central sea n = 0 y
las partculas en los extremos derecho e izquierdo sean, respectivamente, n = N y
n = N. As, la ecuaci on de autovalores se puede escribir de la forma:
_
_
2x
n
x
n1
x
n+1
= (
+
/) x
1
n = 1, N
(1 )x
N
+ x
N1
x
N
=
+
x
N
(1 )x
N
+ x
1N
x
N
=
+
x
N
. (5.9)
Como se ha explicado anteriormente, la inestabilidad que aparece es debida a
la rotura del patr on de vibraci on. Si se observa la forma del autovector correspon-
diente al autovalor
+
(ver gura 5.4), se puede ver que es similar al producido por
un modo de impureza centrado en los extremos de una cadena de osciladores (ver
Captulo 9). As, se puede suponer que las soluciones del sistema (5.9) cumplen las
propiedades siguientes:
x
n
= x
n
, (5.10)
x
n
= r x
n1
si n >> 0; x
n
= r x
n+1
si n << 0. (5.11)
A partir de la segunda propiedad, se obtiene que las ecuaciones segunda y
tercera del sistema (5.9) son equivalentes. En consecuencia, si se introducen las
ecuaciones (5.10) y (5.11) en (5.9), se obtiene:
_
_
_
2 r 1/r =
+
/
+
= 2 + (1/r 1)
. (5.12)
Resolviendo el sistema anterior se llega al resultado:
_
+
=
4
+ 2
r = 1 +
2
. (5.13)
Notese que estos valores son v alidos siempre que el n umero de partculas sea
lo sucientemente grande (aproximadamente, mayor que 9). Sin embargo, aunque
el n umero de partculas fuera menor, el signo de
+
no cambiara.
5.3.2. Dark breathers
En este caso, los sitios en reposo se corresponden con ceros en el codigo. En
la aplicaci on del TEM aparecen nuevos fenomenos dependiendo de las condiciones
de contorno del sistema.
Utilizaremos una notaci on que ayude a presentar los resultados. As, a la matriz
de la perturbaci on de un phonobreather la denominaremos P
v
b
(n), donde n es el
5.3. Aplicaci on del Teorema de Estabilidad de Multibreathers 71
8 6 4 2 0 2 4 6 8
0.8
0.6
0.4
0.2
0
0.2
0.4
0.6
0.8
x
n
n
Fig. 5.4: Autovector de la matriz de perturbaci on Q dada en la ecuaci on (5.8)
para un sistema con un n umero impar de partculas. Este autovector es
responsable de la inestabilidad de paridad.
n umero de partculas del phonobreather, b se corresponde con las condiciones de
contorno del sistema (P si son peri odicas, y L si son de extremos libres) y v indica
el patr on de vibraci on del phonobreather, es decir, v f si deriva del fon on con
q = 0 (vibraci on en fase) y v z si lo hace del fonon con q = (vibraci on en
zigzag). En el caso de dark breathers, utilizaremos la notaci on Q
v
b
(n
), donde n
es el n umero de ceros adyacentes en el c odigo. Con estos preliminares, se puede

proceder al estudio de la estabilidad de dark breathers.
Condiciones de contorno peri odicas
Supongamos un psite dark breather que deriva del fon on con q = 0. En este
caso, la matriz de la perturbaci on es de la forma:
Q
f
P
(p) =
_
_
2 1 0 0 . . . . . . 0 1
1 2 1 0 . . . . . . . . . 0
. . . . . . . . . . . . . . . . . . . . . . . .
0 0 1 1 0 0 . . . 0
0 . . . 0 0 1 1 0 0
. . . . . . . . . . . . . . . . . . . . . . . .
0 . . . . . . . . . 0 1 2 1
1 0 . . . . . . . . . 0 1 2
_
_
. (5.14)
Se puede ver que, utilizando las propiedades de simetra de determinantes bajo
intercambio de las y columnas, el polinomio caracterstico de Q
f
P
(p) es igual que
el de P
f
L
(N p). Es decir, los autovalores de Q
f
P
(p) son de la forma:
m
= 4 sen
2
m
2N
, m = 0, . . . , N p 1, (5.15)
es decir, son todos positivos con independencia de p. Estos resultados est an de
acuerdo con los obtenidos numericamente en la Seccion 5.2.2. Si el dark breather
deriva del fon on con q = , tendremos dos casos. Si
1
=
N
, se obtiene que
Q
z
P
(p) = Q
f
P
(p), de forma que:
m
= 4 sen
2
m
2N
, m = 0, . . . , N p 1. (5.16)
En consecuencia, todos los autovalores son negativos. En cambio, si
1
=
N
,
tenemos un caso similar al de un phonobreather con un n umero de partculas
impar, apareciendo de nuevo las inestabilidades de paridad.
Condiciones de contorno de extremos libres
En el caso de un dark breather que derive de un fon on con q = 0, la matriz de
la perturbaci on es de la forma:
Q
f
L
(p) =
_
_
1 1 0 0 . . . . . . 0 0
1 2 1 0 . . . . . . 0 0
. . . . . . . . . . . . . . . . . . . . . . . .
0 0 1 1 0 0 . . . 0
0 . . . 0 0 1 1 0 0
. . . . . . . . . . . . . . . . . . . . . . . .
0 0 . . . . . . 0 1 2 1
0 0 . . . . . . 0 0 1 1
_
_
, (5.17)
es decir, la matriz Q es diagonal por bloques. Esta matriz puede escribirse utili-
zando una notaci on m as abreviada:
Q
f
L
(p) =
_
_
P
f
L
(n
1
) 0
0 P
f
L
(n
2
)
_
_
, (5.18)
donde n
1
y n
2
son el n umero de valores del c odigo distintos de cero, a un lado y
al otro de los ceros correspondientes al dark breather. Por tanto,
spec Q
f
L
(p) = spec P
f
L
(n
1
) spec P
f
L
(n
2
). (5.19)
Como, por otro lado,
0 spec P
f
L
(n
1
) spec P
f
L
(n
2
), (5.20)
resulta que el autovalor cero es doblemente degenerado en el espectro de Q
f
L
(p),
por lo que el TEM no puede predecir la estabilidad del dark breather en este caso.
5.4. Efectos de tama no nito en la estabilidad de dark breathers 73
El origen fsico de este hecho se debe a la existencia de dos modos de fase globales.
En concreto, el dark breather divide a la matriz Q en dos bloques, de forma que la
ecuacion de autovalores correspondera a una red de osciladores formada por dos
fragmentos desacoplados. As, existen dos modos de fase: uno de ellos corresponde a
todas las componentes del primer fragmento iguales (y distintas de cero), mientras
que las del segundo fragmento seran todas iguales a cero; el otro modo de fase lo
hara con el caso complementario. Este fen omeno no ocurre cuando las condiciones
de contorno son peri odicas, pues en este caso, los dos fragmentos se encuentran
acoplados por los extremos, obligando de esta forma a que todas las componentes
del modo de fase global sean iguales (y distintas de cero).
Sin embargo, el hecho de que no se rompa la degeneraci on del autovalor ce-
ro no impide que se pueda predecir la inestabilidad de ciertas soluciones.

Estas
corresponden a los casos en los que alguno de los fragmentos son inestables por
s solos. Es decir, el TEM predice que un dark breather con potencial blando y
acoplo atractivo o potencial duro y acoplo repulsivo es inestable.
Si el dark breather deriva del fon on q = , la matriz de la perturbaci on es
Q
z
L
(p) = Q
f
L
(p), de forma que se vuelve a tener el autovalor cero degenerado.
Se nalemos, nalmente, que al igual que en el caso de phonobreathers, los re-
sultados para extremos jos son los mismos que para extremos libres.
En resumen, el Teorema de Estabilidad de Multibreathers no siempre puede
predecir la estabilidad de un dark breather con condiciones de contorno de extremos
libres o jos.
5.4. Efectos de tama no nito en la estabilidad de dark
breathers
Tal como muestran Marn y Aubry para bright breathersn [141], existen dos
tipos de bifurcaciones en las que el tama no de la red puede tener un papel impor-
tante. Estas son las bifurcaciones oscilatorias (debidas a colisiones de Krein) y las
subarm onicas, ya que ambos tipos de bifurcaciones implican a modos extendidos,
los cuales se corresponden con oscilaciones de peque na amplitud. En concreto, el
primer tipo implica la colisi on de un modo localizado con otro extendido, mientras
que en el segundo, ambos modos son extendidos. El estudio arroja los siguientes
resultados:
Las bifurcaciones subarm onicas desaparecen cuando el sistema es innito.
Las bifurcaciones oscilatorias no desaparecen cuando el sistema es innito.
Ademas, ambos tipos de bifurcaciones se producen por hibridaci on de las ban-
das, lo que se traduce en la apertura de huecos y en la formaci on de burbujas de
inestabilidad.
En el caso de dark breathers el escenario es algo distinto. En primer lugar, las
bifurcaciones subarm onicas se deben a la perdida de interseccion en = de una
banda correspondiente al fondo. Esta banda est a incluida totalmente en la primera
zona de Brioullin, mientras que en el caso del bright breather, la banda ocupa todas
las zonas de Brioullin, por lo que las unicas perdidas de interseccion se deben
a hibridaciones. Es por ello que, en el caso de dark breathers, las bifurcaciones
subarm onicas no desaparecen cuando el tama no del sistema es innito, ya que no
proceden de la hibridaci on de dos bandas distintas.
El caso de las bifurcaciones de Krein es similar a lo que ocurre con bright brea-
thers, es decir, se producen hibridaciones de bandas que no desaparecen cuando
aumenta el tama no del sistema. Es mas, cuando este es peque no, existen bur-
bujas de inestabilidad que se transforman en una inestabilidad persistente cuando
aumenta el n umero de partculas del sistema. Esto se puede comprobar en la gura
5.5.
En resumen, para un dark breather:
Las bifurcaciones subarm onicas no desaparecen cuando el sistema es innito.
Las bifurcaciones oscilatorias no desaparecen cuando el sistema es innito.
5.5. Estabilidad de dark breathers en sistemas
disipativos
En el estudio de la estabilidad de dark breathers en sistemas disipativos apa-
recen ciertas diferencias sustanciales con respecto al caso hamiltoniano.
Para realizar este estudio partimos del modelo de FrenkelKontorova:
u
n
+ u
n
+ V
(u
n
) + C(2u
n
u
n+1
u
n1
) = F sin
b
t. (5.21)
La estructura de bandas de este sistema (ver Secci on 4.7) muestra que no
existen bandas tangentes al eje E = 0 en el lmite anti-continuo (gura 4.16a).
Por ello, cuando se introduce el acoplamiento, el dark breather ser a estable tanto
si la interacci on es atractiva como si es repulsiva. Por otro lado, se observa que
no existe hibridaci on entre las bandas correspondientes a los modos extendidos,
por lo que las unicas bifurcaciones que pueden aparecer son de tipo arm onico,
las cuales no dependen del tama no del sistema (Esta armaci on es valida en el
modelo considerado por nosotros, cuyo potencial onsite es simetrico; el caso de
un potencial asimetrico, en el que podran aparecer bifurcaciones subarm onicas,
no ha sido estudiado).
Estos hechos se ilustran en la gura 5.6.
5.5. Estabilidad de dark breathers en sistemas disipativos 75
(a) (b)
1.2 1.3 1.4 1.5 1.6 1.7 1.8
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
E
1.4 1.5 1.6 1.7 1.8

0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
E
(c) (d)
0.94 0.96 0.98 1 1.02 1.04 1.06
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
E
||
0.94 0.96 0.98 1 1.02 1.04 1.06
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
E
||
Fig. 5.5: Las guras superiores muestran una ampliaci on de la zona de hibrida-
cion de las bandas correspondientes al oscilador en reposo y al fondo,
mientras que las inferiores muestran la dependencia (E, ||), siendo el
autovalor del operador de Newton correspondiente al modo que bifurca
en la colision de Krein. Las guras de la izquierda corresponden a un
sistema de 21 partculas, mientras que el de la derecha lo hace con uno
de 101. Como se puede observar, cuando aumenta el tama no del siste-
ma, las burbujas de inestabilidad existentes para tama nos peque nos se
transforman en inestabilidades que persisten cuando aumenta el acoplo.
Las guras corresponden a un potencial onsite de Morse e interacci on
repulsiva con C = 0,02 y
b
= 0,8.
(a)
3 2 1 0 1 2 3
2
1.5
1
0.5
0
0.5
1
E
(b) (c)
3 2 1 0 1 2 3
0.1
0.05
0
0.05
0.1
E
3 2 1 0 1 2 3
0.1
0.05
0
0.05
0.1
E
Fig. 5.6: Diagramas de bandas en el lmite anti-continuo (arriba) y con un acoplo

C = 0,005 (abajo) para un dark breather en el modelo de Frenkel
Kontorova cuyos par ametros son:
b
= 0,2, = 0,02 y F = 0,02. (b)
corresponde a un potencial de interacci on atractivo, mientras que (c)
lo hace con un potencial repulsivo. Se observa que el dark breather es
estable en ambos casos.
5.6. Conclusiones 77
5.6. Conclusiones
5.6.1. Sistemas hamiltonianos
Conclusi on 5.1: Al introducir el acoplamiento entre las partculas, el dark brea-
ther es estable para acoplo repulsivo e inestable para acoplo atractivo si el potencial
es blando, mientras que la tendencia es la inversa si el potencial es duro.
Conclusi on 5.2: Cuando el acoplo aumenta, el dark breather estable puede de-
jar de serlo mediante bifurcaciones oscilatorias, arm onicas y subarm onicas. Estas
ultimas no ocurren si el potencial onsite es simetrico.
Conclusi on 5.3: Las bifurcaciones oscilatorias tienen un comportamiento reen-
trante, es decir, se producen burbujas de inestabilidad cuando aumenta el acoplo.
Conclusi on 5.4: Un dark breather con condiciones de contorno peri odicas es in-
estable siempre que derive de un fon on con vector de ondas y las partculas de los
extremos vibren en fase. Este fen omeno lo denominamos inestabilidad de paridad
Conclusi on 5.5: Cuando el tama no del sistema es innito, las bifurcaciones os-
cilatorias persisten (dejan de ser reentrantes) cuando aumenta el acoplo y las bi-
furcaciones subarm onicas siguen existiendo.
5.6.2. Predicciones del Teorema de Estabilidad de Multibreathers
Conclusi on 5.6: El TEM puede predecir la estabilidad de phonobreathers en sis-
temas hamiltonianos. En la siguiente tabla se resumen los resultados para phono-
breathers que derivan de fonones con vector de onda q = 0 y q = , de forma
que el acoplo es atractivo (si fuera repulsivo, los resultados se invierten salvo en
las inestabilidades de paridad, las cuales no dependen del tipo de acoplo y estan
se naladas con un asterisco):
Extremos libres / jos C.c. periodicas
0 Inestable Estable Inestable Estable
(N par) Estable Inestable Estable Inestable
(N impar) Estable Inestable Inestable(*) Inestable(*)
Conclusi on 5.7: El TEM puede predecir la estabilidad de dark multibreathers en
sistemas hamiltonianos, siempre que las condiciones de contorno sean peri odicas.
Si los extremos estuvieran libres o jos, s olo puede predecir los casos en los que los
fragmentos en los que se divide la cadena sean inestables de por s. En la siguiente
tabla se resumen los resultados para psite dark breathers los cuales provienen de
fonones con vector de onda q = 0 y q = , de forma que el acoplo es atractivo (si
fuera repulsivo, los resultados se invierten salvo en las inestabilidades de paridad,
las cuales no dependen del tipo de acoplo y estan se naladas con un asterisco).
Los casos que no pueden ser predichos por el TEM se se nalan con un signo de
interrogaci on.
Extremos libres / jos C.c. periodicas
0 Inestable ? Inestable Estable
(
1
=
N
) ? Inestable Estable Inestable
(
1
=
N
) ? Inestable Inestable(*) Inestable(*)
5.6.3. Sistemas disipativos
Conclusi on 5.8: Al introducirse el acoplo, el dark breather es estable con inde-
pendencia del tipo de potencial onsite y del signo de la constante de acoplo.
Conclusi on 5.9: Si el potencial onsite es de tipo sineGordon, no aparecen bi-
furcaciones oscilatorias.

Floquet

Cargado por

Copyright:

Formatos disponibles

Floquet

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Floquet

Cargado por

Copyright:

Formatos disponibles

4.

(u(t)) + CJ(u(t)), (4.5)

(). Este valor se conserva

(), las cuales no se solapan. Adem as, E

() es una funci on continua, no cons-

() tal que pertenezca a la

(t) con el mismo autovalor E y ar-

(E), las bandas ser an simetricas respecto a

(0). Es decir, para que una soluci on sea estable,

Estas tienen por solucion:

(E) = 0, (t) sera peri odica y reversible en el tiempo, es decir,

Fig. 4.13: Diagramas de bandas en el lmite anti-continuo (izquierda) y para un

de orden N cuyas componentes son iguales a:

, salvo que se han

Fig. 4.16: Diagramas de bandas en el lmite anti-continuo (a) y para un aco-

Fig. 5.2: Diagramas de bandas en el lmite anti-continuo (arriba) y para un aco-

Fig. 5.3: Diagramas de bandas en el lmite anti-continuo (arriba), y para un aco-

es el n umero de ceros adyacentes en el c odigo. Con estos preliminares, se puede

1.4 1.5 1.6 1.7 1.8

Fig. 5.6: Diagramas de bandas en el lmite anti-continuo (arriba) y con un acoplo

También podría gustarte