Entanglement asymmetry in the gapped XYZ spin-$\frac12$ chain

Sintesi Tecnica: Asimmetria dell'Entanglement nella catena XYZ di spin-1/2 con gap

Enunciato del Problema
L'asimmetria dell'entanglement, $\Delta S^{(n)}_A$, quantifica il grado in cui un sottosistema $A$ rompe una simmetria presente nello stato fondamentale globale. Sebbene esistano risultati analitici per questa quantità per le teorie libere e, in modo perturbativo, per la catena critica XXZ, è mancata una caratterizzazione analitica per sistemi interagenti in fasi con gap che presentano simmetrie continue esplicitamente rotte. Nello specifico, il comportamento dell'asimmetria dell'entanglement nella fase gapped, con rottura della simmetria $U(1)$, della catena XYZ di spin-1/2 interagente rimaneva irrisolto. La sfida risiede nel determinare la costante additiva nella crescita logaritmica dell'asimmetria, che codifica la rottura della simmetria, per un modello interagente in cui le tecniche standard dei fermioni liberi falliscono.

Metodologia
Gli autori affrontano questo problema combinando tre distinti ingredienti teorici e numerici per calcolare le asimmetrie di Rényi dell'entanglement per un intervallo ampio nella fase gapped della catena XYZ:

1. Identità del Momento Caricato e Relazione di Suscettibilità: Gli autori dimostrano un'identità del momento caricato per stati gaussiani fermionici e stati a prodotto di matrici (MPS) iniettivi con matrici di trasferimento con gap. Questa identità mette in relazione l'asimmetria dell'entanglement con la suscettibilità statica della carica rotta, $\chi_{zz}$. Nello specifico, stabiliscono che la varianza della carica nella densità ridotta, $V$, è direttamente proporzionale alla varianza della carica ordinaria, $\text{Var}(Q_A)$, permettendo all'asimmetria di essere espressa in termini di $\chi_{zz}$.
2. Somme di Regole di Sine-Gordon: Il limite di scala della catena XXZ perturbata da $\gamma$ viene identificato con il modello sine-Gordon. Gli autori derivano una regola di somma di non conservazione per valutare $\chi_{zz}$ utilizzando i form factor dell'operatore $\sin(\beta\phi)$. Questo comporta il calcolo dei contributi dai canali a due kink e dal singolo breather ($B_1$). Questi contributi forniscono un limite inferiore preciso sull'ampiezza universale $c(\beta^2) = \chi_{zz}/M$ (dove $M$ è la massa del kink).
3. Vincolo della Coda Multi-Particella: Per tenere conto dei contributi oltre i settori a due-kink e un-breather, gli autori impiegano due regole di somma di momenti (un momento di Kubo statico e una regola di somma $f$). Queste regole vincolano la "coda" dei contributi multi-particella (partendo dal termine a quattro-kink) a un corridoio a due lati sul lato attrattivo del punto dei fermioni liberi e forniscono un limite inferiore sul lato repulsivo.
4. Input di Reticolo Esatto e Validazione Numerica: La massa del kink e le costanti di accoppiamento sono fornite dalla soluzione esatta Baxter–Johnson–Krinsky–McCoy (JKM). Questi input analitici vengono combinati con simulazioni di densità-matrix renormalization group a sistema infinito (iDMRG). Le simulazioni validano la formula maestra misurando direttamente la suscettibilità e calcolando l'asimmetria tramite i momenti caricati dell'MPS, confermando le previsioni teoriche attraverso vari accoppiamenti e indici di Rényi.

Contributi Chiave e Risultati

• Formula Maestra: Gli autori derivano una formula maestra per l'asimmetria dell'entanglement nella fase gapped:
  $$\Delta S^{(n)}_A (\ell) = \frac{1}{2} \log (\pi \ell \chi_{zz}) + \frac{\log n}{2(n-1)} + O\left(\frac{1}{M\ell}\right)$$
  Questa formula lega il comportamento asintotico dell'asimmetria alla suscettibilità statica $\chi_{zz}$ e alla massa del kink $M$.
• Limiti dell'Ampiezza Universale: Il documento fornisce un limite inferiore rigoroso per l'ampiezza universale $c(\beta^2)$ composto dai contributi a due-kink e un-breather ($c_{2k} + c_{B1}$).
• Stima della Coda: Il deficit tra il limite inferiore e l'ampiezza reale è attribuito alle code multi-particella. Gli autori mostrano che al punto dei fermioni liberi ($\beta^2 = 1/2$), il termine di coda principale (quattro-kink) svanisce quadraticamente, seguendo una legge tangente; le stime perturbative e i vincoli numerici confermano che questo termine tiene conto del piccolo eccesso osservato nelle simulazioni vicino al punto dei fermioni liberi.
• Conferma Numerica: Le simulazioni iDMRG, estrapolate al limite $\gamma \to 0$, riproducono la formula maestra. La suscettibilità misurata corrisponde ai valori derivati dalle regole di somma entro incertezze sistematiche. Inoltre, il calcolo diretto dell'asimmetria dai densità-matrix degli intervalli conferma la costante sub-leading $c_1$ prevista.
• Comportamento Vicino al Punto dei Fermioni Liberi: L'analisi rivela che l'ampiezza segue una legge tangente, $c(\beta^2) \approx 1 + (\frac{\pi}{2}-1)(\xi_{sG}-1)$, vicino al punto dei fermioni liberi, con la pendenza determinata dal cancellamento tra la soglia del due-kink e il contributo del breather.

Significatività
Il lavoro sostiene di fornire la prima caratterizzazione analitica dell'asimmetria dell'entanglement per un modello interagente con una simmetria $U(1)$ esplicitamente rotta all'equilibrio. Definendo l'ampiezza attraverso l'intera fase — da entrambi i lati del punto dei fermioni liberi per $\xi_{sG} < 1$ e ponendo un limite per $\xi_{sG} > 1$ — il lavoro stabilisce una comprensione completa di come la rottura della simmetria si manifesti nelle misure di entanglement per sistemi interagenti con gap. I risultati evidenziano il ruolo del punto dei fermioni liberi come caso marginale in cui i canali di particelle di ordine superiore si spengono, e offrono una base per comprendere la dinamica post-quench e il comportamento di tipo Mpemba in tali sistemi. Il lavoro unisce le tecniche della teoria dei campi integrabili (form factor di sine-Gordon) con le moderne simulazioni di reti tensoriali (iDMRG/iMPS) per risolvere un problema non banale della fisica della materia quantistica a molti corpi.