Lewis-Ermakov approach for the time-dependent two-level system

Resumen Técnico: Enfoque de Lewis–Ermakov para Sistemas de Dos Niveles Dependientes del Tiempo

Planteamiento del Problema
Los sistemas de dos niveles (TLS, por sus siglas en inglés) dependientes del tiempo son fundamentales para el control cuántico, sin embargo, suelen analizarse como casos distintos del oscilador armónico dependiente del tiempo (THO). Si bien el método de Lewis–Ermakov proporciona un marco robusto para el THO —proporcionando invariantes dinámicos, operadores de evolución en forma cerrada y propagadores exactos para parámetros arbitrarios dependientes del tiempo—, ha faltado un marco de Lewis–Ermakov para los TLS que sea unificado y plenamente desarrollado. Específicamente, no se había establecido una formulación que presentara un invariante dinámico, un Hamiltoniano factorizado dependiente del tiempo y un propagador en forma cerrada válido para cualquier acoplamiento $a(t)$ y desintonía $b(t)$ en el mismo lenguaje unificado utilizado para el THO. Esta brecha es particularmente significativa en entornos abiertos o efectivamente no hermitianos, donde los puntos excepcionales y la ganancia dependiente del estado complican el concepto de seguimiento adiabático.

Metodología
Los autores construyen un invariante dinámico explícito de tipo Lewis–Ermakov para los TLS explotando la correspondencia algebraica entre el TLS y el THO. La metodología se basa en el hecho de que ambos sistemas realizan representaciones de la misma álgebra de Lie complejizada, $sl(2, \mathbb{C})$:

• El THO utiliza la álgebra no compacta $su(1, 1)$ generada por $\{\hat{x}^2, \hat{p}^2, \hat{x}\hat{p} + \hat{p}\hat{x}\}$.
• El TLS utiliza la álgebra compacta $su(2)$ generada por $\{\hat{\sigma}_-, \hat{\sigma}_+, \hat{\sigma}_z\}$.

El procedimiento central consiste en mapear el Hamiltoniano de TLS dependiente del tiempo, $\hat{H}_{TLS}(t) = a(t)(\hat{\sigma}_+ + \hat{\sigma}_-) + b(t)\hat{\sigma}_z$, sobre una forma estacionaria mediante un par de transformaciones dependientes del tiempo, $\hat{G}(t)$ y $\hat{T}(t)$.

1. Transformación: La transformación $\hat{G}(t)$ incorpora la desintonía y el escalamiento, mientras que $\hat{T}(t)$ elimina los términos residuales. Notablemente, $\hat{G}(t)$ no es unitaria debido a la naturaleza nilpotente de $\hat{\sigma}_-$, lo que requiere un espacio de estados intermedio complejizado.
2. Ecuación de Ermakov-Pinney: El mapeo requiere una función de escalamiento auxiliar $\rho(t)$ que satisfaga una ecuación de Ermakov-Pinney compleja:
   $$\ddot{\rho} + \theta^2(t)\rho = \frac{\theta_0^2}{\rho^3}$$
   Aquí, la frecuencia efectiva $\theta^2(t)$ codifica los parámetros físicos $a(t)$, $b(t)$ y sus derivadas temporales:
   $$\theta^2(t) = a^2 + b^2 + i\dot{b} - i\frac{\dot{a}b}{a} + \frac{\ddot{a}}{2a} - \frac{3\dot{a}^2}{4a^2}$$
3. Construcción del Invariante: El invariante dinámico $\hat{I}_{TLS}(t)$ se construye actuando con la transformación compuesta $\hat{S}_{TLS}(t) = \hat{G}(t)\hat{T}(t)$ sobre el operador estacionario $(\hat{\sigma}_+ + \hat{\sigma}_-)$. Este invariante satisface la condición de Lewis $\frac{d\hat{I}}{dt} = \frac{\partial \hat{I}}{\partial t} - i[\hat{I}, \hat{H}] = 0$.
4. Operador de Evolución: Debido a que el Hamiltoniano transformado conmuta consigo mismo en todo momento, el operador de ordenamiento temporal se vuelve trivial. El propagador exacto se deriva en forma cerrada, expresado en términos de la función de escalamiento $\rho(t)$ y una fase acumulada $f(t) = \int dt'/\rho^2(t')$.

Contribuciones Clave y Resultados
El artículo presenta los siguientes resultados específicos:

