Linear Time & Storage Simulation of Non-Clifford Circuits via Symmetric Cartesian Collapse: A Trajectory-Based Solution to the Exponential Bottleneck

Resumen Técnico: Simulación de Tiempo y Almacenamiento Lineal de Circuitos No-Clifford mediante el Colapso Cartesiano Simétrico

Planteamiento del Problema
El artículo identifica dos cuellos de botella exponenciales primarios que asolan la simulación clásica de circuitos cuánticos universales:

1. El Cuello de Botella de la Rotación No-Clifford: Mientras que las operaciones Clifford son rastreables eficientemente (teorema de Gottesman-Knill), las rotaciones no-Clifford (por ejemplo, puertas T, $R_z(\theta)$ arbitrarias) introducen factores de fase complejos. En los marcos de trabajo estándar de vectores de estado, estas operaciones fuerzan representaciones dispersas hacia arreglos densos, escalando los requisitos de memoria para el rastreo de "estados mágicos".
2. El Cuello de Botella del Entrelazamiento Multi-Qubit: Los simuladores de "Matriz Densa" convencionales intentan calcular todos los resultados probabilísticos simultáneamente, lo que el autor argumenta que malinterpreta fundamentalmente la naturaleza física del hardware cuántico al ignorar el potencial de un modelado de trayectoria única y discreta.

Metodología: Colapso Cartesiano Simétrico (SCC)
La solución propuesta, el Colapso Cartesiano Simétrico (SCC), reemplaza las matrices densas de vectores de estado con un enfoque basado en trayectorias utilizando una lista de vectores espaciales 3D. La metodología se basa en tres pilares centrales:

1. Representación Geométrica (Equivalencia de Coordenadas):
   En lugar de almacenar amplitudes complejas en un espacio de Hilbert de $2^n$ dimensiones, el sistema representa cada qubit como un vector de valores reales $V_i = (x_i, y_i, z_i)^T$ dentro de un espacio euclidiano 3D ($\mathbb{R}^3$). Los estados de la base computacional se mapean a los polos del eje Z, mientras que los estados de superposición se mapean a los ejes X e Y. Esto reduce la sobrecarga de memoria global de $O(2^n)$ a $O(n)$.

2. Rotaciones Geométricas Unificadas:
   Las puertas de un solo qubit (tanto Clifford como no-Clifford) se ejecutan no mediante multiplicación de matrices unitarias, sino como rotaciones puras de vectores euclidianos 3D utilizando la Fórmula de Rotación de Rodrigues. Al tratar todas las puertas como rotaciones de vectores geométricos alrededor de ejes específicos, las operaciones no-Clifford heredan la misma huella de memoria lineal y tiempo de ejecución que las puertas Clifford estándar, eliminando la penalización de complejidad de los "estados mágicos".

3. Entrelazamiento Bidireccional y Retroalimentación de Fase (Phase Kickback):
   Las operaciones de entrelazamiento multi-qubit (por ejemplo, CNOT) se modelan como operaciones de trayectoria local de dos vías en lugar de productos tensoriales globales:

   • Acción de la Puerta hacia Adelante (Forward Gate Action): La coordenada Z local del qubit de control determina si la puerta objetivo se activa.
   • Retroalimentación de Fase hacia Atrás (Backward Phase Kickback): La alineación geométrica del qubit objetivo con el eje de la puerta ejerce una rotación de retroalimentación (phase kickback) sobre el qubit de control. Esto preserva las correlaciones sin un rastreo de fase global.

4. Resolución de 3 Ejes Simétricos (El Colapso):
   Para evitar la pérdida de información durante los pasos intermedios, el marco introduce un Postulado de Vértice Cartesiano. En lugar de proyectar los estados sobre un único eje Z (lo que destruiría la información de fase X/Y), el sistema resuelve los vectores continuos simétricamente a través de los tres ejes simultáneamente. El estado colapsa al vértice más cercano de un cubo delimitador $V_{cartesian} = (\pm 1, \pm 1, \pm 1)$.

   • Esta resolución está bloqueada por probabilidad: La probabilidad de colapsar a $+1$ o $-1$ en cualquier eje está determinada por el valor de expectativa previo al colapso de esa coordenada.
   • Este mecanismo preserva la "historia direccional" y la coherencia de fase del estado en un espacio entero discreto.

