A near-term quantum simulation of the transverse field Ising model hints at Glassy Dynamics

Resumo Técnico: Uma Simulação Quântica de Curto Prazo do Modelo de Ising de Campo Transverso Sugere Dinâmicas Vítreas

Declaração do Problema
O artigo aborda o desafio computacional de simular sistemas de vidro de spin em larga escala, que são críticos para compreender fases exóticas da matéria e transições vítreas na física de muitos corpos quânticos. Abordagens computacionais clássicas, como dinâmica molecular e simulações de Monte Carlo, enfrentam dificuldades com o crescimento exponencial das funções de onda eletrônicas em sistemas de muitos elétrons e com a alta complexidade computacional necessária para modelar sistemas desordenados ao longo de escalas de tempo estendidas. Especificamente, sistemas vítreos exibem paisagens de energia acidentadas com inúmeros mínimos locais e tempos de relaxação extremamente longos, tornando-os difíceis de simular com precisão usando computação de alto desempenho (HPC) clássica. Os autores visam explorar se técnicas de simulação quântica digital (DQS), especificamente Algoritmos de Eigensolver Quântico Variacional (VQE), podem investigar efetivamente essas configurações de spin desordenadas em hardware quântico de curto prazo.

Metodologia
Os autores empregaram uma abordagem híbrida quântico-clássica usando o algoritmo Variational Quantum Eigensolver (VQE) para simular o Modelo de Ising de Campo Transverso (TFIM) com a adição de um campo longitudinal. O estudo focou em duas geometrias específicas:

1. Cadeia de Ising 1D: Uma cadeia de comprimento $n=25$ (25 qubits).
2. Rede de Ising 2D: Uma rede quadrada de $6 \times 6$ (36 qubits).

Componentes Técnicos Principais:
*Construção do Hamiltoniano: O Hamiltoniano do sistema inclui interações de vizinhos mais próximos ao longo do eixo z ($J \sum \sigma_i^z \sigma_{i+1}^z$), um campo transverso na direção x ($-B_x \sum \sigma_i^x$) e um campo longitudinal na direção z ($-B_z \sum \sigma_i^z$). Interações diagonais foram excluídas.
*Seleção de Ansatz: Os autores utilizaram um ansatz "TwoLocal" composto por blocos de rotação $R_x$ e $R_y$. Crucialmente, este ansatz foi projetado sem emaranhamento entre os qubits. Esta escolha foi feita para reduzir a profundidade do circuito e mitigar a suscetibilidade ao ruído em dispositivos de curto prazo, enquanto ainda captura as flutuações quânticas necessárias para a estrutura do problema.
*Otimização: O otimizador clássico COBYLA foi utilizado para minimizar o valor de expectativa do Hamiltoniano, com um máximo de 25.000 iterações (embora a maioria das execuções tenha convergido antes).
*Implementação de Software: Os autores desenvolveram um pacote Python customizado, Qiskit Glassy Dynamics (posteriormente referido como QuGlassyIsing), para construir modularmente os Hamiltonianos para redes 1D e 2D. Este pacote lida com a geração de termos de interação ($\sigma^z\sigma^z$) e termos de campo ($\sigma^x, \sigma^z$) de forma eficiente.
*Ambiente de Simulação: As simulações foram conduzidas principalmente utilizando o backend Qiskit Aer GPU. Uma demonstração de prova de conceito também foi realizada no dispositivo real de 7 qubits IBM Q Oslo usando um modelo reduzido de 7 spins.

Resultados Principais
*Paisagens de Energia: As simulações revelaram que o aumento do campo longitudinal ($B_z$) reduz a energia do sistema ao alinhar os spins, enquanto o aumento do campo transverso ($B_x$) induz flutuações quânticas e inversões de spin sem reduzir a energia de forma tão significativa.
*Fases Desordenadas: O estudo identificou um regime específico onde emergem configurações de spin desordenadas. Esta "fase desordenada" é caracterizada por uma mistura de domínios ferromagnéticos e antiferromagnéticos.
* Na ausência de um campo longitudinal ($B_z=0$), o sistema exibiu um ordenamento antiferromagnético perfeito.
* Um estado desordenado foi observado apenas para combinações específicas de $B_x$ e $B_z$ (ex: $B_z=1, B_x=2$), onde a competição entre os campos interrompe o ordenamento uniforme, criando momentos locais de ferromagnetismo dentro de um fundo antiferromagnético.
* Os autores observam que as três configurações de spin mais prováveis nas simulações 2D eram quase idênticas, sugerindo um ordenamento coerente e padronizado mesmo dentro da fase desordenada.
*Demonstração de Hardware: A execução no dispositivo IBM Q Oslo (7 qubits) mostrou um perfil de convergência de energia com flutuações significativas e um valor de energia final ($-4.0139$) superior ao mínimo teórico. Este resultado é atribuído a erros de porta, erros de leitura e decoerência de qubit inerentes aos dispositivos NISQ (Noisy Intermediate-Scale Quantum).

Significância e Alegações
O artigo afirma fornecer uma demonstração de "prova de princípio" de que simuladores quânticos de curto prazo podem detectar fases desordenadas em modelos de Ising que são difíceis de acessar via métodos clássicos. Os autores argumentam que suas descobertas:
*Validam a Utilidade do VQE: Eles demonstram que o VQE, mesmo com ansatzes não emaranhados, pode modelar funções de onda complexas em tempo polinomial, reduzindo os custos computacionais em comparação com métodos clássicos.
*Elucidam Dinâmicas Vítreas: Os resultados sugerem que a interação entre campos transversos e longitudinais cria uma paisagem de energia "acidentada" com fases metaestáveis, espelhando o comportamento de vidros de spin e processos catalíticos onde o desordem leva a mecanismos de reação ineficientes.
*Desenvolvimento de Ferramentas: A introdução do pacote QuGlassyIsing oferece uma interface modular e amigável para a construção de Hamiltonianos para dinâmicas vítreas, facilitando pesquisas futuras em sistemas quânticos complexos.
*Perspectiva Futura: Os autores declaram modestamente que, embora as limitações atuais do NISQ (tempo de coerência, fidelidade de porta) impeçam a escalabilidade para sistemas maiores e mais complexos, seu trabalho estabelece as bases para aplicações futuras de computação quântica tolerante a falhas (FTQC). Eles propõem que, à medida que a correção de erros quânticos melhorar, esses métodos poderão desbloquear simulações de comportamento vítreo em sistemas 3D e geometrias mais complexas, potencialmente informando o design de novos materiais de armazenamento magnético e sistemas catalíticos.

O artigo conclui que, embora a implementação atual seja limitada por restrições de hardware, ela destaca com sucesso o potencial das ferramentas de processamento de informação quântica para trazer uma compreensão matizada da dinâmica e da estrutura de sistemas vítreos.