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

Sintesi Tecnica: Approccio di Lewis–Ermakov per Sistemi a Due Livelli Tempo-Dipendenti

Enunciato del Problema
I sistemi a due livelli (TLS) tempo-dipendenti sono fondamentali per il controllo quantistico, eppure vengono tipicamente analizzati come casi distinti dall'oscillatore armonico tempo-dipendente (THO). Sebbene il metodo di Lewis–Ermakov fornisca un framework robusto per il THO — producendo invarianti dinamici, operatori di evoluzione in forma chiusa e propagatori esatti per parametri tempo-dipendenti arbitrari — un framework di Lewis–Ermakov per i TLS che sia unificato e pienamente sviluppato è stato finora assente. Nello specifico, una formulazione caratterizzata da un invariante dinamico, un hamiltoniano tempo-dipendente fattorizzato e un propagatore in forma chiusa valido per un accoppiamento $a(t)$ e un detuning $b(t)$ arbitrari non è stata stabilita nello stesso linguaggio unificato utilizzato per il THO. Questa lacuna è particolarmente significativa in contesti aperti o effettivamente non-Hermitiani, dove i punti eccezionali e il guadagno dipendente dallo stato complicano il concetto di inseguimento adiabatico.

Metodologia
Gli autori costruiscono un esplicito invariante dinamico di tipo Lewis–Ermakov per i TLS sfruttando la corrispondenza algebrica tra il TLS e il THO. La metodologia si basa sul fatto che entrambi i sistemi realizzano rappresentazioni della stessa algebra complessa, $sl(2, \mathbb{C})$:

• Il THO utilizza l'algebra non compatta $su(1, 1)$ generata da $\{\hat{x}^2, \hat{p}^2, \hat{x}\hat{p} + \hat{p}\hat{x}\}$.
• Il TLS utilizza l'algebra compatta $su(2)$ generata da $\{\hat{\sigma}_-, \hat{\sigma}_+, \hat{\sigma}_z\}$.

La procedura centrale consiste nel mappare l'hamiltoniano tempo-dipendente del TLS, $\hat{H}_{TLS}(t) = a(t)(\hat{\sigma}_+ + \hat{\sigma}_-) + b(t)\hat{\sigma}_z$, su una forma stazionaria tramite una coppia di trasformazioni tempo-dipendenti, $\hat{G}(t)$ e $\hat{T}(t)$.

1. Trasformazione: La trasformazione $\hat{G}(t)$ incorpora il detuning e la scalatura, mentre $\hat{T}(t)$ rimuove i termini residui. Notevolmente, $\hat{G}(t)$ è non-unitaria a causa della natura nilpotente di $\hat{\sigma}_-$, rendendo necessario uno spazio di stato intermedio complesso.
2. Equazione di Ermakov-Pinney: La mappatura richiede una funzione di scalatura ausiliaria $\rho(t)$ che soddisfi un'equazione di Ermakov-Pinney complessa:
   $$\ddot{\rho} + \theta^2(t)\rho = \frac{\theta_0^2}{\rho^3}$$
   Qui, la frequenza efficace $\theta^2(t)$ codifica i parametri fisici $a(t)$, $b(t)$ e le loro derivate temporali:
   $$\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. Costruzione dell'Invariante: L'invariante dinamico $\hat{I}_{TLS}(t)$ è costruito agendo con la trasformazione composta $\hat{S}_{TLS}(t) = \hat{G}(t)\hat{T}(t)$ sull'operatore stazionario $(\hat{\sigma}_+ + \hat{\sigma}_-)$. Questo invariante soddisfa la condizione di Lewis $\frac{d\hat{I}}{dt} = \frac{\partial \hat{I}}{\partial t} - i[\hat{I}, \hat{H}] = 0$.
4. Operatore di Evoluzione: Poiché l'hamiltoniano trasformato commuta con se stesso in ogni istante, l'operatore di ordinamento temporale diventa banale. Il propagatore esatto è derivato in forma chiusa, espresso in termini della funzione di scalatura $\rho(t)$ e di una fase accumulata $f(t) = \int dt'/\rho^2(t')$.

Contributi Chiave e Risultati
Il documento presenta i seguenti risultati specifici:

