Articulo de referencia

teoría de Floquet

En un sistema donde las fuerzas son periódicas —como un péndulo sometido a una fuerza impulsora periódica o un circuito oscilante alimentado por corriente alterna—, el comportam...

En un sistema donde las fuerzas son periódicas —como un péndulo sometido a una fuerza impulsora periódica o un circuito oscilante alimentado por corriente alterna—, el comportamiento general del sistema no es necesariamente completamente periódico. Por ejemplo, consideremos a un niño columpiándose: aunque el movimiento es impulsado por impulsos regulares y periódicos, el columpio puede alcanzar gradualmente mayores alturas mientras continúa oscilando. Esto da como resultado una combinación de periodicidad subyacente y crecimiento.

La teoría de Floquet proporciona una forma de analizar este tipo de sistemas. Su idea fundamental es similar al ejemplo del columpio: la solución se puede descomponer en dos partes: un componente periódico (que refleja el movimiento repetitivo) y un factor exponencial (que refleja el crecimiento, la disminución o la estabilidad neutra). Esta descomposición permite analizar el comportamiento a largo plazo y la estabilidad en sistemas periódicos en el tiempo.

Formalmente, la teoría de Floquet es una rama de las ecuaciones diferenciales ordinarias que se relaciona con la clase de soluciones de ecuaciones diferenciales lineales periódicas de la forma incógnita˙=A(t)incógnita,{\displaystyle {\dot {x}}=A(t)x,} conincógnitaRnorte{\displaystyle x\in {R^{n}}}yA(t)Rnorte×norte{\displaystyle \displaystyle A(t)\in {R^{n\times n}}}siendo una función periódica con períodoT{\displaystyle T}y define el estado de estabilidad de las soluciones.

El teorema principal de la teoría de Floquet, el teorema de Floquet , debido a Gaston Floquet ( 1883 ) , da una forma canónica para cada solución matricial fundamental de este sistema lineal común . Da un cambio de coordenadas. y=Q1(t)incógnita{\displaystyle \displaystyle y=Q^{-1}(t)x}conQ(t+2T)=Q(t){\displaystyle \displaystyle Q(t+2T)=Q(t)}que transforma el sistema periódico en un sistema lineal tradicional con coeficientes reales constantes . Cuando se aplica a sistemas físicos con potenciales periódicos, como los cristales en la física de la materia condensada , el resultado se conoce como el teorema de Bloch .

Nótese que las soluciones de la ecuación diferencial lineal forman un espacio vectorial . Una matrizϕ(t){\displaystyle \phi \,(t)}Se denomina solución matricial fundamental si las columnas forman una base del conjunto solución.Φ(t){\displaystyle \Phi (t)}Se denomina solución de matriz fundamental principal si todas las columnas son soluciones linealmente independientes y existet0{\displaystyle t_{0}}de tal manera queΦ(t0){\displaystyle \Phi (t_{0})}es la identidad. Se puede construir una matriz fundamental principal a partir de una matriz fundamental utilizandoΦ(t)=ϕ(t)ϕ1(t0){\displaystyle \Phi (t)=\phi \,(t){\phi \,}^{-1}(t_{0})}. La solución de la ecuación diferencial lineal con la condición inicialincógnita(0)=incógnita0{\displaystyle x(0)=x_{0}}esincógnita(t)=ϕ(t)ϕ1(0)incógnita0{\displaystyle x(t)=\phi \,(t){\phi \,}^{-1}(0)x_{0}}dóndeϕ(t){\displaystyle \phi \,(t)}es cualquier solución matricial fundamental.

Teorema de Floquet

