Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

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

Definição do Problema
Simular sistemas quânticos é computacionalmente exigente devido ao crescimento exponencial da complexidade com o tamanho do sistema, um desafio exacerbado quando efeitos relativísticos devem ser incorporados para descrever fenômenos físicos com precisão. Embora a teoria de perturbação seja uma abordagem padrão para estudar efeitos relativísticos além do framework não relativístico, implementar essas simulações em dispositivos quânticos de curto prazo permanece difícil. Os métodos existentes frequentemente dependem da segunda quantização, que é padrão para a química quântica, mas pode carecer de um tratamento intuitivo para sistemas com não linearidades. Há uma necessidade de um framework prático de primeira quantização capaz de estimar energias de estado fundamental relativísticas para sistemas sob diferentes condições de contorno em hardware de curto prazo.

Metodologia
Os autores propõem um framework de primeira quantização para simulação quântica relativística que discretiza a função de onda em uma grade finita representada por qubits do sistema. A metodologia central envolve os seguintes componentes:

1. Discretização e Operadores de Translação:
   A função de onda contínua 1D é mapeada para uma grade discreta de $2^L$ pontos usando $L$ qubits, onde a posição $x$ é escalonada para o domínio $0 \le x < 1$. O operador de momento ao quadrado $\hat{p}^2$ é aproximado pelo método de diferença finita baseado em operadores de translação quântica. O operador de translação $\hat{A}_x$ (adicionador quântico) desloca o estado de posição $|j\rangle$ para $|j+\delta x\rangle$, enquanto seu conjugado $\hat{A}_x^\dagger$ (subtrator quântico) desloca na direção oposta. O momento ao quadrado é expresso como $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$, onde $L_m = \lambda_m / \delta x$ é um parâmetro adimensional envolvendo o comprimento de onda de Compton reduzido.

2. Expansão Perturbativa da Energia Cinética:
   O Hamiltoniano de energia cinética relativística $\hat{H}_{re}^K$ é derivado da forma de raiz quadrada da relação energia-momento relativística. Ele é expandido em uma série perturbativa de termos de momento de ordem superior ($\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$) válida no regime onde $\langle \hat{p}^2 \rangle \ll m^2 c^2$. A energia cinética relativística total é aproximada pela truncagem desta série, incorporando coeficientes derivados da equação de Klein-Gordon.

3. Tratamento de Condições de Contorno:
   O framework aborda duas condições de contorno distintas:

   • Condições de Contorno Periódicas (PBC): Os operadores de translação preservam a unitariedade ao envolver a grade ($\hat{A}_x |2^L-1\rangle = |0\rangle$). Os valores esperados de potências pares de momento são expressos como combinações de $\langle \hat{A}^{(l)} \rangle$.
   • Condições de Contorno de Dirichlet (DBC): Aqui, a função de onda não se envolve ($\psi(0) \neq \psi(1)$). Os autores introduzem um operador de correção $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ para contabilizar a fronteira. Os valores esperados para DBC são formulados como os termos de PBC mais termos adicionais envolvendo $\hat{E}_0$ e operadores relacionados ($\hat{E}_1, \hat{E}_2$), que dependem das amplitudes de probabilidade nas fronteiras da grade.

4. Algoritmo de Simulação Quântica:
   Os autores projetam circuitos quânticos para estimar os valores esperados necessários:

   • Energia Cinética: Portas de translação controladas são aplicadas iterativamente a um qubit de controle. Medir o qubit de controle na base Z resulta na parte real de $\langle (\hat{A}_x^\dagger)^l \rangle$, o que corresponde aos termos de momento necessários.
   • Energia Potencial: Uma matriz densidade $\hat{\rho}_V$ representando o potencial é usada em um circuito para estimar $\langle \hat{V} \rangle$ através das probabilidades do estado do sistema.
   • Otimização Variacional: A energia total é estimada minimizando a soma dos valores esperados de energia cinética e potencial sobre um conjunto de parâmetros variacionais $\{\bar{\kappa}\}$ em um estado ansatz $|\psi(\{\bar{\kappa}\}, x)\rangle$.

Principais Contribuições

• Receita de Primeira Quantização: O artigo fornece uma receita específica para simulação relativística usando primeira quantização, oferecendo uma alternativa aos métodos de segunda quantização amplamente utilizados, particularmente benéfica para sistemas com não linearidades.
• Framework de Fronteira Unificado: Apresenta uma abordagem unificada para lidar com condições de contorno Periódicas e de Dirichlet. Crucialmente, demonstra que os dados obtidos para PBC podem ser reutilizados para estimar energias de DBC ao adicionar termos de correção de fronteira específicos, reduzindo o overhead computacional.
• Compatibilidade com Dispositivos de Curto Prazo: O método utiliza primitivas de circuito quântico padrão (adicionadores controlados, portas Hadamard e medições) e otimização variacional, tornando-o adequado para implementação em dispositivos quânticos de curto prazo.
• Designs de Circuitos Explícitos: Os autores fornecem designs detalhados de circuitos para computar os valores esperados dos operadores de translação e os termos de correção de fronteira necessários para a expansão perturbativa.

Resultos e Alegações
O artigo afirma oferecer uma "rota prática" para simular efeitos relativísticos em dispositivos quânticos de curto prazo. Os resultados teóricos mostram que:

• Os termos de energia cinética relativística até uma ordem específica de perturbação podem ser decompostos em combinações lineares de valores esperados de operadores de translação ($\langle \hat{A}^{(l)} \rangle$) e operadores de fronteira ($\langle \hat{E}_d \rangle$).
• Os coeficientes para essas expansões podem ser pré-calculados com base em parâmetros físicos ($m, c, \hbar, \delta x$).
• O método permite a otimização variacional de estados ansatz para estimar energias de estado fundamental tanto não relativísticas quanto relativísticas.

Significância
Os autores posicionam este trabalho como um passo fundamental para apoiar desenvolvimentos futuros em física e química quântica. Ao permitir a simulação de efeitos relativísticos em hardware de curto prazo, o framework aborda uma lacuna nas capacidades atuais de simulação quântica. O artigo enfatiza que, embora aumentar o número de qubits melhore a resolução da grade, deve-se ter cuidado para garantir que o parâmetro adimensional $L_m$ permaneça pequeno o suficiente para que a expansão perturbativa continue válida. O trabalho não pretende resolver a equação de Dirac completa ou lidar com limites ultra-relativísticos, mas foca especificamente no regime onde as correções relativísticas ao framework não relativístico são relevantes. Sugere-se que trabalhos futuros explorem operadores Laplacianos de ordem superior e métodos de construção de ansatz mais precisos.