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

Resumen Técnico: Asimetría de Entrelazamiento en la Cadena de Espín-1/2 XYZ con Brecha

Planteamiento del Problema
La asimetría de entrelazamiento, $\Delta S^{(n)}_A$, cuantifica el grado en que un subsistema $A$ rompe una simetría presente en el estado fundamental global. Si bien existen resultados analíticos para esta cantidad en teorías libres y, de forma perturbativa, para la cadena XXZ crítica, ha habido una falta de caracterización analítica para sistemas interactuantes en fases con brecha (gapped) con simetrías continuas explícitamente rotas. Específicamente, el comportamiento de la asimetría de entrelazamiento en la fase interactuante de la cadena XYZ de espín-1/2 con ruptura de $U(1)$ permanecía sin resolver. El desafío radica en determinar la constante aditiva en el crecimiento logarítmico de la asimetría, la cual codifica la ruptura de la simetría, para un modelo interactuante donde las técnicas estándar de fermiones libres fallan.

Metodología
Los autores abordan este problema combinando tres ingredientes teóricos y numéricos distintos para calcular las asimetrías de entrelazamiento de Rényi para un intervalo grande en la fase con brecha de la cadena XYZ:

1. Identidad de Momentos Cargados y Relación de Susceptibilidad: Los autores demuestran una identidad de momentos cargados para estados gaussianos fermiónicos y estados de producto de matrices (MPS) inyectivos con matrices de transferencia con brecha. Esta identidad relaciona la asimetría de entelazamiento con la susceptibilidad estática de la carga rota, $\chi_{zz}$. Específicamente, establecen que la varianza de la carga en la densidad reducida, $V$, es directamente proporcional a la varianza ordinaria de la carga, $\text{Var}(Q_A)$, permitiendo que la asimetría se exprese en términos de $\chi_{zz}$.
2. Reglas de Suma de Sine-Gordon: El límite de escala de la cadena XXZ perturbada por $\gamma$ se identifica con el modelo sine-Gordon. Los autores derivan una regla de suma de no conservación para evaluar $\chi_{zz}$ utilizando factores de forma del operador $\sin(\beta\phi)$. Esto implica el cálculo de las contribuciones de los canales de dos-kinks y un breather ($B_1$). Estas contribuciones proporcionan un límite inferior preciso sobre la amplitud universal $c(\beta^2) = \chi_{zz}/M$ (donde $M$ es la masa del kink).
3. Acotación de la Cola de Multi-partículas: Para contabilizar las contribuciones más allá de los sectores de dos-kinks y un breather, los autores emplean dos reglas de suma de momentos (un momento de Kubo estático y una regla de suma $f$). Estas reglas restringen la "cola" de las contribuciones de multi-partículas (empezando con el término de cuatro-kinks) a un corredor de dos lados en el lado atractivo del punto de fermiones libres, y proporcionan un límite inferior en el lado repulsivo.
4. Entrada de Red Exacta y Validación Numérica: La masa del kink y las constantes de acoplamiento son suministradas por la solución exacta de Baxter–Johnson–Krinker–McCoy (JKM). Estos insumos analíticos se combinan con simulaciones de grupo de transferencia de densidad de matriz de sistemas infinitos (iDMRG). Las simulaciones validan la fórmula maestra midiendo la susceptibilidad directamente y computando la asimetría mediante momentos cargados del MPS, confirmando las predicciones teóricas a través de varios acoplamientos e índices de Rényi.

Contribuciones Clave y Resultados

• Fórmula Maestra: Los autores derivan una fórmula maestra para la asimetría de entrelazamiento en la fase con brecha:
  $$\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)$$
  Esta fórmula vincula el comportamiento asintótico de la asimetría con la susceptibilidad estática $\chi_{zz}$ y la masa del kink $M$.
• Límites de la Amplitud Universal: El artículo proporciona un límite inferior riguroso para la amplitud universal $c(\beta^2)$ compuesta por las contribuciones de dos-kink y un breather ($c_{2k} + c_{B1}$).
• Estimación de la Cola: El déficit entre el límite inferior y la amplitud real se atribuye a las colas de multi-partículas. Los autores muestran que en el punto de fermiones libres ($\beta^2 = 1/2$), el término de la cola principal (cuatro-kinks) desaparece cuadráticamente, siguiendo una ley de tangente. Las estimaciones perturbativas y los límites numéricos confirman que este término explica el pequeño exceso observado en las simulaciones cerca del punto de fermiones libres.
• Confirmación Numérica: Las simulaciones iDMRG, extrapoladas al límite $\gamma \to 0$, reproducen la fórmula maestra. La susceptibilidad medida coincide con los valores derivados de las reglas de suma dentro de incertidumbres sistemáticas. Además, el cálculo directo de la asimetría a partir de las matrices de densidad de intervalo confirma la constante subdominante $c_1$ predicha.
• Comportamiento Cerca del Punto de Fermiones Libres: El análisis revela que la amplitud sigue una ley de tangente, $c(\beta^2) \approx 1 + (\frac{\pi}{2}-1)(\xi_{sG}-1)$, cerca del punto de fermiones libres, con la pendiente determinada por la cancelación entre el umbral de dos-kinks y la contribución del breather.

Significado
El artículo afirma proporcionar la primera caracterización analítica de la asimetría de entrelazamiento para un modelo interactuante con una simetría $U(1)$ explícitamente rota en equilibrio. Al fijar la amplitud a través de la fase —desde ambos lados del punto de fermiones libres para $\xi_{sG} < 1$ y acotándola para $\xi_{sG} > 1$—, el trabajo establece una comprensión integral de cómo la ruptura de simetría se manifiesta en las medidas de entrelazamiento para sistemas con brecha e interactuantes. Los resultados resaltan el papel del punto de fermiones libres como un caso marginal donde los canales de partículas de orden superior se desactivan, y ofrecen una base para comprender la dinámica post-quench y el comportamiento de tipo Mpemba en tales sistemas. El trabajo une técnicas de teoría de campos integrables (factores de forma de sine-Gordon) con simulaciones modernas de redes de tensores (iDMRG/iMPS) para resolver un problema no trivial en la física de muchos cuerpos cuánticos.