Tomografía de conjuntos de espines: cómo convertir un espectro en un estado cuántico
Contexto
El trabajo va por libre: los vecinos más cercanos en el archivo se dedican a algoritmos cuánticos, corrección de errores e información cuántica, pero no tocan la reconstrucción práctica de estados a partir de espectros. Ya existen herramientas para simular la dinámica de espines (PULSEE) y bibliotecas generales (QuTiP, nmrglue), pero faltaba el eslabón que uniera los datos en bruto del espectrómetro con una matriz de densidad validada. El paquete Tomography-NMR llena ese nicho y se propone como objetivo un procesamiento de datos de RMN transparente y reproducible para la tomografía cuántica.
Métodos
El paquete implementa una cadena de análisis en Python 3.11+: las señales en bruto en el dominio del tiempo se someten a la transformada de Fourier con corrección de fase automática o manual, y después las intensidades de los picos de los dobletes se extraen de tres maneras. La primera es la medición directa de la altura del pico en un solo punto; la segunda, la integración numérica por Simpson en una ventana fija; la tercera, un barrido sistemático de parámetros (desplazamientos de frecuencia, anchura de la ventana, corrección del desdoblamiento) para maximizar la precisión respecto a un estado teórico conocido. La suma y la diferencia de los picos izquierdo y derecho del doblete se asignan a los coeficientes de los operadores producto mediante una tabla de búsqueda, y después la matriz de densidad se reconstruye imponiendo restricciones físicas (hermiticidad, normalización de la traza, semidefinición positiva).
Resultados
En 20 estados experimentales, preparados sobre dos núcleos de fósforo-31 acoplados en ADP, la precisión media de reconstrucción según la métrica de Fortunato fue de 0.975 a 0.995 en función de la categoría del estado. Los estados de la base alcanzaron 0.995, las salidas de la puerta de Hadamard 0.992, las de CNOT 0.983, los estados de Bell 0.980 y la combinación H+T 0.975. La integración optimizada da una precisión superior al 99% para estados concretos (hasta 0.998 para |00⟩), mientras que los métodos que no conocen el estado objetivo se quedan en torno al 98%. Los errores de los elementos de la matriz de densidad se estiman en el orden de 10⁻³–10⁻⁴.
Implicaciones
A diferencia de los trabajos vecinos, centrados en algoritmos o en simulación, este trabajo resuelve el problema inverso: recuperar el estado cuántico a partir de espectros experimentales reales. Aporta la primera cadena abierta y completamente documentada para la tomografía por RMN, elimina la dependencia de programas comerciales cerrados y permite verificar de forma independiente los resultados publicados. Esto cambia la práctica de la informática cuántica experimental y la vuelve más transparente y reproducible.
Desarrollo futuro
Ampliar el método a sistemas de tres o más cúbits exigirá resolver el problema del solapamiento de multipletes y el de la escalabilidad de la optimización. Ya se incluye un prototipo de módulo para tres espines, pero el soporte completo sigue siendo una tarea abierta. También se planean la optimización bayesiana en lugar de la búsqueda en rejilla, el análisis en tiempo real para protocolos adaptativos y nuevas métricas de validación, como los testigos de entrelazamiento.
Impacto
El trabajo influirá en la informática cuántica experimental, sobre todo en la computación cuántica por RMN, al ofrecer una herramienta estandarizada para caracterizar estados. También acelerará la formación de nuevos especialistas y facilitará la comparación de resultados entre grupos.
Próximos pasos
Los próximos pasos incluyen implementar la optimización bayesiana para acelerar la búsqueda de parámetros y desarrollar métodos de deconvolución espectral para sistemas de tres cúbits. También es importante ampliar el soporte de formatos de datos de otros espectrómetros.
Problemas abiertos clave
En información cuántica sigue sin resolverse el problema de la tomografía de estados escalable: el número de mediciones crece exponencialmente con el número de cúbits. Este trabajo no lo resuelve, pero fija una referencia de precisión y transparencia para el caso de dos cúbits, algo necesario para poner a prueba los métodos futuros.
El paquete se llama Tomography-NMR y se distribuye bajo licencia MIT, lo que permite usarlo libremente en proyectos académicos y comerciales. La validación se llevó a cabo con la molécula de ADP, donde dos núcleos de fósforo-31 están acoplados mediante una interacción J de aproximadamente 22 Hz.
Cifras clave
- precisión de reconstrucción (optimizada): hasta 0.998
- precisión de reconstrucción (sin conocer el objetivo): alrededor de 0.980
- número de estados validados: 20
- número de coeficientes no identidad para dos cúbits: 15
- tiempo de optimización para un estado: 2–5 minutos