Un procesador cuántico de IBM recrea una colisión de partículas en una teoría de campo unidimensional

Qué lograron simular

Roland C. Farrell, Nikita A. Zemlevskiy, Marc Illa y John Preskill realizaron una simulación digital, en tiempo real, de una colisión en una teoría de campo de Ising en una dimensión espacial. El cálculo se ejecutó en el procesador cuántico ibm_marrakesh de IBM, con 104 qubits, y siguió la dinámica antes, durante y después de una colisión entre dos partículas. El experimento alcanzó hasta 45 pasos de Trotter y empleó hasta 5.589 puertas de dos qubits. 7

El resultado principal fue evidencia —no una medición con la precisión de un acelerador de partículas— de producción inelástica de partículas: parte de la energía cinética de las partículas ligeras incidentes fue compatible con la aparición, tras el choque, de una excitación pesada junto a una partícula ligera. 7

Dos paquetes de onda que se aproximan y chocan

La simulación comenzó con la preparación de dos paquetes de onda de partículas ligeras, separados espacialmente y con momentos opuestos. Es decir, se configuraron para avanzar uno hacia el otro. Después, el equipo evolucionó el hamiltoniano de la teoría de campo de Ising mediante dinámica en tiempo real trotterizada y midió la densidad local de energía mientras los paquetes se aproximaban, colisionaban y volvían a separarse. 7

Cada paquete de una sola partícula se preparó en dos fases:

1. Se construyó un estado adaptado, similar a ( |W(k_0)\rangle ), que determina la envolvente espacial, el contenido en momento y los números cuánticos del paquete. Esa primera preparación aún puede contener componentes de múltiples partículas.
2. Se aplicaron circuitos variacionales que preservan las simetrías para minimizar la energía y acercar el estado al autoestado deseado de una sola partícula. 7

Antes de llevar los circuitos al procesador de 104 qubits, los autores diseñaron los circuitos de paquetes de onda para la red grande mediante simulaciones clásicas de estados producto de matrices, o MPS por sus siglas en inglés. 7

La señal de una colisión inelástica

En una colisión elástica, las dos partículas ligeras salientes pueden cambiar sus momentos, pero no se crea una especie nueva de partícula. En ese escenario, los perfiles de densidad de energía después del choque deberían ser relativamente simétricos. 7

La situación cambia si aparece una partícula pesada y lenta. Su contribución deja un perfil de energía asimétrico. Para cuantificarlo, los investigadores analizaron la asimetría estadística, o skewness, de la densidad de energía medida. Un valor distinto de cero, coherente con la trayectoria prevista para la partícula pesada, permitió diferenciar el canal inelástico del comportamiento de referencia elástico. 7

Por tanto, la conclusión no se basó únicamente en observar energía desplazándose por el sistema. La asimetría observada era consistente con que la energía cinética inicial de dos partículas ligeras terminara en un estado final formado por una partícula ligera y otra más pesada. 7

Por qué fue decisivo el método basado en estados W

Preparar paquetes de onda extensos con métodos convencionales suele exigir circuitos cada vez más profundos a medida que el paquete ocupa más sitios de la red. El equipo adaptó protocolos de preparación de estados W para producir de forma eficiente el estado inicial parecido a ( |W(k_0)\rangle ). 27

La clave fue usar mediciones durante el circuito y realimentación clásica: una medición intermedia genera un resultado clásico, que se utiliza para decidir condicionalmente operaciones posteriores. Esto permitió un protocolo sin qubits auxiliares para preparar estados W con profundidad constante. Según los autores, la profundidad deja de depender del tamaño del paquete de onda y de la dimensión espacial, salvo por su dependencia de la longitud de correlación física. 7

Los investigadores describen la mejora como superexponencial frente a métodos anteriores, que requerían preprocesamiento clásico exponencialmente costoso o circuitos cuya profundidad crecía polinómicamente con el tamaño del paquete. 7

El vínculo con problemas más ambiciosos, como la QCD

El método no se limitó a la teoría de Ising 1D: los autores también construyeron o demostraron el enfoque de preparación para teoría escalar unidimensional, el modelo de Schwinger y teoría de campo de Ising bidimensional. 7

La motivación de largo plazo es estudiar dinámicas cuánticas de campos fuertemente acoplados y en tiempo real. En las colisiones de alta energía, la energía cinética puede transformarse en nuevas partículas; predecir sus abundancias y su comportamiento colectivo es esencial para comprender la producción de hadrones y, en última instancia, la cromodinámica cuántica o QCD. Las simulaciones clásicas encuentran dificultades particulares en dinámicas genéricas en tiempo real, un terreno donde los ordenadores cuánticos podrían aportar ventajas en el futuro. 7

Qué no demuestra todavía este experimento

Conviene poner el alcance en perspectiva:

• Se trata de una demostración en teoría de campo de Ising en una dimensión, no de una simulación de QCD, de chorros hadrónicos tridimensionales ni de un evento realista de un colisionador. 7
• El hardware era ruidoso y se emplearon técnicas de mitigación de errores, no corrección cuántica de errores plenamente tolerante a fallos. Por ello, el experimento aporta evidencia del proceso modelado, pero no un cálculo de dispersión sistemáticamente preciso y tolerante a fallos. 7
• Llevar este enfoque a teorías gauge en dimensiones superiores, y especialmente a la QCD, exigirá máquinas mucho mayores y más precisas, probablemente tolerantes a fallos. También harán falta métodos escalables para preparar estados, codificar campos gauge, realizar evoluciones largas en tiempo real y extraer observables precisos de muchos cuerpos. 7

El avance es, así, un paso experimental y algorítmico relevante: muestra cómo un procesador cuántico actual puede seguir la firma de una colisión inelástica en un modelo de teoría cuántica de campos simplificado, sin equivaler aún a una simulación realista de física de colisionadores.