• Solución General: Un operador de evolución y un propagador exactos en forma cerrada para cualquier $a(t)$ y $b(t)$, siempre que la ecuación de Ermakov-Pinney compleja para $\rho(t)$ admita una solución.
• Transiciones de Landau-Zener-Majorana-Stückelberg (LZMS): El marco se aplica al modelo estándar de LZMS ($a(t)=g, b(t)=\alpha t$). El método recupera la dinámica no adiabática y las probabilidades de transición conocidas, donde la frecuencia compleja $\theta^2(t)$ captura los efectos de la tasa de barrido mediante un componente imaginario $i\alpha$.
• Dinámica No Hermitiana: El enfoque se extiende a un modelo LZMS disipativo con desintonía imaginaria $b(t) = i\alpha t$. Los autores demuestran que, si bien el punto excepcional estático (donde los autovalores coalescen) está determinado por el espectro instantáneo, el punto excepcional dinámico (donde comienza la amplificación exponencial) se desplaza. Este desplazamiento está gobernado por el cero de la frecuencia efectiva $\theta^2(t)$, resaltando que las tasas de barrido finitas alteran el inicio de la ruptura de simetría.
• Paso Adiabático Rápido (ARP): El método se aplica al modelo de Allen-Eberly (acoplamiento de secante hiperbólica y desintonía de tangente hiperbólica). El enfoque basado en invariantes proporciona una descripción exacta y libre de aproximaciones de la transferencia de población, capturando con precisión las desviaciones de la adiabaticidad perfecta debido a las tasas de barrido finitas.
• Ingeniería Inversa (Atajos hacia la Adiabaticidad): Los autores desarrollan un protocolo para diseñar campos de control $a(t)$ y $b(t)$ que dirijan el sistema a lo largo de una trayectoria específica sin excitaciones residuales (fricción cuántica).
  • Al promover la constante de integración $\theta_0$ a un parámetro geométrico libre $\vartheta_0$, derivan condiciones de contorno para $\rho(t)$ que aseguran que el invariante conmute con el Hamiltoniano al inicio y al final del protocolo.
  • Para una función de escalamiento $\rho(t)$ real, la ingeniería inversa produce campos de control totalmente explícitos y en forma cerrada: $a(t) = \vartheta_0 \cos\phi / \rho^2(t)$ y $b(t) = \vartheta_0 \sin\phi / \rho^2(t)$, donde $\phi$ es un ángulo fijo. Esto restringe el protocolo a una modulación proporcional del acoplamiento y la desintonía, pero garantiza la solvabilidad analítica.
  • Para un $\rho(t)$ complejo, los campos se determinan mediante un sistema de ecuaciones integro-diferenciales no lineales acopladas, lo que ofrece mayor flexibilidad para dar forma a $a(t)$ y $b(t)$ de manera independiente, pero requiere solución numérica.

Significancia y Reivindicaciones
El artículo sostiene que su contribución principal es la construcción explícita y la aplicación práctica del invariante dinámico de Lewis–Ermakov para sistemas de dos niveles dependientes del tiempo. La significancia radica en:

1. Unificación: Establece un vínculo conceptual y matemático entre las realizaciones de variables continuas (oscilador) y variables discretas (dos niveles) de la álgebra $sl(2, \mathbb{C})$, mostrando que la misma secuencia de factorización, ecuaciones auxiliares y construcción de invariantes se aplica a ambos.
2. Utilidad Computacional: Reemplaza las integraciones específicas de cada problema con una única ecuación auxiliar (la ecuación de Ermakov-Pinney compleja), proporcionando propagadores exactos para cualquier campo de excitación, incluyendo regímenes no unitarios.
3. Diseño de Control: Eleva el invariante de una herramienta teórica a un recurso práctico para la ingeniería inversa. Proporciona un método sistemático para diseñar "atajos hacia la adiabaticidad" con campos de control totalmente explícitos (en el caso de $\rho$ real) o campos determinados numéricamente (en el caso de $\rho$ complejo), eliminando la necesidad de restricciones arbitrarias sobre los campos de excitación.
4. Generalizabilidad: Los autores señalan que el marco algebraico se extiende naturalmente a sistemas de $N$ niveles (representaciones de espín-$j$), donde la dinámica sigue gobernada por una ecuación de Ermakov-Pinney con coeficientes modificados, lo que sugiere un enfoque universal para la dinámica cuántica multinivel.

El trabajo no propone nuevos montajes experimentales, sino que proporciona la maquinaria teórica para analizar escenarios existentes (iones atrapados, circuitos superconductores, etc.) y diseñar protocolos de control para ellos.