First-Quantized Relativistic Quantum Simulation with Periodic and Dirichlet Boundary Conditions

Resumo Técnico: Simulação Quântica Relativística de Primeira Quantização com Condições de Contorno Periódicas e de Dirichlet

Enunciado do Problema
A simulação quântica relativística em domínios finitos enfrenta um desafio fundamental: o operador de energia cinética relativística é uma função não linear do momento, exigindo a estimativa de momentos de ordem superior do momento (especificamente $\langle \hat{p}^2 \rangle$ e $\langle \hat{p}^4 \rangle$ no regime fracamente relativístico). Em codificações de grade de primeira quantização usando registradores de qubits, a operação de deslocamento natural é cíclica. Embora este deslocamento cíclico coincida perfeitamente com as Condições de Contorno Periódicas (PBC), ele introduz acoplamentos não físicos de "fechamento de ciclo" (wrap-around) para as Condições de Contorno de Dirichlet (DBC), onde a função de onda deve se anular nas fronteiras. Métodos existentes frequentemente tratam as condições de contorno como uma consideração secundária à discretização; no entanto, este trabalho argumenta que as condições de contorno são intrínsecas à definição do Hamiltoniano de domínio finito, exigindo operadores de momento de rede e estratégias de estimativa distintos para PBC e DBC.

Metodologia
Os autores desenvolvem um framework para simular o operador de partícula única relativística de energia positiva em um intervalo de 1D de comprimento $R$ usando um registrador de posição de $L$ qubits ($N=2^L$ pontos de grade).

1. Formulação do Hamiltoniano:

   • Expansão Fracamente Relativística: O operador cinético relativístico $\hat{T}_{rel}$ é expandido para quarta ordem em momento: $\hat{T}_{rel} \approx \frac{\hat{p}^2}{2m} - \frac{\hat{p}^4}{8m^3c^2}$.
   • Implementação PBC: O operador de momento de rede $\hat{P}_P$ é definido via um operador de translação cíclica unitário $\hat{A}$. Os momentos $\langle \hat{P}_P^2 \rangle$ e $\langle \hat{P}_P^4 \rangle$ são expressos inteiramente em termos das partes reais dos momentos de translação $\langle \hat{A} \rangle$ e $\langle \hat{A}^2 \rangle$.
   • Implementação DBC: O alvo é um operador de diferença finita de cadeia aberta. Os autores constroem este operador tomando o operador cíclico (avaliado com o espaçamento de grade DBC) e adicionando termos de correção locais de fronteira para cancelar o link de fechamento de ciclo não físico. O operador $\hat{P}_D$ é decomposto em uma parte cíclica mais termos que envolvem projetores de fronteira ($\hat{E}_0, \hat{E}_1, \hat{E}_2$) que dependem apenas de probabilidades de extremidade e coerências próximas à extremidade.

2. Estimadores Quânticos:

   • Momentos de Translação: Para PBC e para o bulk de DBC, as partes reais de $\langle \hat{A}^l \rangle$ ($l=1, 2$) são estimadas usando um circuito do tipo teste de Hadamard com uma primitiva de unidade controlada.
   • Sobreposições de Fronteira: Para DBC, os termos de correção são estimados usando uma "identidade de sobreposição de fronteira". Este método reconstrói coerências $\langle |f\rangle\langle g| + |g\rangle\langle f| \rangle$ usando apenas probabilidades da base computacional ($P_f, P_g$) e a probabilidade do estado de superposição $|s^+_{fg}\rangle = (|f\rangle + |g\rangle)/\sqrt{2}$. Isso evita a necessidade de estimação de fase complexa ou tomografia completa do estado.
   • Energia Potencial: Potenciais diagonais são estimados via amostragem direta da base de posição e pós-processamento clássico.

Principais Contribuições

1. Hamiltonianos Resolvidos por Fronteira: A formulação de Hamiltonianos de rede fracamente relativísticos que distinguem explicitamente entre PBC e DBC, garantindo que os operadores de momento sejam consistentes com as restrições físicas de contorno do domínio.
2. Decomposição de Estimadores: A derivação de estimadores de observáveis que reconstroem os momentos necessários de $\langle \hat{P}^2 \rangle$ e $\langle \hat{P}^4 \rangle$. Crucialmente, o estimador DBC é mostrado como uma decomposição do estimador cíclico mais um conjunto finito de termos locais de fronteira, em vez de uma primitiva de medição fundamentalmente diferente.
3. Framework de Validação: Uma estratégia de validação abrangente separando três fontes distintas de erro: discretização de grade finita, truncamento fracamente relativístico e erro estatístico de medição.

Resultados
O framework foi validado analiticamente e numericamente através de vários benchmarks:

• Partícula Livre PBC: Demonstrou que a relação de dispersão da rede converge para a energia relativística contínua conforme o tamanho da grade aumenta, com o estimador perturbativo capturando corretamente o erro de truncamento em relação à energia de raiz quadrada exata da rede.
• Poço Quadrado Infinito DBC: Verificou que o estimador cíclico corrigido pela fronteira reproduz os autovalores exatos do operador de diferença finita de cadeia aberta (e o espectro senoidal analítico) com concordância de precisão de máquina (discrepâncias relativas $< 10^{-11}$).
• Potenciais Suaves: Testado em potenciais de cosseno periódicos (PBC) e harmônicos (DBC). O fluxo de trabalho de energia total (cinética + potencial) mostrou-se visualmente indistinguível da avaliação direta de matriz e da reconstrução do estimador, com diferenças absolutas máximas na ordem de $10^{-7}$.
• Escalonamento Estatístico: Simulações de disparos (shots) finitos confirmaram que o erro quadrático médio (RMSE) dos estimadores de energia cinética escala como $O(M^{-1/2})$ com o número de disparos $M$, consistente com as expectativas estatísticas padrão. O estimador DBC exibiu um preefator ligeiramente maior devido às medições de fronteira adicionais, mas manteve o mesmo escalonamento assintótico.

Significado e Alegações
O artigo afirma estabelecer uma metodologia compacta e sistemática para simulação quântica relativística que respeita as condições de contorno como um componente fundamental da definição do Hamiltoniano. Os autores enfatizam que sua abordagem não requer primitivas de medição qualitativamente novas para a principal correção relativística; ela simplesmente estende o conjunto de observáveis medidos de $\langle \hat{A} \rangle$ para $\langle \hat{A} \rangle, \langle \hat{A}^2 \rangle$ e, para DBC, um pequeno conjunto de termos de sobreposição de fronteira.

O trabalho é explicitamente enquadrado para o regime fracamente relativístico, onde o suporte de momento do estado é pequeno comparado a $mc$. Os autores observam que a validade da simulação depende do estado físico satisfazer esta condição, não das identidades dos estimadores em si, que permanecem exatas para os operadores de rede definidos. O framework é apresentado como um passo prático para simulações relativísticas confiáveis de sistemas contínuos, com o potencial de ser estendido para potenciais e geometrias mais complexas.