Constraint-Aware Quantum Optimization of Defect Configurations in Doped ZrO2: XY-Mixer QAOA and Grover Adaptive Search

Sintesi Tecnica: Ottimizzazione Quantistica Sensibile ai Vincoli delle Configurazioni di Difetto in ZrO2 Drogato

Enunciato del Problema
Il documento affronta la sfida combinatoria dell'ottimizzazione delle configurazioni di difetto nella progettazione dei materiali, specificamente all'interno della zirconia drogata con gadolinio (Gd-ZrO2), un materiale critico per i rivestimenti barriera termica (TBC). L'obiettivo è identificare disposizioni a bassa energia di sostituzioni cationiche (Gd che sostituisce Zr) e vacanze di ossigeno all'interno di una supercella a composizione fissa. Sebbene gli interpotenziali apprenditi tramite machine learning (MLIP) come MACE permettano una rapida valutazione dell'energia, la mappatura di queste configurazioni discrete e vincolate in un flusso di lavoro di ottimizzazione quantistica richiede l'codifica del problema come un modello di Ottimizzazione Binaria Quadratica Non Vincolata (QUBO), gestendo rigorosamente i vincoli di composizione fissa (esattamente due sostituzioni di Gd e una vacanza di ossigeno). Gli autori sostengono che gli approcci quantistici standard "black-box" spesso non riescono a rispettare efficientemente questi vincoli fisici, rendendo necessaria una strategia di progettazione sensibile ai vincoli.

Metodologia
Gli autori hanno sviluppato un flusso di lavoro end-to-end validato rispetto all'enumerazione classica esatta di 448 configurazioni ammissibili (inserite in uno spazio di $2^{24}$). La metodologia consiste in tre componenti primarie:

1. Formulazione QUBO e Modellazione Surrogata:

   • Un modello QUBO a 24 variabili è stato adattato ai dati energetici MACE-MPA-0 (8 variabili cationiche, 16 variabili di vacanza di ossigeno).
   • Il modello include termini lineari per i contributi dipendenti dal sito e termini quadratici per le interazioni a coppie (catione-catione e catione-vacanza). I termini vacanza-vacanza sono stati omessi a causa del vincolo di singola vacanza.
   • La regressione Ridge è stata utilizzata per adattare 180 feature (24 lineari, 156 quadratiche), ottenendo un $R^2$ di 0 even 0.997 sul set di test e un RMSE di 27,6 meV.

