La teoria del funzionale della densità (DFT) si fonda sul principio che la densità elettronica \( n(\mathbf{r}) \), funzione di sole tre coordinate spaziali indipendentemente dal numero di elettroni \( N \), contenga tutta l'informazione necessaria per determinare le proprietà dello stato fondamentale di un sistema quantistico a molti corpi. Questo contrasta con il metodo tradizionale basato sulla funzione d'onda totale
\[
\Psi(\mathbf{r}_1, \mathbf{r}_2, ..., \mathbf{r}_N)
\]
che dipende da \( 3N \) variabili, rendendo il calcolo impraticabile per grandi sistemi dove \( N \) può raggiungere valori dell'ordine di \( 10^{23} \) elettroni [1]. La riduzione dimensionale da uno spazio a molteplici dimensioni a uno tridimensionale è la chiave che permette alla DFT di operare in condizioni computazionali realistiche.
Il meccanismo fondamentale alla base della DFT è costituito dai due teoremi di Hohenberg e Kohn, i quali stabiliscono che:
- Le proprietà dello stato fondamentale di un sistema multielettronico non degenere sono descritte da funzionali della sola densità monoelettronica.
- L’energia totale dello stato fondamentale può essere espressa come un funzionale della densità,
\[
E[n]
\]
e tale energia corrisponde al valore minimo di \( E[n] \) rispetto alle densità elettroniche per cui è definito.
Questo implica che se si conoscesse esattamente la forma funzionale \( E[n] \), sarebbe sufficiente minimizzarla rispetto a tutte le possibili densità ammissibili per trovare quella reale del sistema. La sfida pratica consiste dunque nel costruire approssimazioni affidabili per questo funzionale sconosciuto, mantenendo intatti i principi matematici sottostanti.
Il primo tentativo formale di rappresentare l'energia totale come funzionale della densità risale al modello di Thomas-Fermi sviluppato nel 1927, che forniva una stima dell’energia cinetica degli elettroni come funzione integrale della densità stessa. Il modello combinava questa energia cinetica con termini classici per le interazioni nucleo-elettrone ed elettrone-elettrone, anch’essi espressi in termini di densità elettronica.
La principale limitazione del modello Thomas-Fermi risiedeva nell’omissione dell’energia di scambio, predetta invece dalla teoria di Hartree-Fock. Tale effetto fu introdotto da Dirac nel 1930 mediante un funzionale specifico per lo scambio energetico. Tuttavia, la teoria Thomas-Fermi-Dirac rimaneva inaccurata per molte applicazioni, poiché la maggiore fonte di errore derivava dalla rappresentazione dell'energia cinetica, a cui si aggiungevano gli errori legati all'energia di scambio e alla mancata considerazione delle repulsioni interelettroniche [1].
Nel formalismo più rigoroso, sfruttando l'approssimazione di Born-Oppenheimer applicata ai nuclei, il problema quantistico degli elettroni in un potenziale statico esterno \( V(\mathbf{r}) \) viene descritto dall’Hamiltoniana a molti corpi:
\[
H = \left[ \sum_{i}^{N} -\frac{\hbar^{2}}{2m} \nabla_i^2 + \sum_{i}^{N} V(\vec{r}_i) + \sum_{i<j}^{N} U(\vec{r}_i, \vec{r}_j) \right]
\]
dove il primo termine rappresenta l'energia cinetica, il secondo il potenziale esterno e il terzo l'interazione coulombiana tra coppie di elettroni. Lo stato stazionario elettronico è definito dalla soluzione della equazione di Schrödinger a molti corpi:
\[
H \Psi = E \Psi
\]
con \( \Psi({\vec r}_1,...,{\vec r}_N) \). La complessità intrinseca nella dipendenza da molte variabili rende direttamente questo approccio computazionalmente proibitivo per grandi sistemi.
La DFT aggira questa difficoltà spostando l'attenzione dalla funzione d'onda completa alla sola densità elettronica, ed esprimendo tutte le quantità energetiche come funzionali di essa. Questa trasformazione riduce drasticamente la complessità numerica senza compromettere la capacità descrittiva dello stato fondamentale.
Il funzionale energetico \( E[n] \) può essere idealmente scomposto in varie componenti:
- Energia cinetica degli elettroni espressa come un funzionale della densità;
- Energia dovuta al potenziale esterno;
- Energia di interazione coulombiana tra gli elettroni;
- Termini aggiuntivi come l'energia di scambio e correlazione che incorporano gli effetti quantistici non classici.
L’approssimazione più comune è fornita dal formalismo Kohn-Sham, che introduce orbitali ausiliari non interagenti per calcolare una parte consistente dell’energia cinetica in modo esplicito mentre lascia al funzionale scambio-correlazione il compito più complesso.
La condizione matematica centrale è che il valore minimo del funzionale deve soddisfare:
\[
\frac{\delta E[n]}{\delta n(\mathbf{r})} = 0
\]
condizione stazionaria che porta alle equazioni auto-consistenti tipiche dei calcoli DFT. Questo meccanismo è ciò che consente ai metodi DFT di ottenere risultati affidabili pur con una complessità computazionale contenuta rispetto ai metodi tradizionali basati sulla funzione d’onda completa.
Nonostante il rigore concettuale dei teoremi fondamentali, la pratica applicativa della DFT dipende fortemente dalla qualità delle approssimazioni impiegate per il funzionale energia scambio-correlazione. Questi termini non sono noti esattamente e vengono modellati attraverso formule empiriche o semiempiriche.
Ciò implica che certi fenomeni elettronici fortemente correlati o situazioni fuori equilibrio possono risultare mal descritti se il funzionale adottato non cattura adeguatamente tali effetti. In particolare, sistemi con forti interazioni magnetiche o stati eccitati richiedono spesso trattamenti avanzati o metodi complementari ai classici schemi DFT.
Il passaggio dal calcolo diretto della funzione d'onda a quello della densità consente l’applicazione efficiente ai solidi cristallini e alle molecole complesse tramite software specializzati (ad esempio Quantum Espresso). Questa efficienza ha reso possibile l’esplorazione teorica predittiva “da primi principi” delle proprietà strutturali, delle energie strutturali e delle forze interatomiche senza ricorrere a parametri empirici specifici [2].
La teoria si presta inoltre ad integrazioni con altri metodi numerici come i metodi Monte Carlo quantistici (variazionale, diffusivo, basato sul path-integral) per affrontare problemi dove gli approcci standard risultano insufficienti o troppo semplificati [2].
---
In sintesi la teoria del funzionale della densità si basa sul principio meccanico-fisico secondo cui tutta l’informazione rilevante dello stato fondamentale può essere ricavata dalla sola densità elettronica tridimensionale. I teoremi fondamentali stabiliscono una relazione univoca fra questa funzione scalare e le proprietà energetiche globali del sistema. La complessità computazionale viene così drasticamente ridotta passando da una funzione d’onda ad alta dimensionalità a uno strumento analitico più semplice ma altrettanto completo nel dominio dello stato fondamentale. Le sfide restano legate all’elaborazione accurata dei termini scambio-correlazione e all’estensione verso sistemi più complessi o eccitati.
Sto generando il riassunto…