Many-body post-processing of density functional calculations using the variational quantum eigensolver for Bader charge analysis

Titolo

Post-elaborazione many-body di calcoli DFT utilizzando il Variational Quantum Eigensolver (VQE) per l'analisi delle cariche di Bader.

1. Il Problema

La determinazione accurata della distribuzione della carica elettronica nei materiali è fondamentale per prevedere le loro proprietà industriali e fisiche (es. trasporto di carica, proprietà ottiche). Un metodo comune per quantificare questa distribuzione è l'analisi delle cariche di Bader, che si basano sulla densità elettronica nello spazio reale.

Tuttavia, l'accuratezza delle cariche di Bader dipende direttamente dalla precisione della funzione d'onda o della densità elettronica ottenuta risolvendo l'equazione di Schrödinger.

• Limiti della DFT standard: La Teoria del Funzionale Densità (DFT) basata su Kohn-Sham utilizza approssimazioni per l'interazione di scambio-correlazione (XC). Questo porta a errori di auto-interazione elettronica, causando una sovrade-localizzazione degli elettroni. Di conseguenza, la DFT standard fallisce spesso nel descrivere correttamente sistemi con forti correlazioni elettroniche (es. ossidi di metalli di transizione), dove le proprietà magnetiche e le bande di energia sono critiche.
• Limiti dei metodi many-body classici: Metodi più precisi come Coupled-Cluster o Configuration Interaction sono computazionalmente proibitivi per sistemi periodici di dimensioni realistiche a causa della loro scalatura esponenziale o polinomiale elevata.
• Sfida dei computer quantistici: Gli algoritmi puramente quantistici (es. Quantum Phase Estimation) richiedono risorse hardware attualmente non disponibili. Gli algoritmi ibridi come il Variational Quantum Eigensolver (VQE) sono promettenti ma devono essere integrati efficacemente nei flussi di lavoro esistenti della chimica computazionale.

2. Metodologia

Gli autori propongono un approccio di post-elaborazione che combina calcoli DFT classici con simulazioni quantistiche many-body tramite VQE. Il flusso di lavoro si articola come segue:

1. Calcolo DFT Iniziale: Vengono eseguiti calcoli DFT (utilizzando il codice Quantum ESPRESSO) per ottenere gli orbitali di Kohn-Sham e la struttura elettronica di base del sistema periodico.
2. Costruzione dell'Hamiltoniana Many-Body:
   • Gli orbitali di Kohn-Sham vengono utilizzati per costruire un'Hamiltoniana many-body in un spazio attivo selezionato (orbitale intorno all'energia di Fermi).
   • Vengono calcolati gli elementi di matrice a uno e due elettroni (cinetica, pseudopotenziali, interazione Coulombiana) trasformando le funzioni d'onda da basi di onde piane alla base degli orbitali di Kohn-Sham.
   • Viene applicata l'approssimazione del nucleo congelato (frozen core) per ridurre il numero di elettroni attivi, trattando gli elettroni di core come uno sfondo che modifica i potenziali effettivi.
3. Risoluzione con VQE:
   • L'Hamiltoniana ridotta viene mappata su qubit (tramite la mappatura di Jordan-Wigner).
   • Viene utilizzato l'algoritmo VQE con un ansatz UCCSD (Unitary Coupled Cluster con singoli e doppi eccitamenti) per trovare lo stato fondamentale dell'Hamiltoniana many-body.
   • La densità elettronica viene ricostruita a partire dallo stato fondamentale ottimizzato dal VQE.
4. Calcolo delle Cariche di Bader:
   • Sulla densità elettronica ottenuta dal VQE (o dalla DFT standard per confronto) viene eseguita un'analisi di Bader per assegnare le cariche atomiche.
5. Software: L'intero framework è implementato nel pacchetto open-source Dopyqo.

3. Contributi Chiave