Contribuciones Clave y Resultados
El artículo presenta un marco teórico y pruebas de rendimiento empíricas que afirman abordar el problema del escalamiento exponencial:

• Verificación de Escalamiento Lineal: El autor demuestra que la huella de memoria escala estrictamente de forma lineal ($O(n)$) con el número de qubits, independientemente de la profundidad de las puertas o la presencia de puertas no-Clifford.
• La Prueba de Doble Hadamard: Se presenta una prueba de reversibilidad donde un estado $|0\rangle$ experimenta una puerta Hadamard, un colapso estocástico de 3 ejes y una segunda puerta Hadamard. Los resultados muestran que el sistema regresa exitosamente al estado $|0\rangle$, demostiendo que el colapso estocástico intermedio preserva la coherencia de fase si se retiene la información de los 3 ejes.
• Prueba de Bell Paralela de 1000 Qubits: El simulador fue ejecutado en una PC de consumo estándar (8GB de RAM) para simular 1,000 qubits (500 pares entrelazados). El experimento consistió en pasar este circuito de 1,000 qubits a través de 1,000 trayectorias independientes para recolectar las frecuencias finales.
  • Preservación del Entrelazamiento: El simulador mantuvo una fuerte correlación entre los pares entrelazados, registrando la frecuencia empírica de estados mixtos inválidos (por ejemplo, $|01\rangle$ o $|10\rangle$) como 0.00% (un valor redondeado) en los datos reportados.
  • Precisión: Las frecuencias empíricas para varios ángulos de rotación ($\theta = 30^\circ, 45^\circ, 60^\circ$) coincidieron estrechamente con las predicciones teóricas de la mecánica cuántica con desviaciones menores (por ejemplo, 84.9% frente a 85.36% para $\theta=45^\circ$).
• Eficiencia No-Clifford: Las pruebas de rendimiento confirmaron que los circuitos saturados con puertas T y rotaciones arbitrarias escalan de manera idéntica a los circuitos de solo Clifford, mostrando que no hay un pico exponencial en el costo computacional.

Significado y Reivindicaciones
El artículo afirma que el marco SCC cambia fundamentalmente el paradigma de la simulación cuántica de "calcular todas las posibilidades" a "rastrear trayectorias únicas bloqueadas por probabilidad".

• Rompiendo la Barrera de los 50 Qubits: El autor asegura que este método permite la simulación de más de 1,000 qubits en una máquina de 8GB de RAM, señalando que este enfoque puede lograrlo en "menos de diez segundos", aunque el tiempo de ejecución específico para la prueba de 1,000 qubits no se detalla explícitamente en el texto.
• Reversión del Cuello de Botella: El trabajo argumenta que el cuello de botella exponencial de memoria es un artefacto del método de simulación, no una limitación física. Al imponer un colapso simétrico de 3 ejes, el marco preserva la coherencia cuántica sin matrices de densidad globales.
• Potencial Híbrido: El modelo basado en trayectorias sugiere una nueva arquitectura donde los algoritmos clásicos pueden interactuar directamente con las puertas cuánticas, lo que podría simplificar operaciones complejas como la aritmética modular de Shor al tratar las mismas como funciones clásicas integradas en la trayectoria cuántica.

Limitaciones Reconocidas
El artículo reconoce modestamente que, si bien el SCC elimina los requisitos de memoria exponenciales, transfiere el desafío computacional a la construcción de puertas. Las operaciones que no corresponden naturalmente a rotaciones geométricas simples (como la Transformada de Fourier Cuántica o la aritmética modular) deben descomponerse en secuencias de puertas de rotación fundamentales, lo que requiere métodos de diseño eficientes para puertas compuestas. El marco no elimina todas las dificultades, sino que las reubica del almacenamiento a la construcción.