Pauli spectrum and nonstabilizerness of random fermionic Gaussian states

Resumen Técnico: Espectro de Pauli y no-estabilizerness de estados gaussianos fermiónicos aleatorios

Planteamiento del Problema
El artículo aborda la caracterización de la "magia" (no-estabilizerness) en estados gaussianos fermiónicos aleatorios, una clase de estados cuánticos que incluye los autoestados de Hamiltonianos de Sachdev–Ye–Kitaev (SYK2) cuadráticos y el resultado de circuitos de matchgate profundos. Aunque estos estados son clásicamente simulables, estudios recientes sugieren que suelen estar cerca de ser máximamente mágicos desde una perspectiva de Clifford; sin embargo, la estructura detallada de su magia —específicamente el espectro de Pauli completo, el comportamiento de las entropías de Rényi de estabilizadores (SREs) a través de diferentes índices de Rényi $q$, el papel de la simetría de número de partículas y la fiabilidad de las estimaciones basadas en muestreo— permanecía elusiva. Los autores pretenden derivar expresiones exactas de tamaño finito para las puridades de estabilizadores, reconstruir el espectro de Pauli exacto y analizar el comportamiento termodinámico de la densidad de la magia, particularmente en regímenes donde los métodos de muestreo estándar pueden fallar.

Metodología
Los autores analizan tres conjuntos aleatorios distintos de estados gaussianos fermiónicos puros:

1. Conjunto de Majorana Haar: Estados con pares permitidos con matrices de covarianza $\Gamma = OJO^T$, donde $O$ es Haar-distribuido en $SO(2N)$.
2. Conjunto de Slater Complejo: Estados que conservan el número con un llenado fijo $\nu = r/N$, correspondientes a rotaciones unitarias Haar-aleatorias en el Grassmanniano complejo $Gr_{\mathbb{C}}(r, N)$.
3. Conjunto de Slater Real: Estados que conservan el número con saltos reales, correspondientes a rotaciones ortogonales Haar-aleatorias en $Gr_{\mathbb{R}}(r, N)$.

El análisis se basa en las siguientes herramientas teóricas:

• Teorema de Wick y Pfaffianos: Expresar los valores esperados de Pauli al cuadrado $x_P^2$ como determinantes de submatrices de covarianza ($\det \Gamma_A$).
• Integrales de Jacobi y Selberg: Utilizar la invarianza rotacional de los conjuntos para mapear los momentos del espectro de Pauli a integrales de Selberg sobre conjuntos de Jacobi.
• Teoría de Réplicas y Dualidad de Skew Howe: Para los conjuntos que conservan el número, expresar los momentos enteros como integrales de coeficientes sobre $SU(2q)$, reduciendo el problema a integrales de grupo independientes del tamaño del sistema $N$.
• Teoría de Grandes Desviaciones: Analizar la competencia entre la entropía combinatoria de los sectores de peso de Majorana y el costo energético de los valores esperados de Pauli para determinar las contribuciones dominantes a las SREs.
• Análisis de Muestreo: Comparar promedios de tipo annealed (promediar la pureza antes de tomar el logaritmo) con estimaciones directas de muestreo de la característica para identificar sesgos en regímenes de muestras finitas.

Contribuciones Clave y Resultados

1. Momentos Exactos y Espectro de Pauli (Conjunto de Majorana):

   • Los autores derivan expresiones cerradas de tamaño finito $N$ para las puridades de estabilizadores promediadas por el conjunto $\mathbb{E}[\zeta_q]$ como sumas finitas de productos de Selberg.
   • Reconstruyen el espectro de Pauli promediado por el conjunto como una mezcla de productos de variables aleatorias beta-distribuidas independientes, resuelto por peso de Majorana.
   • El espectro se organiza en una jerarquía: un "bulk" exponencialmente grande de cadenas de alto peso que es asintóticamente lognormal, y colas de bajo peso raras donde los valores esperados se suprimen solo polinómicamente.

2. Transición de Congelación en $q_c = 2$:

   • Ocurre una transición de fase en el índice de Rényi $q_c = 2$.
   • Para $q \le 2$, la densidad de la magia está gobernada por los sectores centrales (de alto peso).
   • Para $q > 2$, la densidad de la magia está dominada por los sectores de frontera (cadenas de bajo peso, específicamente cadenas de Majorana de peso dos y sus complementos).
   • Consecuentemente, la densidad de la magia annealed $D_q = \lim_{N\to\infty} M_q/N$ se "congela" en $1/(q-1)$ para $q > 2$. Esto persiste incluso después de filtrar los operadores de identidad y paridad, distinguiendo estos estados de los estados típicos de Haar donde el filtrado restaura la densidad unitaria.

