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

Sintesi Tecnica di GroupFunctions.jl

Enunciato del Problema
Le rappresentazioni del gruppo unitario $U(d)$ sono fondamentali per vari sottocampi della fisica e della matematica, in particolare nell'ottica quantistica e nella combinatoria algebrica. Un compito computazionale centrale in questi campi è la valutazione dei singoli elementi matriciali di rappresentazioni irriducibili (irreps), noti come "funzioni di gruppo". Questi elementi, denotati come $D^{(\lambda)}_{\text{out,init}}(U) = \langle \text{out} | D^{(\lambda)}(U) | \text{init} \rangle$, sono necessari per compiti che spaziano dal calcolo delle ampiezze di transizione dei fotoni in reti ottiche lineari alla caratterizzazione di dispositivi quantistici e all'analisi delle proprietà di simmetria.

Sebbene la stima numerica sia spesso sufficiente per alcune applicazioni, altri contesti — come la derivazione di espressioni simboliche esatte per ampiezze di transizione o caratteri — richiedono una computazione simbolica precisa. Le soluzioni software esistenti affrontano problemi correlati ma distinti: alcune calcolano i coefficienti di Clebsch-Gordan (ad esempio, SUNRepresentations.jl), altre forniscono un accesso numerico indiretto alle funzioni di gruppo (ad esempio, RepLAB), e altre ancora si concentrano sull'integrazione simbolica su gruppi (ad esempio, IntegrateUnitary.jl). Le librerie dedicate all'ottica quantistica (ad esempio, BosonSampling.jl, The Walrus) tipicamente modellano i dispositivi numericamente o calcolano ampiezze specifiche come i permanenti, ma nessuna è progettata per calcolare simbolicamente le singole voci delle matrici di rappresentazione. Inoltre, l'assemblaggio di un'intera matrice di rappresentazione è spesso computazionalmente inefficiente quando è richiesta solo una singola voce.

Metodologia
Il documento introduce GroupFunctions.jl, una libreria Julia progettata per calcolare queste singole voci matriciali in modo simbolico o numerico. La libreria opera all'interno di una base dello spazio portante enumerata tramite pattern di Gelfand-Tsetlin (GT).

L'approccio algoritmico centrale si basa sulla formula di Grabmeier-Kerber, che esprime un elemento matriciale come una somma di monomi nelle entrate della matrice di input $U$. Questi monomi sono pesati dalla rappresentazione irriducibile e dagli specifici stati della base coinvolti. L'implementazione ottimizza questo processo tramite:

1. Enumerazione dei doppi cosetti: Invece di iterare su tutte le permutazioni nel gruppo simmetrico $S_N$ (dove $N$ è il numero di caselle nel diagramma di Young), l'algoritmo raggruppa le permutazioni che contribuiscono a monomi identici.
2. Sommatoria Efficiente: Costruisce doppi cosetti e itera sui loro elementi di permutazione rappresentativi, riducendo significativamente l'overhead computazionale rispetto a un approccio ingenuo permutazione per permutazione.
3. Rappresentazione Interna: Gli stati della base sono rappresentati internamente come tavole di Young semistandard, mentre l'interfaccia utente espone pattern GT equivalenti tramite la struttura dati GTPattern.

La libreria supporta matrici di input simboliche (utilizzando SymEngine.Basic) e fornisce funzioni di utilità per tradurre i pattern GT in ket di numero di occupazione, costruire unitarie di input da comuni parametrizzazioni dell'ottica quantistica e calcolare le associate funzioni di Schur.

Contributi Chiave

• Computazione Simbolica di Singole Voci: GroupFunctions.jl è il primo pacchetto, a conoscenza degli autori, capace di calcolare singole voci di matrici di rappresentazione di $U(d)$ in modo simbolico.
• Framework Unificato: Fornisce un metodo unificato per calcolare gli elementi matriciali per irreps specificate da partizioni intere di lunghezza al massimo $d$, gestendo sia input simbolici che numerici.
• Specializzazione per l'Ottica Quantistica: La libreria include strumenti specifici per scenari di ottica quantistica, come la valutazione delle ampiezze di transizione per il boson sampling e la gestione di irreps simmetriche (che modellano bosoni completamente indistinguibili).
• Capacità di Esportazione: I risultati possono essere esportati in un formato compatibile con Mathematica, facilitando ulteriori manipolazioni simboliche.

Risultati e Prestazioni
Il documento dimostra le capacità della libreria attraverso esempi, inclusa la computazione simbolica di funzioni di gruppo $SU(2)$ (che si riducono alle funzioni D di Wigner) e una irrep simmetrica di $U(10)$ coinvolgente 7 bosoni.

• Prestazioni: Su una CPU AMD Ryzen 7 PRO 4750U, la computazione simbolica di un elemento matriciale specifico di $U(10)$ ha richiesto circa 4 secondi dopo il riscaldamento della compilazione (compilation warmup), senza riutilizzo di risultati memorizzati nella cache.
• Scalabilità: Il costo computazionale è governato da tre fattori: la dimensione del gruppo $m$, la dimensione dello spazio di Gelfand-Tsetlin $d_\lambda(m)$ e la dimensione della rappresentazione del gruppo simmetrico $f^\lambda$ (il numero di tavole di Young standard). Gli autori osservano che, sebbene il lavoro sottostante comporti permutazioni che possono avvicinarsi a $N!$, il collo di bottiglia pratico è solitamente la crescita di $N$ o della dimensione $f^\lambda$, piuttosto che l'aumento di $m$ tramite zeri residui.
• Applicazione: La libreria è già stata utilizzata per valutare numericamente i caratteri di gruppo $SU(d)$ nel benchmarking randomizzato.

Significatività
La significatività di GroupFunctions.jl risiede nella sua capacità di fornire soluzioni simboliche esatte per le singole voci delle matrici di rappresentazione, una capacità precedentemente non disponibile negli ecosistemi software esistenti. Evitando la necessità di assemblare intere matrici di rappresentazione, la libreria offre una via più efficiente per calcoli specifici richiesti nella combinatoria algebrica e nella fisica quantistica. Gli autori pongono lo strumento come una risorsa specializzata per ricercatori che necessitano di espressioni analitiche esatte per ampiezze di transizione o valutazioni di caratteri, colmando il divario tra i simulatori ottici puramente numerici e le librerie di teoria dei gruppi a scopo generale che mancano di una diretta computazione simbolica delle singole voci.