Semiclassical Langevin dynamics of long-range dissipative time crystals

Sintesi Tecnica: Dinamica di Langevin Semiclassica per Cristalli Temporali Dissipativi a Lungo Raggio

Definizione del Problema
I cristalli temporali di confine (BTC) dissipativi rappresentano una fase della materia in cui la simmetria di traslazione temporale continua o discreta viene spontaneamente rotta nel limite termodinamico, dando origine a oscillazioni collettive persistenti. Mentre i modelli idealizzati con interazioni a lungo raggio (totalmente collettive) supportano un comportamento stabile di ciclo limite descritto da equazioni di campo medio, i sistemi reali spesso coinvolgono interazioni a corto raggio e dissipazione locale. Una questione critica aperta è la robustezza della fase del cristallo temporale quando la struttura collettiva ideale viene indebolita. Nello specifico, è necessario comprendere come gli effetti di dimensione finita influenzino la durata di queste oscillazioni e identificare le soglie dell'intervallo di interazione (parametrizzato da un esponente di legge di potenza $\alpha$) oltre le quali la fase del cristallo temporale si dissolve in uno stato stazionario. Le simulazioni numeriche esatte utilizzando l'equazione maestra di Lindblad sono limitate a dimensioni di sistema ridotte, rendendo necessari metodi alternativi per indagare lo scaling di dimensione finita nel limite termodinamico.

Metodologia
Gli autori sviluppano un approccio di Langevin semiclassico radicato nei framework dell'approssimazione di Wigner troncata (TWA) e della DTWA (discrete TWA). Invece di risolvere la completa equazione maestra di Lindblad, la dinamica è formulata tramite equazioni di Langevin quantistiche per operatori di spin locali.

