Efficient classical simulation of two-dimensional long-range systems: Rydberg arrays and beyond

Resumo Técnico: Simulação Clássica Eficiente de Sistemas de Longo Alcance Bidimensionais

Enunciado do Problema
O Método de Monte Carlo Variacional (VMC) é uma ferramenta primária para realizar o benchmarking de simuladores quânticos, particularmente para sistemas bidimensionais (2D) com interações de longo alcance (ex: arranjos de átomos de Rydberg com acoplamentos $1/r^3$ ou $1/r^6$). No entanto, existe um gargalo computacional significativo na avaliação da "energia local" ($H_{\text{loc}}$) para sistemas com interações de dois corpos de todos-para-todos arbitrários.

• No VMC padrão, enquanto os custos de amostragem são $O(N)$ para Estados de Rede de Tensores (TNS) e $O(N^2)$ para Estados Quânticos Neurais (NQS), a avaliação de $H_{\text{loc}}$ para acoplamentos não locais fora da diagonal de longo alcance requer a soma de $O(N^2)$ termos não locais.
• Como cada termo requer uma avaliação completa da rede levando $O(N)$ de tempo, o custo total por varredura de Monte Carlo (MC) escala como $O(N^3)$. Esse escalonamento cúbico torna a simulação clássica precisa de grandes sistemas 2D de longo alcance (ex: redes $10 \times 10$) intratável.
• Adicionalmente, uma patologia técnica surge ao evoluir TNS a partir de estados de produto (comum em experimentos de Rydberg): em estados de produto exatos, o Tensor de Geometria Quântica (QGT) desenvolve direções nulas se a dimensão de ligação for maior que um, fazendo com que as equações do Princípio Variacional Dependente do Tempo (TDVP) tornem-se singulares.

Metodologia
Os autores introduzem uma estrutura tVMC de rede de tensores de base adaptativa projetada para reduzir a complexidade computacional da avaliação da energia local de $O(N^3)$ para $O(N)$, mantendo a capacidade de lidar com acoplamentos desiguais gerais e geometrias realistas.

1. Amostragem de Base Adaptativa:

   • Em vez de amostrar em uma única base fixa (ex: a base $z$), o Hamiltoniano é decomposto em componentes que são diagonais em bases específicas.
   • Para o modelo XY dipolar ($H = H_{XY} + H_Z$), os autores dividem a evolução em $H_X$ (diagonal na base $x$) e $H_Y + H_Z$ (diagonal na base $y$).
   • Duas cadeias de Markov independentes são mantidas: uma amostrando de $|\langle x|\psi_\theta\rangle|^2$ e outra de $|\langle y|U|\psi_\theta\rangle|^2$, onde $U$ é o operador de mudança de base.
   • Como os acoplamentos de longo alcance fora da diagonal desaparecem dentro de cada base específica, a avaliação da energia local para cada componente torna-se $O(N)$.
   • As atualizações de parâmetros são computadas separadamente para cada base ($\dot{\theta}_X$ e $\dot{\theta}_Y$) e combinadas via $\dot{\theta} = \dot{\theta}_X + U^\dagger \dot{\theta}_Y$.

2. Tratamento de Condições Iniciais de Estado de Produto:

   • Para resolver a singularidade no QGT ao iniciar de um estado de produto (onde a dimensão de ligação $D > 1$), os autores empregam uma abordagem de tVMC Projetado (p-tVMC) para o passo de tempo inicial.
   • Em vez de resolver a equação TDVP singular, o primeiro passo ($\delta t$) é determinado pela maximização da fidelidade $|\langle \psi_\theta | V | \psi(0) \rangle|^2$, onde $V$ é a expansão de segunda ordem do operador de evolução temporal.
   • Esta otimização é eficiente porque o estado inicial é um estado de produto, tornando a distribuição de amostragem $|\langle s|V|\psi(0)\rangle|^2$ computacionalmente tratável. Após este passo inicial, o tVMC padrão é retomado.

3. Exploração de Simetria:

   • A estrutura impõe explicitamente a simetria de reflexão ao longo das diagonais da rede no ansatz TNS. Isso reduz o número de parâmetros variacionais independentes por um fator de quatro, aliviando significativamente os requisitos de memória da GPU e reduzindo a variância das atualizações de parâmetros.

Resultados Principais
O método foi aplicado para simular o protocolo de tempo real adiabático de um modelo XY dipolar $10 \times 10$ realizado em um experimento recente de Rydberg (Chen et al., Nature 2023), um tamanho de sistema anteriormente fora do alcance da simulação clássica.

• Benchmarking: Em um sistema $4 \times 4$, os resultados do TNS-tVMC (com dimensões de ligação $D=4, 5$) convergiram sistematicamente para os resultados de diagonalização exata, reproduzindo com precisão a dinâmica de relaxação induzida e oscilações rápidas.
• Simulação $10 \times 10$:
  • A simulação rastreou a evolução da energia, mostrando o sistema relaxando para um platô de quase estado estacionário ligeiramente acima da energia do estado fundamental.
  • Correlações Antiferromagnéticas (AFM): O estudo abordou uma discrepância no experimento original, onde o estado preparado dinamicamente mostrou correlações AFM de longo alcance suprimidas em comparação com cálculos de Estado de Matriz de Produto (MPS) do estado fundamental.
  • A simulação de dinâmica ideal dos autores (partindo de um estado de Néel perfeito com um detuning maior $\delta_0 = 19.48$) demonstrou que o rampa adiabática prescrita de fato produz correlações AFM espacialmente estendidas consistentes com o estado fundamental variacional.
  • Isso sugere que a supressão observada no experimento é provavelmente devida a imperfeições experimentais (como potência limitada do laser impedindo o uso do $\delta_0$ ideal na rede $10 \times 10$) em vez de uma falha intrínseca do protocolo adiabático.
• Convergência: Os resultados mostraram convergência robusta em relação ao tamanho da amostra ($N_s = 81.920$ vs. $163.840$), confirmando a confiabilidade da abordagem de base adaptativa.

Significado e Alegações
O artigo afirma fornecer uma ferramenta de benchmarking clássico escalável para plataformas quânticas de longo alcance, especificamente simuladores de Rydberg e armadilhas de íons, que atualmente necessitam urgentemente de tais ferramentas.

• Conquista Técnica: O trabalho resolve o gargalo $O(N^3)$, permitindo a simulação precisa de dinâmicas de tempo real e estados fundamentais para sistemas 2D de longo alcance com custo de $O(N)$ por varredura de MC.
• Resolução de Patologia: Oferece uma solução robusta para evoluir TNS a partir de estados de produto, um cenário comum em protocolos experimentais.
• Impacto Científico: Ao fornecer uma referência clássica direta para o experimento de Rydberg $10 \times 10$, o trabalho esclarece a natureza da discrepância entre os resultados experimentais e as previsões teóricas do estado fundamental. Demonstra que o protocolo adiabático ideal é capaz de gerar a ordem de longo alcance esperada, deslocando o foco da otimização experimental para as imperfeições de realização.
• Perspectiva Futura: Os autores posicionam este método como uma ponte para a colaboração clássico-quântica, onde simulações clássicas podem guiar a identificação de configurações físicas interessantes antes que o crescimento do emaranhamento se torne severo demais para redes de tensores, momento em que a dinâmica pode ser entregue aos simuladores quânticos. O método também é notado como um novo caminho para estudar sistemas 2D de longo alcance com problemas de sinal, onde outros métodos numéricos têm dificuldade.