Dejarincógnita˙=A(t)incógnita{\displaystyle {\dot {x}}=A(t)x}sea ​​una ecuación diferencial lineal de primer orden, dondeincógnita(t){\displaystyle x(t)}es un vector columna de longitudnorte{\displaystyle n}yA(t){\displaystyle A(t)}unnorte×norte{\displaystyle n\times n}matriz periódica con períodoT{\displaystyle T}(eso esA(t+T)=A(t){\displaystyle A(t+T)=A(t)}para todos los valores reales det{\displaystyle t}). Dejarϕ(t){\displaystyle \phi \,(t)}Sea una solución matricial fundamental de esta ecuación diferencial. Entonces, para todotR{\displaystyle t\in \mathbb {R} }, ϕ(t+T)=ϕ(t)ϕ1(0)ϕ(T).{\displaystyle \phi (t+T)=\phi (t)\phi ^{-1}(0)\phi (T).} Aquí ϕ1(0)ϕ(T){\displaystyle \phi ^{-1}(0)\phi (T)} se conoce como la matriz de monodromía . Además, para cada elección de matrizB{\displaystyle B}(posiblemente complejo) de tal manera que miTB=ϕ1(0)ϕ(T),{\displaystyle e^{TB}=\phi ^{-1}(0)\phi (T),} Hay un período (con período)T{\displaystyle T}) función matricialtPAG(t){\displaystyle t\mapsto P(t)}de tal manera que ϕ(t)=PAG(t)mitB a pesar de tR.{\displaystyle \phi (t)=P(t)e^{tB}{\text{ para todo }}t\in \mathbb {R} .} Esta representación se denomina forma normal de Floquet para la solución matricial fundamental.ϕ(t){\displaystyle \phi \,(t)}. [ 1 ]

Además, siϕ(t){\displaystyle \phi (t)}es una matriz real para cada valor det{\displaystyle t}, entonces también existe al menos una matriz realR{\displaystyle R}lo cual satisface ϕ1(0)ϕ(2T)=(ϕ1(0)ϕ(T))2=mi2TR.{\displaystyle \phi ^{-1}(0)\phi (2T)=(\phi ^{-1}(0)\phi (T))^{2}=e^{2TR}.} Para cualquier elección de talR{\displaystyle R}, hay un período real (con período-2T{\displaystyle 2T}) función matricialtQ(t){\displaystyle t\mapsto Q(t)}de tal manera que ϕ(t)=Q(t)mitR a pesar de tR.{\displaystyle \phi (t)=Q(t)e^{tR}{\text{ para todo }}t\in \mathbb {R} .} En lo anteriorB{\displaystyle B},PAG{\displaystyle P},Q{\displaystyle Q}yR{\displaystyle R}sonnorte×norte{\displaystyle n\times n}matrices. [ 2 ]

Ecuaciones periódicas de segundo orden

La teoría de Floquet se aplica al caso especial de una ecuación de segundo orden. y+V(incógnita)y=miy{\displaystyle -y''+V(x)y=Ey} dóndeV(incógnita){\displaystyle V(x)}es un potencial periódico,V(incógnita+T)=V(incógnita){\displaystyle V(x+T)=V(x)}. [ 3 ] [ 4 ] Estas ecuaciones surgen en el estudio de la ecuación de Hill , la ecuación de Lamé y la ecuación de Mathieu . El parámetro espectralmi{\displaystyle E}corresponde a los niveles de energía del operador de SchrödingerL=d2dincógnita2+V(incógnita){\displaystyle L=-{\frac {d^{2}}{dx^{2}}}+V(x)}. [ 5 ]

Para un valor dado del parámetromi{\displaystyle E}, dejar(incógnita,mi){\displaystyle u(x,E)}yv(incógnita,mi){\displaystyle v(x,E)}sea ​​el par de soluciones tales que (0,mi)=1,(0,mi)=0,v(0,mi)=0,v(0,mi)=1.{\displaystyle u(0,E)=1,\quad u'(0,E)=0,\quad v(0,E)=0,\quad v'(0,E)=1.} Debido a que el potencial es periódico, la superposición implica que existen constantesA(mi),B(mi),do(mi),D(mi){\displaystyle A(E),B(E),C(E),D(E)}con (incógnita+T,mi)=A(mi)(incógnita,mi)+B(mi)v(incógnita,mi),v(incógnita+T,mi)=do(mi)(incógnita,mi)+D(mi)v(incógnita,mi).{\displaystyle u(x+T,E)=A(E)u(x,E)+B(E)v(x,E),\quad v(x+T,E)=C(E)u(x,E)+D(E)v(x,E).} En forma matricial, ((incógnita+T,mi)v(incógnita+T,mi))=(A(mi)B(mi)do(mi)D(mi))((incógnita,mi)v(incógnita,mi)){\displaystyle {\begin{pmatrix}u(x+T,E)\\v(x+T,E)\end{pmatrix}}={\begin{pmatrix}A(E)&B(E)\\C(E)&D(E)\end{pmatrix}}{\begin{pmatrix}u(x,E)\\v(x,E)\end{pmatrix}}} expresando la monodromía de la solución fundamental. (La matriz aquí se llama matriz de monodromía)METRO(mi){\displaystyle M(E)}.) [ 3 ] [ 4 ]