• Formalismo: Il metodo sostituisce i prodotti di operatori con prodotti dei loro valori di aspettazione (approssimazione di campo medio), pur mantenendo termini di rumore stocastico per rendere conto delle fluttuazioni quantistiche di ordine superiore generate dai canali dissipativi.
• Modello Spin-1/2: Per una catena di spin-1/2 con operatori di Lindblad a legge di potenza (la dissipazione è a lungo raggio, l'Hamiltoniana è locale), la dinamica è mappata in equazioni differenziali stocastiche per le componenti di spin locali $s_i^\mu$. Il rumore è rappresentato da rumori bianchi gaussiani reali.
• Modello Spin-1: Per un sistema spin-1 con dissipazione locale e interazioni Hamiltoniane a lungo raggio, la dinamica è formulata utilizzando una base di matrici di Gell-Mann (i generatori di $SU(3)$). Le variabili locali sono i valori di aspettazione di queste matrici, $\lambda_i^\mu$.
• Implementazione Numerica: Le equazioni stocastiche sono risolte utilizzando un'interpretazione di Itô con uno schema Runge-Kutta per il drift deterministico e incrementi di Wiener gaussiani per il rumore. Fondamentalmente, gli autori riscontrano che, per i modelli studiati, il campionamento su diverse condizioni iniziali (standard nella DTWA) non è necessario; una singola condizione iniziale è sufficiente perché la memoria dello stato iniziale viene persa rapidamente rispetto al periodo di oscillazione.
• Diagnosi: L'ordine del cristallo temporale viene diagnosticato analizzando lo scaling di dimensione finita della stabilità delle oscillazioni. Per lo spin-1/2, l'attenzione è rivolta al tasso di decadimento dell'inviluppo delle oscillazioni. Per lo spin-1, vengono utilizzati due strumenti: il "tempo di deviazione" (quanto tempo la traiettoria a dimensione finita rimane vicina al ciclo limite di campo medio) e lo scaling della potenza del picco di Fourier dominante.

Risultati Chiave

1. Modello Spin-1/2 (Dissipazione a Legge di Potenza):

   • Nel limite totalmente collettivo ($\alpha=0$), il sistema esibisce oscillazioni robuste dove lo smorzamento dovuto alla dimensione finita è soppresso all'aumentare della dimensione del sistema $L$.
   • Per esponenti di legge di potenza finiti $\alpha$, l'inviluppo delle oscillazioni decade esponenzialmente nel tempo, ma il tasso di decadimento $D$ scala algebricamente con la dimensione del sistema: $D(L) \sim L^{-\beta}$.
   • Un esponente di scaling positivo $\beta$ indica che il tasso di decadimento svanisce nel limite termodinamico, segnalando un cristallo temporale stabile.
   • L'esponente $\beta$ diminuisce monotonicamente all'aumentare di $\alpha$. Il regime del cristallo temporale persiste per $\alpha \lesssim 1.2$. Notevolmente, questa soglia è leggermente superiore ad $\alpha=1$. Mentre la teoria di campo medio prevede il rilassamento verso un punto fisso per $\alpha > 1$ (dove il campo medio non è più esatto nel limite termodinamico), l'approccio di Langevin semiclassico prevede oscillazioni persistenti nell'intervallo $1 < \alpha \lesssim 1.2$. Ciò suggerisce che le fluttuazioni quantistiche, catturate dal metodo, stabilizzano il cristallo temporale oltre il regime di campo medio.

2. Modello Spin-1 (Dissipazione Locale, Hamiltoniana a Lungo Raggio):

   • Il sistema esibisce un comportamento da cristallo temporale quando l'interazione coerente è sufficientemente a lungo raggio.
   • Tempo di Deviazione: Per $\alpha \leq 1$, il tempo $t^*$ in cui la traiettoria a dimensione finita devia dalla soluzione di campo medio scala come $t^*(L) \sim L^b$. L'esponente $b$ è positivo per piccoli $\alpha$ (circa 0.37 per $\alpha=0$) e diminuisce man mano che $\alpha$ si avvicina a 1, svanendo vicino a $\alpha \approx 1.2$. Questa analisi è strettamente valida solo per $\alpha \leq 1$, dove la teoria di campo medio descrive il limite termodinamico.
   • Scaling del Picco di Fourier: La potenza del picco di Fourier dominante scala come $P(L) \sim L^s$. L'esponente $s$ è positivo per $\alpha \lesssim 1.2$, indicando un comportamento da cristallo temporale. A differenza dell'analisi del tempo di deviazione, questo strumento non dipende dall'assunzione che il limite termodinamico sia di campo medio, rendendolo robusto per $\alpha > 1$.
   • Entrambi i metodi identificano una soglia critica per la scomparsa del cristallo temporale nell'intervallo $\alpha \in [1, 1.5]$, in accordo con i precedenti risultati delle espansioni dei cumulanti.

Significato e Rivendicazioni
Il documento afferma che la dinamica di Langevin semiclassica fornisce una sonda di dimensione finita utile per i cristalli temporali dissipativi in regimi in cui le simulazioni esatte di Lindblad sono computazionalmente intrattabili.

• Oltre il Campo Medio: L'approccio cattura con successo la robustezza dell'ordine del cristallo temporale in regimi in cui la teoria di campo medio fallisce (specificamente per $\alpha > 1$ nel modello spin-1/2). Dimostra che la fase del cristallo temporale si estende oltre il limite stretto di campo medio, a condizione che la dissipazione o l'interazione rimangano sufficientemente a lungo raggio.
• Strumenti Diagnostici: Il lavoro stabilisce gli esponenti di scaling di dimensione finita (tasso di decadimento $\beta$, esponente del tempo di deviazione $b$ e esponente del picco di Fourier $s$) come diagnostici diretti della robustezza della fase del cristallo temporale.
• Validazione Metodologica: I risultati per il modello spin-1 sono in linea con le precedenti scoperte delle espansioni dei cumulanti, convalidando l'approccio di Langevin semiclassico come uno strumento affidabile per studiare sistemi quantistici aperti con interazioni a lungo raggio.
• Intuizione Fisica: Lo studio evidenzia che, sebbene il limite termodinamico per $\alpha \leq 1$ sia descritto dalla dinamica di campo medio, lo scaling di dimensione finita verso tale limite è non banale e dipende da $\alpha$. Inoltre, lo scaling algebrico del tempo di deviazione (anziché logaritmico) suggerisce una dinamica regolare e non caotica in queste fasi di cristallo temporale.