A QSVT-Based Quantum Jacobi Algorithm for Linear Systems with Application to the Poisson Equation

Prevedere come si muovono i fluidi è un pilastro dell'ingegneria moderna, dalla progettazione di ali di aeromobili efficienti alla modellazione dei modelli meteorologici. Per farlo, gli scienziati scompongono le complesse, vorticose equazioni che governano il flusso dei fluidi in una massiccia griglia di punti. In ogni punto, il computer deve risolvere una semplice relazione tra i suoi vicini. Quando queste relazioni vengono combinate, formano un gigantesco puzzle di equazioni lineari. Risolvere questo puzzle è la parte più costosa della simulazione, e spesso consuma la maggior parte del tempo di un supercomputer. Per decenni, gli ingegneri si sono affidati a metodi iterativi per risolvere questi puzzle, un processo simile al raffinamento di uno schizzo grezzo in un disegno dettagliato, regolando ripetutamente le linee finché non si assestano nella forma corretta. Sebbene i computer classici gestiscano bene questo compito, la dimensione enorme delle griglie richieste per simulazioni realistiche spinge ai limiti anche le macchine più potenti.

L'informatica quantistica offre una strada diversa, promettendo di gestire queste griglie massicce codificando le informazioni negli strani stati sovrapposti delle particelle quantistiche. Tuttavia, la maggior parte degli approcci quantistici esistenti cerca di risolvere questi puzzle dei fluidi invertendo l'intera struttura matematica in un colpo solo, un metodo che è spesso troppo rigido o dispendioso in termini di risorse per la natura specifica e passo dopo passo delle simulazioni dei fluidi. Un team di ricercatori ha sviluppato un nuovo algoritmo quantistico che rispetta lo spirito iterativo dell'ingegneria classica. Invece di costringere un computer quantistico a invertire l'intero sistema, hanno adattato il classico metodo di Jacobi — una tecnica semplice di regolazione ripetuta — per farlo girare su hardware quantistico. Il loro lavoro dimostra che, utilizzando uno strumento matematico specifico chiamato trasformazione del valore singolare quantistico, possono eseguire queste regolazioni ripetute in modo efficiente, mantenendo costante il numero di bit quantistici extra necessari indipendentemente dal numero di passi della fase di calcolo.

I ricercatori, lavorando con istituzioni tra cui Volkswagen e il Centro Aerospaziale Tedesco, si sono concentrati su un tipo specifico di problema comune nella dinamica dei fluidi: l'equazione di Poisson. Questa equazione appare ogni volta che gli scienziati devono garantire che il flusso di un fluido rimanga coerente, come ad esempio quando si calcolano i cambiamenti di pressione in una simulazione del movimento dell'aria sopra un'auto o all'interno di una scatola. In una simulazione standard su computer, la risoluzione di questa equazione comporta un ciclo in cui il computer ipotizza una soluzione, controlla l'errore e corregge l'ipotesi, ripetendo questo ciclo centinaia o migliaia di volte. Il team ha tradotto questo ciclo in un circuito quantistico. Si sono resi conto che ogni passo della correzione poteva essere visto come una trasformazione polinomiale, un'operazione matematica che rimodella i valori nel sistema. Utilizzando la trasformazione del valore singolare quantistico, potevano applicare questa rimodellazione direttamente allo stato quantistico che rappresenta il fluido, eseguendo efficacemente molti passi del processo di correzione in un unico, coerente' operazione quantistica.

Un ostacolo importante nei tentativi precedenti di eseguire metodi iterativi su computer quantistici era l'accumulo di risorse. Ogni volta che un computer quantistico moltiplica due matrici per eseguire un passo del calcolo, richiede tipicamente bit di memoria extra, noti come qubit ancillari, per contenere i risultati intermedi. Nelle vecchie versioni quantistiche del metodo di Jacobi, il numero di questi bit extra cresceva linearmente con il numero di passi, esaurendo rapidamente la capacità limitata degli attuali dispositivi quantistici. Il nuovo algoritmo risolve questo problema riformulando l'iterazione in modo che i bit extra rimangano costanti. I ricercatori hanno dimostrato che il loro metodo richiede solo un numero fisso e piccolo di qubit extra, indipendentemente dal fatto che la simulazione duri dieci o cento passi. Questa efficienza avviene a costo di una profondità del circuito che cresce linearmente con il numero di iterazioni, ma questo compromesso è favorevole perché evita l'esplosione esponenziale delle risorse che ha afflitto i tentativi precedenti.

Per testare la loro idea, il team ha simulato l'algoritmo su un computer classico per vedere come si sarebbe comportato su una vera macchina quantistica. Hanno applicato prima l'algoritmo a un semplice problema monodimensionale, un modello base della pressione di un fluido, e hanno scoperto che l'algoritmo quantistico riproduceva esattamente gli stessi risultati del metodo classico, con errori così piccoli da essere praticamente nulli. Sono poi passati a uno scenario bidimensionale più complesso: il flusso d'aria all'interno di una scatola quadrata dove il coperchio superiore si muove, trascinando con sé l'aria. Questo è un caso di test standard in ingegneria noto come "lid-driven cavity flow" (flusso di cavità guidato dal coperchio). In questa simulazione, l'equazione della pressione deve essere risolta ripetutamente mentre il fluido si muove. I ricercatori hanno dimostrato che il loro approccio quantistico poteva gestire i calcoli della pressione per questo fluido in movimento, eguagliando la soluzione classica con alta precisione. Hanno anche affrontato una difficoltà tecnica in cui la struttura matematica delle equazioni dei fluidi non si adattava perfettamente ai requisiti del loro strumento quantistico. Hanno scoperto che, modificando leggermente le condizioni al contorno del problema, potevano rendere le equazioni compatibili senza influenzare significativamente il risultato fisico finale, un compromesso che ha retto bene nelle loro simulazioni.

Lo studio conferma che questo approccio basato sui polinomi è un mattone vitale per le future simulazioni quantistiche dei fluidi. Esso allinea l'informatica quantistica più da vicino alle strategie effettivamente utilizzate oggi dagli ingegneri, piuttosto che forzare i problemi dei fluidi in uno stampo progettato per la pura matematica. Sebbene l'attuale lavoro sia una simulazione e si basi sull'assunto che le equazioni dei fluidi possano essere strutturate in un modo specifico, i risultati suggeriscono una chiara strada da seguire. Gli autori osservano che questo metodo potrebbe eventualmente servire come componente in flussi di lavoro quantistici più ampi e avanzati, come quelli utilizzati per i metodi multigrid, che rappresentano il gold standard per velocizzare le simulazioni dei fluidi sui computer classici. Dimostrando che un semplice processo di regolazione ripetuta può essere codificato efficientemente in un circuito quantistico, i ricercatori hanno fornito un passo pratico verso l'integrazione dell'informatica quantistica nel kit di strumenti quotidiano della fluidodinamica computazionale.