GroupFunctions.jl: computing individual entries of the irreducible representations of the unitary group U(d)

Resumen Técnico de GroupFunctions.jl

Planteamiento del Problema
Las representaciones del grupo unitario $U(d)$ son fundamentales para diversos subcampos de la física y las matemáticas, particularmente en la óptica cuántica y la combinatoria algebraica. Una tarea computacional central en estos campos es la evaluación de elementos matriciales individuales de representaciones irreducibles (irreps), conocidos como "funciones de grupo". Estos elementos, denotados como $D^{(\lambda)}_{\text{out,init}}(U) = \langle \text{out} | D^{(\lambda)}(U) | \text{init} \rangle$, son necesarios para tareas que van desde el cálculo de amplitudes de transición de fotones en redes ópticas lineales hasta la caracterización de dispositivos cuánticos y el análisis de propiedades de simetría.

Si bien la estimación numérica suele ser suficiente para algunas aplicaciones, otros contextos —como la derivación de expresiones simbólicas exactas para amplitudes de transición o caracteres— requieren una computación simbólica precisa. Las soluciones de software existentes abordan problemas relacionados pero distintos: algunas computan coeficientes de Clebsch-Gordan (por ejemplo, SUNRepresentations.jl), otras proporcionan acceso numérico indirecto a funciones de grupo (por ejemplo, RepLAB), y algunas se centran en la integración simbólica sobre grupos (por ejemplo, IntegrateUnitary.jl). Las librerías dedicadas a la óptica cuántica (por ejemplo, BosonSampling.jl, The Walrus) suelen modelar dispositivos numéricamente o computar amplitudes específicas como permanentes, pero ninguna está diseñada para computar individualmente las entradas de la matriz de representación de forma simbólica. Además, ensamblar una matriz de representación completa suele ser computacionalmente ineficiente cuando solo se requiere una única entrada.

Metodología
El artículo presenta GroupFunctions.jl, una librería de Julia diseñada para computar estos elementos de matriz individuales de forma simbólica o numérica. La librería opera dentro de una base de espacio portador enumerada mediante patrones de Gelfand-Tsetlin (GT).

El enfoque algorítmico central se basa en la fórmula de Grabmeier-Kerber, que expresa un elemento de matriz como una suma de monomios en las entradas de la matriz de entrada $U$. Estos monomios están ponderados por la representación irreducible y los estados de la base específicos involucrados. La implementación optimiza este proceso mediante:

1. Enumeración de Dobles Cocets: En lugar de iterar sobre todas las permutaciones en el grupo simétrico $S_N$ (donde $N$ es el número de cajas en el diagrama de Young), el algoritmo agrupa las permutaciones que contribuyen con monomios idénticos.
2. Sumatoria Eficiente: Construye dobles cocets e itera sobre sus elementos de permutación representativos, reduciendo significativamente la carga computacional en comparación con un enfoque ingenuo de permutación por permutación.
3. Representación Interna: Los estados de la base se representan internamente como tablas de Young semiestándar, mientras que la interfaz de usuario expone patrones GT equivalentes a través de la estructura de datos GTPattern.

La librería soporta matrices de entrada simbólicas (usando SymEngine.Basic) y proporciona funciones de utilidad para traducir patrones GT en kets de número de ocupación, construir unitarias de entrada a partir de parametrizaciones comunes de la óptica cuántica y computar funciones de Schur asociadas.

Contribuciones Clave

• Computación Simbólica de Entradas Individuales: GroupFunctions.jl es la primera librería, según conocen los autores, capaz de computar entradas individuales de matrices de representación de $U(d)$ de forma simbólica.
• Marco Unificado: Proporciona un método unificado para computar elementos de matriz para irreps especificadas por particiones enteras de longitud a lo sumo $d$, manejando tanto entradas simbólicas como numéricas.
• Especialización para Óptica Cuántica: La librería incluye herramientas específicas para escenarios de la óptica cuántica, tales como la evaluación de amplitudes de transición para el muestreo de bosones (boson sampling) y el manejo de irreps simétricas (que modelan bosones completamente indistinguibles).
• Capacidades de Exportación: Los resultados pueden exportarse en un formato compatible con Mathematica, facilitando la manipulación simbólica posterior.

Resultados y Rendimiento
El artículo demuestra las capacidades de la librería a través de ejemplos, incluyendo la computación simbólica de funciones de grupo de $SU(2)$ (que se reducen a funciones D de Wigner) y una irrep simétrica de $U(10)$ que involucra 7 bosones.

• Rendimiento: En un CPU AMD Ryzen 7 PRO 4750U, la computación simbólica de un elemento de matriz específico de $U(10)$ tomó aproximadamente 4 segundos tras la calentamiento de la compilación, sin reutilizar resultados almacenados en caché.
• Escalabilidad: El costo computacional está gobernado por tres factores: la dimensión del grupo $m$, la dimensión del espacio de Gelfand-Tsetlin $d_\lambda(m)$ y la dimensión de la irrep del grupo simétrico $f^\lambda$ (el número de tablas de Young estándar). Los autores señalan que, aunque el trabajo subyacente involucra permutaciones que pueden acercarse a $N!$, el cuello de botella práctico suele ser el crecimiento de $N$ o de la dimensión $f^\lambda$, más que el aumento de $m$ mediante ceros finales.
• Aplicación: La librería ya ha sido utilizada para evaluar numéricamente los caracteres de grupos $SU(d)$ en pruebas de benchmarking aleatorio (randomized benchmarking).

Significado
El significado de GroupFunctions.jl radica en su capacidad para proporcionar soluciones simbólicas exactas para las entradas individuales de las matrices de representación, una capacidad que no estaba disponible previamente en los ecosistemas de software existentes. Al evitar la necesidad de ensamblar matrices de representación completas, la librería ofrece una ruta más eficiente para cálculos específicos requeridos en la combinatoria algebraica y la física cuántica. Los autores posicionan la herramienta como un recurso especializado para investigadores que necesitan expresiones analíticas exactas para amplitudes de transición o evaluaciones de caracteres, llenando un vacío entre los simuladores ópticos puramente numéricos y las librerías de teoría de grupos de propósito general que carecen de la computación directa de entradas simbólicas.