2. Approccio Variazionale Near-Term (QAOA):

   • Due strategie QAOA sono state messe a confronto: una Penalty QAOA standard (aggiunta di penalità di vincolo all'Hamiltoniana del costo) e una QAOA con Mixer XY-Preserving dei Vincoli.
   • L'approccio con mixer XY utilizza mixer a peso di Hamming fisso che operano separatamente sui settori cationici e di vacanza. Questa costruzione confina matematicamente l'evoluzione quantistica all'interno del sottospazio ammissibile a 448 dimensioni, garantendo la fattibilità senza necessità di calibrazione delle penalità.

3. Ricerca della Soglia Fault-Tolerant (Grover Adaptive Search - GAS):

   • Un oracle GAS vincolato è stato esplicitamente costruito e validato. A differenza delle astrazioni di alto livello, ciò ha comportato la costruzione di circuiti aritmetici reversibili per la valutazione dell'energia a punto fisso.
   • L'oracle include: mappatura dei flag di termine, accumulo pesato positivo/negativo, confronto sicuro tra rami (gestione di interi con segno tramite logica senza segno), controllo di fattibilità e phase kickback.
   • Gli autori hanno evitato l'aritmetica in virgola mobile, convertendo i coefficienti in interi a punto fisso a 18 bit per garantire un confronto reversibile esatto.

Risultati Chiave

• Fedeltà del Surrogato: Il surrogato QUBO riproduce accuratamente il panorama energetico MACE ($R^2 = 0.999$ sui dati completi), sebbene gli autori notino che il manifold a bassa energia è quasi degenere. Il minimo globale del QUBO differisce leggermente dallo stato fondamentale MACE, suggerendo un flusso di lavoro ibrido in cui il campionamento quantistico è seguito da un re-ranking tramite MACE.
• Prestazioni QAOA:
  • Penalty QAOA: Non è riuscita a produrre campioni ammissibili in modo affidabile; in 6 casi su 12 testati (variando profondità e forza della penalità), la probabilità di trovare uno stato ammissibile era zero. Il miglior caso ammissibile ha prodotto solo il 7% di probabilità.
  • XY-Mixer QAOA: Ha raggiunto il 100% di fattibilità per costruzione. Alla profondità $p=3$, ha concentrato l'85,9% della massa di probabilità entro 1 meV dall'optimum MACE e l'86,1% entro le prime 50 configurazioni MACE.
  • Analisi del Rumore: Un'analisi del fattore di sopravvivenza del primo ordine ha indicato che il conteggio delle porte a due qubit è la fonte dominante di errore. Con una fedeltà a due qubit del 99,99%, il circuito $p=2$ mantiene un'alta probabilità di successo, mentre al 99,9%, il circuito più superficiale $p=1$ è preferibile.
• Stima delle Risorse GAS:
  • L'oracle di fase vincolato richiede 324 qubit logici di alto livello (352–358 considerando un conteggio conservativo degli ancilla puliti).
  • Il costo delle porte Toffoli per iterazione varia da 35.799 a 43.175, a seconda del caso di soglia (numero di stati marcati).
  • Il costo è guidato dalla precisione dell'aritmetica a punto fisso e dal numero di termini QUBO, non semplicemente dal numero di variabili.

Significatività e Reclami
L'articolo sostiene che la consapevolezza dei vincoli sia la scelta di progettazione decisiva per l'ottimizzazione quantistica nella scienza dei materiali, applicabile sia ai regimi NISQ che a quelli fault-tolerant (FTQC).

• Per NISQ: Lo studio dimostra che codificare i vincoli direttamente nel mixer (XY-mixer) è vastamente superiore ai metodi basati su penalità, che spesso falliscono nel trovare stati fattibili. Ciò consente all'algoritmo di concentrare interamente le risorse computazionali sullo spazio fisicamente rilevante.
• Per il Calcolo Fault-Tolerant: Il lavoro va oltre il trattare la ricerca adattiva di Grover come una black-box. Costruendo e validando i livelli aritmetici reversibili, gli autori forniscono una stima delle risorse concreta e verificabile. Dimostrano che il costo dell'oracle è governato dalla complessità dell'aritmetica (precisione a punto fisso, accumulo dei termini) piuttosto che dal solo conteggio delle variabili.
• Amplificazione dello Spazio Ammissibile: Gli autori presentano un limite superiore idealizzato mostrando che limitare l'amplificazione di Grover allo spazio ammissibile di 448 stati (invece del pieno spazio di $2^{24}$) potrebbe produrre una riduzione di circa $\sim 240\times$ del costo totale delle porte Toffoli. Dichiarano esplicitamente che questo è un obiettivo teorico che richiede un diffusore che preservi i vincoli, il quale non è ancora stato implementato.

Conclusione
L'articolo stabilisce un ponte validato tra la modellazione QUBO informata dai materiali e l'ottimizzazione quantistica. Fornisce una stima calibrata delle risorse per un problema specifico della scienza dei materiali e sostiene che i futuri sforzi di scalabilità debbano dare priorità ad algoritmi che preservino i vincoli e a modelli surrogati sparsi e locali per mantenere l'efficienza. Il flusso di lavoro è progettato per far parte di un ciclo ibrido in cui il campionamento quantistico propone candidati, che vengono poi raffinati da MLIP o calcoli primi-principi.