Particle-preserving fermionic shadows with mode-independent sample complexity

Resumo Técnico: Sombras fermiônicas preservadoras de partículas com complexidade de amostragem independente de modos

Definição do Problema
O artigo aborda o desafio de aprender eficientemente valores de esperança de operadores específicos em relação a um estado fermiônico de $n$ modos e $\eta$ partículas desconhecido, utilizando o framework de sombras clássicas (classical shadows). O foco principal recai sobre duas classes de observáveis:

1. Sobreposições de determinantes de Slater: Estimar a sobreposição entre um estado desconhecido e um determinante de Slater arbitrário.
2. Observáveis quadráticos preservadores de partículas: Estimar os valores de esperança de operadores fermiônicos quadráticos gerais que conservam o número de partículas.

Embora trabalhos anteriores tenham estabelecido que a complexidade de amostragem de caso médio para sobreposições de Slater é constante, os melhores limites conhecidos para o pior caso eram de $O(\sqrt{n} \log n)$ para sombras que não preservam partículas (Ref. [6]) e $O(n)$ para esquemas preservadores de partículas (Ref. [5]). Os autores visam derivar limites de pior caso rigorosos que sejam independentes do número total de modos $n$, dependendo, em vez disso, do número de partículas $\eta$.

Metodologia
Os autores empregam o protocolo de sombras clássicas onde cópias do estado desconhecido $\rho$ são evoluídas sob unitárias aleatórias amostradas do grupo unitário $U(n)$ (especificamente, a representação da potência externa $U_\eta$) e medidas na base de ocupação (determinantes de Slater). A principal inovação técnica reside na análise rigorosa da variância dos estimadores resultantes.

1. Análise de Espaço Simétrico: A prova reduz o problema de limitar a variância extrema da sombra à análise harmônica no espaço simétrico do tipo AIII, $G/K = U(n)/(U(\eta) \times U(n-\eta))$.
2. Funções Esféricas Zonais: Ao alavancar a comutação do operador de segundo momento $\Gamma_R$ com o subgrupo $K$, os autores decompõem o espaço de Hilbert em camadas de Johnson (Johnson shells) (representações irredutíveis de $K$). Os autovalores de $\Gamma_R$ são expressos como integrais sobre o espaço simétrico envolvendo funções esféricas zonais.
3. Teoria de Matrizes Aleatórias e Polinômios Ortogonais: A avaliação do autovalor líder (que dita a complexidade de amostragem) é transformada em uma integral sobre o ensemble unitário de Jacobi. Os autores utilizam técnicas da teoria de polinômios ortogonais (especificamente polinômios de Jacobi deslocados) e matrizes de momentos para derivar expressões exatas de soma finita para esses autovalores.
4. Análise de Profundidade de Circuito: O artigo analisa o custo de implementação das unitárias aleatórias necessárias tanto em primeira quanto em segunda quantização, comparando as profundidades de circuito necessárias para alcançar 3-designs em cada regime.

Contribuições Principais e Resultados

• Limite de Pior Caso Estrito para Sobreposições de Slater:
  Os autores provam que a complexidade de amostragem para estimar a sobreposição com um determinante de Slater arbitrário é $O(\eta \log \eta)$. Isso representa uma melhoria significativa em relação ao limite anterior de pior caso de $O(\sqrt{n} \log n)$ (Ref. [6]) e à escala de $O(n)$ implicada por análises anteriores de preservação de partículas.

  • Teorema 1: Estabelece um limite superior na variância do estimador: $\text{Var}[\hat{\Pi}_\phi] \leq (\eta + 1)H_{\eta+1}$, onde $H_m$ é o $m$-ésimo número harmônico.
  • Teorema 2: Demonstra que este limite é estrito (assintoticamente) para estados com uma razão de preenchimento fixa, mostrando que a variância escala como $\sim \frac{c-1}{c}(\eta+1)H_{\eta+1}$ quando $n=c\eta$.
  • Custo Computacional: O pós-processamento clássico para computar uma estimativa de sombra de tiro único (single-shot) para uma sobreposição de Slater escala como $O(n\eta^2)$.

• Complexidade de Amostragem para Observáveis Quadráticos:
  Para observáveis fermiônicos quadráticos preservadores de partículas arbitrários $h$, os autores derivam um limite de complexidade de amostragem de $O(\eta \|h_0\|_2^2)$, onde $h_0$ é a componente traço-nula de $h$.

  • Teorema 4: Prova que $\text{Var}[\hat{h}] \leq (2\eta + 1)\|h_0\|_2^2$.
  • Teorema 5: Mostra que esta escala é ótima via um contraexemplo explícito.
  • Custo Computacional: O custo de pós-processamento para estes observáveis escala como $O(n^2\eta)$.
  • Distinção Técnica: Os autores observam que, enquanto o caso do determinante de Slater requer teoria avançada de espaços simétricos, o caso do observável quadrático depende de argumentos de teoria de representação mais simples, pois operadores quadráticos formam uma única representação irredutível sob a ação relevante, enquanto determinantes de Slater não o fazem.

• Implementação de Profundidade de Circuito:
  O artigo distingue a profundidade de circuitos necessária para implementar as unitárias aleatórias em primeira e segunda quantização:

  • Primeira Quantização: Usando designs unitários aproximados, a profundidade de circuito necessária é polilogarítmica em $n$, especificamente $O(\log \eta \log \log (\eta \log n))$.
  • Segunda Quantização: Para matchgates de vizinho próximo que preservam partículas, os autores argumentam (baseado em argumentos de cone de luz) que designs de profundidade sublinear são impossíveis, exigindo uma profundidade de $\Theta(n)$.

Significância e Alegações
O artigo afirma fornecer os primeiros limites de complexidade de amostragem de pior caso independentes de modo para sombras fermiônicas preservadoras de partículas. Ao alcançar a escala de $O(\eta \log \eta)$, o método torna-se significativamente mais eficiente do que abordagens anteriores quando o número de partículas $\eta$ é muito menor que a raiz quadrada do número de modos ($\eta \ll \sqrt{n}$).

Os autores enfatizam que sua contribuição técnica reside na avaliação exata das integrais de variância, que são tipicamente difíceis de computar na teoria de sombras clássicas devido ao envolvimento de momentos de terceira ordem (cubos tensoriais). Ao utilizar a teoria de funções esféricas zonais em espaços simétricos e polinômios ortogonais, eles contornam os limites frouxos que costumam ser derivados de estimativas simples de norma.

Finalmente, o trabalho destaca um trade-off na implementação: enquanto a primeira quantização oferece profundidades de circuito polilogarítmicas favoráveis para a randomização necessária, a segunda quantização (frequentemente preferida para simulações fermiônicas) exige atualmente circuitos de profundidade linear para garantias estatísticas equivalentes, sugerindo um potencial gargalo para implementações de hardware de curto prazo nesse tipo de codificação.