3. Fallo del Muestreo Directo para $q > 2$:

   • El artículo demuestra que el muestreo de la característica directa (dibujar cadenas de Pauli con una probabilidad proporcional a $x_P^2$) falla al estimar la SRE filtrada para $q > 2$ con recursos polinomiales.
   • Aunque los sectores de peso dos dominan la pureza promediada, su probabilidad combinada bajo la distribución de muestreo es exponencialmente pequeña ($\sim 2N/2^N$).
   • El muestreo directo pierde estas cadenas raras, lo que conduce a una densidad de la magia aparente de 1 (el valor del estado típico) en lugar del $1/(q-1)$ correcto. Esto explica las discrepancias en informes numéricos previos (por ejemplo, en $q=3$).

4. Conjuntos que Conservan el Número y Conexión con SYK2:

   • Para un llenado fijo, la pureza del estabilizador promediada en $q$ entero se expresa como una integral sobre $SU(2q)$.
   • En $q=2$, el conjunto de Slater complejo se evalúa exactamente, recuperando resultados conocidos.
   • Los autores establecen un diccionario entre estos conjuntos y los estados fundamentales de modelos SYK2 cuadráticos (GUE para el complejo, GOE para el real).
   • Si bien la densidad extensa de la magia es idéntica entre los conjuntos, las correcciones de orden inferior (términos algebraicos en $N$) distinguen a los conjuntos de Majorana, Slater complejo y Slater real.

5. Magia bajo Traza Parcial:

   • Los autores investigan la robustez de la magia cuando se descarta una fracción $\kappa$ de los modos.
   • Derivan un umbral $\kappa_G^\star \approx 0.7613$ por debajo del cual el estado reducido está certificado como mágico (fuera del politopo del estabilizador) con alta probabilidad.
   • Este umbral es estrictamente mayor que la transición $\kappa = 2/3$ para los estados de Haar aleatorios. Así, a pesar de ser clásicamente simulables, los estados gaussianos aleatorios retienen una magia certificable sobre un rango de pérdida de información más amplio que los estados de Haar aleatorios genéricos.

6. Promedios Annealed vs. Quenched:

   • El artículo prueba que para los conjuntos de Majorana y de Slater complejo a medio llenado, las SREs annealed y quenched comparten el mismo coeficiente extenso ($N/(q-1)$) para $q \ge 2$. La diferencia es subextensa ($O(\log N)$), confirmando que los resultados annealed describen estados típicos.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera caracterización exacta del espectro de Pauli y las entropías de estabilizador para estados fermiónicos gaussianos aleatorios, resolviendo la estructura de la magia en estos sistemas. Las reivindicaciones clave son:

• Exactitud: La derivación de fórmulas exactas de tamaño finito $N$ para los momentos y el espectro completo, evitando aproximaciones comunes en estudios numéricos previos.
• Mecanismo de Congelación: Identificar que la congelación de las SREs para $q > 2$ es impulsada por cadenas de Majorana de bajo peso raras, un mecanismo distinto de la dominancia de la cadena de identidad en los estados de Haar en el espacio de Hilbert completo.
• Limitaciones de Muestreo: Demostrar rigurosamente que las técnicas de muestreo estándar son insuficientes para estimar las SREs de orden superior en estos conjuntos debido a la probabilidad exponencialmente pequeña de las cadenas raras dominantes.
• Robustez de la Magia: Establecer que los estados gaussianos aleatorios poseen una forma de magia que es más resiliente a la traza parcial (pérdida de información) que los estados de Haar aleatorios genéricos, desafiando la intuición de que los estados clásicamente simulables son menos complejos o robustos.
• Diccionario SYK2: Proporcionar un marco analítico para distinguir los estados fundamentales de diferentes modelos SYK2 cuadráticos basándose en las correcciones de orden inferior a sus propiedades de la magia.

Los autores concluyen que, si bien la densidad de la magia extensa es universal entre estos conjuntos, la estructura espectral detallada y las correcciones de tamaño finito ofrecen una clasificación refinada de la no-estabilizerness en sistemas de fermiones libres.