Phase-sensitive representation of Majorana stabilizer states

Resumen Técnico: Representación Sensible a la Fase de Estados Estabilizadores de Majorana

1. El Problema

La simulación clásica de sistemas cuánticos fermiónicos (como los utilizados en química cuántica) es computacionalmente costosa debido a la complejidad de representar estados altamente entrelazados.

• Limitaciones actuales: Los estados estabilizadores (estados generados aplicando puertas Clifford a un estado base) son una clase importante de estados que pueden ser altamente entrelazados pero eficientemente simulados en sistemas de qubits. Sin embargo, en sistemas fermiónicos, la representación tradicional a menudo requiere transformaciones a qubits (como la transformación de Jordan-Wigner), lo que pierde la localidad de los operadores fermiónicos y la interpretabilidad física.
• La brecha: Aunque se ha definido el grupo Clifford de Majorana (operadores que transforman operadores de Majorana en otros operadores de Majorana), faltaba una representación sensible a la fase completa para los estados estabilizadores de Majorana. Sin esta fase, no es posible representar estados fermiónicos generales como superposiciones de estados estabilizadores, ni calcular amplitudes de probabilidad o productos internos de manera exacta, lo cual es crucial para métodos de simulación de bajo rango (low-rank).

2. Metodología

Los autores proponen una extensión del formalismo CH (nombre derivado de los autores Bravyi, Browne, et al., que lo introdujeron para qubits) al contexto fermiónico utilizando operadores de Majorana.

• Definición del Estado: Un estado estabilizador de Majorana se define como $|\psi\rangle = U|0\rangle$, donde $U$ es un producto de operadores Clifford de Majorana y $|0\rangle$ es el estado de vacío fermiónico.
• Representación Propuesta: El estado se representa mediante la tupla $(\phi, U_C, b, s)$, donde:
  • $\phi$: Un entero que codifica la fase global ($e^{i\pi\phi/4}$).
  • $U_C$: Un operador Clifford de tipo "control" (que deja invariante el estado de vacío), representado mediante una tabla de estabilizadores (matrices binarias $E, F, G$ y un vector de fases $\omega$).
  • $b$: Un vector binario que define un producto de operadores de trenzado (braiding operators) $B_k = e^{-i\frac{\pi}{4}c_k\tilde{c}_{k+1}}$.
  • $s$: Una cadena de bits que define el estado de la base computacional.
• Operadores de Majorana: Se utiliza una notación específica $(\varphi, z, x)$ para representar cadenas de operadores de Majorana, análoga a la representación simpléctica de los operadores de Pauli, permitiendo verificar conmutaciones y productos eficientemente.

3. Contribuciones Clave

El artículo presenta un conjunto completo de algoritmos para manipular estos estados:

1. Reglas de Actualización (Update Rules): Se derivan fórmulas explícitas para actualizar la representación $(\phi, U_C, b, s)$ bajo la acción de:

   • Puertas Clifford de Majorana específicas ($\eta_j, W_{jk}, \eta_{jk}$).
   • Rotaciones Clifford generales de Majorana.
   • Operadores de Majorana generales.
   • Nota: Las actualizaciones para $\eta_{jk}$ son las más complejas, siguiendo una estrategia similar a la aplicación de puertas Hadamard en el formalismo CH para qubits, implicando la manipulación de la tabla de estabilizadores y la reestructuración de la superposición de estados base.

2. Cálculo de Amplitudes de Probabilidad: Se proporciona un algoritmo para calcular $\langle x | \psi \rangle$ (la amplitud de un estado base computacional $|x\rangle$ en el estado $|\psi\rangle$) con una complejidad de $O((|x| + c)n)$. Esto es esencial para muestrear distribuciones de probabilidad.

3. Cálculo de Productos Internos: Se desarrolla un método para calcular el producto interno $\langle \phi | \psi \rangle$ entre dos estados estabilizadores de Majorana con una complejidad de $O(n^3)$. Esto permite evaluar superposiciones y solapamientos necesarios para métodos variacionales.

4. Valores Esperados: Se establece cómo calcular el valor esperado de un operador de Majorana en un estado estabilizador, lo cual es fundamental para evaluar energías en simulaciones de química cuántica.

4. Resultados y Complejidad Computacional

Los algoritmos propuestos mantienen la eficiencia característica de los métodos de estabilizadores, permitiendo simulaciones escalables:

• Actualización de Puertas:
  • Puertas diagonales ($\eta_j, W_{jk}$): $O(n)$.
  • Puertas de trenzado ($\eta_{jk}$) y rotaciones generales: $O(n^2)$.
• Amplitudes: $O(n^2)$ (dependiendo del peso de la cadena de bits).
• Producto Interno: $O(n^3)$, dominado por la inversión de matrices y multiplicaciones necesarias para actualizar las tablas de estabilizadores.
• Validación: Los autores proporcionan una implementación en Python (disponible en GitHub) que valida la corrección de los algoritmos comparándolos con el simulador de vectores de estado de Qiskit (usando la representación de Jordan-Wigner), confirmando la exactitud de los resultados.

5. Significado e Impacto

Este trabajo es fundamental para el avance de la simulación cuántica fermiónica y la química cuántica:

• Simulación Directa en Picture Fermiónico: Permite realizar simulaciones de circuitos cuánticos fermiónicos directamente en el lenguaje de operadores de Majorana, evitando la sobrecarga y la pérdida de localidad inherente a las transformaciones a qubits.
• Métodos de Bajo Rango (Low-Rank): Facilita la representación de estados fermiónicos generales como superposiciones de estados estabilizadores de Majorana ($|\Psi\rangle = \sum \alpha_i |\phi_i\rangle$). Esto es crucial para métodos de descomposición de estabilizadores que buscan reducir la complejidad de la simulación de sistemas fuertemente correlacionados.
• Herramienta para Nuevos Algoritmos: Proporciona la infraestructura matemática necesaria para desarrollar nuevos algoritmos de corrección de errores, estimación de fase y optimización variacional en el dominio fermiónico, aprovechando la estructura del grupo Clifford de Majorana.

En resumen, el artículo cierra una brecha teórica y práctica al proporcionar las herramientas necesarias para manipular, actualizar y medir estados fermiónicos complejos de manera eficiente, manteniendo la localidad y la estructura física intrínseca de los sistemas de muchos cuerpos.