Una solución se denomina solución de Floquet si y(incógnita+T)=ξy(incógnita){\displaystyle y(x+T)=\xi \,y(x)} por alguna constanteξ{\displaystyle \xi }llamado multiplicador de Floquet . Escribir una solución de Floquet comoy(incógnita)=a(incógnita,mi)+bv(incógnita,mi){\displaystyle y(x)=a\,u(x,E)+b\,v(x,E)}muestra que(a,b){\displaystyle (a,b)}debe ser un vector propio (izquierdo) deMETRO(mi){\displaystyle M(E)}, entoncesξ{\displaystyle \xi }Satisface ξ2(A(mi)+D(mi))ξ+detMETRO(mi)=0.{\displaystyle \xi ^{2}-(A(E)+D(E))\xi +\det M(E)=0.} Porque el Wronskiano de las soluciones(incógnita,mi){\displaystyle u(x,E)}yv(incógnita,mi){\displaystyle v(x,E)}es constante,detMETRO(mi)=1{\displaystyle \det M(E)=1}por lo que esta ecuación se simplifica a ξ2(A(mi)+D(mi))ξ+1=0{\displaystyle \xi ^{2}-(A(E)+D(E))\xi +1=0} o equivalentemente ξ+ξ1=A(mi)+D(mi).{\displaystyle \xi +\xi ^{-1}=A(E)+D(E).} Por lo tanto, existen genéricamente dos multiplicadores de Floquet para cada energía.mi{\displaystyle E}. [ 3 ] [ 4 ]

Los valores demi{\displaystyle E}para los cuales el multiplicador de Floquet satisfaceξ=±1{\displaystyle \xi =\pm 1}se denominan autovalores periódicos y antiperiódicos . Forman los extremos de las bandas espectrales del operador de Schrödinger periódico. Para valores reales demi{\displaystyle E}El comportamiento de las soluciones está determinado por la magnitud de los multiplicadores de Floquet. Una solución no trivial permanece acotada en la recta real si y solo si los multiplicadores correspondientes satisfacen|ξ|=1{\displaystyle |\xi |=1}. [ 4 ] [ 5 ] El conjunto de tales energías consiste en una unión de intervalos [mi1,mi2],[mi3,mi4],,{\displaystyle [E_{1},E_{2}],\;[E_{3},E_{4}],\;\ldots ,} llamadas bandas espectrales , separadas por intervalos ( brechas de banda ) donde las soluciones no están acotadas. Comomi{\displaystyle E}varía a lo largo de una banda espectral, los multiplicadores de Floquet se mueven a lo largo del círculo unitario en el plano complejo, y los puntos finales de cada banda ocurren precisamente cuandoξ=±1{\displaystyle \xi =\pm 1}.

