Operator delocalization in disordered spin chains via exact MPO marginals

Título: Delocalização de Operadores em Cadeias de Spin Desordenadas via Marginais Exatos de MPO

1. O Problema

O artigo aborda a questão fundamental da termalização quântica e do espalhamento de informação (scrambling) em sistemas de muitos corpos desordenados. Em sistemas caóticos, os operadores locais evoluem para strings de operadores cada vez mais extensas, espalhando a informação por graus de liberdade não locais. No entanto, em sistemas desordenados, observa-se o fenômeno de Localização de Muitos Corpos (MBL - Many-Body Localization), onde a termalização é suprimida.

O desafio central é quantificar como os operadores se delocalizam e como a informação se espalha nesses regimes. Métricas tradicionais, como o Correlacionador de Ordem Fora do Tempo (OTOC), fornecem informações sobre a frente de onda do espalhamento, mas muitas vezes exigem amostragem estocástica ou não capturam a estrutura interna da complexidade do operador. Além disso, distinguir entre a localização de Anderson (sistemas não interagentes) e a MBL (sistemas interagentes) requer medidas sensíveis à dinâmica lenta e logarítmica característica da MBL.

2. Metodologia

Os autores propõem uma abordagem baseada em Redes de Tensor (Tensor Networks), especificamente utilizando o formalismo de Operadores Produto Matricial (MPO).

• Novas Métricas: Introduzem duas grandezas complementares para medir a complexidade do operador, derivadas da expansão do operador no tempo em uma base de strings de Pauli:
  1. Massa do Operador ($m$): O número médio de sítios onde o operador atua de forma não trivial (contagem de matrizes de Pauli que não são a identidade).
  2. Comprimento do Operador ($h$): A extensão espacial média do suporte do operador, definida pela posição do sítio mais à direita onde o operador atua não trivialmente.
• Formulação MPO: Em vez de realizar tomografia completa do operador (que seria exponencialmente cara), os autores mapeiam o operador $\hat{O}$ para um vetor de estado em um espaço de Hilbert duplo (representação de Choi-Liouville). Este vetor é representado como um MPO.
• Cálculo Exato e Determinístico: Demonstram que a massa e o comprimento podem ser calculados exatamente e de forma determinística dentro do framework MPO, sem a necessidade de amostragem estocástica. Eles derivam fórmulas recursivas para obter as distribuições de probabilidade completas dessas grandezas a partir das entropias de Rényi de segunda ordem do MPO.
• Modelo Estudado: Dinâmica unitária da cadeia de spin-1/2 XXZ desordenada com campos aleatórios no sítio ($h_j$) e interações ($\Delta$). O operador inicial é local ($\hat{\sigma}^z_1$).
• Simulações: Utilizam Diagonalização Exata (ED) para sistemas pequenos ($L=12$) e evolução temporal de MPO (tDMRG/TEBD) para sistemas maiores ($L$ até 32), permitindo tempos de evolução longos.

3. Principais Contribuições

• Definição e Cálulo Eficiente: Estabelecem o "comprimento do operador" como uma métrica crucial e fornecem um algoritmo eficiente para calculá-lo e calcular sua distribuição de probabilidade exata via MPO.
• Distinção entre Regimes: Demonstram uma diferença qualitativa e quantitativa clara entre o regime de localização de Anderson (não interagentes, $\Delta=0$) e o regime MBL (interagentes, $\Delta \neq 0$).
• Validação do Modelo $\ell$-bit: Confirmam numericamente que a dinâmica lenta observada é quantitativamente descrita por um modelo efetivo de $\ell$-bits (integrais de movimento localizadas), onde as interações decaem exponencialmente com a distância.
• Viabilidade Experimental: Propõem um protocolo experimental viável para medir essas grandezas em simuladores quânticos, utilizando o isomorfismo Choi-Jamiołkowski, evolução temporal controlada (frente e trás) e o método de "sombras clássicas" (classical shadows) com medições de pares de Bell.

4. Resultados Chave

• Caso Não Interagente (Localização de Anderson, $\Delta=0$):

  • A massa, o comprimento e a entropia de emaranhamento do operador saturam rapidamente para valores constantes.
  • Não há crescimento logarítmico; o operador permanece estritamente localizado, refletindo a ausência de espalhamento de informação.

• Caso Interagente (MBL, $\Delta \neq 0$):

  • Crescimento Logarítmico: Para qualquer interação fraca não nula, a massa, o comprimento e a entropia de emaranhamento do operador exibem um crescimento logarítmico no tempo ($h(t) \sim \ln t$).
  • Cones de Luz Logarítmicos: Este comportamento é consistente com o "cone de luz logarítmico" conhecido para a propagação de correlações quânticas e entropia de emaranhamento em MBL.
  • Robustez: O crescimento logarítmico persiste em diferentes tamanhos de sistema e é observado mesmo para interações muito fracas, indicando que a localização de Anderson é singular em $\Delta=0$; qualquer interação induz a dinâmica lenta da MBL.
  • Eficiência Computacional: O crescimento logarítmico da entropia de emaranhamento implica que a dimensão de ligação (bond dimension) necessária para representar o operador MPO cresce apenas linearmente com o tempo. Isso torna as simulações computacionalmente viáveis para tempos longos, ao contrário de sistemas caóticos onde o crescimento é exponencial.

• Análise de Médias Temporais:

  • Ao analisar médias temporais do comprimento e massa, observa-se um colapso de dados quando escalados por $\ln L$, confirmando a previsão teórica $h(t) \sim \xi(W) \ln t$.
  • A dependência da desordem mostra um comportamento não monótono, com um pico em desordens intermediárias, possivelmente relacionado à transição entre regimes integráveis e localizados.

5. Significado e Impacto

• Ferramenta Teórica: O trabalho fornece uma nova lente para estudar a dinâmica de operadores, complementando o OTOC e a entropia de emaranhamento. A capacidade de acessar a distribuição de probabilidade completa (não apenas a média) sem amostragem estocástica é um avanço metodológico significativo.
• Compreensão da MBL: Os resultados reforçam a visão de que a MBL é caracterizada por uma dinâmica extremamente lenta governada por interações efetivas exponencialmente fracas entre $\ell$-bits, distinguindo-a claramente da localização de Anderson.
• Aplicabilidade Experimental: Ao fornecer um protocolo baseado em sombras clássicas e medições locais, o artigo torna acessível a verificação experimental dessas propriedades teóricas em plataformas de simulação quântica atuais (átomos neutros, íons presos, qubits supercondutores).
• Eficiência Computacional: Demonstra que a estrutura de emaranhamento específica dos sistemas MBL permite simulações de longo prazo que seriam impossíveis em sistemas genéricos, validando o uso de MPOs como a ferramenta ideal para estudar a fase MBL.

Em resumo, o artigo estabelece o "comprimento do operador" como uma métrica robusta e calculável para caracterizar a delocalização e o scrambling em sistemas desordenados, confirmando a natureza logarítmica da dinâmica na fase MBL e oferecendo caminhos tanto para simulações numéricas eficientes quanto para verificações experimentais.