Exact chiral symmetry with quantum signal processing

Resumen Técnico: Simetría Quiral Exacta con Procesamiento de Señal Cuántico

Planteamiento del Problema
La simulación cuántica de observables no perturbativos y fuera del equilibrio en la cromodinámica cuántica (QCD) en el retículo enfrenta desafíos significativos relacionados con el doble de fermiones y la preservación de la simetría quiral en las formulaciones hamiltonianas. Si bien la relación de Ginsparg-Wilson (GW) y el operador de overlap proporcionan soluciones elegantes en el espacio-tiempo euclídeo, ha faltado un análogo hamiltoniano canónico para los fermiones de overlap. Los enfoques existentes, como los fermiones de pared de dominio (domain-wall fermions), se mapean bien a los entornos hamiltonianos pero requieren una dimensión extra explícita, lo que aumenta los costos de cúbits. Por el contrario, los fermiones de overlap ofrecen una formulación en las dimensiones espaciales físicas, pero implican interacciones de todos contra todos altamente no locales, específicamente la función signo del Hamiltoniano de Wilson-Dirac, $\varepsilon(h_W)$, la cual es difícil de implementar eficientemente en hardware cuántico. El problema central abordado es cómo construir algoritmos cuánticos eficientes para fermiones de retículo con simetría quiral que equilibren los costos de memoria (cúbits) frente a la complejidad de compuertas, manteniendo una simetría quiral exacta (o controladamente rota).

Metodología
Los autores proponen un algoritmo de Procesamiento de Señal Cuántico (QSP) para simular el Hamiltoniano de fermiones de overlap. La metodología procede a través de pasos clave:

1. Formulación Hamiltoniana: El trabajo utiliza el Hamiltoniano de overlap $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, donde $h_W$ es el Hamiltoniano de partícula única de Wilson-Dirac. La función signo $\varepsilon(X) = X/\sqrt{X^\dagger X}$ se aproxima mediante un polinomio $p_M(x)$ de grado $M$.
2. Codificación de Bloque (Block Encoding): Los autores construyen una codificación de bloque del Hamiltoniano de Wilson de partícula única $\hat{h}_W$ (incluyendo campos de gauge) utilizando un operador "preparar" $P$ y un operador "seleccionar" $S$. Esto codifica el Hamiltoniano en un operador unitario $U_{\hat{h}_W}$ que actúa en un espacio de Hilbert extendido con cúbits ancila.
3. Procesamiento de Señal Cuántico (QSP): Para implementar la aproximación de la función signo, los autores emplean QSP (específicamente Transformación de Valores Singulares Cuánticos). Esto permite la aplicación de un polinomio de grado $M$, $p_M(U_{\hat{h}_W})$, al operador codificado en bloque. El polinomio se elige para aproximar la función signo dentro de un brecha espectral $[\kappa, 1]$ con un error $\epsilon_e$.
4. Evolución Temporal: Una vez que el Hamiltoniano de overlap está codificado en bloque, el operador de evolución temporal $e^{-iHt}$ se aproxima utilizando QSP, requiriendo un número de compuertas que escala con el costo de la codificación de bloque y el tiempo de evolución.

Contribuciones Clave y Resultados

• Construcción Algorítmica: El artículo proporciona un algoritmo concreto basado en QSP para el Hamiltoniano de overlap que preserva la relación GW hasta un error controlable $\epsilon_e$. Se demuestra que el operador quiral modificado $\hat{\gamma}_5$ conmuta con el Hamiltoniano aproximado hasta un error de orden $\epsilon_e$.
• Análisis de Complejidad:
  • Complejidad de Compuertas: El costo de codificar en bloque el Hamiltoniano de overlap escala como $O(Q \kappa^{-1} \log(1/\epsilon_e))$, donde $Q$ es el número de sitios del retículo (veces los grados de libertad internos) y $\kappa$ está relacionado con la brecha espectral. La complejidad total de compuertas para la evolución temporal escala como $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Costo de Cúbits: El algoritmo requiere $O(Q + \log Q)$ cúbits. Esta es una reducción significativa comparada con los fermiones de pared de dominio, que requieren $O(Q L_5)$ cúbits, donde $L_5$ es la extensión de la dimensión extra.
  • Comparación: Mientras que los fermiones de pared de dominio se benefician de la localidad geométrica (permitiendo un costo casi lineal $O(Q t \text{polylog})$ vía Suzuki-Trotter o métodos similares), la formulación de overlap incurre en una mayor profundidad de compuertas debido a la naturaleza no local de la aproximación de la función signo. Sin embargo, el enfoque de overlap ofrece una escala asintótica superior en memoria.
• Interpretación Física de QSP: Los autores demuestran que el grado del polinomio $M$ requerido para lograr un error $\epsilon_e$ escala como $M = O(\kappa^{-1} \log(1/\epsilon_e))$. Identifican este escalamiento con el tamaño de la dimensión extra $L_5$ en los fermiones de pared de dominio, donde el error escala como $e^{-c L_5}$. Así, la implementación de QSP efectivamente "construye" una dimensión extra a través de la profundidad del circuito, reflejando la correspondencia física entre el operador de overlap y la teoría de frontera de los fermiones de pared de dominio.

Significancia y Reivindicaciones
El artículo afirma que las simulaciones cuánticas de fermiones de Dirac con simetría quiral exacta son "casi gratuitas" en el sentido de que aplicar el Hamiltoniano de overlap cuesta solo un factor logarítmico más (en términos de tolerancia al error) que el Hamiltoniano de Wilson-Dirac.

La idea central es que QSP proporciona una realización algorítmica cuántica de la conocida correspondencia entre el operador de overlap y la dimensión extra de los fermiones de pared de dominio. El compromiso (trade-off) es explícito:

• Fermiones de pared de dominio: Mayor costo de cúbits ($O(Q L_5)$) pero menor profundidad de compuertas debido a la localidad.
• Fermiones de overlap (vía QSP): Menor costo de cúbits ($O(Q)$) pero mayor profundidad de compuertas ($O(Q^2)$ en el escalamiento de peor caso para la evolución temporal debido a la no localidad).

Los autores concluyen que la elección entre formulaciones depende de las restricciones del hardware (conteo de cúbits frente a profundidad de circuito). Señalan que, si bien la QCD clásica a menudo utiliza aproximaciones racionales (por ejemplo, Zolotarev) para la función signo, las cuales son más eficientes, estas no tienen análogos directos en el QSP estándar, lo que sugiere una dirección para trabajos futuros. El artículo también aclara que los resultados de escalamiento son robustos frente a la codificación específica de los campos de gauge, aunque los factores constantes pueden variar.