Introduciendo una nueva variableν=ξξ1{\displaystyle \nu =\xi -\xi ^{-1}}, uno tiene ν2=(ξ+ξ1)24=(A(mi)+D(mi))24.{\displaystyle \nu ^{2}=(\xi +\xi ^{-1})^{2}-4=(A(E)+D(E))^{2}-4.} Esto define una cubierta de dos hojas de la línea de energía, llamada curva espectral . Para un potencial analítico general, esta cubierta es analítica en lugar de algebraica. Sin embargo, en situaciones especiales, se convierte en una curva algebraica . Por ejemplo, para la ecuación de Lamé.y(incógnita)gramo(gramo+1)(incógnita)y(incógnita)=miy(incógnita){\displaystyle y''(x)-g(g+1)\wp (x)y(x)=Ey(x)}con grado enterogramo{\displaystyle g}, la curva correspondiente es hiperelíptica y se puede escribir ν2=R2gramo+1(mi){\displaystyle \nu ^{2}=R_{2g+1}(E)} dóndeR2gramo+1{\displaystyle R_{2g+1}}es el polinomio espectral de Lamé. [ 6 ]

Estabilidad

Este mapeoϕ(t)=Q(t)mitR{\displaystyle \phi \,(t)=Q(t)e^{tR}}da lugar a un cambio de coordenadas dependiente del tiempo (y=Q1(t)incógnita{\displaystyle y=Q^{-1}(t)x}), bajo el cual nuestro sistema original se convierte en un sistema lineal con coeficientes constantes realesy˙=Ry{\displaystyle {\dot {y}}=Ry}. DesdeQ(t){\displaystyle Q(t)}Si es continua y periódica, debe estar acotada. Por lo tanto, la estabilidad de las solucionesy(t){\displaystyle y(t)}yincógnita(t){\displaystyle x(t)}están determinados por los valores propios deR{\displaystyle R}.

Los valores propios demiTB{\displaystyle e^{TB}}se denominan multiplicadores característicos del sistema. Mientras que la matrizB{\displaystyle B}no es único, los valores propios demiTB{\displaystyle e^{TB}}son lo mismo para las elecciones deB{\displaystyle B}También son los valores propios de los mapas de Poincaré (lineales).incógnita(t)incógnita(t+T){\displaystyle x(t)\to x(t+T)}Un exponente de Floquet (a veces llamado exponente característico) es un exponente complejo .μ{\displaystyle \mu }de tal manera quemiμT{\displaystyle e^{\mu T}}es un multiplicador característico del sistema. Nótese que los exponentes de Floquet no son únicos, ya quemi(μ+2πikT)T=miμT{\displaystyle e^{\left(\mu +{\frac {2\pi ik}{T}}\right)T}=e^{\mu T}}, dóndek{\displaystyle k}es un número entero. Las partes reales de los exponentes de Floquet son exponentes de Lyapunov . La solución cero es asintóticamente estable si todos los exponentes de Floquet tienen parte real negativa. Es estable en el sentido de Lyapunov si todos los exponentes de Floquet tienen parte real no positiva y, además, los multiplicadores de Floquet de módulo 1 son semisimples (equivalentemente, sus multiplicidades algebraica y geométrica coinciden); de lo contrario, es inestable. [ 1 ]

Aplicaciones

  • La teoría de Floquet se utiliza en el estudio de sistemas dinámicos , como la ecuación de Mathieu (que recibe su nombre de Émile Léonard Mathieu ).
  • La teoría de Floquet también se puede aplicar a sistemas dinámicos discretos y ecuaciones en diferencias .
  • La teoría de Floquet demuestra la estabilidad de las ecuaciones diferenciales de Hill (introducidas por George William Hill ). Dichas ecuaciones surgen, por ejemplo, al aproximar el movimiento de la luna como un oscilador armónico en un campo gravitatorio periódico .
  • El ablandamiento y el endurecimiento de los enlaces en campos láser intensos pueden describirse en términos de soluciones obtenidas a partir del teorema de Floquet.
  • La dinámica de los sistemas cuánticos fuertemente impulsados ​​se suele estudiar mediante la teoría de Floquet. En los circuitos superconductores, el marco teórico de Floquet se ha utilizado para dilucidar la electrodinámica cuántica de las interacciones entre múltiples cúbits inducidas por la fuerza impulsora y las transiciones de estado inducidas por dicha fuerza.
  • La teoría de Floquet también se utiliza en sistemas cuánticos de materia condensada impulsados ​​que son topológicamente no triviales. La topología del sistema impulsado se analiza mediante el estudio del hamiltoniano de Floquet.

