Research
Hamiltoniano efectivo de un circuito LC con carga discreta bajo la aproximación semiclásica
E. Mamani1
M. Calcina-Nogales1
D. Sanjinés1
1Carrera de Física, Universidad Mayor de San Andrés, c. 27 Cota-Cota, Casilla de Correos 8635, La Paz, Bolivia.
Resumen
Se utiliza un modelo de un circuito LC con carga discreta y una fuente AC de frecuencia ω elevada para determinar su Hamiltoniano efectivo aplicando la aproximación semiclásica y el método de promediación temporal de Kapitza hasta 0ω-2. El resultado corresponde al de la propagación de una partícula en una red de enlace fuerte por el mecanismo de hopping bajo la acción simultánea de un campo homogéneo rápidamente oscilante y un campo estático lineal [1], donde la coordenada de posición de la partícula equivale a la carga eléctrica del condensador en el circuito LC. De manera análoga a resultados previos para el oscilador armónico en la red [2], en este caso también se encuentra una condición de bifurcación dada por la carga inicial del condensador que determina la aparición súbita de un voltaje DC. Otra predicción relevante es la supresión controlada de la corriente efectiva al elegir los parámetros de la fuente AC.
Palabras clave: Sistemas mesoscópicos; aproximación semiclásica; Hamiltoniano efectivo
Abstract
We use a model for a LC circuit with discrete charge and an AC source with high frequency ω to determine its effective Hamiltonian by the application of the semiclasical approximation and the Kapitza’s time-average method up to 0ω-2. The result corresponds to that of the propagation of a hopping particle in a lattice in the presence of a combined homogeneous rapidly oscillating field and a static linear field [1], where the particle’s position coordinate corresponds to the capacitor’s electrical charge in the LC circuit. We find results which are analogous to the previous ones concerning the bifurcation condition for the harmonic oscillator in the lattice [2]; in the LC circuit such a condition is given by the initial capacitor’s charge which determines the sudden onset of a DC voltage. Another relevant prediction is the controlled suppression of the effective current by choosing the parameters of the AC source.
Keywords: Mesoscopic systems; semiclassical approximation; effective Hamiltonian
PACS: 73.23.-b; 03.65.Sq; 05.45.-a; 72.30.+q
1.Introducción
El estudio de Hamiltonianos efectivos en la física del estado sólido, derivados a partir de técnicas de promediacion temporal, ha adquirido reciente interés desde el punto de vista de la factibilidad de manipulación de los elementos de salto hopping en el modelo de enlace fuerte 1,3,4. Dicha manipulación supone la supresión controlada de la interacción del electrón con sus primeros o segundos vecinos en la red, lo que da lugar a fenómenos interesantes tales como la oscilación de Bloch efectiva 1. El método de promediación temporal que se usa en este trabajo es el que desarrolló originalmente P. L. Kapitza 5,6 en el que la solución de las ecuaciones de movimiento se separa en una parte rápida con promedio nulo y una parte lenta cuyo promedio evoluciona según ecuaciones efectivas de movimiento. Una de las motivaciones para la deducción de Hamiltonianos efectivos es el hecho de que los algortimos numéricos para resolver las ecuaciones diferenciales que modelan los sistema físicos suelen presentar inestabilidades cuando la funciones oscilan rápidamente.
En el presente trabajo se aprovecha la analogía formal entre dos sistemas físicos diferentes, una red de estado sólido y un circuito cuántico con carga discreta, para aplicar el método de promediación temporal de Kapitza al segundo sistema, que consiste de un circuito LC en serie con una fuente rápidamente oscilante Fig. 1a. Asumiremos el marco conceptual de la teoría cuántica de circuitos con carga discreta 7 donde se supone que los autovalores del operador carga qˆ toman valores discretos nqe con n entero y qe la carga del electrón. El operador conjugado correspondiente a qˆ es el operador flujo ϕˆ=-iℏ∂/∂q tal que vale el conmutador fundamental [qˆ,ϕˆ&]=iℏ. El operador
ϕˆ
se puede reemplazar por un operador de diferencias finitas cf., por ejemplo, 8; sin embargo, la forma directa usual de tomar en cuenta los efectos de la carga discreta es a través del cambio formal ϕˆ→2ℏ/qesinqeϕˆ/2ℏ tal como se hace, por ejemplo, en 9,10. Los fenómenos físicos investigados en este tipo de sistemas han sido ampliamente reportados en la literatura; tales fenómenos comprenden, por ejemplo: bloqueo de Coulomb 11, corrientes persistentes 7, magnificación de la corriente cuántica 12. En estos casos se ha demostrado exitosamente la factibilidad de modelar la física de los dispositivos mesoscópicos a través del marco conceptual referido.
Para el propósito anunciado, este trabajo se organiza de la siguiente forma: en la Sec. 2 derivamos el Hamiltoniano efectivo independiente del tiempo del circuito LC, donde ya no aparece explícitamente la fuente fwt Fig. 1b. En la Sec. 3 describimos brevemente la analogía del transporte electrónico en una red de enlace fuerte y en un circuito LC; dicha analogía permitirá entender mejor la física asociada al concepto de “punto de bifurcación”. En la Sec. 4 se deducen dos fenómenos físicos asociados a la elección de los parámetros de control de la fuente AC, lo que permitiría, en principio, manipular la carga y la corriente efectivas del circuito. En la Sec. 5 se concluyen los aspectos más relevantes de este trabajo y se ofrece algunas perspectivas para futuras investigaciones.
2.Derivación del Hamiltoniano efectivo
Consideremos el circuito LC de la Fig. 1a con una fuente oscilatoria fωt de alta frecuencia. El Hamiltoniano clásico o función Hamiltoniana correspondiente es:
H0(ϕ,q;t)=ϕ22L+q22C+qf(ωt)
(1)
donde q es la carga eléctrica y ϕ=Lq̇ es el flujo magnético; ω0=1/LC es la frecuencia natural del circuito LC cuando fωt=0. La sustitución ϕ→2ℏ/qesinqeϕ/2ℏ referida en la Sec. 1 permite obtener el Hamiltoniano con carga discreta
H(ϕ,q;t)=-ℏ2qe2Lcosqeℏϕ+q22C+qf(ωt)+ℏ2qe2L,
(2)
donde se usó sin2x=1-cos2x/2. A continuación aplicaremos el formalismo semiclásico al Hamiltoniano a través de las ecuaciones clásicas de Hamilton
q˙=∂H∂ϕ, ϕ˙=-∂H∂Q
(3)
de donde se obtiene el sistema no-lineal de ecuaciones diferenciales acopladas
q̇=ℏqeLsinqeℏϕ, ϕ̇=-qC-f(ωt).
(4)
A continuación se invoca la transformación canónica ϕt=Φt-gt, con
gt≡∫0tfωtdt,
de forma que el sistema de Ecs. (4) se escribe ahora de manera más conveniente:
q̇=ℏqeLsinqeℏΦ-g, Φ̇=-qC.
(5)
Para este sistema de ecuaciones se propone la solución
q(t)=Q(t)+ξ(τ), Φ(t)=Ψ(t)+η(τ),
(6)
donde se definieron: t y τ≡ωt como los tiempos “lento” y “rápido”, Qt y Ψt como las “coordenadas lentas”, ξτ y ητ como las “coordenadas rápidas”. Reemplazando (6) en (5) se tiene
Q̇(t)+ωdξdτ=ℏqeLsinqeℏΨ+η-g,Ψ̇(t)+ωdηdτ=-1CQ+ξ.
(7)
Utilicemos a continuación la definición usual del promedio temporal
⟨⋅⟩≡1/2π∫02π⋅dτ,
tal que las coordenadas rápidas tengan un promedio nulo en el intervalo [0,2π]: ⟨ξt⟩=0,⟨ηt⟩=0, y las coordenadas lentas varíen muy poco: ⟨Ψt⟩=Φt,⟨Qt⟩=Qt. Así, aplicando dicho promedio temporal al sistema de Ecs. (7), se tiene
Q̇(t)=ℏqeLsinqeℏΨ+η-g,Ψ̇(t)=-QC.
(8)
Este es el sistema de ecuaciones de movimiento efectivas para las coordenadas lentas. De manera similar, (8) restando de (7), se obtiene el sistema de ecuaciones de movimiento efectivas para las coordenadas rápidas,
ωdξdτ=ℏqeLsinqeℏΨ+η-g-sinqeℏΨ+η-g,ωdηdτ=-ξC.
(9)
A continuación se introducen las series de potencias
ξ=∑i=0∞ϵiξi=ξ0+ϵξ1+ϵ2ξ2+⋯,η=∑i=0∞ϵiηi=η0+ϵη1+ϵ2η2+⋯,
(10)
donde ϵ≡t/τ=1/ω es el parámetro de pequeñez. Las expansiones subsecuentes se hacen de manera consistente hasta ϵ2 en el Apéndice A; allí se obtienen las ecuaciones efectivas de movimiento para Q̇t y Ψ̇t (cf. (A.17)),
Q̇(t)=ℏqeLF0sinqeℏΨ+ℏϵ22qeL2CF˜0sin2qeℏΨ,
(11)
Ψ̇(t)=-QC.
(12)
Con este resultado se calculará el correspondiente Hamiltoniano efectivo HQ,Ψ. De (11) y (12) se obtienen el cociente Q̇/Ψ̇=dQ/dΨ y la ecuación diferencial -Ψ̇dQ+Q̇dΨ=0 que resulta ser exacta pues -∂Ψ̇/∂Ψ=∂Q̇/∂Q. Así, la función constante HQ,Ψ, cuya diferencial exacta es dHQ,Ψ=-Ψ̇dQ+Q̇dΨ=0, se deduce a partir de la teoría de ecuaciones diferenciales como
H(Q,Ψ)=∫-Ψ̇dQ+∫Q̇-∂∂Ψ∫-Ψ̇dQdΨ.
(13)
Sustituyendo (11) y (12) en (13) se tiene
H(Q,Ψ)=-ℏ2qe2LF0cosqeℏΨ-ℏ2qe2Lω02ω2F˜0cos2qeℏΨ+Q22C,
(14)
que es finalmente la expresión buscada para el Hamiltoniano efectivo HQ,Ψ, a menos de una constante aditiva, con la misma estructura canónica de: Q̇=∂HQ,Ψ/∂Ψ, Ψ̇=-∂HQ,Ψ/∂Q; omitimos dicha constante ya que no contiene a las variables dinámicas Q,Ψ y por lo tanto resulta inmaterial en (11) y (12). Nótese que HQ,Ψ ya no contiene explícitamente el tiempo t ni la fuente externa que depende de t, pero su efecto se traduce en la aparición de una inductancia efectiva Lef. Así, HQ,Ψ representa al circuito LC “efectivo” para la carga lenta Qt Fig. 1 b. En el límite clásico se obtiene
limqe→0H(Q,Ψ)=Ψ22L+Q22C,
(15)
pues F0=1+Oqe2 y F˜0=Oqe2 cf. Apéndice B, así que adopta la forma del Hamiltoniano clásico para el circuito LC sin fuente fωt=0, como debe ser.
3. Analogía del transporte electrónico en dos sistemas físicos: red de enlace fuerte y circuito LC cuántico
La dinámica de una partícula cargada que se traslada en una red unidimensional con constante de red a por el mecanismo de hopping en presencia de un campo estático arbitrario U(x) y de un campo uniforme rápidamente oscilante fwt, se puede describir aproximadamente por el Hamiltoniano de enlace fuerte
Hx,p=-2Acosap
-2Bcos2ap+Ux+xfωt,
(16)
donde A, B son los elementos de hopping que caracterizan a las interacciones de la partícula con sus primeros y segundos vecinos en la red las interacciones con vecinos más distantes se suprimieron en debido a que sus elementos de hopping son pequeños comparados con 𝐴 y 𝐵; la frecuencia 𝜔 se considera mucho mayor que la del sistema físico en ausencia de la fuente fwt 13, por ejemplo, para el caso de la oscilación de Bloch cuando UX es lineal. La aplicación del método de Kapitza permite obtener el Hamiltoniano efectivo independiente del tiempo 1
H(X,Y)=-2A˜(X)cosaY-2B˜(X)cos2aY-2C˜(X)cos3aY-2D˜(X)cos4aY+U(X),
(17)
donde las coordenadas canónicas X,Y corresponden a las partes lentas de la posición x y del momentum desplazado y=p+g, con ġ=fωt. Para el caso particular de un potencial estático 𝑈𝑋 cuadrático, los coeficientes A˜,B˜,C˜,D˜ corresponden a los elementos efectivos de hopping que dependen de los parámetros de los campos externos estático y oscilante, tales que A˜→A, B˜→B, C˜→0, D˜→0, cuando ω→∞. En este caso, los Hamiltonianos HX,Y y HQ,Ψ en (14) son formalmente equivalentes después de establecer las identificaciones: X↔Q, Y↔qe/aℏΨ y B˜/A˜↔L˜/Lˆ.
La equivalencia formal de los problemas físicos descritos por sus correspondientes Hamiltonianos HX,Y en y HQ,Ψ en (17), permite la interesante posibilidad de asociar fenómenos en contextos físicos diferentes. En el caso que exploramos en este trabajo, invocaremos al “oscilador armónico en la red” descrito previamente por Gallinar y Chalbaud 2, donde, entre otros resultados, se dedujo la existencia de dos regímenes dinámicos referidos por “régimen del continuo” y “régimen de Bloch”, en los que la partícula cargada e.g., el electrón oscila con periodos característicos diferentes, dependiendo de la posición inicial de la partícula X0 en un potencial cuadrático desplazado de la forma UX=αX-β2+γ. La transición entre los regímenes dinámicos mencionados ocurre de manera súbita en un valor específico de X0 que se caracteriza como un “punto de bifurcación”, pues la dinámica de la partícula es diferente en ambos regímenes. Este mismo comportamiento (desarrollado en detalle en la Sec. 4.1) debería verificarse experimentalmente en el caso de un circuito LC cuántico en vista de la analogía formal entre los dos sistemas referidos en esta sección.
4. Ingeniería de la carga y de la corriente efectivas
En esta sección consideramos la expresión (11) para la variación temporal de la carga efectiva Qt. Los coeficientes F0 y F˜0 dependen de los parámetros de la fuente cf. Apéndice B, así que en principio es posible “manipular” los valores de la corriente y carga efectivas:Q̇t y Qt. Dicha “manipulación” corresponde al significado atribuido al término “ingeniería” en este trabajo. En particular, veremos que es posible lograr dos fenómenos que ocurren como consecuencia directa de la naturaleza discreta de la carga eléctrica: la generación súbita de un voltaje DC a través del condensador (Sec. 4.1) y la supresión de la corriente efectiva en la malla (Sec. 4.2).
4.1. Generación súbita de un voltaje DC
Consideremos la variación temporal de la corriente efectiva Q̇t en (11),
Q¨=1LF0cos(qeℏΨ)+(ω0ω)2F~0cos(2qeℏΨ)Ψ˙
(18)
De (14) se tiene la energía total del sistema,
E=Q22C-ℏ2qe2L[F0cosqeℏΨ+ω02ω2F˜0cos2qeℏΨ].
(19)
Para frecuencias ω altas, los términos entre corchetes en (18) y (19) son iguales entre sí hasta ω-2
Apéndice B, de tal forma que
Q̈=-qeℏ2Q32C2-QCE;
(20)
esta ecuación se integra utilizando Q̈
=Q̇dQ̇/dQ con las condiciones iniciales Q̇0=0 y Q0≠0,
Q̇2=1Cqeℏ2Q02-Q2Q02+Q24C-E.
(21)
Nótese que en el límite qe→0 con E=Eqe en la Ec. (20) se reduce a
Q̈+ω02Q=0,
(22)
que describe a un oscilador armónico simple, como debe ser. La Ec. (21) se puede representar como un diagrama en el espacio de fases Q,Q̇ que se muestra en la Fig. 2a para el caso del campo externo fωt=2f1cosωt
Apéndice B. Una característica peculiar de este tipo de diagramas es la transición súbita entre dos tipos diferentes de curvas: aquellas que oscilan en torno a Q=0 y las que oscilan en torno a Q≠0; la curva separatriz entre ambos tipos de curvas comprende al punto Q,Q̇=0,0 que en corresponde al valor inicial de la carga Q0=0. El “punto de bifurcación” que separa esos tipos de curvas es el valor especial de la carga inicial del condensador
Qˆ0=2ℏqeCLF0+ω02ω2F˜0,
(23)
que se calcula estableciendo Q̇=0 en (21) para los casos Q02<Qˆ02 y Q02>Qˆ02; en el primer caso se obtiene que Q∈[-Q0,Q0], mientras que en el segundo caso se obtiene que
Q∈-Q0,-Q02-Qˆ02∪Q02-Qˆ02,Q0.
(24)
Una simulación numérica de la solución de (21) para los mismos valores de la Fig. 2a se muestra en la Fig. 2b. Se observa los dos regímenes dinámicos característicos de la transición a través del punto de bifurcación: oscilaciones centradas en Q=0 Q0=±1,±2 y oscilaciones centradas en Q≠0
Q0=±3. Estas últimas son oscilaciones desplazadas que se pueden asociar al fenómeno físico de aparición súbita de un voltaje DC VDC diferente de cero, como una consecuencia de la naturaleza discreta de la carga eléctrica qe, pues se sabe que en el límite clásico qe→0, el valor de dicho voltaje a través del capacitor es VDC=0 caso de la Fig. 1a
i
.
La definición formal de VDC corresponde al valor RMS de Vt = Qt del capacitor para un valor unitario de la capacitancia: VDC≡〈Q2⟩T', donde el promedio temporal ⟨⋅⟩T' se toma sobre el periodo T´ de Qt. Debe notarse, sin embargo, que dicha definición se aplica sólo a una componente de Fourier de la función periódica Qt. El caso que corresponde a Qt en la Fig. 2b comprende, ciertamente, varias componentes de Fourier, por lo que la definición del voltaje DC se debe extender como VDC≡∑i⟨Qi2⟩T'i, donde el índice i corresponde a cada componente de Fourier. Un procedimiento más simple para calcular un valor numérico del voltaje DC es proponer la definición alternativa
VDC≡∫0T'dt QQ̇2∫0T'dt Q̇2,
(25)
que comprende de manera global a Q(t) y, por lo tanto, a todas sus componentes de Fourier. Esta re-definición de VDC no es más que el centroide de las trayectorias en el diagrama del espacio fase. En el límite clásico qe→0 se conoce Qt de acuerdo a (25), así que se evalúa las integrales de VDC en , de donde resultaVDC=0; esto se verifica sobre la base de la simetría de inversión inferida de (21) para la trayectoria en el espacio fase, pues limqe→0Qˆ0=∞ en (23) implica que cualquier valor finito de Q0 será menor que Qˆ0, entonces la oscilación simétrica de Q(t) en torno a Q=0 necesariamente corresponde a VDC=0. Por la misma razón de simetría, cuando qe>0 la evaluación de (25) para las trayectorias centradas en Q=0 es VDC=0. La ventaja de utilizar esta definición de VDC es que los valores de Q(t) y Q̇t que resultan de la solución numérica de (11) y (12) se incorporan de manera directa en la evaluación de (25) a través de una rutina numérica muy simple.
4.2. Supresión de la corriente efectiva en la malla
Éste es otro fenómeno físico interesante cuya eventual observación se podría sugerir en un intervalo adecuado de valores del campo eléctrico fωt=2f1cosωt+2f2cos2ωt. Notemos que la expresión para Q̇ en (11) es cero si F0 = 0 y F˜0=0 simultáneamente de acuerdo a (B.2) y (B.3). De manera equivalente, a partir de se demuestra que |Q̇max|=Qˆ02 qe/4ℏC para Q0≥Qˆ0, donde se verifica que |Q̇max|=0 cuando F0=F˜0=0. En la Fig. 3a se obtuvieron numéricamente las regiones en la que F0=0 curva continua y F˜0=0 curva segmentada, así que las intersecciones de dichas curvas corresponden a los valores del campo eléctrico que permitirían suprimir la corriente efectiva Q̇ de la malla. La condición F0=F˜0=0 implica que Qˆ0=0 y esto corresponde, según (21), al sitio geométrico en el espacio fase consistente de los puntos Q,Q̇=Q0,0 que yacen sobre el eje Q. En la Fig. 3b se aprecia la dependencia de |Q̇max| vs. f1,f2 de tal forma que se vea, por ejemplo, cuál sería la “ruta” óptima para suprimir la corriente de la manera menos abrupta posible. En el límite clásico del circuito LC se obtiene
limqe→0Q̇=ΨL
(26)
en (11), de manera consistente con la relación entre el flujo magnético y la corriente, ϕ=Lq̇, en el Hamiltoniano clásico Fig. 1a donde no es posible suprimir la corriente.
5. Conclusiones
En este trabajo se utilizó el modelo semiclásico de la dinámica electrónica y la técnica de promediación temporal de Kapitza para obtener las ecuaciones de movimiento efectivas y de las variables dinámicas carga y flujo de un circuito LC de carga discreta y fuente arbitraria rápidamente oscilante, de donde se derivó un Hamiltoniano efectivo independiente del tiempo. A partir de dichas ecuaciones de movimiento fue posible inferir dos fenómenos interesantes como consecuencia directa de la naturaleza discreta de la carga eléctrica. El primero (Sec. 4.1) se refiere al régimen dinámico especial en el que la carga del condensador oscila en torno a un valor Q≠0; dicha oscilación desplazada Fig. 2b se puede asociar al fenómeno físico de la aparición súbita de un voltaje DC en el condensador. El segundo fenómeno (sección 4.2) es la supresión de la corriente efectiva de la malla, que ocurriría para valores adecuados de las amplitudes de la fuente Fig. 3a.
La analogía formal del transporte electrónico en dos sistemas físicos diferentes, la red de enlace fuerte y el circuito LC cuántico Sec. 3, nos lleva a sugerir que aquellos fenómenos dinámicos que serían difícilmente observables en la red por varias limitaciones experimentales, podrían ser observados de manera equivalente en un circuito LC cuántico, así como el fenómeno de “localización asintótica” predicho para la dinámica electrónica en una red de enlace fuerte 14 fue observado de manera equivalente en un arreglo de fibras ópticas15. Por ejemplo, en el caso del “oscilador armónico en la red” (Sec. 3) estudiado por Gallinar y Chalbaud 2 se demostró la existencia de dos regímenes, “el régimen del continuo” y el “régimen de Bloch”, que están separados por un punto de bifurcación; en el caso del circuito LC cuántico, dicho punto de bifurcación separa los regímenes de carga continua VDC=0 y carga discreta VDC>0.
References
1. E. Mamani, M. Calcina-Nogales, y D. Sanjinés, Int. J. Mod. Phys. B 31 2017 1750116,.
[ Links ]
2. J.P. Gallinar y E. Chalbaud, Phys. Rev. B, 43 1991 2322.
[ Links ]
3. A.P. Itin y A. Neishtadt, Phys. Lett. A 378 2014 822.
[ Links ]
4. A. P. Itin y M. I. Katsnelson, Phys. Rev. Lett. 115 2015 075301.
[ Links ]
5. P. L. Kapitza, Zh. Eksp. Teor. Fiz. 21 1951 588.
[ Links ]
6. P. L. Kapitza, Collected papers of P.L. Kapitza. ed. D. ter Haar, Pergamon Press, Oxford, 1965.
[ Links ]
7. B. Chen y Y. Q. Li Phys. Rev. B 53 1996 4027.
[ Links ]
8. C. A. Utreras y D. Laroze, Mod. Phys. Lett. B 26 2012 1250138.
[ Links ]
9. J. C. Flores y E. Lazo, IEEE. 4 2005 688.
[ Links ]
10. M. Calcina-Nogales, Mod. Phys. Lett. B 27 2013 1350138.
[ Links ]
11. B. Chen, Phys. Lett. A 335 2005 103.
[ Links ]
12. J. C. Flores, Phys. Rev. B 66 2002 153410.
[ Links ]
13. L. D. Landau, Mecánica. Reverté, Barcelona, 1985.
[ Links ]
14. D. Sanjinés y J. P. Gallinar, Phys. Rev. B 64 2001 054301.
[ Links ]
15. F. Dreisow, O. Bahat-Treidel, D. Welt, y A. Szameit Opt. Lett. 36 2011 2065-2067.
[ Links ]
16. E. Mamani, Tesis Doctoral. Universidad Mayor de San Andrés, La Paz, Bolivia, 2017.
[ Links ]
17. M. Abramowitz y I. Stegun, Handjournal of Mathematical Functions. Dover, New York, 1965.
[ Links ]
Apéndice
A. Derivación de las ecuaciones de movimiento efectivas
En este apéndice derivaremos (11) y (12) como resultado de suponer las expansiones en series (10). Reemplazando en se tiene (10) en (9) se tiene
dξ0dτ+ϵdξ1dτ+ϵ2dξ2dτ=ℏqeLϵsinα-sinα,
(A.1)
dη0dτ+ϵdη1dτ+ϵ2dη2dτ=-1Cϵξ0+ϵξ1+ϵ2ξ2
(A.2)
donde αϵ=qe/ℏΨ+η0+ϵη1+ϵ2η2-g; luego se expande las funciones sin(⋅) en el lado derecho de (A.1) alrededor de ϵ=0 hasta ϵ2,
sinα=sinα0+ϵqeℏη1cosα0+qeℏη2cosα0-qe2ℏη12sinα0ϵ2,
(A.3)
con α0≡α0=qe/ℏΨ+η0-g. Sustituyendo (A.3) en (A.1), y comparando los términos de las mismas potencias de ϵ de las ecuaciones resultantes en (A.1) y (A.2), se obtienen:
(i) para ϵ0:
dξ0dτ=0, dη0dτ=0,
(A.4)
nótese que al integrar las Ecs. (A.4) se obtienen constantes, que se eligen iguales a cero a fin de que las funciones periódicas ξτ y ητ tengan promedios temporales nulos, por lo que queda ξ0=η0=0;
(ii) para ϵ1:
dξ1dτ=ℏqeLsinα1-⟨sinα1⟩,dη1dτ=0 ⇒ η1=0;
(A.5)
donde α1≡α0=qe/ℏΨ-g se evaluó en η0=0;
(iii) para ϵ2:
dξ2dτ=0 ⇒ ξ2=0,dη2dτ=-1Cξ1.
(A.6)
Con estos resultados, las ecuaciones de movimiento “efectivas” para Qt y Ψt dadas 8 en con ayuda de (A.3) quedan como
Q̇(t)=ℏqeLsinα1+qeℏϵ2η2cosα1,
(A.7)
Ψ̇(t)=-QC.
(A.8)
A continuación, se busca integrar (A.7) y (A.8) para encontrar un Hamiltoniano efectivo HQ,Ψ. Para comenzar expresemos sinα1 en forma compleja,
sinα1=12iei(qeℏ)Ψe-iqeℏg-e-i(qeℏ)Ψeiqeℏg;
(A.9)
ya que gτ es una función periódica, entonces Ft≡expiqe/ℏg también es una función periódica cuya expansión en serie de Fourier es Ft=∑nFnexpinτ. Así, (A.9) queda como
sinα1=12i∑nei(qeℏ)ΨFn*-e-i(qeℏ)ΨF-ne-inτ
(A.10)
Realizando el promedio temporal de y utilizando expin±mτ=δn,∓m se obtiene
sinα1=12iei(qe/ℏ)ΨF0*-e-i(qe/ℏ)ΨF0.
(A.11)
A continuación (A.10) y (A.11) se reemplazan en la primera ecuación de (A.5) y se integra en τ,
ξ1τ=ℏ2qeL∑n≠1n×ei(qeℏ)ΨFn*-e-i(qeℏ)ΨF-ne-inτ.
(A.12)
Por otra parte, reemplazando (A.12) en la segunda ecuación de (A.6) e integrando en τ se obtiene
η2τ=ℏ2qeLC∑n≠0e-inτin2×ei(qeℏ)ΨFn*-e-i(qeℏ)ΨF-n.
(A.13)
Con un procedimiento análogo al que se utilizó para obtener (A.10), se obtiene
cosα1=12∑n×ei(qeℏ)ΨF~n*-e-i(qeℏ)ΨF-ne-inτ
(A.14)
A continuación se calcula la expresión η2cosα1 según lo requerido en (A.7). Así, utilizando (A.13) y (A.14) resulta
η2cosα1=ℏ4iqeLC×e2i(qeℏ)ΨF~0*-e-2i(qeℏ)ΨF~0,
(A.15)
donde se definió F˜m≡∑n≠0FnFm-n/n2. Sustituyendo (A.11) y (A.15) en (A.7) se tiene que
Q̇(t)=ℏ2iqeLei(qe/ℏ)ΨF0*-e-i(qe/ℏ)ΨF0+ϵ22LCe2i(qe/ℏ)ΨF˜0*-e-2i(qe/ℏ)ΨF˜0.
(A.16)
Ya que fωt de la fuente es una función real, entonces se expande como fωt=∑nfnexpinωt, donde los coeficientes f
n
son reales y pares, i.e., fn=f-n=fn*; asimismo se supone que el promedio temporal de fwt en un periodo 2π/ω es nulo, fωt=0, por lo que f0=0. Así, los coeficientes Fn y F˜n son reales. Finalmente, las ecuaciones de movimiento efectivas para la carga Qt y el flujo magnético desplazado Ψt son (cf. (11) y (12))
Q̇(t)=ℏqeLF0sin(qe/ℏ)Ψ+ℏϵ22qeL2CF˜0sin2qeℏΨ,Ψ̇(t)=-QC.
(A.17)
B. Deducción de los coeficientes Fn
Para la función de la fuente fωt considerada en el Apéndice A, se sabe que los coeficientes Fn están dados por 16
Fn=δn,0+qeℏω∑m≠0fmmδn,m+qe22ℏ2ω2∑(m,p)≠0fmfpmpδn,m+p.
(B.1)
Por otra parte, de la definición de F˜n usada en (A.15), se tiene que
F˜0=∑n≠0(-1)nn2Fn2.
(B.2)
Luego, en el límite ω→∞ se tiene F0→1 y F˜0→0, lo que se usa en (15).
Consideremos a continuación el caso del campo fωt=2f1cosωt usado en la (Sec. 4.1). En este caso se obtienen
g(t)≡∫0tf(ωt)dt=(2f1/ω)sin(ωt)
(B.3)
y
F(t)=exp2iqesin(ωt)f1/(ℏω),
(B.4)
así que los coeficientes de Fourier en (A.11) son
Fn=(1/2π)∫02πF(τ)e-inτdτ=Jn(2qef1/(ℏω)),
(B.5)
donde se usó la función generatriz de las funciones de Bessel de primera clase y orden entero
eiξsin(ωt)=∑neinωtJn(ξ),
(B.6)
y
∫02πei(m-n)ωτdτ=2πδm,n.
(B.7)
Para valores suficientemente grandes de ω≡1/ϵ en el argumento de Jn en usamos las expansiones de Jn en series de potencias 17:
J0(ϵ)=1-(ϵ/2)2+O(ϵ4),Jn(ϵ)=(ϵ/2)/n!+O(ϵ3),
(B.8)
para n≥1. De (B.5) y (B.2) se obtiene las expansiones
F0=1-(qef1/(ℏω))2+O(ω-4),F˜0=-(7/8)(qef1/(ℏω))2+O(ω-4),
(B.9)
que se usan en la Sec. 4.1 junto con (ω0/(2ω))2F˜0/F0=0.05, f1/ω=0.22 y valores unitarios de L,C,qe,ℏ, de tal forma que Qˆ0=2 en (23)
Ahora consideremos el caso del campo f(ωt)=2f1cos(ωt)+2f2cos(2ωt) usado en la Sec. 4.2. Seguimos el mismo procedimiento de la primera parte de este Apéndice (Ecs. (B.3)-(B.8)) y obtenemos
Fn=∑mJn-2m2qef1ℏωJmqef2ℏω.
(B.10)
De aquí se obtienen los valores de F0 y F˜0 por evaluación numérica. En particular, se pueden calcular los valores de 2f1/ω y f2/ω tales que F0=0 y F˜0=0, como se muestra en la Fig. 3a.