Come IBM ha simulato la produzione di particelle su un processore da 104 qubit

Non una collisione nel chip, ma una simulazione quantistica

Roland C. Farrell, Nikita A. Zemlevskiy, Marc Illa e John Preskill hanno usato il processore ibm_marrakesh di IBM, con 104 qubit, per eseguire una simulazione digitale in tempo reale di una collisione nella teoria di campo di Ising in una dimensione. Hanno preparato due pacchetti d’onda di particelle leggere, inizialmente distanti e diretti uno contro l’altro, quindi hanno seguito l’evoluzione del sistema anche dopo l’urto. Per raggiungere questa fase hanno impiegato fino a 5.589 porte quantistiche a due qubit. 7

Non significa che particelle fisiche abbiano urtato all’interno del processore. Il chip rappresenta invece, qubit per qubit, la dinamica quantistica di un modello semplificato di teoria dei campi. Lo scopo è affrontare un problema difficile per i calcolatori classici: descrivere l’evoluzione in tempo reale di molti gradi di libertà quantistici durante un processo di scattering.

Come è stata ricostruita la collisione

La simulazione si è svolta in tre passaggi principali:

1. Preparazione delle particelle in arrivo. I ricercatori hanno creato due pacchetti d’onda localizzati di particelle leggere, con momenti opposti e quindi in avvicinamento.
2. Evoluzione temporale. La dinamica della teoria è stata approssimata con circuiti quantistici trotterizzati: una successione di piccoli passi temporali che descrive l’avvicinamento, la sovrapposizione e la successiva separazione dei pacchetti.
3. Misura dello stato finale. Dopo la collisione è stata misurata la densità locale di energia, cioè come l’energia risultava distribuita lungo il reticolo simulato. 7

Il punto decisivo è che i dati non sono compatibili soltanto con due particelle leggere che si allontanano dopo un urto elastico. Indicano invece evidenza di produzione anelastica di particelle: un canale finale composto da una particella leggera e da un’eccitazione più pesante e più lenta. 7

Il segnale della particella pesante: l’asimmetria dell’energia

La particella pesante non è stata identificata tramite un’immagine di rivelatore nel senso tradizionale della fisica delle alte energie. Il suo indizio è arrivato dalla forma del profilo di densità di energia misurato dopo la collisione.

In uno scenario relativamente elastico, con due particelle leggere in uscita, ci si aspetta una distribuzione più simmetrica. Se invece una delle particelle finali è pesante, si muove più lentamente e lascia una firma asimmetrica nel profilo energetico. 7

Per quantificare tale asimmetria, il gruppo ha usato la skewness — in italiano, coefficiente di asimmetria — una misura statistica che segnala se una distribuzione presenta una coda più marcata su uno dei due lati. Una skewness diversa da zero, coerente con la traiettoria attesa per la particella pesante, è stata interpretata come la firma del processo anelastico. 7

In questa lettura, una parte dell’energia cinetica iniziale delle due particelle leggere è stata convertita in massa-energia di una nuova eccitazione pesante, mentre l’altra particella in uscita rimane leggera. È un risultato presentato come evidenza del canale anelastico, non come una misura di precisione paragonabile a quelle di un acceleratore. 7

Perché gli stati W rendono più praticabile la simulazione

Preparare correttamente le particelle iniziali è uno degli ostacoli centrali nelle simulazioni di scattering. Per ciascun pacchetto d’onda, gli autori hanno prima costruito uno stato simile a ( |W(k_0)\rangle ), progettato per fissarne inviluppo spaziale, contenuto in quantità di moto e numeri quantici. Poiché lo stato iniziale può contenere componenti indesiderate a più particelle, è stato poi raffinato con circuiti variazionali che preservano le simmetrie e ne abbassano l’energia verso l’autostato a una particella desiderato. 7

I circuiti per i pacchetti d’onda della teoria di Ising su reticoli grandi sono stati progettati inizialmente tramite simulazioni classiche con stati a prodotto di matrici, prima dell’esecuzione sull’hardware a 104 qubit. 7

La novità algoritmica è l’adattamento di tecniche di preparazione degli stati W mediante:

• misure a metà circuito, che producono risultati classici prima della fine del programma quantistico;
• feed-forward classico, che usa quei risultati per scegliere condizionatamente le operazioni successive.

Questa combinazione consente un protocollo senza ancille in cui la profondità del circuito non dipende dalla dimensione del pacchetto d’onda né dalla dimensionalità spaziale, salvo la dipendenza dalla lunghezza di correlazione fisica. Secondo gli autori, ciò rappresenta un miglioramento superesponenziale rispetto a metodi precedenti che richiedevano un pre-processing classico molto oneroso oppure circuiti la cui profondità cresceva polynomialmente con la dimensione del pacchetto. 27

Dall’Ising alla QCD: perché il risultato interessa

La produzione di nuove particelle è un elemento fondamentale delle collisioni ad alta energia. L’obiettivo di lungo periodo è simulare teorie di campo fortemente accoppiate in tempo reale, dove le collisioni producono particelle e comportamento collettivo. È il tipo di problema che, su scale realistiche, è legato alla produzione di adroni e alla cromodinamica quantistica, o QCD, la teoria delle interazioni forti.

La dimostrazione non riguarda ancora la QCD, ma il metodo di preparazione dei pacchetti d’onda non è limitato al caso mostrato su hardware. Gli autori lo hanno costruito o dimostrato anche per una teoria di campo scalare unidimensionale, per il modello di Schwinger e per la teoria di campo di Ising bidimensionale. 27

I limiti da tenere presenti

Il risultato è un passo concreto, ma non equivale alla simulazione di un evento in un grande collisore.

• Il modello studiato è la teoria di campo di Ising in una dimensione, non la QCD, non collisioni adroniche tridimensionali e non getti di particelle realistici. 7
• L’esperimento ha fornito evidenza di un canale anelastico, non una previsione di scattering sistematicamente precisa e protetta da correzione d’errore completa. 7
• Il processore era rumoroso e il lavoro ha fatto ricorso a tecniche di mitigazione degli errori; questo è diverso da un calcolo pienamente fault-tolerant, cioè eseguito con correzione quantistica degli errori su larga scala. 78
• Per arrivare a teorie di gauge realistiche serviranno codifiche in dimensioni superiori, evoluzioni temporali più lunghe, preparazione scalabile degli stati iniziali e misure precise di osservabili con molte particelle.

Il traguardo immediato, quindi, è circoscritto ma significativo: un processore quantistico ha seguito uno stato a due particelle attraverso una collisione e ha prodotto un segnale misurabile compatibile con la comparsa di una particella più pesante. Questo rende più credibile l’idea di studiare in futuro, con hardware quantistico più potente, processi di scattering in tempo reale sempre più complessi. 7