• Framework Ibrido Dopyqo: Sviluppo e rilascio di un software open-source che collega i codici DFT classici (Quantum ESPRESSO) ai simulatori quantistici (TenCirChem/PySCF) per calcolare proprietà many-body in sistemi periodici.
• Metodologia di Post-Elaborazione: Dimostrazione che è possibile calcolare cariche di Bader accurate partendo da orbitali DFT e correggendoli con metodi many-body su uno spazio attivo, senza dover riscrivere l'intero codice DFT per i computer quantistici.
• Gestione della Correlazione Forte: Applicazione riuscita del VQE a sistemi di ossidi di metalli di transizione, dove la DFT standard fallisce, mostrando un miglioramento significativo rispetto ai risultati DFT+U (che richiedono la calibrazione manuale del parametro U).
• Validazione su Sistemi Reali: Test estesi su supercelle di $MgH_2$ (sistemi debolmente correlati) e vari ossidi di metalli di transizione ($TiO_2$, $CrO_2$, $RuO_2$, ecc.) e disolfuri.

4. Risultati

Lo studio confronta le cariche di Bader eccessive (BEC) ottenute con tre approcci: DFT standard (NC), DFT+U (riferimento per sistemi correlati) e il metodo proposto (Dopyqo/VQE).

• Sistemi Debolmente Correlati ($MgH_2$):
  • Per le supercelle di $MgH_2$ drogate con metalli di transizione (Ti, Cr, Fe, Ni, Zn), sia la DFT standard che il metodo Dopyqo producono risultati molto vicini ai dati di riferimento.
  • Questo conferma che il metodo VQE è accurato anche quando non è strettamente necessario (sistemi dove la DFT funziona bene), validando l'implementazione del software.
• Sistemi Fortemente Correlati (Ossidi di Metalli di Transizione):
  • Per sistemi come $CrO_2$, $RuO_2$, $RhO_2$ e $TiS_2$, la DFT standard sottostima significativamente le cariche di Bader rispetto ai dati di riferimento (DFT+U).
  • Il metodo Dopyqo/VQE riesce a recuperare gran parte di questa discrepanza, spostando i valori delle cariche verso quelli ottenuti con DFT+U.
  • Ad esempio, per $CrO_2$, la differenza tra DFT standard e riferimento è di 0.49, mentre con Dopyqo si riduce a 0.11.
  • Il metodo riesce a catturare gli effetti di correlazione forte senza la necessità di calibrare il parametro di Hubbard $U$, che è spesso difficile da determinare e non trasferibile tra diversi sistemi.

5. Significato e Implicazioni

• Superamento dei Limiti della DFT: Il lavoro dimostra che l'uso di computer quantistici (o simulatori quantistici) tramite VQE può correggere i difetti intrinseci della DFT (come la sovrade-localizzazione) in sistemi complessi, offrendo una descrizione più fedele della distribuzione di carica.
• Efficienza Operativa: Rispetto al DFT+U, il flusso di lavoro proposto è più diretto: non richiede la sintonizzazione manuale del parametro $U$ né la scelta arbitraria degli orbitali su cui applicare la correzione. Si basa su un calcolo DFT standard seguito da un post-processing many-body in uno spazio attivo selezionato.
• Scalabilità e Futuro: Sebbene i calcoli many-body siano attualmente costosi, questo approccio dimostra la fattibilità di utilizzare risorse quantistiche ibride per migliorare le proprietà predittive dei materiali. Il framework è pronto per essere esteso a calcoli DFT polarizzati in spin e mesh di k-point più densi.
• Accessibilità: La disponibilità del codice Dopyqo su GitHub permette alla comunità scientifica di replicare e estendere questi studi, accelerando lo sviluppo di materiali per applicazioni industriali (es. catalisi, batterie, materiali elettronici).

In sintesi, il paper rappresenta un passo significativo verso l'integrazione pratica degli algoritmi quantistici nella scienza dei materiali, offrendo un metodo robusto per analizzare le proprietà elettroniche di sistemi fortemente correlati che sfidano i metodi computazionali classici tradizionali.