A Quantum Computing Approach to Track Reconstruction in Strip-Type Detectors

Resumen Técnico: Un enfoque de computación cuántica para la reconstrucción de trazas en detectores de tipo tira

Declaración del problema
El seguimiento de partículas en detectores gaseosos de tipo tira, particularmente en entornos de alta multiplicidad como los experimentos de base corta para la búsqueda de sectores oscuros, enfrenta desafíos combinatorios significativos. La presencia de "hits fantasma" (señales espurias derivadas de ambigüedades del detector) y múltiples combinaciones de hits transforma el reconocimiento de patrones en un problema de optimización complejo. Los métodos tradicionales tienen dificultades para seleccionar eficientemente las combinaciones de hits consistentes que describen con precisión las trayectorias de las partículas en medio de un alto nivel de pileup y fondos de radiación. Este estudio investiga si el recocido cuántico (quantum annealing) puede resolver eficazmente los problemas de optimización combinatoria inherentes a la reconstrucción de trazas para tales detectores.

Metodología
Los autores formulan dos subproblemas distintos de la reconstrucción de trazas como problemas de Optimización Binaria Cuadrática sin Restricciones (QUBO, por sus siglas en inglés), diseñados para ser resueltos en un recocedor cuántico D-Wave. El estudio utiliza eventos simulados del detector DAMSA, una configuración compacta de blanco fijo que cuenta con un apilamiento de 10 capas de centelladores de CsI intercalados con detectores de tiras $\mu$RWELL.

1. Selección de Hits de Traz Única:

   • Objetivo: Seleccionar un hit geométricamente consistente de cada una de tres capas de detectores dentro de un candidato a cúmulo de fotones localizado, rechazando los hits fantasma.
   • Formulación: Las variables binarias representan pares de hits (segmentos) entre dos capas e hits individuales en la tercera capa. La función de energía QUBO incluye:
     • Restricciones: Términos de penalización cuadrática que imponen la selección de exactamente un segmento y exactamente un hit de la tercera capa.
     • Objetivo: Un término de "recompensa de rectitud" que maximiza la linealidad del triplete seleccionado al minimizar el residuo geométrico entre el segmento extrapolado y el tercer hit.
   • Implementación: Para gestionar los recursos de la QPU, el problema se descompone en bloques independientes de 3 capas (trozos o chunks), que se ensamblan en una matriz QUBO diagonal por bloques y se envían como un lote.

2. Asociación Simultánea de Candidatos de Trazas:

   • Objetivo: Asociar candidatos a cúmulos de diferentes capas de detectores en candidatos a trazas 3D, permitiendo la selección simultánea de múltiples candidatos a trazas (por ejemplo, dos fotones de una desintegración de ALP) dentro de un solo envío a la QPU.
   • Formulación: Las variables binarias representan aristas entre cúmulos en capas adyacentes. La función de energía QUBO incluye:
     • Recompensas: Términos que favorecen la selección de aristas y la rectitud geométrica (minimizando la curvatura).
     • Restricciones: Penalizaciones de grado que suprimen la ramificación/fusión (múltiples aristas hacia/desde el mismo cúmulo) y restricciones de flujo que aseguran la continuidad a través de la capa intermedia (las aristas entrantes son iguales a las salientes).
   • Post-procesamiento: Los tripletes de cúmulos seleccionados se decodifican y se ensamblan en candidatos a trazas más largos utilizando una regla de solapamiento de grafos, donde los tripletes que comparten al menos dos asignaciones de capa-cúmulo idénticas se conectan. No se realiza ningún ajuste de Kalman (Kalman fitting) ni optimización de residuos adicional en este paso.

Resultados Clave
Las formulaciones propuestas fueron probadas contra eventos simulados de DAMSA dominados por señales de dos fotones provenientes de la desintegración de un Axion-Like Particle (ALP) con bajo pileup.

• Desempeño de la Selección de Hits de Traz Única:

  • La reconstrucción basada en QPU se comparó con una reconstrucción estándar de filtro de Kalman.
  • Resoluciones: El método QPU logró resoluciones de posición y angulares cercanas a la línea base de Kalman. Para los residuos de posición, la media de la QPU fue de 0.4204 cm (frente a 0.3925 cm para Kalman) con una anchura del 68% de 0.2951 cm (frente a 0.2548 cm). Para los residuos angulares, la media de la QPU fue de 0.0679 rad (frente a 0.0635 rad para Kalman) con una anchura del 68% de 0.0446 rad (frente a 0.0425 rad).
  • Forma de la Distribución: Las comparaciones estadísticas ( $\chi^2$ reducido $< 1$, baja divergencia KL) indicaron que el método QUBO reprodujo con éxito la forma general de las distribuciones de residuos de Kalman, a pesar de la naturaleza discreta de la selección por bloques.

• Desempeño de la Asociación Simultánea:

  • Recuperación de Tripletes: La QPU identificó con éxito los tripletes de cúmulos correctos con una fracción de recuperación mediana cercana a la unidad a través de varias longitudes de referencia de candidatos.
  • Ensamblaje de Trazas: El paso de ensamblaje de grafos posterior recuperó los candidatos a trazas de referencia con alta eficiencia, demostando que los tripletes superpuestos pueden conectarse eficazmente sin necesidad de ajustes adicionales.

• Recursos y Latencia:

  • Escalabilidad: El costo de qubits físicos por bloque escaló suavemente con el número de variables lógicas, lo que indica que la estrategia de matriz de bloques mantiene controlados los costos de incrustación (embedding).
  • Latencia: Para la selección de traza única, la latencia media total de extremo a extremo fue de aproximadamente 96.5 ms por evento. Para la asociación simultánea, la latencia fue de aproximadamente 1137.5 ms por evento, impulsada principalmente por las sobrecargas de programación e incrustación más que por el tiempo de muestreo en sí.

Significancia y Reivindicaciones
Los autores posicionan este trabajo como una demostración de prueba de concepto más que como una solución definitiva para entornos de alta ocupación. El estudio afirma que:

1. Las formulaciones QUBO pueden reproducir eficazmente decisiones de reconstrucción local (selección de hits) comparables con los filtros de Kalman clásicos en una topología de detector específica.
2. El recocido cuántico puede manejar la selección simultánea de múltiples candidatos a trazas dentro de un solo envío, separando la optimización cuántica de los tripletes locales del ensamblaje de grafos clásico de las trazas globales.
3. La estrategia de envío de diagonal por bloques es un método práctico para enviar múltiples problemas de reconstrucción independientes al hardware cuántico actual.

El documento señala explícitamente que los resultados actuales se basan en un entorno de bajo pileup y disperso, optimizado para la identificación de señales. Los autores declaran que se requieren estudios futuros para probar estas formulaciones bajo condiciones de pileup más fuertes, mayores multiplicidades de cúmulos y ruido de detector realista para evaluar plenamente la robustez del método para topologías de seguimiento complejas. El trabajo proporciona una base práctica para el desarrollo de métodos de reconstrucción híbridos cuántico-clásicos para futuros experimentos.