Notas

Referencias

  • Chicone, Carmen Charles (1999). Ecuaciones diferenciales ordinarias con aplicaciones . Nueva York: Springer. ISBN 0-387-98535-2.
  • Deconinck, Bernard; Kutz, J. Nathan. "Cálculo de espectros de operadores lineales mediante el método de Hill" (PDF) . Universidad de Washington . Consultado el 10 de abril de 2026 .
  • Eastham, MSP (1973). La teoría espectral de las ecuaciones diferenciales periódicas . Edimburgo: Scottish Academic Press [distribuido por Chatto & Windus, Londres. ISBN 978-0-7011-1936-2.
  • Ekeland, Ivar (1990). "Uno". Métodos de convexidad en mecánica hamiltoniana . Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Resultados en Matemáticas y Áreas Afines (3)]. vol.  19. Berlín: Springer-Verlag. págs.x  +247. ISBN 3-540-50613-6. MR 1051888 . 
  • Floquet, Gaston (1883), "Sur les équations différentielles linéaires à coeficientes périodiques" (PDF) , Annales Scientifiques de l'École Normale Supérieure , 12 : 47– 88, doi : 10.24033/asens.220
  • Krasnosel'skii, MA (1968), El operador de traslación a lo largo de las trayectorias de ecuaciones diferenciales , Providence : American Mathematical Society, Traducción de Monografías Matemáticas, 19, 294 págs.
  • Magnus, Wilhelm; Winkler, Stanley (2013). La ecuación de Hill . Mineola, NY: Dover Publications. ISBN 978-0-486-15029-1.
  • Maier, Robert S (28-03-2008). "Polinomios de Lamé, reducciones hiperelípticas y estructura de banda de Lamé". Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences . 366 (1867): 1115– 1153. arXiv : math-ph/0309005 . doi : 10.1098/rsta.2007.2063 . ISSN 1364-503X . 
  • NW McLachlan, Teoría y aplicación de las funciones de Mathieu , Nueva York: Dover, 1964.
  • Montagnier, Pierre; Paige, Christopher C.; Spiteri, Raymond J. (2003). "Factores de Floquet reales de sistemas lineales periódicos en el tiempo" . Systems & Control Letters . 50 (4): 251– 262. doi : 10.1016/S0167-6911(03)00158-0 .
  • Teschl, Gerald (2012). Ecuaciones diferenciales ordinarias y sistemas dinámicos . Providence : American Mathematical Society . ISBN 978-0-8218-8328-0.
  • Deng, Chunqing; Shen, Feiruo; Ashhab, Sahel; Lupascu, Adrian (27-09-2016). "Dinámica de un sistema de dos niveles bajo fuerte excitación: optimización de puerta cuántica basada en la teoría de Floquet". Physical Review A. 94 ( 3) 032323. arXiv : 1605.08826 . Bibcode : 2016PhRvA..94c2323D . doi : 10.1103/PhysRevA.94.032323 . ISSN 2469-9926 . 
  • Huang, Ziwen; Mundada, Pranav S.; Gyenis, András; Schuster, David I.; Houck, Andrew A.; Koch, Jens (22 de marzo de 2021). "Ingeniería de puntos óptimos dinámicos para proteger los cúbits del ruido 1/f". Physical Review Applied . 15 (3) 034065. arXiv : 2004.12458 . Bibcode : 2021PhRvP..15c4065H . doi : 10.1103/PhysRevApplied.15.034065 . ISSN 2331-7019 . 
  • Nguyen, LB; Kim, Y.; Hashim, A.; Goss, N.; Marinelli, B.; Bhandari, B.; Das, D.; Naik, RK; Kreikebaum, JM; Jordan, A.; Santiago, DI; Siddiqi, I. (16 de enero de 2024). "Interacciones de Heisenberg programables entre cúbits de Floquet" . Nature Physics . 20 (1): 240– 246. arXiv : 2211.10383 . Bibcode : 2024NatPh..20..240N . doi : 10.1038/s41567-023-02326-7 .