Symmetry-Protected Basin Localization in Variational Quantum Eigensolvers

Resumo Técnico: Localização de Bacia Protegida por Simetria em Solucionadores Variacionais de Autovalores Quânticos

Enunciado do Problema
Os Solucionadores Variacionais de Autovalores Quânticos (VQEs) enfrentam um obstáculo fundamental em regimes moleculares fortemente correlacionados: a paisagem de energia não convexa é particionada pela simetria molecular em bacias de atração concorrentes. A principal modalidade de falha não ocorre durante o loop de otimização, mas na inicialização.

1. Inicialização Aleatória: Aproxima um conjunto unitário do tipo Haar, levando a "planícies estéreis" onde a variância do gradiente desaparece exponencialmente ($O(1/2^N)$) devido à concentração de medida, tornando a descida local não informativa.
2. Inicialização de Campo Médio (Hartree–Fock): Falha em regimes fortemente correlacionados (além do ponto de Coulson–Fischer) porque quebra espontaneamente a simetria de spin. Isso viés o estado variacional para variedades de simetria quebrada com sobreposição negligenciável com o setor do estado fundamental correlacionado, desviando a otimização para um mínimo local concorrente.

Em ambos os casos, a medida de inicialização é deslocada da bacia do estado fundamental antes que a primeira etapa de gradiente seja tomada, fazendo com que o algoritmo falhe antes que a otimização comece.

Metodologia
Os autores introduzem um pré-condicionador condicionado à geometria, denotado $P_{eq}: \mathbb{R} \mapsto \theta_0$, que mapeia a geometria nuclear diretamente para os parâmetros do circuito. Este mapa é restrito pela covariância $SE(3)$ do Hamiltoniano molecular.

• Estrutura Teórica: O pré-condicionador é definido como uma restrição holonômica que comuta com o grupo de simetria física do Hamiltoniano. Ao exigir $| \psi(P_{eq}(g \cdot R)) \rangle \simeq \hat{U}(g) | \psi(P_{eq}(R)) \rangle$ para todo $g \in SE(3)$, o mapa remove graus de liberdade redundantes de rotação e translação. Isso restringe a medida de inicialização induzida ao espaço quociente $\mathbb{R}^{3M}/SE(3)$ e força o estado inicial $\theta_0$ para o componente consistente com a simetria da bacia do estado fundamental.
• Implementação: $P_{eq}$ é realizado usando uma rede de passagem de mensagens equivariante a $SE(3)$ (tipo MACE), construída a partir de princípios de acoplamento de tensores usados em modelos atômicos. A arquitetura da rede respeita as restrições de localidade do circuito quântico (correspondendo o campo receptivo ao cone causal do circuito).
• Treinamento: A rede é treinada em benchmarks de vetores de estado (obtidos via busca por saltos em bacias e refinamento L-BFGS-B) para prever representantes de parâmetros atrativos à bacia. Uma função objetivo consciente de calibre é usada para lidar com redundâncias de parâmetros, combinando ancoragem de parâmetros com termos de fidelidade de estado.
• Protocolo Híbrido: Para sistemas desordenados, os autores propõem um método híbrido combinando direcionamento determinístico de bacia equivariante com fuga estocástica para resolver armadilhas em escala fina dentro da bacia localizada.

Principais Contribuições

1. Mudança do Controle de Concentração para Controle de Curvatura: O artigo estabelece uma mudança teórica nas estatísticas de gradiente. VQEs padrão sofrem com gradientes controlados por concentração (variância suprimida pelo volume do espaço de Hilbert). A localização de bacia via $P_{eq}$ muda o regime para gradientes controlados por curvatura, onde a variância é determinada pela curvatura local do Hessiano ($\kappa_j^2 \sigma_j$) dentro de uma vizinhança estritamente convexa, independente da dimensão total do espaço de Hilbert.
2. Seleção de Bacia Pré-Tiro: O procedimento realiza a identificação da bacia antes do loop quântico limitado por tiros. O mapa clássico seleciona a bacia compatível com a simetria, e o circuito quântico é reservado exclusivamente para refinar a correlação dentro dessa bacia selecionada.
3. Derivação Formal: Os autores derivam a Proposição 1 (e o Teorema S1 no Material Suplementar), provando formalmente que, se a medida de inicialização for suportada dentro de uma bacia quadrática estritamente convexa, a variância do gradiente escala com a curvatura local em vez de exponencialmente com o tamanho do sistema.

Resultados
O método foi avaliado em seis benchmarks moleculares esticados ($N_2$, $H_2O$, $BeH_2$, $CO$, $LiH$, $H_8$) e cadeias desordenadas de $H_{10}$ usando simulações de vetores de estado:

• Redução do Erro de Inicialização: $P_{eq}$ reduziu os erros de inicialização de Hartree–Fock por fatores de 38$\times$ a 6250$\times$.
• Precisão:
  • Alcançou precisão de inicialização sub-milihartree (sub-mHa) para $CO$, $LiH$e$H_8$.
  • Colocou $N_2$, $H_2O$ e $BeH_2$ dentro da bacia correlacionada na escala de mHa.
  • Em $N_2$ esticado ($R=2.5$ Å), a inicialização de Hartree–Fock caiu em uma bacia de simetria quebrada com $\Delta E > 0.2$ Ha, enquanto $P_{eq}$ colocou o estado na bacia do estado fundamental correlacionado.
• Robustez à Desordem: Em cadeias desordenadas de $H_{10}$, o protocolo híbrido (direcionamento equivariante + fuga estocástica) alcançou probabilidade de sucesso unitária com um orçamento de otimização fixo, enquanto a inicialização aleatória sofreu de distribuições de erro de cauda pesada e identificação frequente de bacias equivocada.
• Transferibilidade: O método demonstrou transferibilidade para sistemas poliatômicos (ex.: estiramento O–H em $H_2O$), mantendo precisão na escala de mHa em geometrias onde a inicialização aleatória falhou por ordens de grandeza.

Significado e Alegações
O artigo alega que a principal barreira ao VQE em sistemas fortemente correlacionados é a seleção de bacia, não a ineficiência do otimizador. Ao impor equivariância $SE(3)$ no mapa de inicialização, os autores demonstram que:

• A descoberta de bacia deslocalizada pode ser substituída por descida local controlada por curvatura.
• O fenômeno da "planície estéril" pode ser visto como uma propriedade de medidas de inicialização deslocalizadas; restringir a medida a um suporte consistente com a simetria e localizado na bacia contorna a supressão exponencial dos gradientes.
• A abordagem desacopla o custo da busca geométrica (tratada classicamente e de forma barata) do custo do refinamento de correlação quântica. Isso desloca o ônus computacional do loop quântico limitado por tiros para um mapa clássico protegido por simetria, oferecendo uma aceleração de tempo real de $10^4$–$10^6$ para a fase de inicialização.

Os autores permanecem modestos quanto ao escopo, observando que a construção é específica para Hamiltonianos de estrutura eletrônica com simetria $SE(3)$ e ansatz que preservam a localidade, e que a extensão para sistemas periódicos ou Hamiltonianos com acoplamento spin-órbita requer a adaptação dos mapas equivariantes aos grupos de simetria relevantes.