Il metodo degli orbitali molecolari di Hückel è un metodo quanto meccanico semi-empirico usato per calcolare le proprietà degli orbitali molecolari degli elettroni π in sistemi d’idrocarburici coniugati, come, per esempio, l’etilene, il benzene o il butadiene. In questa serie di articoli, ho riassunto gli aspetti principali della teoria con esempi pratici di applicazioni e programmazione del metodo.

Il metodo fu ideato da Erich Hückel nel 1930 e, successivamente ampliato e migliorato, da altri studiosi. Il metodo ha aiutato a ottenere importanti risultati teorici nell’ambito dello studio delle proprietà di molecole organiche, e dei meccanismi delle loro reazioni chimiche, in periodo precedente allo sviluppo di calcolatori elettronici. Per esempio, ha fornito le basi teoriche della regola di Hückel per l’aromaticità di (4n + 2) π sistemi planari ciclici di elettroni.

Hückel ha ipotizzato che per in molecole organiche insature gli elettroni nei legami possano essere trattati separatamente da quelli coinvolti nei legami . La simmetria dei due tipi di orbitali molecolare in cui sono localizzati gli elettroni per i due tipi di legame giustifica questa assunzione. Infatti, un orbitale π è antisimmetrico per riflessione sul piano della molecola mentre un è simmetrico.

Con tale ipotesi, possiamo scrivere l’energia totale del sistema molecolare come . Ciò significa che la funzione d’onda della molecola è data dal prodotto della funzione d’onda che descrive gli e elettroni:

Possiamo anche supporre che la funzione d’onda molecolare dell’intero sistema possa essere approssimata da orbitali descritti come prodotto di funzioni d’onda con un elettrone

L’energia di un sistema può essere calcolata utilizzando l’operatore quantistico di Hamilton (Hamiltoniano) definito come

con

• : l’energia cinetica del nucleo.
• : l’energia cinetica degli elettroni.
• : l’energia potenziale di attrazione protone-elettrone.
• : l’energia potenziale di repulsione protone-protone.
• : l’energia di repulsione elettrone-elettrone.

Nel caso di una molecola composta da atomi M, l’hamiltoniano totale può essere scritto come

Per un sistema molecolare con elettroni π, possiamo ulteriormente distinguere questi elettroni da quelli come

Per gli elettroni non d’interazione, possiamo anche supporre che l’energia totale del sistema possa essere calcolata da un operatore Hamiltoniano totale (vedi il mio blog sulla meccanica classica) può essere espressa come

e gli autostati di ciascun elettrone possono essere calcolati dalle equazioni di Schrödingen

Come ulteriore ipotesi, la funzione d’onda degli elettroni può essere fattorizzata come segue

Usando queste approssimazioni, l’energia elettronica totale degli elettroni può essere espressa come

Ora considereremo un semplice sistema molecolare (la molecola dell’allile) per mostrare come assegnare gli orbitali atomici degli atomi per costruire il determinante di Hückel. Una molecola allilica ha la formula strutturale con un gruppo metilico () attaccato ad un gruppo vinilico (). Il nome della molecola deriva dalla parola latina per aglio, Allium sativum, e fu scelto da Theodor Wertheim nel 1844 quando riuscì a isolarne un derivato dall’olio d’aglio.

Consideriamo i tre orbitali di ciascun carbonio descritti dalle funzioni d’onda . Poiché sono antisimmetrici per riflessione lungo il piano molecolare possiedono lobi (ellissi colorate) dello stesso segno nelle due parti del piano.

L’orbitale molecolare totale può essere considerato come una combinazione lineare di questi orbitali atomici.

Per trovare l’energia del sistema molecolare, utilizzeremo il metodo basato sul risultato del teorema variazionale.

Il teorema di variazione afferma che dato un sistema con un operatore hamiltoniano allora se  è una qualsiasi funzione liscia che soddisfa correttamente le condizioni al contorno del sistema, allora vale la relazione

dove  è il vero valore dell’autovalore di energia più basso di Questa espressione è anche chiamata quoziente di Rayleigh.

Il principio variazionale ci consente di calcolare un limite superiore per l’energia dello stato fondamentale trovando la funzione d’onda di prova  per la quale l’integrale è minimizzato. Il principio ci assicura che solo se la funzione è esatta e in tutti gli altri casi .

Il metodo variazionale lineare

Utilizzeremo il cosiddetto metodo della variazione lineare (o metodo di Rayleigh-Ritz) in cui la funzione di variazione lineare è una combinazione lineare di funzioni orbitali atomiche linearmente indipendenti che soddisfano le condizioni al contorno del problema ($\lattice \Psi_{trial}$). I coefficienti () della funzione d’onda di prova sono parametri da determinare minimizzando l’integrale variazionale.

Inserendo la funzione della molecola allile nel quoziente di Rayleigh, il numeratore diventa

dove è l’elemento della matrice hamiltoniana, e il denominatore

dove è l’elemento della matrice di sovrapposizione.

Possiamo quindi scrivere

o nella forma riorganizzata come

(1)

Ora dobbiamo trovare il valore dei coefficienti che minimizzano la funzione usando il metodo variazionale. Questo può essere ottenuto utilizzando la definizione di un minimo di una funzione, ovvero calcolando le equazioni

per tutte le funzioni i.

Differenziando entrambi i membri della (1) e utilizzando le proprietà di e e dei coefficienti , otteniamo infine la cosiddetta equazione secolare

Abbiamo quindi un insieme di k equazioni omogenee simultanee (equazioni secolari) in k incognite. Nel caso della molecola d’allile, il sistema è il seguente

Per una soluzione non banale (cioè per tutti i valori di i), il determinante della matrice secolare deve essere uguale a zero

L’integrale è chiamato integrale di Coulomb e fornisce l’energia media del elettrone nell’orbitale atomico i sperimentando un campo potenziale dovuto all’intero molecola. Si presume che siano uguali alla stessa quantità .

L’integrale è chiamato integrale di risonanza (o integrale di legame) e fornisce l’energia di sovrapposizione tra i diversi orbitali. La sovrapposizione di atomi collegati da legame contribuiscono alla sovrapposizione dei loro orbitali della stessa quantità mentre atomi non primi vicini danno un valore zero.

Pertanto possiamo impostare i valori dell’integrale nel determinante di Hückel come costanti e sulla base della connettività molecolare con valori

e poiché gli orbitali sono normalizzati, gli integrali di sovrapposizione sono dati dalla espressione

Con queste semplificazioni diventa il determinante di Hückel per la molecola dell’ allile è il seguente

Nella figura che segue, sono mostrati altri esempi di determinante di Hückel per semplici molecole insature

Nella seconda parte di questo articolo, impareremo come risolvere l’equazione caratteristica del determinate secolare per trovare l’energia degli orbitali molecolari.

Se hai trovato utile questo articolo, condividilo con i tuoi amici e non dimenticare di aggiungere un Like e d’iscriviti per ricevere aggiornamenti!

REFERENCES

• J.P. Lowe. Quantum Chemistry. 1993, Academic Press.