• Soluzione Generale: Un operatore di evoluzione e un propagatore esatti in forma chiusa per arbitrari $a(t)$ e $b(t)$, a condizione che l'equazione di Ermakov-Pinney complessa per $\rho(t)$ ammetta una soluzione.
• Transizioni di Landau-Zener-Majorana-Stückelberg (LZMS): Il framework viene applicato al modello standard LZMS ($a(t)=g, b(t)=\alpha t$). Il metodo recupera la dinamica non-adiabatica e le probabilità di transizione note, con la frequenza complessa $\theta^2(t)$ che cattura gli effetti del tasso di scansione tramite una componente immaginaria $i\alpha$.
• Dinamica Non-Hermitiana: L'approccio è esteso a un modello LZMS dissipativo con detuning immaginario $b(t) = i\alpha t$. Gli autori dimostrano che, mentre il punto eccezionale statico (dove gli autovalori coalescono) è determinato dallo spettro istantaneo, il punto eccezionale dinamico (dove inizia l'amplificazione esponenziale) è traslato. Questa traslazione è governata dallo zero della frequenza efficace $\theta^2(t)$, evidenziando come i tassi di scansione finiti alterino l'inizio della rottura della simmetria.
• Passaggio Adiabatico Rapido (ARP): Il metodo viene applicato al modello Allen-Eberly (accoppiamento secante iperbolica e detuning tangente iperbolica). L'approccio basato sull'invariante fornisce una descrizione esatta e priva di approssimazioni del trasferimento di popolazione, catturando accuratamente le deviazioni dalla perfetta adiabaticità dovute ai tassi di scansione finiti.
• Ingegneria Inversa (Shortcuts to Adiabaticity): Gli autori sviluppano un protocollo per progettare campi di controllo $a(t)$ e $b(t)$ che guidino il sistema lungo una traiettoria specificata senza eccitazioni residue (attrito quantistico).
  • Promuovendo la costante di integrazione $\theta_0$ a un parametro geometrico libero $\vartheta_0$, derivano condizioni al contorno per $\rho(t)$ che garantiscono che l'invariante commuti con l'hamiltoniano all'inizio e alla fine del protocollo.
  • Per una funzione di scalatura $\rho(t)$ reale, l'ingegneria inversa produce campi di controllo completamente espliciti e in forma chiusa: $a(t) = \vartheta_0 \cos\phi / \rho^2(t)$ e $b(t) = \vartheta_0 \sin\phi / \rho^2(t)$, dove $\phi$ è un angolo fisso. Ciò restringe il protocollo a una modulazione proporzionale di accoppiamento e detuning, ma garantisce la risolvibilità analitica.
  • Per una $\rho(t)$ complessa, i campi sono determinati da un sistema di equazioni differenziali integrali non lineari accoppiate, offrendo maggiore flessibilità nella modellazione indipendente di $a(t)$ e $b(t)$, ma richiedendo una soluzione numerica.

Significatività e Rivendicazioni
Il documento sostiene che il suo contributo primario sia la costruzione esplicita e l'applicazione pratica dell'invariante dinamico di Lewis–Ermakov per i sistemi a due livelli tempo-dipendenti. La significatività risiede in:

1. Unificazione: Stabilisce un legame concettuale e matematico tra le realizzazioni di variabili continue (oscillatore) e variabili discrete (due livelli) dell'algebra $sl(2, \mathbb{C})$, mostrando che la stessa sequenza di fattorizzazione, equazioni ausiliarie e costruzione dell'invariante si applica a entrambi.
2. Utilità Computazionale: Sostituisce le integrazioni specifiche per problema con una singola equazione ausiliaria (l'equazione di Ermakov-Pinney complessa), fornendo propagatori esatti per arbitrari campi di guida, inclusi i regimi non-unitari.
3. Design del Controllo: Eleva l'invariante da strumento teorico a risorsa pratica per l'ingegneria inversa. Fornisce un metodo sistematico per progettare "scorciatoie verso l'adiabaticità" con campi di controllo completamente espliciti (nel caso $\rho$-reale) o determinati numericamente (nel caso $\rho$-complesso), eliminando la necessità di vincoli arbitrari sui campi di guida.
4. Generalizzabilità: Gli autori osservano che il framework algebrico si estende naturalmente ad arbitrarie sistemi a $N$ livelli (rappresentazioni spin-$j$), dove la dinamica rimane governata da un'equazione di Ermakov-Pinney con coefficienti modificati, suggerendo un approccio universale per la dinamica quantistica multi-livello.

Il lavoro non propone nuovi setup sperimentali, ma fornisce piuttosto il apparato teorico per analizzare scenari esistenti (ioni intrappolati, circuiti superconduttori, ecc.) e progettare protocolli di controllo per essi.