Capitolo 7, "Stati, livelli energetici e transizioni"

Capitolo 7
Stati, livelli energetici e
transizioni.
7.1
Separazione dei moti molecolari.
Come in meccanica classica, anche in meccanica quantistica lo studio della dinamica
di un insieme di particelle interagenti, quale una molecola, si semplifica se `e possibile trovare un sistema di riferimento mobile, tale da evidenziare l’indipendenza
dei moti traslazionale e rotazionale complessivi da quelli interni. La nozione stessa
di molecola, come insieme di elettroni in moto intorno a nuclei con posizioni relative quasi fisse, richiede di separare le coordinate interne da quelle di posizione ed
orientazione. Questo argomento precede logicamente lo studio dell’interazione dei
moti interni nucleari e di quelli elettronici in una molecola. Per un approfondimento
indichiamo alcuni testi fondamentali: Wilson, Decius e Cross (1955), Allen e Cross
(1973), Califano (1976), Howard e Moss (1970 e 1971).
Per quanto riguarda il moto traslazionale complessivo, `e possibile mostrare che esso
si separa esattamente dalle altre coordinate. Ad esempio, per due sole particelle,
possiamo far uso delle coordinate del baricentro:
~ CM =
R
m2
m1
~r1 +
~r2
m1 + m2
m1 + m2
(7.1)
e del vettore che esprime la posizione relativa,
~r12 = ~r2 − ~r1
(7.2)
In assenza di campi esterni, l’energia potenziale dipende solo da ~r12 . L’energia
cinetica, nel nuovo sistema di coordinate, risulta una somma di due termini, ciascuno
159
~ CM oppure di ~r12 :
dei quali contiene derivate seconde rispetto alle componenti di R
!
!
2
2
2
2
2
2
2
2
h
¯
∂
∂
∂
∂
h
¯
∂
∂
Tˆ = −
+
+
−
+ 2 + 2
(7.3)
2
2
2
2
2M ∂RCM
2µ ∂r12x
∂RCM
∂RCM
∂r12y
∂r12z
x
y
z
dove M = m1 + m2 `e la massa totale e µ = m1 m2 /(m1 + m2 ) `e la massa ridotta. Ne
segue che l’hamiltoniano si pu`
o separare in due termini, ciascuno dei quali dipende
da (ed agisce su) un certo set di coordinate; le autofunzioni sono quindi prodotti
di un fattore traslazionale e di uno rotazionale-interno; i livelli energetici sono somma di una energia traslazionale (spettro continuo, da zero a +∞), e di un’energia
rotazionale-interna. Per sistemi di N particelle (N > 2), analogamente, possiamo
separare il moto del baricentro da quello di altre 3N − 3 coordinate relative.
La separazione del moto rotazionale `e questione pi`
u complessa. Per le molecole
biatomiche `e possibile definire un sistema ruotante, con un asse coincidente con la
congiungente dei nuclei, di modo che il momento angolare nucleare in quel sistema
si annulla. Nel caso di molecole poliatomiche, non siamo in grado di definire rigorosamente un sistema di assi ruotanti “ancorato” alla molecola, perch´e i nuclei si
muovono l’uno rispetto all’altro (per non dire degli elettroni!). Si ottiene una separazione del moto rotazionale soltanto mediante l’approssimazione di considerare i
nuclei connessi rigidamente uno all’altro, e trascurando la massa degli elettroni. Se
i nuclei compiono piccole oscillazioni intorno alle posizioni di equilibrio, l’hamiltoniano risulta somma di un termine rotazionale, uno interno (che dipende da 3N − 6
coordinate relative) e termini misti non troppo grandi che accoppiano il moto elettronico e quello nucleare interno con la rotazione complessiva. Avremo allora, con
buona approssimazione:
Ψ(X) = Ψtrasl (RCM ) Ψrot (α, β, γ) Ψint (q1 , q2 . . . q3N −6 )
(7.4)
dove le q sono coordinate interne. Inoltre:
Etot = Etrasl + Erot + Eint
(7.5)
Per molecole non rigide, oppure nella trattazione di reazioni chimiche, `e poco
giustificato separare la rotazione dalle deformazioni della molecola.
La separazione pi`
u o meno rigorosa dei moti traslazionali e rotazionali da quelli interni `e fondata su una trasformazione di coordinate, che permetta di scrivere l’operatore
di energia cinetica come somma di termini, dipendenti da variabili diverse. In ogni
caso, l’energia potenziale non dipende dalle coordinate di traslazione o rotazione
complessiva, quindi `e a priori “separabile”. Il passo successivo nella semplificazione e interpretazione delle funzioni d’onda molecolari consiste nella separazione del
moto elettronico da quelli nucleari. In questo caso, per`
o, abbiamo un contributo
importante all’energia potenziale, l’attrazione coulombiana nucleo-elettrone, che dipende contemporaneamente da coordinate di nuclei ed elettroni. Il procedimento da
seguire `e perci`
o sostanzialmente diverso.
160
D’ora in poi indichiamo con Q l’insieme delle coordinate dei nuclei, che potrebbero
essere coordinate cartesiane riferite ad assi “del laboratorio”, oppure coordinate
interne, definite dopo aver separato i moti traslazionali e rotazionali. Le coordinate
degli elettroni sono indicate con q. La funzione d’onda molecolare Ψ(q, Q) pu`
o
essere quindi intesa come una funzione d’onda complessiva, oppure come il fattore
Ψint nell’eq. (7.4). L’approssimazione di Born-Oppenheimer o approssimazione
adiabatica (Born e Oppenheimer 1927) consiste nell’assumere che Ψ(q, Q) si possa
scrivere come il prodotto di una funzione d’onda per gli elettroni ed una per i nuclei,
e che ciascuna delle due sia soluzione di una speciale equazione di Schr¨
odinger, come
vedremo nel seguito:
Ψ(q, Q) = ψ(q; Q) χ(Q)
(7.6)
Qui si intende che |χ(Q)|2 `e la densit`a di probabilit`a di trovare i nuclei nelle posizioni corrispondenti a Q1 . . . Qr . . .; ci`o implica che la funzione d’onda elettronica sia
normalizzata per ogni scelta delle Q, nel senso che:
Z
X
|ψ(q; Q)|2 dq1x dq1 dq2 . . . dq3n = hψ(Q) |ψ(Q) i = 1
(7.7)
R
Qui il simbolo Σ indica al contempo somma sulle coordinate di spin ed integrale
su quelle spaziali; d’ora in poi, per chiarezza, i simboli h e i, che indicano la stessa
operazione, saranno riservati alle sole coordinate degli elettroni. Bench´e ψ(q; Q)
dipenda anch’essa dalle Q, non ha nulla a che vedere con la probabilit`a di localizzare i nuclei: |ψ(q; Q)|2 rappresenta solo la probabilit`a di trovare gli n elettroni in
q1 , q2 , . . . q3n , una volta stabilito che i nuclei si trovano in Q1 . . . Qr . . .. ψ(q; Q) `e la
funzione d’onda elettronica “a nuclei fissi”.
Tutti gli integrali che coinvolgono grandezze fisiche dipendenti dalle coordinate Q e
q, si valutano integrando prima sulle coordinate elettroniche e poi su quelle nucleari:
Z
X
ˆ ψl (q; Q)χv (Q) . . . dqi . . . dQr . . . =
(7.8)
ψk∗ (q; Q)χ∗u (Q) O
=
dove
Z
X
χ∗u (Q) Okl (Q) χv (Q) . . . dQr . . .
Okl (Q) =
Z
X
D
E
ˆ ψl (q; Q) . . . dqi . . . = ψk (Q) O
ˆ ψl (Q)
ψk∗ (q; Q) O
(7.9)
`e una funzione delle sole Q o un operatore che agisce solo su χv (Q).
L’hamiltoniano molecolare ha la forma:
ˆ=−
H
X h
h2 X ∂ 2
¯
¯ 2 ∂2
−
+ V (q, Q)
2µr ∂Q2r
2m
∂qi2
r
i
161
(7.10)
Se intendiamo le Q come coordinate cartesiane, µr `e la massa del nucleo associato
con la coordinata Qr ; se invece le Qr sono coordinate interne, µr `e una massa ridotta; nel secondo caso, l’hamiltoniano (7.10) `e un’espressione semplificata, priva dei
termini di accoppiamento tra rotazione e moti interni. Rispetto alla funzione d’onda ψ(q; Q)χ(Q), il termine di energia potenziale V (q, Q) `e puramente moltiplicativo,
∂2
∂2
mentre ∂q
2 agisce solo su ψ; ∂Q2 invece agisce su entrambi i fattori:
r
i
∂2
∂2ψ
∂ψ ∂χ
∂2χ
ψ(q;
Q)χ(Q)
=
χ
+
2
+
ψ
∂Q2r
∂Q2r
∂Qr ∂Qr
∂Q2r
(7.11)
ˆ alla funzione d’onda ψ(q; Q)χ(Q) si ottiene:
Cos`ı, applicando H
ˆ ψ(q; Q)χ(Q) =
H
(7.12)
#
"
#
X ¯h2 ∂ 2
X h
¯ 2 ∂2
ψ(q; Q) − ψ(q; Q)
= χ(Q) V (q, Q) −
χ(Q)−
2m ∂qi2
2µr ∂Q2r
r
i
X ¯
∂2ψ
∂ψ ∂χ
h2
+χ 2
2
−
2µr ∂Qr ∂Qr
∂Qr
r
"
ˆ induce
La richiesta che il prodotto ψχ approssimi al meglio un autostato di H
a definire ψ e χ mediante due condizioni. In prima istanza si impone che ψ sia
ˆ l’hamiltoniano elettronico H
ˆ el :
autofunzione di una parte di H,
¯ 2 X ∂2
ˆ el = V (q, Q) − h
H
2m
∂qi2
(7.13)
ˆ el ψk = Uk (Q)ψk
H
(7.14)
i
ˆ el dipende dalle coordinate nucleari, anche se `e un operatore puramente moltiplicaH
tivo rispetto ad esse; perci`
o gli autovalori Uk e le autofunzioni ψk devono dipendere
ˆ el e quindi le energie
a loro volta dalle Qr . Notiamo che l’hamiltoniano elettronico H
Uk convenzionalmente includono anche i termini di repulsione nucleo-nucleo. Coerentemente col significato attribuito alla funzione d’onda elettronica, la (7.14) si
pu`
o interpretare come un’equazione di Schr¨
odinger per gli elettroni, “a nuclei fissi”. L’indice k viene introdotto perch´e dovremo distinguere molte diverse soluzioni
dell’equazione, corrispondenti agli stati elettronici fondamentale ed eccitati.
L’equazione (7.12) si semplifica in:
ˆ ψk (q; Q)χ(Q) =
H
(7.15)
#
X ¯h2 ∂ψk ∂χ
X h
∂ 2 ψk
¯ 2 ∂2
χ(Q) −
+χ
2
= ψk (q; Q) Uk (Q) −
2µr ∂Q2r
2µr ∂Qr ∂Qr
∂Q2r
r
r
"
162
La seconda condizione che viene imposta riguarda la χ(Q); si definisce un
hamiltoniano per il moto dei nuclei, il cosiddetto hamiltoniano “vibrazionale”:
ˆ vib = Uk (Q) −
H
X h
¯ 2 ∂2
2µr ∂Q2r
r
(7.16)
dove Uk (Q), “energia elettronica”, gioca il ruolo di un’energia potenziale efficace per
il moto nucleare. Si impone che le funzioni d’onda nucleari o “vibrazionali” siano
ˆ vib . Poich´e il potenziale Uk (Q) dipende dallo stato elettronico,
autofunzioni di H
anche le funzioni vibrazionali dovranno portare un indice k e saranno diverse a
seconda dello stato elettronico cui appartengono:
ˆ vib χku = Eku χku
H
(7.17)
A questo punto, l’equazione (7.12) assume una forma ancora pi`
u semplice:
ˆ ψk (q; Q)χku (Q) =
H
(7.18)
X h
¯2
∂ 2 ψk
∂ψk ∂χku
= Eku ψk (q; Q)χku (Q) −
+ χku
2
2µr ∂Qr ∂Qr
∂Q2r
r
Fin qui, non `e stata introdotta alcuna approssimazione; tuttavia, affermare che
ψk χku `e autofunzione dell’hamiltoniano, significa trascurare alcuni termini nell’ultima equazione, precisamente quelli che implicano derivate di ψk (q; Q) rispetto alle
coordinate dei nuclei. Rimandiamo l’esame della validit`
a dell’approssimazione di
B.-O. ad un’altra sezione e vediamo ora quale modello della molecola e dei fenomeni
chimici ne consegue.
Per ogni stato elettronico abbiamo una diversa superficie di energia potenziale,
Uk (Q), nella quale si potranno identificare uno o pi`
u punti di minimo relativo. Consideriamo il moto dei nuclei nelle vicinanze di uno dei minimi del potenziale, Umin .
La situazione `e del tutto analoga a quella degli stati legati in una dimensione; avremo
quindi uno spettro discreto di autovalori Eku , derivanti dalla soluzione dell’eq.(7.17).
Tutti gli autovalori saranno maggiori di Umin , perci`
o si potr`
a formalmente scrivere:
E = Eel + Evib + Erot + Etrasl
(7.19)
dove Eel = Umin e Evib = Eku − Umin . Le energie elettroniche dipendono dall’equazione di Schr¨
odinger a nuclei fissi (7.14) e quelle vibrazionali dalla (7.17), che a sua
volta contiene il potenziale elettronico Uk (Q). Finch`e i nuclei compiono solo piccole
oscillazioni intorno alle posizioni di equilibrio, i minimi relativi di Uk (Q) determinano la forma delle molecole, cio`e gran parte della stereochimica. Ad autovalori
(energie) vibrazionali pi`
u alti, corrispondono autofunzioni con pi`
u nodi, nelle diverse dimensioni; inoltre, la densit`a di probabilit`a |χ(Q)|2 diventa non trascurabile in
regioni sempre pi`
u ampie dello spazio delle configurazioni nucleari: in pratica tutte
quelle regioni dove Eku > Uk (Q) e anche limitate zone contigue dove Eku < Uk (Q).
163
Supponiamo di connettere due minimi nella superficie di energia potenziale con
un percorso di minima energia, che definisce la coordinata di reazione: si tratta
di andare da un minimo all’altro limitando il pi`
u possibile l’aumento di energia
potenziale. Il punto pi`
u alto del percorso sar`
a in genere un punto di sella, la cui
energia U ∗ determina l’energia di attivazione ∆E ∗ = U ∗ − Umin . Quando l’energia
vibrazionale Eku supera U ∗ , possiamo avere stati vibrazionali delocalizzati, ossia
probabilit`a non trascurabili di trovare il sistema nelle vicinanze di entrambi i minimi,
oltre che in ampie regioni tra i due e tutto intorno. In uno stato del genere il sistema
pu`
o passare liberamente da un minimo (reagenti) all’altro (prodotti). Una reazione
termica, in questa rappresentazione della realt`
a molecolare, consiste di due passaggi:
A) eccitazione vibrazionale, dallo stato fondamentale k = 0, u = 0, ad uno stato
con u sufficientemente alto da oltrepassare la barriera di attivazione; l’energia `e
fornita dagli urti con altre molecole, che sono molto efficienti nel trasformare energia
cinetica traslazionale in energia rotazionale e vibrazionale, o viceversa; naturalmente,
nell’urto tra una molecola vibrazionalmente eccitata ed una “fredda”, `e pi`
u facile che
la prima ceda energia alla seconda ed al moto traslazionale complessivo, piuttosto
che viceversa; in situazioni di equilibrio termodinamico, la popolazione degli stati
vibrazionali obbedir`
a alla statistica di Boltzmann.
B) diseccitazione, con cessione di energia al mezzo, cadendo in uno degli stati legati
di bassa energia corrispondenti ai prodotti, oppure ancora ai reagenti.
Tutto il processo si svolge su una sola superficie di energia potenziale, U0 (Q).
Implicitamente, quando abbiamo parlato dei minimi di Uk (Q) come di punti isolati
sulla superficie, ci siamo riferiti a reagenti e prodotti rappresentati da una singola molecola, quindi a reazioni di isomerizzazione. In una reazione di dissociazione,
lungo la coordinata di reazione, che si identifica essenzialmente con la lunghezza del
legame che si rompe, l’energia potenziale sale fino ad un valore asintotico Uk (∞),
eventualmente passando per un massimo (barriera di attivazione). La coordinata di
reazione, che era una delle coordinate vibrazionali dei reagenti, diventa una coordinata traslazionale relativa nella zona dei prodotti. Altre due coordinate vibrazionali
diventano coordinate di traslazione, e altre tre diventano coordinate rotazionali per
uno dei due frammenti. Alla fine avremo 6 coordinate traslazionali e 6 rotazionali,
invece delle 3+3 iniziali (eccezion fatta per i casi di molecole lineari, biatomiche o
atomi) e corrispondentemente 3N − 12 coordinate vibrazionali invece che 3N − 6.
Gli stati vibrazionali con Eku > Uk (∞) sono dissociativi, con uno spettro continuo.
La reazione consiste nell’accumulo di energia vibrazionale fino a sorpassare Uk (∞)
oppure l’eventuale energia di attivazione: a quel punto avviene la dissociazione e i
due frammenti seguono ciascuno una propria storia, che non `e pi`
u determinata dal
solo hamiltoniano molecolare, ma anche dalle interazioni con altre molecole.
E’ difficile arrivare ad una descrizione pi`
u dettagliata degli stati vibro-rotazionali
per sistemi di pi`
u di tre atomi, in presenza di moti interni di grande ampiezza, come
quelli che caratterizzano le reazioni chimiche. In alternativa ai metodi quantistici, si pu`
o ricorrere ad una trattazione classica del moto dei nuclei: ci`o si giustifica
164
considerando che i nuclei sono particelle abbastanza pesanti da ricadere, in prima
approssimazione, sotto le leggi della meccanica classica. Naturalmente, alcuni aspetti tipicamente quantistici del moto nucleare (livelli energetici discreti, penetrazione
nelle zone classicamente proibite, esistenza di un’energia minima di punto zero) si
perdono nella descrizione classica. Tuttavia, ne deriva una visione qualitativamente
corretta delle molecole e delle reazioni chimiche, capace anche di fornire previsioni
semiquantitative mediante il calcolo. Una molecola `e descritta da 3N coordinate che
specificano le posizioni dei nuclei e dalla corrispondente energia elettronica. Questo
punto materiale si sposta sulla superficie di energia potenziale Uk (Q); quando oscilla
intorno ad un minimo, abbiamo un composto stabile; quando ha sufficiente energia
per spostarsi da un minimo all’altro, passando pi`
u o meno vicino ad un punto di
sella, abbiamo la reazione chimica.
In ogni caso, finch´e non si induce un’eccitazione elettronica, tutta la dinamica si
svolge sulla pi`
u bassa ipersuperficie di potenziale, U0 (Q): si pu`
o dire che il 90% dei
fenomeni chimici e chimico-fisici siano semplicemente legati alla forma della U0 (Q).
In pratica solo in campi quali la spettroscopia elettronica o la fotochimica `e necessario
considerare esplicitamente le funzioni d’onda elettroniche, oltre che le energie.
7.2
Vibrazioni molecolari.
Supponiamo di conoscere l’energia elettronica per lo stato k, in funzione delle coordinate cartesiane dei nuclei, X1 . . . X3N . Intorno ad un minimo possiamo espandere
Uk (X) in serie di Taylor, limitandoci ai termini di secondo grado (approssimazione
armonica):
1X
Kij ∆Xi ∆Xj + . . .
(7.20)
Uk (X) = Umin +
2
i,j
Qui i ∆Xi indicano spostamenti dal punto di minimo. Inoltre:
2
∂ Uk
Kij =
∂Xi ∂Xj ∆X=0
(7.21)
Operiamo un cambiamento di variabili, definendo le nuove coordinate come
combinazioni lineari degli spostamenti cartesiani:
X
√
(7.22)
Qr =
Lir mi ∆Xi
i
dove mi `e la massa della particella associata alla coordinata Xi . Si pu`
o dimostrare
che, per un’opportuna scelta delle coordinate Qr , il potenziale Uk contiene solo
termini che dipendono da una singola variabile:
1X 2 2
Uk (X) = Umin +
ω Q
(7.23)
2 r r r
165
Delle 3N coordinate, 3 rappresentano traslazioni e 3 rotazioni complessive, quindi il
potenziale non pu`
o dipendere da queste 6 coordinate (5 nel caso di molecole lineari)
e le loro frequenze ωr risulteranno nulle. Al tempo stesso l’energia cinetica mantiene
una semplice espressione:
2
h
¯
Tˆ = −
2
X ∂2
∂Q2r
r
(7.24)
In questo sistema di coordinate, dette coordinate normali, l’hamiltoniano
vibrazionale `e quindi:
3N −6 2
1 X
2 2
2 ∂
ˆ
+ ω r Qr
Hvib =
−¯
h
2 r
∂Q2r
(7.25)
Abbiamo qui una somma di termini, ciascuno dipendente da una sola coordinata
Qr , perci`
o la funzione d’onda vibrazionale si fattorizza in un prodotto:
χu (Q) = χ1 (Q1 ) χ2 (Q2 ) . . . χ3N −6 (Q3N −6 )
(7.26)
Ogni χr `e soluzione di un’equazione di Schr¨
odinger monodimensionale con potenziale
armonico:
2
1
2 2
2 ∂
+ ωr Qr χr (Qr ) = Er χr (Qr )
−¯
h
(7.27)
2
∂Q2r
Per ogni coordinata normale Qr avremo quindi un numero quantico ur che identifica
il particolare autostato ed autovalore: in questo senso, l’indice u della χu (Q) diventa
una collezione di numeri quantici ur . Gli autovalori dell’energia vibrazionale per ogni
singola coordinata, in questa approssimazione, sono equispaziati:
Er = (ur + 1/2) h
¯ ωr
(7.28)
dove ωr `e la frequenza caratteristica di quella vibrazione. L’energia vibrazionale
complessiva `e una somma di termini, ciascuno relativo ad una coordinata normale:
Evib = h
¯
3N
−6
X
(ur + 1/2)ωr
(7.29)
r
Le frequenze vibrazionali variano da qualche decina di cm−1 a circa 3000: frequenze
pi`
u basse si hanno per moti che coinvolgono nuclei pi`
u pesanti e causano minori
variazioni del potenziale; frequenze pi`
u alte per nuclei leggeri e ripida salita del
potenziale (ossia legami pi`
u rigidi).
Con l’aumentare del numero quantico ur aumenta l’ampiezza della vibrazione, ossia
si allarga l’intervallo di Qr dove `e probabile trovare il sistema. Il troncamento al
secondo ordine della serie (7.20) diventa allora un’approssimazione pi`
u scadente; in
166
genere, il potenziale sale pi`
u lentamente di quanto indichi l’approssimazione armonica, ad esempio per l’allungamento di un legame o per l’allargamento di un angolo.
Ne consegue che, per un singolo modo vibrazionale, i livelli si infittiscono verso l’alto.
Un’altra conseguenza dell’anarmonicit`a `e pi`
u importante: per un potenziale qualsiasi
non vale la separazione dei moti vibrazionali secondo le coordinate normali, perci`
o
gli stati (7.26) non sono stati stazionari esatti; al pi`
u, sono buone approssimazioni
degli stati stazionari.
Consideriamo ad esempio una molecola poliatomica in cui diversi modi normali siano
moderatamente eccitati (ur piccoli ma non tutti nulli). Non `e difficile che l’energia vibrazionale totale (7.29) superi l’energia di dissociazione di qualche legame: in
questo caso ci troviamo nella regione di spettro continuo e gli stati sono dissociativi.
Se valesse una perfetta separazione dei moti vibrazionali, avremmo autostati legati
con l’eccitazione vibrazionale distribuita su vari modi ed autostati dissociativi con
l’eccitazione concentrata in una coordinata di stretch. In realt`
a, l’anarmonicit`a del
potenziale accoppia stati legati e stati dissociativi; perci`
o, se la molecola si trova
inizialmente in uno stato legato, essa decade sotto l’influenza di un perturbatore costante. In effetti, vari processi di eccitazione tendono a creare stati vibrazionalmente
eccitati con pochi quanti vibrazionali di alta frequenza. Per esempio, l’assorbimento di radiazione infrarossa in genere provoca salti vibrazionali di un solo quanto,
con frequenza uguale o vicina a quella della radiazione stessa; pu`
o accadere, in casi particolari, che una singola frequenza vibrazionale sia pi`
u grande dell’energia di
dissociazione di un legame, come nel caso di complessi legati da forze di van der
Waals; l’eccitazione di una delle frequenze pi`
u alte non provoca immediatamente la
dissociazione, ma piuttosto un decadimento del tipo visto sopra, che prende il nome di “predissociazione vibrazionale”. L’assorbimento di radiazione visibile o UV,
provoca la transizione ad un altro stato elettronico, accompagnata da un’eccitazione
vibrazionale: anche qui, come vedremo in seguito, pu`
o essere favorito uno stato di
arrivo quasi-legato, rispetto ad uno dissociativo, a seconda dei rispettivi fattori di
Franck-Condon. Un altro processo che pu`
o creare stati quasi-legati `e l’impatto con
elettroni e conseguente ionizzazione della molecola, sfruttato nella spettrometria di
massa; in alcuni casi gli ioni si frammentano in tempi comparabili con quelli di volo,
producendo picchi che non corrispondono a masse intere.
Pi`
u in generale, lo studio delle reazioni unimolecolari ruota intorno alla seguente
domanda: qual `e la probabilit`a che l’energia vibrazionale disponibile si concentri nel
moto lungo la coordinata di reazione (dissociazione o isomerizzazione)? A questo
proposito, bisogna tener conto che le coordinate normali stesse non sono in genere identificabili con moti localizzati di stretch, torsione o altro, bens`ı implicano
spostamenti di tutti i nuclei che compongono la molecola. Questo fatto, insieme
all’anarmonicit`a del potenziale, rende molto rapida la redistribuzione dell’energia
vibrazionale tra i moti interni (IVR, Intramolecular Vibrational Energy Redistribution): se, per esempio, si trova il modo di eccitare una vibrazione di stretch che
corrisponde alla dissociazione di un atomo posto ad un’estremit`a della molecola, pre-
167
sto tutti gli atomi, anche i pi`
u lontani, entreranno in vibrazione, assorbendo energia
dalla coordinata di reazione: in un certo senso, il fenomeno inverso della predissociazione vibrazionale. Una tipica vibrazione molecolare ha un periodo di 10−14 ÷ 10−13
sec; il tempo necessario per la “randomizzazione” dell’energia vibrazionale in una
molecola poliatomica in genere non `e molto pi`
u lungo.
Densit`
a di stati in funzione dell’energia e del numero
di atomi, poste tutte le frequenze pari a 1000 cm−1
numero energia N
q
ρ
−1
di atomi ( cm )
(stati/ cm−1 )
3
30000
3 30
0.5
4
30000
6 30
330
5
30000
9 30
5 × 104
6
30000
12 30
3 × 106
8
30000
18 30
3 × 109
10
3000
24 3
3
10
5000
24 5
100
10
10000
24 10
105
10
20000
24 20
109
10
30000
24 30
1012
Data l’energia vibrazionale totale, la probabilit`a di trovarsi in un particolare stato
dipende dal numero totale di stati che hanno la stessa energia, o quasi: questo in
assenza di agenti esterni che popolino preferenzialmente certi stati vibrazionali. In
una molecola poliatomica, i primi stati vibrazionali di ogni superficie elettronica sono
ben separati in energia: i ∆Eku,kv sono dell’ordine di una frequenza vibrazionale, o
di una differenza ωr − ωs (10−2 ÷ 10−3 a.u., salvo accidentali degenerazioni). Salendo
in energia, per`
o, abbiamo molti modi di distribuire l’energia vibrazionale tra le varie
coordinate, perci`
o il numero di livelli in un dato intervallo di energie pu`
o diventare
molto grande. Consideriamo per semplicit`a N modi normali di eguale frequenza, ω;
se l’energia totale `e E, si tratta di distribuire q = E
ω quanti vibrazionali (bosoni) in
(N +q−1)!
N oscillatori: la degenerazione `e q!(N −1)! . Poich´e i livelli energetici sono distanziati
di ω l’uno dall’altro, la densit`a di stati (numero di stati per unit`
a di energia), sar`
a:
s
1 N + q N + q N −1 N + q q
1 (N + q − 1)!
(7.30)
'
ρ(E) =
ω q!(N − 1)!
ω 2πN q
N
q
√
dove abbiamo applicato la formula di Stirling: n! ' 2πn nn e−n . ρ(E) `e una funzione rapidamente crescente dell’energia e del numero di atomi nella molecola, come
si vede dalla tabella.
In realt`
a, tenendo presente che in ogni molecola esiste di norma qualche frequenza
minore di 1000 cm−1 , e che l’anarmonicit`a del potenziale implica un infittimento
dei livelli al crescere dell’energia, la densit`a di stati aumenta ancor pi`
u rapidamente
168
di quanto indicato in tabella. Ne segue che, per una molecola grande (pi`
u di 4 o 5
atomi), la probabilit`a che l’energia si concentri in un particolare modo di vibrazione
`e molto bassa. Questo tipo di considerazioni sono alla base delle teorie delle velocit`a
di reazioni unimolecolari, tipo RRKM (vedi ad es. Gardiner, sez.5.2, oppure Smith,
sez.4.5).
7.3
Gli accoppiamenti nonadiabatici.
Gli stati ψk χku , autostati approssimati dell’hamiltoniano molecolare, sono accoppiati
ˆ cio`e quelli che contengono derivate delle funzioni
dai termini “nonadiabatici” di H,
d’onda elettroniche rispetto alle coordinate nucleari (vedi eq.(7.18)). L’entit`a della
perturbazione che accoppia due stati con indici k, u ed l, v `e data dall’integrale:
Z +∞
X ¯
h2
∂
(r)
(r)
Vku,lv =
χ∗ku tkl (Q) + 2gkl (Q)
χlv . . . dQs . . .
(7.31)
2µ
∂Q
r
r
−∞
r
dove:
(r)
gkl
(r)
tkl
∂ ψl
= ψk ∂Qr (7.32)
2 ∂ = ψk 2 ψl
∂Qr
(7.33)
(r)
Quando l’indice k `e uguale all’indice l (stesso stato elettronico), gkl = 0; infatti, le
ˆ el , ψk , si possono sempre considerare reali, e d’altra parte:
autofunzioni di H
(r)
(r)
gkl + glk =
=
=
Z
X
(7.34)
ψk∗ (q; Q)
Z
X
∂ψl
∂ψk
. . . dqi . . . +
. . . dqi . . . =
ψl∗ (q; Q)
∂Qr
∂Qr
Z
∂
∂ X ∗
ψk (q; Q) ψl (q; Q) . . . dqi . . . =
δkl = 0
∂Qr
∂Qr
(r)
Perci`
o, gli integrali gkl non contribuiscono all’accoppiamento di stati vibrazionali
diversi appartenenti allo stesso stato elettronico. Per quanto riguarda gli integrali
(r)
tkl , con k = l, essi giocano il ruolo di una piccola correzione del potenziale Uk ,
dipendente dalle masse nucleari: alterano quindi le energie Eku , ma non inficiano la
separazione dei moti elettronici e nucleari.
Per giustificare l’approssimazione di Born-Oppenheimer occorre anzitutto valutare
l’entit`
a degli accoppiamenti nonadiabatici (7.31). Naturalmente, i valori degli accoppiamenti variano a seconda della molecola e dello stato quantico, ed in alcuni
169
casi non sono affatto trascurabili. In genere per`
o, Vku,lv `e molto pi`
u piccolo dell’energia cinetica dei nuclei (qualche centinaio o migliaio di cm−1 ); infatti, di solito
la variazione delle funzioni d’onda elettroniche rispetto alle coordinate nucleari `e
molto pi`
u lenta di quelle vibrazionali. Per esempio, lungo la coordinata di stretch
o dissociazione di un legame, la χku compie diverse oscillazioni, o almeno passa da
zero al massimo e di nuovo a zero, per uno spostamento di ≈ 0.1 a.u.; la funzione
d’onda elettronica a distanza di equilibrio rappresenta un legame chimico, mentre
alla dissociazione rappresenta due radicali o due ioni: la variazione tra questi due
estremi implica uno spostamento di diverse a.u.
Il valore di Vku,lv deve essere confrontato con la differenza di energia degli stati che
accoppia. Se ψk `e lo stato elettronico fondamentale (k = 0) e χ0u uno degli stati
vibrazionali pi`
u bassi, esiste un accoppiamento solo con stati elettronicamente eccitati (l ≥ 1); la differenza di energia Elv − E0u `e normalmente dell’ordine di qualche
decina di migliaia di cm−1 . In questo caso, gli accoppiamenti nonadiabatici possono
causare solo piccole correzioni nei livelli energetici e nelle funzioni d’onda stazionarie: l’approssimazione di Born-Oppenheimer `e ben giustificata. Ad esempio, per la
molecola di BH, che `e stata oggetto di un’accurata indagine, la correzione `e inferiore
a 1 cm−1 ; prevedibilmente, si troveranno valori ancora pi`
u piccoli quando sono in
gioco nuclei pi`
u pesanti di H. Notiamo che i rispettivi ordini di grandezza delle
differenze di energia elettroniche, di quelle vibrazionali, e degli accoppiamenti nonadiabatici, sono legati rispettivamente all’inverso delle masse di elettroni e nuclei,
che compaiono nelle equazioni (7.13), (7.17) e (7.31): perci`
o `e valida l’idea intuitiva
che si possono considerare i nuclei pesanti e lenti, mentre gli elettroni aggiustano
quasi istantaneamente il loro moto in conseguenza delle variazioni di geometria della
molecola. Al crescere del numero quantico vibrazionale u diminuisce ∆E0u,lv e perci`o la correzione nonadiabatica diventa pi`
u importante; ora, alti numeri quantici u
sono necessari per raggiungere quelle zone dello spazio delle configurazioni nucleari
dove il potenziale U0 (Q) `e sensibilmente maggiore che nel minimo. Ci`
o equivale a
dire che per geometrie molecolari fortemente distorte, ad esempio per uno stato di
transizione, `e facile che l’approssimazione di Born-Oppenheimer sia meno valida che
in vicinanza della geometria di equilibrio.
Per gli stati elettronici eccitati la situazione `e differente. I ∆Eku,lv , con k ≥ 1 e l ≥ 1
sono in genere pi`
u piccoli di ∆E0u,lv , perch´e gli stati elettronici si infittiscono salendo
in energia; ma soprattutto, ogni stato vibronico elettronicamente eccitato ψk χku
`e degenere o quasi con qualche stato ψl χlv appartenente a superfici elettroniche
pi`
u basse. Di norma la funzione d’onda ψl χlv dovr`a essere vibrazionalmente molto
eccitata e conterr`
a molti nodi; perci`
o l’elemento di matrice (7.31) che la accoppia
con ψk χku sar`
a molto piccolo. Avremo comunque transizioni non radiative tra i due
stati e conversione di energia elettronica in energia vibrazionale. L’accoppiamento
diventa rapidamente meno efficace con l’aumentare della differenza di energia tra
i due stati elettronici; perci`
o, si pu`
o dire che, quanto pi`
u grande `e la quantit`a di
energia elettronica da trasformare in energia cinetica dei nuclei, tanto meno probabile
170
`e l’evento.
Figura 7.1: Stati vibronici legati interagenti.
0.16
U1
0.14
χ10
0.12
χ01
0.1
0.08
U0
0.06
0.04
χ00
0.02
0
2
2.5
3
3.5
4
Qr
4.5
5
5.5
6
Figura 7.2: Predissociazione.
0.16
U1
0.14
χ10
0.12
χ01
0.1
0.08
U0
0.06
0.04
χ00
0.02
0
2
3
4
5
Qr
6
7
8
Nel caso in cui χlv appartenga ad uno spettro continuo di stati dissociativi, siamo
nelle tipiche condizioni in cui `e valida la regola d’oro di Fermi. Se eccitiamo la
molecola nello stato legato ψk χku , avremo un decadimento di forma esponenziale,
con tempo di vita determinato dall’accoppiamento nonadiabatico e dissociazione
finale. Questo processo `e detto predissociazione. In molecole poliatomiche con pi`
u
171
di 5 o 6 atomi, abbiamo visto che gli stati con energia di eccitazione vibrazionale di
qualche decina di migliaia di cm−1 , pari cio`e ad un’eccitazione elettronica, danno
luogo ad un insieme fittissimo di livelli, quasi continuo. La situazione `e analoga
a quella descritta sopra per la predissociazione di piccole molecole. Anche qui,
l’accoppiamento nonadiabatico `e causa di transizioni non radiative, che vanno sotto
il nome di Internal Conversion (IC).
Di norma, radiazione di lunghezza d’onda adeguata pu`
o causare transizioni dallo
stato fondamentale ψ0 χ00 a ψk χku , con una probabilit`a che dipende dal dipolo di
transizione, un integrale del tipo di eq. (6.65). L’eccitazione diretta a stati dissociativi `e molto meno probabile, perch´e la funzione d’onda vibrazionale d’arrivo, χlv ,
contiene molti nodi e quindi d`
a luogo a un piccolo valore dell’integrale. L’accoppiamento di ψk χku con gli stati dissociativi crea per`
o un insieme continuo di livelli
che corrispondono a stati misti, con un contributo dello stato legato decrescente via
via che ci si allontana dall’energia Eku . L’intervallo di energie interessato `e proporzionale al modulo quadro dell’accoppiamento, perci`
o `e inversamente proporzionale
al tempo di vita (regola d’oro di Fermi). La predissociazione in molecole piccole `e
quindi associata a spettri di assorbimento discreti, ma con righe allargate o “diffuse”.
Naturalmente, le molecole che compiono la transizione allo stato sottostante dissociano. Perci`
o, sempre parlando di piccole molecole in fase vapore, abbiamo le
seguenti corrispondenze:
spettro di righe ⇐⇒ stati legati ⇐⇒ l’eccitazione non causa dissociazione
spettro continuo ⇐⇒ stati dissociativi ⇐⇒ dissociazione nello stato d’arrivo
spettro di righe diffuse ⇐⇒ stati metastabili ⇐⇒ predissociazione.
Troviamo esempi dei tre casi nello spettro e nella fotochimica dell’ossigeno molecolare (vedi figura 7.3). Transizioni agli stati 1 ∆g e 1 Σ+
g non possono dissociare O2 ;
lo spettro `e discreto, con intensit`
a di assorbimento molto debole, poich´e si tratta di
transizioni proibite di spin; il decadimento allo stato fondamentale, in gas rarefatti,
quindi in assenza di urti, `e unicamente radiativo (emissione IR/visibile). Un’altra
serie di righe inizia intorno a λ = 280 nm (stato 3 Σ+
u ) e per λ < 250 nm si trasforma
in uno spettro continuo: si ha quindi dissociazione in due atomi di ossigeno, entrambi nello stato fondamentale 3 P . Una serie di righe pi`
u intense si sovrappone al
continuo, a partire da λ = 200 nm: `e la prima transizione permessa per spin e simmetria, allo stato B 3 Σ−
u . Anche qui, con λ < 175 nm si passa allo spettro continuo,
ma questa volta uno dei prodotti di dissociazione `e un atomo di O eccitato, nello
stato 1 D. Pi`
u su ancora, una serie di righe molto diffuse indica la presenza dello
3 −
0
stato B Σu , che interagisce fortemente col sottostante stato B, tramite l’accoppiamento nonadiabatico. La predissociazione genera gli stessi prodotti che si ottengono
eccitando direttamente in B. Poich´e l’allargamento delle righe `e di qualche decina
di cm−1 , il tempo di vita dello stato B 0 `e dell’ordine di 10−13 secondi.
172
Figura 7.3: Curve di potenziale e spettro di O2 .
7.4
Principio di Franck-Condon
L’assorbimento od emissione di luce visibile o UV implica transizioni tra stati elettronici. Considerando le funzioni d’onda iniziale e finale nell’approssimazione di
Born-Oppenheimer, ψk χku e ψl χlv , il dipolo di transizione si potr`
a scrivere:
Z
Z
µku,lv = dQ χ∗ku (Q) χlv (Q) dq ψk∗ (q; Q) ~µ ψl (q; Q)
~
(7.35)
Ricordiamo che il momento di dipolo `e un vettore che ha per componenti:
X
X
µλ =
Zα Rαλ −
riλ
α
i
173
(7.36)
dove Zα , Rαλ e riλ sono cariche e coordinate cartesiane di nuclei ed elettroni
rispettivamente. Per ogni componente λ, l’integrale:
Z
µklλ (Q) = ψk∗ (q; Q) µλ ψl (q; Q) dq
(7.37)
`e una funzione delle coordinate nucleari Q, in quanto dipendono da Q sia le funzioni
d’onda ψk e ψl , sia la componente µλ . Notiamo che i termini contenenti cariche
nucleari in µλ non contribuiscono all’integrale; infatti, poich´e questi
R termini non
dipendono dalle q elettroniche, i rispettivi integrali si riconducono a ψk ψl dq, che `e
nullo per k 6= l. Una volta definito il dipolo di transizione elettronico ~µkl , il dipolo di
transizione tra stati di Born-Oppenheimer si ottiene integrando sulle sole coordinate
nucleari:
Z
µkl (Q) χlv (Q) dQ
(7.38)
µku,lv = χ∗ku (Q) ~
~
Spesso `e lecito trascurare la dipendenza di ~µkl dalle Q, considerandolo approssimativamente costante e uguale al suo valore nella geometria di equilibrio per lo stato
elettronico di partenza (approssimazione di Condon): ~µkl (Q) ' ~µkl (Qeq ). Allora:
Z
(7.39)
µku,lv ' ~
~
µkl (Qeq )
χ∗ku (Q) χlv (Q) dQ
Le probabilit`a di transizione sono proporzionali al modulo quadro di ~µku,lv :
Bku,lv =
2π |µku,lv |2
3¯
h2
(7.40)
R
2
e quindi ai quadrati degli integrali di sovrapposizione, χ∗ku χlv dQ , detti fattori
di Franck-Condon.
Le funzioni vibrazionali, quindi i fattori di Franck-Condon, dipendono dalla forma
dei due potenziali elettronici, Uk (Q) e Ul (Q). Alle temperature ordinarie, abbiamo
visto che lo stato vibrazionale pi`
u basso dello stato elettronico fondamentale, χ00 , `e
di gran lunga il pi`
u popolato. Anche negli stati elettronici eccitati, l’energia vibrazionale eventualmente presente, dovuta al processo di eccitazione, viene rapidamente
dispersa dagli urti con altre molecole, eccetto il caso di gas molto rarefatti. Cos`ı, in
gas a pressioni ordinarie o in soluzione, in assorbimento come in emissione, lo stato
vibrazionale di partenza `e il pi`
u basso.
7.5
Processi di eccitazione e decadimento.
Se il risultato di un processo di eccitazione (assorbimento di energia dal campo
eletromagnetico, collisione, o altro) fosse quello di popolare un autostato dell’hamiltoniano molecolare diverso da quello di partenza, non ci sarebbe alcuna evoluzione
174
temporale, perch´e avremmo creato un altro stato stazionario. In particolare, non
avremmo i processi di decadimento (IC, predissociazione) descritti in sezione 7.3.
Teoricamente questo `e possibile, se la radiazione eccitante `e sufficientemente monocromatica, in modo da poter selezionare un solo livello energetico. In caso contrario,
lo stato eccitato risultante sar`
a una combinazione lineare di molti autostati, pi`
uo
meno vicini in energia, e perci`
o uno stato non stazionario. Quale tipo di stato non
stazionario venga prodotto dipende dalle caratteristiche della radiazione usata.
Tratteremo due casi ideali di grande interesse. Il primo permette di chiarire perch´e,
e in quali condizioni, lo stato che si crea al momento dell’eccitazione `e un puro
prodotto di Born-Oppenheimer, soggetto a decadimento, piuttosto che un autostato
molecolare. Nel secondo caso vedremo come si crea una funzione d’onda vibrazionale
non stazionaria (nel senso del moto dei nuclei) in concomitanza con una transizione
elettronica.
Consideriamo uno stato vibrazionale legato χ1u nel potenziale di uno stato elettronico eccitato ψ1 ; lo stato ψ1 χ1u interagisce, per esempio a causa di accoppiamenti
nonadiabatici, con un continuo di stati vibrazionali dissociativi appartenenti ad uno
stato elettronico sottostante ψ0 (il fondamentale o un altro stato), come in figura
7.1; oppure, con un quasi-continuo di stati legati molto fitti, χ0v , come in figura 7.2.
Un autostato esatto di energia Ek vicina a E1u `e essenzialmente una combinazione
lineare di ψ1 χ1u e degli ψ0 χ0v :
X
Ψk = ak ψ1 χ1u +
bvk ψ0 χ0v
(7.41)
v
Come gi`
a osservato, le χ0v sono funzioni d’onda con molti nodi, quindi i dipoli di transizione che le collegano allo stato di partenza ψ0 χ00 sono piccolissimi:
hψ0 χ00 |µx | ψ0 χ0v i ' 0. Perci`
o, il dipolo di transizione all’autostato Ψk dipende
solo dal coefficiente ak :
µxk0 = hΨk |µx | ψ0 χ00 i = a∗k hψ1 χ1u |µx | ψ0 χ00 i
(7.42)
Lo stato ψ1 χ1u , unico responsabile dell’esistenza di una banda di assorbimento in
questo intervallo di frequenze, `e detto “stato porta” o “doorway state”. A loro
volta, i coefficienti dipendono da energie relative ed accoppiamenti; l’equazione di
Schr¨
odinger per la funzione d’onda Ψk `e:
ˆ k = Ek Ψk
HΨ
(7.43)
ossia, moltiplicando per ψ1∗ χ∗1u e integrando:
X
ak (E1u − Ek ) =
bvk V1u,0v
(7.44)
v
dove V1u,0v pu`
o essere un accoppiamento nonadiabatico tipo (7.31) o avere un’altra
origine fisica. Poich´e la forza dell’oscillatore `e legata a |ak |2 , l’eq. (7.44) mostra che
175
lo spettro di assorbimento `e composto da un insieme di righe (in pratica, un’unica
banda) con intensit`
a decrescente man mano che ω0k = (Ek − E00 )/¯h si allontana da
ωmax = (E1u − E00 )/¯
h; supponendo che l’accoppiamento V1u,0v non dipenda troppo
irregolarmente dallo stato χ0v , l’intensit`
a varia come (E1u − Ek )−2 = ¯h2 (ω0k −
−2
ωmax ) ; la larghezza di banda `e proporzionale a |V1u,0v |2 .
Figura 7.4: Processo di eccitazione.
ω
impulso breve
impulso lungo
ψ0 χ0v
Ψk
V1u,0v
ψ1 χ1u
impulso lungo
E0
imp. breve
t
2
˜ E0 Se consideriamo un singolo impulso di radiazione, la teoria delle perturbazioni (sezione 6.6) dice quali sono i coefficienti degli stati Ψk dopo l’eccitazione; omettendo
il contributo relativo allo stato di partenza (che normalmente non interferisce con
l’evoluzione temporale dello stato eccitato), dopo il passaggio dell’impulso eccitante
la funzione d’onda sar`
a:
r
i πX
µxk0 E˜0 (ω − ωk0 ) Ψk
(7.45)
Ψ=
h 2
¯
k
dove E˜0 `e la trasformata
di Fourier della funzione di modulazione dell’impulso di
˜ 2
radiazione; E0 `e lo spettro della luce eccitante, e rappresenta una banda centrata sulla frequenza ω. Le relazioni tra livelli energetici, spettro della radiazione e
lunghezza dell’impulso sono rappresentate in figura 7.4. I livelli corrispondenti agli
stati Ψk sono in parte a tratto pieno e in parte tratteggiati, per visualizzare il peso
|ak |2 dello stato porta, ovvero la forza dell’oscillatore.
Se la luce `e altamente monocromatica (impulso lungo), cio`e E˜0 (∆ω) ha un picco
molto alto e stretto in corrispondenza di ∆ω = 0, lo stato eccitato Ψ coincide con
l’autostato Ψk , quindi `e uno stato stazionario. Questo `e praticamente possibile solo
quando la densit`a degli stati `e abbastanza bassa (piccole molecole). Viceversa, se
176
l’impulso `e di breve durata, la larghezza di banda della radiazione eccitante supera di
molto quella della transizione 00 → kv; con ω ' ωmax , possiamo allora considerare
la funzione E˜0 praticamente costante all’interno dell’intervallo di frequenze per le
quali si hanno dipoli di transizione µxk0 non trascurabili. In questo caso l’eq. (7.45)
diventa
r
X
i π ˜
Ψ'
E0 (ω − ωmax )
µxk0 Ψk =
(7.46)
h 2
¯
k
i
=
h
¯
r
X
π ˜
E0 (ω − ωmax ) hψ1 χ1u |µx | ψ0 χ00 i
a∗k
2
k
"
ak ψ1 χ1u +
X
v
bvk ψ0 χ0v
#
Ora, per l’unitariet`a della matrice dei coefficienti
che esprimono
P
P ∗ gli autostati esatti
2
Ψk nell’insieme {ψ1 χ1u , ψ0 χ0v }, abbiamo:
|a
|
=
1
e
k k
k ak bvk = 0. Ne segue
che lo stato eccitato creato da un impulso abbastanza breve `e una combinazione di
autostati Ψk tale da “ricostruire” lo stato porta ψ1 χ1u , senza alcun contributo degli
stati privi di forza dell’oscillatore ψ0 χ0v :
r
i π ˜
E0 (ω − ωmax ) hψ1 χ1u |µx | ψ0 χ00 i ψ1 χ1u
(7.47)
Ψ'
h 2
¯
La successiva evoluzione temporale in assenza di radiazione si pu`
o descrivere come il
decadimento dello stato non stazionario ψ1 χ1u , “causato” dal suo accoppiamento con
gli stati ψ0 χ0v . Applicando la regola d’oro di Fermi, e ricordando la relazione vista
sopra tra accoppiamento ed allargamento della riga spettrale, vediamo che tempo di
vita e larghezza di riga sono inversamente proporzionali. Queste considerazioni si
applicano a rigore solo a treni d’onda coerenti, quali sono prodotti dai laser; la luce
prodotta dalle pi`
u comuni sorgenti si pu`
o descrivere come una sovrapposizione di
brevi treni d’onda, non correlati l’uno con l’altro; tuttavia, considerando separatamente l’azione di ciascun impulso, si pu`
o estendere la conclusione contenuta in eq.
(7.47) a luce non altamente monocromatica prodotta da una sorgente qualsiasi.
Come secondo processo ideale di eccitazione, consideriamo il caso di un impulso
ancora pi`
u breve e quindi di maggiore estensione nel dominio delle frequenze, tale
da abbracciare diversi livelli vibrazionali E1u . In base a quanto visto sopra, possiamo
senz’altro trascurare il contributo di eventuali stati vibrazionalmente molto eccitati,
appartenenti ad altri stati elettronici, debolmente accoppiati con gli ψ1 χ1u . La
funzione d’onda eccitata sar`
a quindi una combinazione di diversi stati vibrazionali
tutti appartenenti a ψ1 , cio`e un ”pacchetto d’onda” non stazionario che si muover`
a
sulla superficie di potenziale U1 (Q):
r
i π X˜
E0 (ω − ω0u ) hψ1 χ1u |µx | ψ0 χ00 i ψ1 χ1u
(7.48)
Ψ=
h 2 u
¯
Qui ω0u = (E1u − E00 )/¯
h `e la frequenza di risonanza per arrivare al livello E1u .
Se l’impulso di radiazione `e ultrabreve, dell’ordine di pochi femtosecondi, ancora
177
una volta larghezza di banda della radiazione `e tanto grande da poter considerare
E˜0 (ω−ω0u ) quasi indipendente dall’indice corrente u, sostituendola con una costante;
inoltre, introduciamo l’approssimazione di Condon, eq. (7.39), per valutare il dipolo
di transizione:
r
X
i π ˜
hχ1u |χ00 i χ1u
(7.49)
Ψ=
E0 (ω − ωvert ) hψ1 |µx | ψ0 i ψ1
h 2
¯
u
Qui ωvert corrisponde alla cosiddetta energia di “eccitazione verticale”, cio`e U1 − U0
valutata alla geometria di equilibrio dello stato di partenza; ωvert di norma approssima bene la frequenza della sottobanda vibrazionale pi`
u intensa nello spettro di
assorbimento. La somma in eq. (7.49) altro non `e che lo sviluppo della funzione
d’onda χ00 , in termini dell’insieme {χ1u }. Perci`
o, la funzione d’onda eccitata si pu`
o
scrivere pi`
u semplicemente:
r
i π ˜
Ψ=
E0 (ω − ωvert ) hψ1 |µx | ψ0 i ψ1 χ00
(7.50)
h 2
¯
Dunque la radiazione, sommando insieme contributi di tanti stati vibrazionali appartenenti allo stato ψ1 con opportuni coefficienti, “ricostruisce” lo stato vibrazionale
di partenza χ00 , cambiando soltanto la funzione d’onda elettronica. Naturalmente
χ00 , portato in un potenziale che non `e “il suo”, non `e pi`
u uno stato stazionario, e
si muover`
a in discesa sulla superficie di energia dello stato eccitato.
` bene sottolineare che questa forma di eccitazione (traslazione della funzione d’onda
E
vibrazionale da uno stato elettronico all’altro) `e una conseguenza del principio di
Franck-Condon, e avviene solo se l’impulso di luce eccitante `e estremamente breve.
Luce moderatamente monocromatica pu`
o selezionare un livello vibrazionale nello
stato elettronico d’arrivo, creando perci`
o uno stato stazionario dal punto di vista del
moto nucleare. La descrizione classica dell’eccitazione tipo Franck-Condon implica
una transizione elettronica senza variazione nelle posizioni e velocit`a dei nuclei; a
questo punto `e evidente che tale descrizione somiglia alla creazione di un pacchetto
d’onda mediante un impulso di radiazione brevissimo, piuttosto che alla popolazione
di un autostato vibrazionale con luce approssimativamente monocromatica.
7.6
Un modello per
poliatomiche.
il
decadimento
in
molecole
In una molecola poliatomica, gli elementi di matrice di accoppiamento nonadiabatico
o spin-orbita possono subire ampie variazioni, praticamente casuali, in dipendenza
dai numeri quantici vibrazionali. Tuttavia, nell’impossibilit`a di esaminare in dettaglio tali variazioni, Bixon e Jortner (J. Chem. Phys. 48, 715 (1968)) hanno proposto
un modello semplificato di interazione tra uno stato |Ai ed un quasi-continuo di stati
178
|Ki, che potremmo identificare rispettivamente con uno dei primi stati vibrazionali
di un singoletto o multipletto eccitato di una poliatomica e con gli stati vibrazionalmente molto eccitati del fondamentale elettronico. In questo modello gli stati
|Ki sono tutti equidistanti in energia ed hanno lo stesso elemento di matrice di
accoppiamento con |Ai.
Figura 7.5: Il modello di Bixon-Jortner.
Eα
εK
K
4 }
ε
3
2
1
0
-1
-2
-3
-4
εA
V
}∆
Definiamo un hamiltoniano di ordine zero che ha per autostati |Ai e l’insieme dei
|Ki:
X
ˆ (0) = |Ai εA hA| +
H
|Ki εK hK|
(7.51)
K
Le energie εK sono equidistanziate e quella immediatamente sottostante εA
corrisponde a K = 0, con εA − εK = ∆:
εK = εA − ∆ + Kε
(7.52)
ˆ 1 accoppia
Il parametro ε `e l’inverso delle densit`a di stati: ε = ρ−1 . Il perturbatore H
gli stati |Ki solo con |Ai e non tra di loro:
X
ˆ (1) = V
H
(|Ai hK| + |Ki hA|)
(7.53)
K
Indichiamo con |αi gli autostati dell’hamiltoniano completo:
ˆ |αi = (H
ˆ (0) + H
ˆ (1) ) |αi = Eα |αi
H
(7.54)
Supponiamo che il sistema si trovi inizialmente nello stato |Ai, in seguito all’assorbimento di radiazione o per altra causa. Come discusso nella sezione 7.5, questa
179
assunzione `e realistica se si verificano le seguenti condizioni:
A) gli stati |Ki non sono accoppiati con lo stato di partenza, di solito il fondamentale, tramite la radiazione (“non portano forza dell’oscillatore”); questo normalmente
avviene perch´e gli elementi di matrice del dipolo e degli altri operatori di accoppiamento radiazione-molecola sono estremamente piccoli tra uno stato con bassi numeri
quantici vibrazionali ed un altro molto pi`
u eccitato; B) la larghezza spettrale della
radiazione `e tale da impedire di eccitare precisamente un solo stato quantico |αi; ci`o
`e probabile, data l’enorme densit`a degli stati |Ki;
C) la larghezza spettrale della radiazione non `e cos´ı grande da eccitare
contemporaneamente altri stati di tipo |Ai, appartenenti allo stesso stato elettronico.
L’evoluzione temporale a partire dallo stato |Ai `e determinata applicando
l’equazione di Schr¨
odinger dipendente dal tempo alla funzione d’onda |Ψ(t)i:
X
|Ψ(t)i = CA (t)e−iωA t |Ai +
CK (t)e−iωK t |Ki
(7.55)
K
dove ωA = εA /¯
h e ωK = εK /¯
h. Si ottiene un set di equazioni differenziali accoppiate
per i coefficienti CA e CK :
iX
CK eiωAK t V
(7.56)
C˙ A (t) = −
h
¯
K
i
C˙ K (t) = − CA e−iωAK t V
h
¯
CA (0) = 1
e
(7.57)
CK (0) = 0 ∀K
dove ωAK = (∆ − Kε)/¯
h. La (7.57) fornisce:
Z t
0
i
CA (t0 )e−iωAK t dt0
CK (t) = − V
h
¯
0
(7.58)
(7.59)
Sostituendo la (7.59) nella (7.56) si ottiene un’equazione per il solo coefficiente CA :
Z
V2 X t
0
˙
CA (t) = − 2
CA (t0 )e−i(Kε−∆)(t−t )/¯h dt0
(7.60)
h K 0
¯
Se la densit`a di stati `e molto elevata (in particolare, se εt/¯h 1), la somma in
equazione (7.60) pu`
o essere sostituita con un integrale. Si ha allora:
Z +∞
Z
0
V2 t
eix(t −t) dxdt0
(7.61)
CA (t0 )
C˙ A (t) ' −
hε 0
¯
−∞
dove ho posto x = (Kε − ∆)/¯
h e quindi dx = εdK/¯h. Come gi`a sappiamo dalla
(6.22) la trasformata di Fourier di una costante `e proporzionale ad una delta di
Dirac:
Z +∞
1
ei(x−x0 )ω dω
(7.62)
δ(x − x0 ) =
2π −∞
180
La (7.61) pu`
o dunque essere riscritta in questo modo:
2πV
C˙ A (t) = −
hε
¯
2
t
Z
0
CA (t0 )δ(t0 − t)dt0
(7.63)
A questo punto `e possibile risolvere l’integrale a secondo membro della (7.63)
notando che:
Z b
f (b)
δ(x − b)f (x)dx =
(b > a)
(7.64)
2
a
Si ottiene in questo modo un’equazione differenziale del primo ordine lineare e
omogenea per il coefficiante CA :
2
πV
CA (t)
C˙ A (t) = −sgn(t)
hε
¯
(7.65)
dove sgn(t) = 1 per t > 0 e sgn(t) = −1 per t < 0. Con la condizione iniziale
CA (0) = 1 si ha:
CA (t) = e−πV
2 |t|/¯
hε
(7.66)
La probabilit`a di trovare la molecola nello stato |Ai dopo un tempo t > 0 `e data
da:
PA (t) = |CA (t)|2 = e−2πV
2 t/¯
hε
= e−t/τ
(7.67)
con τ = h
¯ ε/2πV 2 = h
¯ /2πV 2 ρ. L’equazione (7.67) mostra che lo stato |Ai “decade”
nel tempo esponenzialmente, con un tempo di vita inversamente proporzionale a
V 2 , forza dell’accoppiamento, e alla densit`a dei livelli. La probabilit`a di trovare la
molecola nello stato |Ai diminuisce, mentre aumenta quella di trovarlo in uno dei
tanti stati |Ki, con minore energia elettronica e maggiore energia vibrazionale.
Possiamo ora determinare la forma di riga, ovvero la distribuzione dello stato |Ai
popolato inizialmente negli autostati esatti. A tale scopo scriviamo:
X
X
|Ψ(t = 0)i = |Ai =
|αi hα |A i =
cαA |αi
(7.68)
α
|Ψ(t)i =
X
α
α
cαA e−iωα t |αi
(7.69)
dove ωα = Eα /¯
h. Naturalmente, gli stati |αi sono a loro volta combinazioni lineari
di |Ai e |Ki:
!
X
X
|αi = |Ai hA| +
|Ki hK| |αi = cAα |Ai +
cKα |Ki
(7.70)
K
K
181
Sostituendo nella (7.69) abbiamo:
X
X
|Ψ(t)i =
|cAα |2 e−iωα t |Ai +
cαA cKα e−iωα t |Ki
(7.71)
Confrontando la (7.71) con la (7.55) si ha:
X
CA (t)e−iωA t =
|cAα |2 e−iωα t
(7.72)
Ma CA (t) `e stato gi`
a ricavato ed `e dato dalla (7.66). Si ottiene quindi:
X
2
|cAα |2 e−iωαA t = e−πV |t|/¯hε
(7.73)
α
α,K
α
α
con ωαA = (Eα − εA )/¯
h. Supponiamo a questo punto che la densit`a degli autostati
esatti sia uguale a quella degli stati di ordine zero, ovvero a ε−1 . Ci`
o `e ragionevole
se si considera che ogni stato |Ki interagisce solo con |Ai. D’altra parte `e possibile
dimostrare rigorosamente (ma noi non lo faremo) che fra ogni coppia di stati |Ki
` quindi lecito approssimare Eα ' εA + αε con
`e contenuto un solo stato |αi. E
α = 0, ±1, . . . ± ∞. Inoltre, data l’enorme densit`a di stati sostituiamo la somma in
equazione (7.73) con un integrale:
Z +∞
2
ε
(7.74)
|cAx |2 e−ixt dx = e−πV |t|/¯hε
h
¯
−∞
con x = (Eα − εA )/¯
h = αε/¯
h. La (7.74) ci dice che i coefficienti di espansione
di |Ai negli autostati esatti sono correlati a CA (t) da una trasformata di Fourier.
Un’espressione esplicita per la forma di riga si pu`
o dunque ottenere in questo modo:
moltiplico ambo i membri della (7.74) per eix0 t e integro da −∞ a +∞.
Z
Z +∞
Z +∞
ε +∞ ix0 t −πV 2 |t|/¯hε
e e
dt
(7.75)
|cAx |2
e−i(x−x0 )t dtdx =
¯h −∞
−∞
−∞
Dalla (7.62) si ottiene:
Z +∞
ε
2
|cAx0 |2 =
eix0 t e−πV |t|/¯hε dt
2π¯
h −∞
(7.76)
L’integrale che compare nella (7.76) si pu`
o valutare facilmente. Infatti, se a e b sono
due numeri reali con a > 0 si ha:
Z +∞
Z 0
Z +∞
2a
−a|t| ibt
at ibt
e
e dt =
e e dt +
e−at eibt dt = 2
(7.77)
a + b2
−∞
0
−∞
Quindi la (7.76) diventa, rimpiazzando x0 con (Eα − εA )/¯h:
"
2 #−1
E
−
ε
α
A
2
|cAα | = π 2 V 2 ε−2 +
V
182
(7.78)
Questa `e una lorenziana, funzione di Eα , con massimo in εA e semilarghezza a met`
a
altezza ∆E = πV 2 /ε. Ci`
o significa che lo stato |Ai si distribuisce sugli autostati
esatti con peso decrescente quando Eα si allontana da εA . Se |Ai “porta forza
dell’oscillatore” e i |Ki no, la riga spettrale di assorbimento da uno stato qualsiasi
verso |Ai avr`a forma lorenziana con ∆ν = V 2 /2ε¯h.
Notiamo infine che per tempi molto lunghi (t ≈ ¯h/ε = h
¯ ρ), non sarebbe lecito
trasformare la somma in eq. (7.60) nell’integrale di eq. (7.61), perch´e il fattore esponenziale complesso varierebbe sensibilmente passando da un K al successivo; questo
significa che, in luogo di esibire un decadimento esponenziale, la probabilit`a PA (t)
dovrebbe oscillare nel tempo, come `e caratteristico dei sistemi quantistici con livelli
discreti. Tuttavia, `e di norma impossibile considerare il sistema veramente isolato per tempi cos´ı lunghi: una molecola in uno stato vibrazionalmente eccitato |Ki
perde facilmente energia negli urti con altre molecole. L’interazione con l’ambiente
circostante fa quindi decadere anche gli stati |Ki e garantisce un decadimento esponenziale o di forma simile anche al di l`a dei limiti di validit`
a delle approssimazioni
fatte.
7.7
La regola di non intersezione e gli incroci “evitati”.
L’accoppiamento nonadiabatico `e particolarmente efficace nella cosiddetta situazione
di “incrocio evitato”. Consideriamo le curve di energia potenziale per due stati
elettronici di una molecola biatomica. Supponiamo che i due stati, η1 ed η2 , siano
rappresentati con buona approssimazione da due configurazioni, ossia due diverse
disposizioni degli elettroni negli orbitali, con certi spin. Alternativamente, η1 ed η2
potrebbero essere approssimati da due strutture tipo Valence Bond, ad esempio una
ionica ed una covalente. Funzioni d’onda elettroniche che conservano a geometrie
diverse lo stesso carattere, cio`e la stessa distribuzione degli elettroni nello spazio e
negli orbitali, sono dette diabatiche. Pu`o accadere che le curve dell’energia di η1
ed η2 tendano ad incrociarsi, perch´e a certe distanze `e pi`
u stabile (ha energia pi`
u
bassa) η1 , ad altre distanze `e pi`
u stabile η2 . Sorge quindi il problema: `e possibile
che per qualche distanza R gli autostati elettronici siano degeneri, ossia che le curve
Uk (R) si incrocino? Che effetto pu`
o avere l’incrocio sulla probabilit`a di transizione
da uno stato all’altro? Indichiamo con Hij gli integrali:
D E
ˆ (7.79)
Hij (R) = ηi H
el ηj
H11 e H22 sono le curve di energia potenziale che si incrociano, H12 `e un’interazione
elettronica tra i due stati η1 ed η2 . Se H12 6= 0, gli stati diabatici non sono autostati
ˆ el . I veri autostati saranno, con buona approssimazione, combinazioni lineari
di H
degli ηi :
ψ1 = cosθ η1 + sinθ η2
(7.80)
183
ψ2 = −sinθ η1 + cosθ η2
(7.81)
Figura 7.6: Incrocio evitato in NaCl.
energia (a.u.)
0.1
Na+ +Cl−
E2 ' H11 (ionico)
0.05
E2 ' H22 (neutro)
0
E1 ' H22 (neutro)
-0.05
E1 ' H11 (ionico)
-0.1
-0.15
accopp. nonadiab. (a.u.)
5
10
15
20
R(Na − Cl), bohr
25
30
g12
4
2
0
E2 ' H11
H12 (x3)
energia (a.u.)
0.0005
E2 ' H22
0
E1 ' H22
-0.0005
-0.0010
E1 ' H11
17
17.5
18
R(Na − Cl), bohr
184
18.5
19
I coefficienti cosθ e ±sinθ assicurano la normalizzazione. Gli autovalori sono:
s
H11 + H22
H11 − H22 2
2
U1,2 (R) =
±
+ H12
(7.82)
2
2
(il segno + vale per U2 e il − per U1 ).
Per avere U1 = U2 occorre che siano soddisfatte contemporaneamente le due condizioni: H11 = H22 e H12 = 0. Avendo a disposizione una sola incognita, R, le
soluzioni delle due equazioni in generale non coincidono. Ne segue la regola di non
intersezione (non-crossing rule): curve dell’energia di stati con la stessa simmetria
spaziale e di spin non si intersecano. Se gli stati hanno diversa simmetria, H12 = 0
`e un’identit`
a, indipendentemente da R: allora U1 e U2 coincidono con H11 e H22 e
possono intersecarsi, come accade, per esempio, tra singoletti e tripletti.
Un esempio classico ci `e dato dai primi due stati di singoletto degli alogenuri alcalini,
come NaCl (vedi figura 7.6). A brevi distanze internucleari lo stato fondamentale `e
essenzialmente ionico (Na+ Cl− ). A grandi distanze invece lo stato a distribuzione di
carica neutra (Na2 S + Cl2 P) `e il pi`
u basso, perch´e l’affinit`
a elettronica di Cl (0.133
a.u.) `e inferiore al potenziale di ionizzazione di Na (0.189 a.u.). La differenza tra
le due situazioni, a piccole ed a grandi distanze, sta principalmente nell’attrazione
coulombiana tra Na+ e Cl− , che abbassa notevolmente l’energia dello stato ionico.
A distanze medie o grandi (R > 10 bohr), se trascuriamo la reciproca polarizzazione
di Na+ e Cl− , lo stato ionico ha un’energia pari a: H11 = IP (Na) − EA(Cl) − R1 ,
mentre, nella stessa scala l’energia dello stato covalente `e H22 = 0, praticamente
costante. Le due curve dovrebbero incrociarsi per R = 18.2 bohr, ma la regola
di non intersezione lo impedisce. In un caso tipico come questo, l’incrocio evitato
si riconosce facilmente esaminando le curve di potenziale oppure gli autostati di
ˆ el , detti anche stati adiabatici: questi ultimi cambiano rapidamente di natura,
H
il fondamentale da ionico a covalente e l’eccitato viceversa, passando attraverso la
zona dell’incrocio da sinistra a destra. Un altro esempio `e l’incrocio evitato tra le
curve B e B 0 di O2 ; le curve diabatiche sono tratteggiate in figura 7.3.
In prossimit`
a di un incrocio evitato, la differenza di energia tra le due curve di
potenziale diventa minima: se, per R = R∗ , H11 = H22 , allora |U2 − U1 | = 2H12 .
Le funzioni d’onda adiabatiche, ψk , cambiano rapidamente, perch´e il parametro θ
passa da circa 0 a π/2: per esempio, se ψ1 coincide approssimativamente con η1
per R < R∗ , diventa η2 per R > R∗ . Gli accoppiamenti nonadiabatici saranno
importanti, perch´e la derivata rispetto ad R delle funzioni d’onda `e pi`
u grande del
solito. Questo fatto, associato con la piccola differenza di energia, rende molto pi`
u
probabili le transizioni non radiative tra stati elettronici, in presenza di un incrocio
evitato.
Landau e Zener hanno fornito una prima trattazione di questo problema nel 1932,
basata sulla rappresentazione classica del moto dei nuclei. Consideriamo un’equazione di Schr¨
odinger per i soli elettroni, dove l’hamiltoniano elettronico dipende dal
185
tempo attraverso le coordinate dei nuclei:
i¯
h
d
ˆ el (t)ψel (t)
ψel (t) = H
dt
(7.83)
Le coordinate dei nuclei, Qr (t), sono funzioni del tempo determinate dalle leggi della meccanica classica. ψel (t) `e la funzione d’onda elettronica; anch’essa dipende dal
tempo, per due ragioni: indirettamente, attraverso la geometria molecolare, e direttamente, perch´e varia la probabilit`a di trovarsi nei vari stati elettronici. Possiamo
ˆ el ,
espandere ψel per mezzo delle funzioni adiabatiche (ψl ), cio`e le autofunzioni di H
oppure di quelle diabatiche (ηj ); nel primo caso scriviamo:
ψel (t) =
X
i
cl (t) e− h¯
Rt
0
Ul (t0 )dt0
ψl
(7.84)
l
Ricordiamo che gli integrali HIJ e gli autovalori Ul dipendono dal tempo a causa del moto
R t dei nuclei. Il fattore di fase dipendente dal potenziale elettronico
exp(−i/¯
h 0 Ul (t0 )dt0 ) sarebbe contenuto in ψel (t) anche se ci fosse un solo stato
elettronico, nel qual caso l’unico coefficiente cl non dipenderebbe dal tempo; si tratta di un fattore inessenziale, nel senso che non influisce sulla probabilit`a di trovare
il sistema nello stato elettronico ψl , data da |cl |2 ; la sua introduzione semplifica i
calcoli successivi. Sostituendo nella (7.83) si ottiene:
Rt
X dcl
i
i
dψl
0
0
(7.85)
ψl − cl (t)Ul (t)ψl + cl
e− h¯ 0 Ul (t )dt =
i¯
h
dt
h
¯
dt
l
=
X
i
cl e− h¯
Rt
0
Ul (t0 )dt0
l
ossia
X dcl
l
dt
i
ψl e− h¯
Rt
0
ˆ el ψl
H
Ul (t0 )dt0
=−
X
i
cl e− h¯
Rt
0
Ul (t0 )dt0
l
X dQr ∂ψl
dt ∂Qr
r
(7.86)
Premoltiplicando per ψk e integrando, otteniamo la derivata del singolo coefficiente:
Rt
X dQr (r)
X
i
0
dck
=−
cl (t)e− h¯ 0 (Ul −Uk )dt
g
dt
dt kl
r
(7.87)
l(6=k)
(r)
dove gkl `e l’accoppiamento nonadiabatico gi`a definito in eq.(7.32). Vediamo quindi
che la probabilit`a di saltare da uno stato adiabatico ad un altro dipende dalla velocit`a
Rt
(r)
del moto nucleare oltre che dall’accoppiamento gkl . La funzione exp(−i/¯h 0 (Ul −
Uk )dt0 ) compie rapide oscillazioni nel tempo se Ul − Uk `e grande; quando si integra
la (7.87), queste oscillazioni tendono ad annullare la probabilit`a di transizione se le
due superfici di potenziale sono lontane.
186
Se si espande ψel usando le funzioni d’onda diabatiche, gli integrali Hjj prendono il
posto dei potenziali Uk :
Rt
X
i
0
0
(7.88)
ψel (t) =
cj (t) e− h¯ 0 Hjj (t )dt ηj
j
Sostituendo la (7.88) nella (7.83) otteniamo:
Rt
X dcj
i
i
0
0
i¯
h
− cj (t)Hjj (t) e− h¯ 0 Hjj (t )dt ηj =
dt
h
¯
(7.89)
j
=
X
i
cj e− h¯
Rt
0
Hjj (t0 )dt0
j
ˆ el ηj
H
Qui non occorre considerare la derivata rispetto a t di ηj , poich´e si tratta di una
funzione d’onda quasi indipendente dalle coordinate nucleari e quindi dal tempo.
Premoltiplicando per ηi e integrando, otteniamo:
dci
i X − i R t (Hjj −Hii )dt0
Hij (t)cj (t)
e h¯ 0
=−
dt
h
(7.90)
j(6=i)
Le transizioni tra stati diabatici sono dovute ad accoppiamenti di tipo elettronico,
(r)
mentre quelle tra stati adiabatici sono dovute agli accoppiamenti “dinamici” gkl .
Anche qui, per`
o, una grande differenza di energia Hjj − Hii ha come conseguenza
una piccola probabilit`a di transizione.
Vediamo ora la soluzione approssimata del problema a due stati, usando la rappresentazione diabatica che porta ad equazioni pi`
u semplici. Le approssimazioni
sono valide se l’incrocio `e strettamente evitato, e quindi la zona in cui pu`
o avvenire la transizione viene attraversata in un tempo molto breve. Chiamiamo Q la
coordinata di avanzamento lungo la traiettoria e sia Q = Q∗ il punto di incrocio,
ˆ el in un intorno di Q∗ , ponendo:
dove H11 = H22 . Approssimiamo gli integrali di H
∗
H22 − H11 = F (Q − Q ) e H12 = costante. Assumiamo inoltre che la velocit`a rimanga costante durante l’attraversamento dell’incrocio, in modo che Q(t) = Q∗ + vt
(la scala dei tempi viene scelta arbitrariamente in modo che Q = Q∗ per t = 0).
Allora:
Rt
i
i
i
dc2
2
0 0
= − H12 e h¯ 0 F vt dt c1 (t) = − H12 eiF vt /2¯h c1 (t)
dt
h
¯
¯h
(7.91)
Questa equazione si potrebbe integrare esattamente, ma le approssimazioni gi`a introdotte per arrivare a questo punto si accordano con quelle della teoria perturbativa,
che facilita di molto il calcolo; infatti, un incrocio strettamente evitato presuppone
un piccolo accoppiamento H12 , quindi una piccola probabilit`a di transizione. Se il
sistema si trovava nello stato η1 al tempo iniziale, potremo porre approssimativamente c1 (t) = 1 a tutti i tempi. Per comodit`
a matematica, il tempo iniziale (molto
187
prima dell’attraversamento dell’incrocio) sar`
a posto a t = −∞, e quello finale a
t = +∞; allora:
Z
2
H12 +∞ iF vt2 /2¯h 2 2πH12
|c2 (f inale)| = e
dt =
h −∞
¯
¯ |vF |
h
2
(7.92)
L’estensione dei limiti d’integrazione a ±∞ `e giustificata se pensiamo che, lontano dall’incrocio, la probabilit`a di transizione cade rapidamente a zero; ci`o si riflette sul fatto che l’esponenziale rapidamente oscillante contribuisce sempre meno
all’integrale quando t2 cresce.
2
In conclusione, per piccoli accoppiamenti (H12
¯h |vF |), la probabilit`a di
transizione all’incrocio tra due stati diabatici `e:
P12 =
2
2πH12
h |vF |
¯
(7.93)
Essa `e proporzionale all’intensit`
a dell’accoppiamento elettronico ed inversamente
proporzionale alla differenza di pendenza tra le due curve ed alla velocit`a, due
parametri che determinano il tempo necessario ad attraversare l’incrocio. La generalizzazione non perturbativa della soluzione di eq.(7.91) porta alla regola di
Landau-Zener:
2
P12 = 1 − e−2πH12 /¯h|vF |
(7.94)
Le equazioni (7.93) e (7.94) coincidono fino al prim’ordine di uno sviluppo in
2 /¯
h |vF |.
H12
Ricordiamo che la probabilit`a di saltare da uno stato diabatico all’altro `e uguale
alla probabilit`a di rimanere nello stesso stato adiabatico, e viceversa. Perci`
o, in
rappresentazione adiabatica la probabilit`a di transizione `e:
2
P12 = e−2πH12 /¯h|vF |
(7.95)
I parametri essenziali che determinano la probabilit`a di transizione variano a seconda del sistema e del tipo di processo considerati. Negli incroci tra stati ionici
e neutri, H12 `e una funzione rapidamente decrescente della distanza internucleare:
infatti, le due funzioni d’onda coinvolte nell’integrale (7.79) differiscono per la collocazione di un elettrone in orbitali appartenenti all’uno o all’altro centro, e il valore
dell’integrale dipende dalla sovrapposizione dei due orbitali. Quel che conta `e il
valore di H12 nel punto di incrocio R∗ ; a sua volta la distanza R∗ dipende dal potenziale di ionizzazione dell’elettron donatore e dall’affinit`a elettronica dell’accettore:
R∗ ' 1/(IP − EA). Cos`ı, diverse coppie alogeno-metallo presentano incroci pi`
uo
meno evitati, con diverse probabilit`a di transizione. La velocit`a di attraversamento
dell’incrocio dipende dal tipo di processo ed eventualmente dalle curve di energia
potenziale. In tabella abbiamo i parametri per quattro diversi processi, coinvolgenti
188
due molecole per le quali sono state calcolate energie elettroniche e funzioni d’onda,
LiF e NaCl.
Le collisioni tra atomi neutri avvengono a velocit`a relative che obbediscono alla distribuzione di Maxwell, ossia con energia cinetica media pari a KB T /2. A bassissime
temperature (10−3 K), oggi raggiungibili in particolari apparati, la velocit`a `e cos`ı
bassa che la funzione d’onda elettronica ha tutto il tempo di adattarsi gradualmente
ai cambiamenti di posizione dei nuclei, anche attraversando un incrocio strettamente
evitato come quello di NaCl; la descrizione adiabatica (o di Born-Oppenheimer) `e
perfettamente valida. A temperatura ambiente, invece, una frazione importante dei
sistemi collidenti salta da una curva di potenziale all’altra. L’eccitazione verticale, secondo la versione classica del principio di Franck-Condon, consiste nel passare
dallo stato ψ1 allo stato ψ2 con i nuclei fissi alla distanza di equilibrio per U1 ; trovandosi in un potenziale repulsivo, i nuclei si allontanano acquistando energia cinetica,
e attraversano l’incrocio con una velocit`a notevolmente superiore a quella termica;
la funzione d’onda elettronica non fa in tempo ad adattarsi e rimane neutra, con
un’alta probabilit`a di saltare da una curva adiabatica all’altra; questo comportamento si pu`
o classificare come “prevalentemente diabatico”. Un risultato simile si
ottiene nella collisione di due ioni, che vengono accelerati dalla reciproca attrazione
coulombiana (in tutti i casi si `e trascurato per semplicit`a l’effetto del potenziale
centrifugo, considerando solo collisioni “frontali”).
189
Probabilit`a di transizione P12
per un incrocio evitato ionico/neutro.
Parametri in unit`
a atomiche.
LiF
IP (metallo)-EA(alogeno)
0.0732
13.65
distanza dell’incrocio, R∗
differenza di pendenza, F = (R∗ )−2
0.00537
0.00065
interazione, H12 (R∗ )
9340
massa ridotta, M
collisione fredda, T = 10−3 K
energia cinetica,
1.58 · 10−9
p EK = KB T /2
5.8 · 10−7
velocit`
a, v = 2EK /M
1.000
P12 neutro→ionico
P12 tra stati adiabatici
0.000
collisione termica, T = 300 K
energia cinetica,
4.75 · 10−4
p EK = KB T /2
3.2 · 10−4
velocit`
a, v = 2EK /M
0.789
P12 neutro→ionico
P12 tra stati adiabatici
0.211
eccitazione verticale
∗
energia cinetica,
0.0475
p EK = U2 (Re ) − U2 (R )
0.00319
velocit`
a, v = 2EK /M
P12 neutro→ionico
0.144
0.856
P12 tra stati adiabatici
collisione tra ioni
energia cinetica,
0.0732
p EK = IP − EA
velocit`
a, v = 2EK /M
0.00396
0.118
P12 ionico→neutro
0.882
P12 tra stati adiabatici
7.8
NaCl
0.0559
17.88
0.00313
0.00015
25285
1.58 · 10−9
3.5 · 10−7
1.000
0.000
4.75 · 10−4
1.9 · 10−7
0.217
0.783
0.0410
0.00180
0.026
0.974
0.0559
0.00210
0.022
0.978
Intersezioni coniche.
La situazione `e un po’ pi`
u complessa quando si passa a considerare molecole poliatomiche, quindi superfici di energia in molte dimensioni; in questo caso la regola di non
intersezione non `e strettamente valida. Consideriamo le due condizioni necessarie
per avere degenerazione di due stati, H11 = H22 e H12 = 0, e due coordinate interne,
X e Y . Ciascuna delle due condizioni `e un’equazione che lega X e Y : le soluzioni, se
esistono, definiscono una o pi`
u curve nel piano X, Y . Una curva tale che H11 = H22
pu`
o incrociarne un’altra tale che H12 = 0, in uno o pi`
u punti del piano: in ciascuno
di questi punti abbiamo degenerazione (U1 = U2 ), ma allontandosi in qualsiasi direzione dal punto di contatto le due superfici si separano una dall’altra, producendo
due coni che si toccano per i vertici; questa situazione `e detta intersezione conica.
190
Quando `e presente simmetria, la visualizzazione di un’intersezione conica pu`
o essere
pi`
u facile. Consideriamo un tipo di deformazione della geometria che lasci intatta la simmetria molecolare (coordinata X) ed un’altra coordinata (Y ) che invece
rimuove qualche elemento di simmetria. Per esempio, nel cammino di isomerizzazione cis-trans dell’azometano, CH3 − N = N − CH3 , possiamo conservare un asse
binario, variando solo l’angolo di torsione CNNC; altre coordinate, come la combinazione antisimmetrica dei due angoli di legame CNN, rompono la simmetria. Lo
stato elettronico fondamentale `e simmetrico rispetto all’asse binario, mentre il primo
singoletto eccitato (n → π ∗ ) `e antisimmetrico. Di conseguenza, finch´e non si rompe
la simmetria, H12 `e identicamente nullo: cio`e, H12 = 0 per Y = 0. Se, variando la
coordinata X, troviamo un punto in cui H11 = H22 , questo rappresenta un’intersezione conica. In effetti, poich`e l’energia dello stato fondamentale ha un massimo per
CNNC' 90◦ e quella dello stato eccitato ha invece un minimo, le curve si incrociano
due volte e danno luogo a due intersezioni coniche (figura 7.7). Se la molecola percorre una traiettoria lungo la coordinata X, con Y costante, un’intersezione conica
appare come un incrocio evitato; pi`
u piccolo `e Y , pi`
u vicino si passa al punto di vero
incrocio, e pi`
u alta `e la probabilit`a di transizione. In genere, se la molecola arriva
nelle vicinanze di un’intersezione conica, il decadimento `e molto rapido. Il tempo
di vita dello stato superiore `e circa pari al tempo necessario perch´e si produca la
deformazione richiesta per passare dalla geometria iniziale (solitamente, quella di
equilibrio nello stato fondamentale) a quella dell’intersezione conica. Ne consegue
che altri processi pi`
u lenti, come l’emissione di fotoni, risultano del tutto inibiti o
quasi: gli azocomposti non sono fluorescenti.
Un’altra possibilit`a `e che due stati elettronici siano degeneri per simmetria, come avviene nel caso di gruppi del punto non abeliani; questo caso `e frequente nei
complessi di metalli di transizione. Abbiamo allora necessariamente un’intersezione
conica, dato che uno spostamento lungo una coordinata non totalsimmetrica rimuove la degenerazione; il teorema di Jahn-Teller ci assicura che il minimo di energia per
lo stato pi`
u stabile non coincide con l’intersezione conica (vedi Atkins, sez. 10.11).
In generale, potremo classificare le intersezioni coniche in tre classi, come suggerito
da D. R. Yarkony in Acc. Chem. Res. 31, 511 (1998): quelle che coinvolgono due
stati della stessa simmetria; quelle permesse dalla presenza di simmetria e quelle
obbligate per simmetria (caso Jahn-Teller)
Considerando N coordinate interne, senza condizioni particolari imposte dalla simmetria, potremo avere degenerazione degli stati in un insieme di punti di dimensione
N − 2: uno o alcuni punti in due dimensioni, una o pi`
u curve in tre dimensioni, ecc.
L’insieme dei punti ove si ha degenerazione `e detto “crossing seam”. Come abbiamo gi`
a visto, l’occorrenza di un’intersezione conica determina la presenza di un
“funnel”, cio`e di una zona delle superfici di energia potenziale dove l’approssimazione di Born-Oppenheimer non `e valida, e le transizioni tra i due stati elettronici
coinvolti sono molto rapide e probabili. L’importanza del funnel nel meccanismo
di reazione fotochimico dipende dalla sua accessibilit`a. Spesso la punta del doppio
191
Figura 7.7: Intersezioni coniche nell’azometano.
kcal/mol
S
1
140
120
100
80
cis
60
S
0
trans
40
20
0
0
30
60
90
120 150 180
20
1015
0 5
-5
-10
bend asimmetrico
torsione
cono rappresenta un minimo per l’energia potenziale dello stato soprastante; questo
avviene quando, lungo tutte le coordinate (o almeno quelle totalsimmetriche, se c’`e
simmetria), le curve di potenziale dei due stati hanno pendenze opposte (peaked
crossing, nella classificazione di Ruedenberg e collaboratori, J. Chem. Phys. 95,
1862 (1991)). In questo caso, se la molecola si trova nello stato superiore, `e attratta
verso il funnel, e dopo la transizione allo stato sottostante, se ne allontana. In altri
casi, si trova qualche coordinata lungo la quale le due curve di potenziale hanno pendenza dello stesso segno; perci`
o il punto di intersezione non `e un minimo per nessuna
delle due superficie e pu`
o trovarsi ad energia pi`
u o meno alta rispetto ai rispettivi
minimi. Ne consegue una maggiore o minore difficolt`
a a raggiungere il funnel, con
una possibile dipendenza dall’energia disponibile, cio`e dalla lunghezza d’onda della
luce eccitante ed eventualmente dalla temperatura. Per esempio, nell’azometano,
secondo i calcoli pi`
u accurati finora eseguiti, l’intersezione conica S0 − S1 `e di tipo
“sloped” e giace ad energia un po’ pi`
u alta del minimo di S1 ; questo per`
o non ha
grande influenza sul tempo di vita dello stato eccitato, perch´e l’eccesso di energia
vibrazionale posseduto dalla molecola dopo l’eccitazione a S1 in questo caso `e pi`
u
192
` da notare che
che sufficiente per raggiungere le vicinanze dell’intersezione conica. E
in entrambi i tipi di intersezione conica, “peaked” o “sloped”, la superficie di energia potenziale sottostante ha sempre una pendenza tale da allontanare il sistema
dal punto di intersezione; questo fatto tende a rendere irreversibile la transizione in
molecole poliatomiche.
Un esempio un po’ pi`
u complesso `e fornito dalle olefine; qui, come negli azocomposti,
lo stato fondamentale ha un massimo di energia quando l’angolo di torsione del doppio legame raggiunge circa 90◦ , mentre il primo stato eccitato ha un minimo; le due
curve di energia per`
o non arrivano a incrociarsi e la differenza minima, nelle monoolefine, rimane dell’ordine di 60 kcal/mol (`e circa 170 kcal/mol quando la molecola
`e planare); tuttavia, un’ulteriore deformazione della geometria molecolare, come la
piramidalizzazione di un atomo di carbonio, accompagnata da una forte polarizzazione degli elettroni π, pu`
o abbassare ancora l’energia di S1 e soprattutto innalzare
quella di S0 : si giunge cos`ı in prossimit`a di un’intersezione conica; a seconda dei
sostituenti, questa pu`
o essere di tipo “peaked” oppure “sloped”. Di conseguenza,
la conversione da S1 a S0 `e molto rapida e neppure le monolefine sono fluorescenti.
L’esito finale pu`
o essere l’isomerizzazione cis-trans oppure il ritorno all’isomero di
partenza.
Possiamo dimostrare due propriet`a formali delle intersezioni coniche, con importanti
conseguenze per la dinamica e per la pratica computazionale. La prima propriet`a
consiste nell’esistenza di due sole coordinate che “levano la degenerazione” di un’intersezione conica; ossia, due sole direzioni lungo le quali le superfici si separano,
partendo da un punto di degenerazione. La seconda propriet`
a riguarda la fase delle
funzioni d’onda adiabatiche nelle vicinanze di un’intersezione conica.
Preventivamente, ricaviamo due formule che forniscono la derivata dell’energia elet(r)
tronica rispetto ad una coordinata Qr e l’accoppiamento gkl . In realt`
a, le espressioni
che troveremo potrebbero riferirsi ad un qualunque altro parametro contenuto nell’hamiltoniano elettronico, per esempio ad una componente di un campo esterno.
Partiamo dall’equazione di Schr¨
odinger per gli elettroni,
ˆ el |ψl i = Ul (Q) |ψl i
H
(7.96)
e deriviamo ambo i membri rispetto a Qr :
ˆ el
∂H
ˆ el ∂ |ψl i = ∂Ul |ψl i + Ul ∂ |ψl i
|ψl i + H
∂Qr
∂Qr
∂Qr
∂Qr
(7.97)
Moltiplicando a destra per hψl | si ottiene la formula per la derivata dell’energia
(teorema di Hellmann-Feynman):
+
* ∂H
∂Ul
ˆ el = ψl (7.98)
ψ
∂Qr l
∂Qr
193
Invece, moltiplicando per hψk |, con k 6= l, abbiamo l’espressione dell’accoppiamento
(r)
nonadiabatico gkl :
D ˆ E
Hel ψk ∂∂Q
ψl
∂ r
(r)
gkl = ψk ψ
=
(7.99)
l
∂Qr U −U
l
k
Queste due espressioni sono importanti dal punto di vista teorico; sono anche utilizzate in ambito computazionale, ma occorre tener conto che, nella dimostrazione,
abbiamo supposto di conoscere le autofunzioni esatte ψk e ψl , mentre di norma abbiamo solo autofunzioni approssimate. Il teorema di Hellmann-Feynman, eq. (7.98),
vale
per
D anche
E la derivata dell’energia relativa ad una funzione d’onda variazionale,
ˆ ∂ ψl Hel ψl /∂Qr , purch´e la ψl sia ottimizzata rispetto a tutti i parametri che
possono variare in funzione di Qr . L’espressione (7.99) mostra che l’accoppiamento
(r)
gkl ha un massimo in corrispondenza di un incrocio evitato e diverge nelle vicinanze
di un’intersezione conica.
Possiamo ora considerare il caso di due stati ψ1 e ψ2 , espressi nella base {η1 , η2 }
dagli autovettori
cosθ
(7.100)
C1 =
sinθ
e
C2 =
−sinθ
cosθ
(7.101)
~ x un vettore che individua il punto di
in accordo con le equazioni (7.80-7.81). Sia Q
~ x ) = U2 (Q
~ x ). In tale caso
intersezione nello spazio delle coordinate nucleari: U1 (Q
~ x ) − H11 (Q
~ x ) = 0 e H12 (Q
~ x ) = 0. Possiamo espandere H22 − H11 e
avremo H22 (Q
H12 in serie di Taylor in un intorno del punto di intersezione. Se ∆H = H22 − H11
~ −Q
~ x:
si ha, al primo ordine nello spostamento Q
~
~ x ) · (Q
~ −Q
~ x)
∆H = ∇∆H(
Q
~ 12 (Q
~ x ) · (Q
~ −Q
~ x)
H12 = ∇H
(7.102)
(7.103)
~1 e V
~2 i gradienti di H22 − H11 e di H12 nello spazio delle coordinate
Chiamiamo V
nucleari, valutati nel punto di intersezione:
~1 = ∇∆H(
~
~ x)
V
Q
(7.104)
~2 = ∇H
~ 12 (Q
~ x)
V
(7.105)
~1 e V
~2 definiscono un piano e solo i vettori ortogonali ad esso garantiscono
I vettori V
il mantenimento della degenerazione. In altre parole, partendo da un punto di in~1 , V
~2 )
tersezione conica, solo la componente dello spostamento che giace sul piano (V
194
~1 e V
~2 sono legati agli elemen`e efficace nel rimuovere la degenerazione. I vettori V
ˆ
ti di matrice di Hel nella base diabatica, che a loro volta non sono direttamente
determinabili con i metodi pi`
u in uso nella chimica teorica. Infatti, le funzioni d’onda diabatiche ηI sono in genere definite a partire da quelle adiabatiche ψK , ossia
ˆ el . Sono state avanzate diverse definiziogli autostati (di solito approssimati) di H
ni operative delle ηI , che richiedono direttamente o indirettamente di minimizzare
la loro variazione in funzione delle coordinate interne Qr ovvero gli accoppiamenti
hηI |∂/∂Qr | ηJ i (vedi per esempio M. Persico, voce “Electronic diabatic states” nella
“Encyclopedia of Computational Chemistry”). Allo scopo di individuare il piano
~1 , V
~2 ), `e pi`
(V
u semplice e conveniente metterlo in relazione con vettori legati direttamente alle energie ed agli stati adiabatici. Uno di essi `e il gradiente di U2 − U1 .
Infatti, dal teorema di H.-F. abbiamo:
∂H
∂UK
= C†K
CK
∂Qr
∂Qr
(7.106)
Ne consegue:
∂(H22 − H11 )
∂H12
∂(U2 − U1 )
= cos(2θ)
− 2sin(2θ)
∂Qr
∂Qr
∂Qr
(7.107)
Il gradiente di U2 − U1 `e dunque:
~ 1 = ∇(U
~ 2 − U1 ) = cos(2θ) V
~1 − 2sin(2θ) V
~2
W
(7.108)
(r)
L’altro `e dato dal vettore degli accoppiamenti nonadiabatici g12 , moltiplicati per
(U2 − U1 ) in modo da eliminare la divergenza nel punto di intersezione conica.
Infatti, dall’eq. (7.99) abbiamo:
(r)
(U2 − U1 ) g12 = C†1
∂H
∂(H22 − H11 )
∂H12
1
C2 = 2 sin(2θ)
+ cos(2θ)
∂Qr
∂Qr
∂Qr
(7.109)
Qui abbiamo trascurato l’elemento di matrice hηI |∂/∂Qr | ηJ i, che oltre ad essere
piccolo se le ηI sono funzioni diabatiche, `e moltiplicato per il fattore (U2 − U1 )
tendente a zero. Cos`ı possiamo definire il vettore:
~1 + cos(2θ) V
~2
~ 2 = (U2 − U1 ) ~g(r) = 1 sin(2θ) V
W
12
2
(7.110)
` evidente che i vettori W
~1 e W
~ 2 individuano lo stesso piano di V
~1 e V
~2 . Questa
E
propriet`a permette di implementare metodi per l’ottimizzazione condizionata della
geometria di un’intersezione conica, ossia per individuare il minimo dell’energia con
il vincolo della degenerazione, U1 = U2 . Infatti, dopo aver trovato un qualsiasi
punto sul “crossing seam”, per minimizzare l’energia mantenendo la degenerazione
~1 e W
~ 2 (naturalmente
basta variare le coordinate in una direzione ortogonale a W
servono tecniche efficienti per il calcolo delle derivate delle energie elettroniche e degli
(r)
accoppiamenti gkl ). Trovare il punto pi`
u basso di un crossing seam `e importante per
195
valutare la sua accessibilit`a e costituisce un’indicazione del cammino pi`
u probabile
che la molecola segue per passare da uno stato eccitato a quello sottostante. In
questo senso, il minimo di energia di un crossing seam svolge un ruolo analogo a
quello dello stato di transizione nelle reazioni termiche; tuttavia, non va dimenticato
che l’eccesso di energia normalmente disponibile in fotochimica rende i percorsi di
reazione meno “obbligati” e meno prevedibili sulla base delle sole propriet`a statiche
delle superfici di energia potenziale.
L’altra propriet`a notevole delle intersezioni coniche `e legata al modo in cui varia
una funzione d’onda adiabatica in funzione delle coordinate interne. Il teorema che
stiamo per dimostrare si deve ad un chimico teorico (H. C. Longuet-Higgins, Adv.
Spectrosc. 2, 429 (1961); vedi anche C. Alden Mead, Rev. Mod. Phys. 64, 51
(1992)). Pi`
u tardi la sua generalizzazione al caso di un hamiltoniano dipendente da
un qualsiasi set di parametri esterni (M. V. Berry, Proc. R. Soc. London A 392,
45 (1984)) ha trovato applicazioni in diversi campi della fisica. Consideriamo una
qualsiasi delle due funzioni d’onda (7.80-7.81), per esempio ψ1 . Il parametro θ `e
dato da:
p
2
∆H − ∆H 2 + 4H12
(7.111)
tgθ =
2H12
dove ∆H = H22 − H11 . Questa equazione assicura che i coefficienti cosθ e sinθ sono
definiti in qualsiasi punto dello spazio delle coordinate, tranne che nei punti in cui
∆H = H12 = 0. Supponiamo che il sistema descriva un percorso chiuso Γ nello
spazio delle coordinate.
Figura 7.8: Fase geometrica: curve chiuse Γ
Y
Y
γ
X
γ
X
Se Γ racchiude un’intersezione conica, possiamo deformare il cammino fino a ridurlo
ad un percorso di dimensioni infinitesime intorno al punto di intersezione, proiettato
196
~1 , V
~2 ). Infatti, qualunque spostamento infinitesimo al di fuori di questo
sul piano (V
piano non altera ∆H n´e H12 , cio`e le due variabili da cui dipende θ. Supponiamo per
~1 e V
~2 siano fra loro ortogonali (cosa che non altera la generalit`
semplicit`a che V
a della
dimostrazione). ∆H e H12 possono essere espansi in serie di Taylor degli spostamenti
x, y dal punto di intersezione. Per spostamenti infinitesimi si pu`
o troncare lo sviluppo
al primo ordine, ottenendo:
~ (7.112)
∆H = V
1 x
~ = V
2 y
(7.113)
~ Conviene a questo punto definire un sistema di assi (X, Y ) in modo che X = V
1 x e
~ Y = 2 V2 y e passare a coordinate polari (ρ, γ) definite da X = ρ cosγ, Y = ρ sinγ
(vedere la figura 7.8). Si ha dunque:
H12
tgθ =
cosγ − 1
= −tg(γ/2)
sinγ
(7.114)
Ne segue che seno e coseno di θ sono in modulo uguali a quelli di γ/2, e data
la parziale arbitrariet`a dei segni dei coefficienti di una funzione d’onda, possiamo
assumere cosθ = cos(γ/2) e sinθ = −sin(γ/2), ovvero θ = −γ/2. Si ha dunque:
ψ1 = cos(γ/2) η1 − sin(γ/2) η2
(7.115)
Se il cammino chiuso Γ non racchiude un punto di intersezione, l’angolo γ semplicemente ritorna al suo valore iniziale, e quindi la funzione d’onda non cambia. Se, al
contrario, la curva chiusa Γ racchiude un punto di intersezione, l’angolo γ passa dal
suo valore iniziale γi a γi ± 2π, e la funzione d’onda definita dalla (7.115) cambia di
segno.
Vediamo ad esempio il caso di H3 . Nella conformazione simmetrica D3h lo stato
elettronico fondamentale 2 E0 `e degenere; per il teorema di Jahn-Teller questa non
pu`
o essere la sua geometria di equilibrio: per trovarla si deve avere una deformazione
che abbassa la simmetria e rimuove la degenerazione. Si ha dunque una intersezione conica obbligata dalla simmetria. Consideriamo un cammino che gira intorno
all’intersezione in cui, ciclicamente, due dei tre atomi vengono avvicinati (vedere
figura 7.9). Supponiamo per semplificare che i tre atomi si trovino, per tutte le
geometrie considerate, a distanza relativamente grande l’uno dall’altro in modo tale
che la sovrapposizione fra i tre orbitali 1s dei tre H sia trascurabile, cos`ı come il contributo delle configurazioni “ioniche” (con due elettroni su uno stesso orbitale 1s)
alla funzione d’onda dello stato fondamentale. In questo modo possiamo limitarci a
considerare 3 determinanti di Slater:
Φ1 = k φa φb φc k
(7.116)
Φ3 = k φa φb φc k
(7.118)
Φ2 = k φa φb φc k
(7.117)
197
Ha
Hc
A
Hb
Ha
Hc
Hb
Ha
Ha
C
Hb
Hb
Hc
B
Hc
Figura 7.9: Cammino attorno al crossing seam in H3
dove φa (rispettivamente, φa ) indica uno spinorbitale 1s sull’atomo Ha con parte
di spin α (rispettivamente, β). Alla geometria D3h (indicata con “O”) i tre determinanti di Slater Φ1 , Φ2 e Φ3 fanno ottenere uno stato di quartetto, che non ci
interessa, e due doppietti degeneri, di simmetria 2 E0 :
ψ1 (O) =
ψ2 (O) =
1
√ (Φ1 − Φ2 )
2
1
√ (2Φ3 − Φ1 − Φ2 ) .
6
(7.119)
(7.120)
Alla geometria C2v indicata con “A” in figura 7.9 lo stato fondamentale si ottiene
appaiando gli spin degli elettroni su Ha e Hb in modo da ottenere un singoletto (che
198
naturalmente diventa un doppietto quando si aggiunge il terzo elettrone):
ψ1 (A) =
1
√ (Φ1 − Φ2 ) .
2
(7.121)
Vogliamo ora passare da “A” a “B” con continuit`
a, in modo che hψ1 (A) |ψ1 (B) i > 0.
Si avr`a allora:
ψ1 (B) =
1
√ (Φ1 − Φ3 ) .
2
(7.122)
Passiamo poi da “B” a “C”:
ψ1 (C) =
1
√ (Φ2 − Φ3 )
2
(7.123)
e torniamo da “C” ad ”A”:
ψ10 (A) =
1
√ (Φ2 − Φ1 ) .
2
(7.124)
In un giro completo attorno al punto di intersezione la funzione d’onda elettronica
adiabatica reale `e cambiata di segno: ψ10 (A) = −ψ1 (A).
Si dimostra, cio`e, che percorrendo con continuit`
a un percorso chiuso che racchiude
l’intersezione conica si ottiene il cambiamento di segno della funzione d’onda elettronica adiabatica (“fase geometrica” o “fase di Berry”). Ci`
o ha un’influenza non
trascurabile sul moto nucleare: poich´e ψ1 `e polidroma, anche la funzione d’onda
nucleare deve esserlo, dato che la funzione d’onda totale deve naturalmente essere
monodroma. Si noti anche che tale influenza si pu`
o manifestare anche se il punto di
intersezione non `e energeticamente accessibile.
La funzione d’onda elettronica adiabatica ψ1 definita dalla (7.115) pu`
o facilmente
essere resa monodroma, basta moltiplicarla per un fattore di fase dipendente dalle
coordinate nucleari: ψ˜1 = ψ1 eiγ/2 `e monodroma. Il prezzo che si paga in questa
(r)
operazione `e dato dal fatto che ψ˜1 `e una funzione complessa, con g˜11 6= 0. Infatti:
∂ (r)
˜
ψ˜1 = ψ˜1 ∂ ψ˜1 ∂γ = i ∂γ
(7.125)
g˜11 = ψ1 ∂γ ∂Qr ∂Qr
2 ∂Qr
Si arriva dunque alla stessa conclusione: la presenza dell’intersezione conica causa
un’interazione fra il moto elettronico e quello nucleare, perfino in assenza di effetti
nonadiabatici.
Per esempio, un pacchetto d’onda senza nodi, che incontri sulla sua strada un’intersezione conica, in genere si divider`
a in due componenti, una che rimane nello stato
iniziale ed una che transisce all’altro stato elettronico. Quest’ultima pu`
o anche essere piccola o trascurabile, se si tratta di un’intersezione conica di tipo “peaked” e il
pacchetto si muove nello stato pi`
u basso, aggirando l’intersezione conica su entrambi
i lati da una certa distanza. In ogni caso, le due frazioni del pacchetto che passano
199
Figura 7.10: Pacchetto d’onda che “aggira” un’intersezione conica. I simboli
WP1, WP2 e WP3 indicano il pacchetto, rappresentato come densit`a di probabilit`a
riportata sulla superficie di energia potenziale, a tre tempi successivi.
intersezione conica
WP1
WP2
WP2
WP3
Q_2
Q_1
da lati opposti rispetto all’intersezione conica restando sulla superficie di partenza
subiranno il fenomeno della fase geometrica; perci`
o, se torneranno a incontrarsi “a
valle” dell’intersezione conica, avranno segni opposti e produrranno interferenza distruttiva; nel suo complesso, la funzione d’onda acquister`
a un nodo, cio`e un certo
grado di eccitazione vibrazionale (vedi figura 7.10).
200
7.9
Regola di non intersezione e spazio di rimozione
della degenerazione: un approfondimento.
Nella sezione 7.7 abbiamo visto che per avere U1 = U2 occorre che siano verificate
le seguenti due condizioni:
H22 − H11 = 0
(7.126)
H12 = 0
(7.127)
Come si deduce immediatamente dalle equazioni (7.102) e (7.103), questo si traduce
nel fatto che, per una intersezione conica, lo spazio di rimozione della degenerazione
(branching space) `e bidimensionale. Notiamo per`
o che gli elementi di matrice fuori
ˆ
diagonale di Hel non sono necessariamente reali: le due condizioni sopra riportate
diventano cos`ı, in generale, tre condizioni, dato che le due componenti (reale e
immaginaria) di H12 devono essere poste separatamente uguali a zero. Ci sono
inoltre altri fattori da considerare, che derivano dalle propriet`a di simmetria per
inversione temporale delle funzioni d’onda elettroniche.
Indichiamo con Tˆ l’operatore di inversione temporale. Tˆ inverte lo stato di moto
di un sistema, lasciando inalterate le coordinate spaziali e cambiando di segno al
momento lineare ed angolare, spin compreso. Vediamo quali devono essere le propriet`
a di tale operatore. Consideriamo un sistema molecolare caratterizzato da un
ˆ Se H
ˆ `e l’hamiltoniano elettrostatico (6.26), ci aspettiamo che Tˆ
hamiltoniano H.
ˆ
ˆ sia simmetrico per inversione temporale), in quanto
commuti con H (ovvero che H
l’inversione del segno del momento lo lascia inalterato. Pi`
u in generale, questo sar`
a
vero per qualunque hamiltoniano molecolare descrivente un sistema isolato oppure
interagente con un campo conservativo. La simmetria rispetto all’inversione temporale viene invece a cadere in presenza, ad esempio, dell’interazione con un campo
magnetico. Consideriamo dunque un sistema molecolare simmetrico rispetto all’inversione temporale. Di conseguenza, si dovr`a ottenere lo stesso risultato se si inverte
lo stato di moto Ψ del sistema e poi si fa evolvere per un tempo ∆t, oppure se si
fa evolvere Ψ per un tempo −∆t e poi si applica l’inversione temporale. Ovvero,
tenendo conto che l’operatore di evoluzione temporale per un tempo infinitesimo ∆t
ˆ h (si veda la (6.32)):
`e dato da 1 − i∆tH/¯
(1 −
i∆t ˆ
i∆t ˆ ˆ
H)T Ψ = Tˆ(1 +
H)Ψ.
h
¯
¯h
(7.128)
Poich´e l’equazione soprascritta deve valere per uno stato di moto Ψ qualsiasi, abbiaˆ Tˆ = TˆiH.
ˆ Se Tˆ fosse un operatore lineare otterremmo dunque, eliminando
mo: −iH
ˆ Tˆ = TˆH,
ˆ contro l’assunto che il sistema sia simmetrico
l’unit`
a immaginaria, −H
ˆ Di conseguenza
rispetto all’inversione temporale (e quindi che Tˆ commuti con H).
ˆ
ˆ
ˆ
dobbiamo supporre che T i = −iT . Quindi T `e un operatore antilineare:
Tˆ(a |Ψ1 i + b |Ψ2 i) = a∗ Tˆ |Ψ1 i + b∗ Tˆ |Ψ2 i .
201
(7.129)
L’operazione di inversione temporale `e un’operazione di simmetria e, come tale,
deve conservare la norma. Tuttavia, a causa dell’antilinearit`
a, Tˆ non pu`
o essere un
operatore unitario; abbiamo infatti:
D
E
TˆΨ TˆΦ = hΨ |Φ i∗ .
(7.130)
Le due relazioni (7.129) and (7.130) definiscono Tˆ come operatore antiunitario. A
ˆ=H
ˆ Tˆ.
questo punto la (7.128) fornisce, correttamente, TˆH
Sia ψ una funzione d’onda elettronica; abbiamo il seguente risultato importante
(vedi per esempio M.J. Klein, Am. J. Phys. 20, 65 (1952)):
Tˆ2 ψ = (−1)N ψ
(7.131)
dove N `e il numero di elettroni. Sia ψ 0 = Tˆψ e supponiamo che ψ e ψ 0 siano
linearmente dipendenti, ovvero ψ 0 = aψ. Dunque:
Tˆ2 ψ = Tˆaψ = |a|2 ψ
(7.132)
Dalla (7.131) otteniamo |a|2 = (−1)N che pu`
o essere vero solo se il numero di
elettroni N `e pari. Dunque, con un numero dispari di elettroni, ψ e Tˆψ sono necesˆ el e TˆH
ˆ el = H
ˆ el Tˆ,
sariamente linearmente indipendenti. Ora, se ψ `e autostato di H
0
ˆ
ˆ
allora anche ψ = T ψ `e autostato di Hel , con lo stesso autovalore. Di conseguenza
abbiamo che, per un sistema simmetrico rispetto all’inversione temporale e con un
ˆ el sono doppiamente degeneri
numero dispari di elettroni, tutti gli autostati di H
(degenerazione di Kramers).
Possiamo a questo punto considerare tre differenti casi.
Caso A: sistema simmetrico rispetto all’inversione temporale con numero pari di elettroni. Abbiamo Tˆ2 = 1 e possiamo sempre scegliere una base completa ortonormale
di funzioni elettroniche {η} tale che Tˆηi = ηi . Si ha:
D E D
E D
E D E∗
ˆ ˆ ˆ
ˆ ˆ ˆ
ˆ
ˆ
=
ηi H
η
T
η
H
=
T
η
(7.133)
T
η
T
H
η
i
j
i
el j
el
el j = ηi Hel ηj
dove l’ultima uguaglianza segue dalla (7.130). Dunque, tutti gli elementi di matrice
ˆ el sono reali; le condizioni di degenerazione (7.126) e (7.127) portano a due sole
di H
equazioni. Questo `e il caso che `e stato considerato nelle sezioni 7.7 e 7.8.
Caso B : sistema simmetrico rispetto all’inversione temporale con numero dispari di
elettroni. Abbiamo Tˆ2 = −1; a causa della degenerazione di Kramers, un incrocio
fra due superfici di energia potenziale U1 e U2 coinvolge in realt`
a 4 stati anzich´e
ˆ
ˆ
due: ψ1 , ψ2 , T ψ1 e T ψ2 . Dobbiamo dunque considerare, oltre a η1 e η2 , anche Tˆη1
e Tˆη2 . Cominciamo con l’osservare che, tenendo conto della (7.130), del fatto che
ˆ el e Tˆ commutano e che gli elementi di matrice diagonali dell’operatore hermitiano
H
202
ˆ el sono reali:
H
E
E∗ D
E D
D
ˆ ˆ
2 ˆ ˆ2
2 ˆ ˆ ˆ
ˆ
ˆ
T
η
H
=
T
η
T
H
T
η
T
η
=
T
η
Tˆηi H
i
i el i
i
i
el
el D E
ˆ = ηi H
el ηi = Hii .
Procedendo nello stesso modo abbiamo anche:
E
D
ˆ ˆ
∗
= Hij
T
η
Tˆηi H
j
el
E
D
ˆ Tˆηi H
el ηi = 0
E
D
ˆ ∗
Tˆηi H
el ηj = −HiT j
(7.134)
(7.135)
(7.136)
(7.137)
E
D ˆ ˆ
ˆ
dove abbiamo definito HiT j = ηi H
el T ηj . La matrice di Hel nella base costituita
dai quattro stati η1 , η2 , Tˆη1 e Tˆη2 risulta quindi:


H11
H12
0
H1T 2
∗
 H12
H22
−H1T 2
0 


(7.138)
∗
∗ 
 0
−H1T 2
H11
H12
∗
H1T
0
H12
H22
2
La diagonalizzazione di tale matrice fornisce i due autovalori (doppiamente degeneri)
U1 e U2 . L’espressione che si ottiene `e analoga alla (7.82):
r
H11 + H22
∆H 2
U1,2 =
±
+ |H12 |2 + |H1T 2 |2
(7.139)
2
4
Dove ∆H = H22 − H11 . Per poter avere U1 = U2 dobbiamo dunque imporre le
seguenti 5 condizioni:
H11 = H22
Re{H12 } = 0
Re{H1T 2 } = 0
(7.140)
Im{H12 } = 0
(7.141)
Im{H1T 2 } = 0
(7.142)
dove Re e Im indicano, rispettivamente, la parte reale e immaginaria di un numero
complesso. Lo spazio di rimozione della degenerazione risulta dunque di dimensione
ˆ el contiene solo termini elettrostatici, allora H12
5. Notiamo d’altra parte che, se H
`e reale e H1T 2 = 0: ritorniamo allora al caso A (branching space di dimensione 2).
` questo il caso,
Caso C : sistema non simmetrico rispetto all’inversione temporale. E
ad esempio, di una molecola in un campo magnetico esterno. La degenerazione
ˆ el sono, in
di Kramers `e rimossa ma gli elementi di matrice fuori diagonale di H
generale, complessi. Le condizioni da imporre per avere degenerazione sono dunque
3: H11 = H22 , Re{H12 } = 0 e Im{H12 } = 0 (branching space tridimensionale).
203
7.10
Interazioni a lungo raggio tra singoletti eccitati.
In questa sezione esamineremo il meccanismo del trasferimento di energia tra singoletti eccitati. In tabella sono riportate alcune velocit`a di trasferimento di energia da
S1 a S1 . In tutti gli esempi l’energia dello stato eccitato del donatore `e superiore a
quella dell’accettore. Nei primi quattro la costante di velocit`a `e all’incirca quella che
ci si aspetta per un processo controllato dalla diffusione nel solvente usato (esano a
28 ◦ C, Kdiff ' 2.4 · 1010 l · mol−1 · s−1 ). Gli altri tre mostrano costanti di velocit`a pi`
u
grandi di almeno un ordine di grandezza.
Evidentemente devono agire due meccanismi di trasferimento dell’energia ben diversi. Nel primo caso le molecole devono arrivare “a contatto”; a distanza di interazione
non si pu`
o trascurare il fatto che gli elettroni sono indistinguibili, e che la funzione
d’onda complessiva `e antisimmetrica rispetto allo scambio di tutti gli elettroni di X
e di Y; cos`ı, questo tipo di interazione `e detta “di scambio”, e implica una modifica
pi`
u o meno rilevante delle funzioni d’onda dei due partner.
Le costanti di velocit`
a pi`
u alte corrispondono ad un’interazione a lunga distanza, tale
che le molecole, al di sopra di una certa concentrazione, non hanno neppure bisogno
di diffondere per scambiare energia: infatti la costante di velocit`a `e largamente
indipendente dalla viscosit`
a del solvente. Questo secondo tipo di interazione `e detta
“coulombiana”, perch´e `e dovuta essenzialmente alle propriet`a elettrostatiche degli
stati quantici di X e di Y, in particolare ai loro momenti di dipolo di transizione.
Il termine di hamiltoniano che descrive l’interazione tra X e Y a grandi distanze si
pu`
o sviluppare secondo i momenti multipolari delle due distribuzioni di carica; per
molecole neutre, il termine pi`
u importante `e quello dipolo-dipolo; l’accoppiamento
tra stato iniziale e stato finale sar`
a allora:
Z
µ⊥
~
µ⊥
µkX · ~µkY
X ·~
Y − 2~
∗
∗
ΨX (S0 )ΨY (S1 ) dqX dqY =
(7.143)
HXY ' ΨX (S1 )ΨY (S0 )
R3
=R
−3
Z
Ψ∗X (S1 )(~
µ⊥
X
Z
∗
ΨY (S0 )~µY ΨY (S1 ) dqY
− 2~
µX )ΨX (S0 ) dqX ·
k
Qui R `e la distanza tra X e Y, ossia la lunghezza di un vettore congiungente i “centri”
dei due cromofori; µ
~⊥ e ~
µk sono componenti del dipolo molecolare, rispettivamente
~
perpendicolare e parallela a R.
204
Trasferimenti di energia tra stati di singoletto.
Costanti di velocit`a in l · mol−1 · s−1 .
donatore
accettore
K
"b
n
b"
benzene
biacetile
3.3 · 1010
"b
"
b
toluene
biacetile
3.7 · 1010
bifenile
biacetile
2.2 · 1010
naftalene
biacetile
2.2 · 1010
perilene
2.0 · 1011
rubrene
2.0 · 1011
rubrene
3.0 · 1011
b
"
b"
CH3
O
@
@
C
TT
TT
T
T
"b "b
n n
b" b"
C
@
@
O
CH3
"b
"b
"
b"
b
b
"
b
"
"b
b"
b"
1-cloroantracene
b
"
"b
"b
b"
"
b
"b
"
b
b"
b"
Cl
"b
"b
"b
"
b"
b"
b
"b
"b
"
b
b"
b"
b"
1-cloroantracene
Φ
Φ
"b
"b
"b
"b
"
b"
b"
b"
b
b
"b
"b
"b
"
b"
b"
b"
b"
Φ
Φ
CN
"b
"b
"b
"
b"
b"
b
9-cianoantracene
b
"b
"b
"
b"
b"
b"
Secondo la regola d’oro di Fermi, la velocit`a del processo dipende dalla somma degli
2
HXY
per tutti gli stati vibronici finali possibili, con energia uguale o quasi a quella
dello stato iniziale. Gli stati finali differiscono tra loro per la quantit`a di energia
vibrazionale che viene attribuita alla molecola X piuttosto che ad Y, fermo restando
che la somma delle due `e costante:
EX (S1 , 0) = EX (S0 , v) + EY (S1 , v 0 )
(7.144)
Ora, hν = EX (S1 , 0) − EX (S0 , v) `e l’energia ceduta da X; essa deve coincidere con
EY (S1 , v 0 ), che `e l’energia ricevuta da Y. D’altra parte, gli integrali in eq. (7.143),
elevati al quadrato e sottoposti ad un’operazione di media orientazionale (basta
immaginare di far ruotare la molecola Y), sono proporzionali alle forze dell’oscillatore
per le transizioni che caratterizzano rispettivamente lo spettro di fluorescenza di X
205
e quello di assorbimento di Y; perci`
o la costante di velocit`a `e proporzionale a:
X
ν −2 fX (S1 , 0 → S0 , v) fY (S0 , 0 → S1 , v 0 )
(7.145)
R−6
v
ovvero a:
R−6
Z
ν −4 FX (ν)εY (ν) dν
(7.146)
dove FX `e l’intensit`
a di fluorescenza di X e εY il coefficiente di estinzione molare di
Y. In pratica l’interazione coulombiana `e efficace nel trasferimento di energia quando
entrambi gli spettri sono intensi ed hanno una buona sovrapposizione.
7.11
Accoppiamento eccitonico.
Quando una molecola contiene due o pi`
u cromofori che interagiscono, gli stati eccitati
(adiabatici) si possono descrivere come combinazioni lineari di eccitazioni localizzate,
cio`e di funzioni d’onda che descrivono stati eccitati del singolo cromoforo; fanno eccezione eventuali stati a trasferimento di carica, che sono presenti quando due cromofori hanno rispettivamente caratteristiche di elettron-accettore ed elettron-donatore
(vedi sezione 4.8). Non considereremo il mescolamento con lo stato fondamentale,
che normalmente `e abbastanza lontano in energia. Se anche gli stati eccitati sono
ben separati energeticamente, il mescolamento tra loro `e minimo, cio`e ogni autostato
esatto somiglia molto ad uno stato eccitato localizzato di un certo cromoforo. In
questo caso, ogni cromoforo pu`
o essere eccitato separatamente; se lo stato che vien
popolato non `e il pi`
u basso in energia, potr`
a aver luogo un trasferimento di eccitazione da un cromoforo all’altro; questo si pu`
o vedere come un normale fenomeno di
Internal Conversion Sn → S1 , ma con un ben definito spostamento dell’eccitazione
nello spazio (vedi sezione 5.3).
Qui ci occuperemo del caso in cui la differenza di energia tra gli stati eccitati di due
cromofori `e piccola, ossia paragonabile all’interazione tra gli stati stessi. Ci`
o pu`
o
avvenire, ad esempio, perch´e i due cromofori sono identici o simili, per struttura
ed intorno chimico. Adottando il linguaggio di sezione 7.7, chiamiamo η1 ed η2 i
due stati con eccitazione localizzata rispettivamente sui cromofori X e Y; infatti,
gli stati eccitati localizzati (“eccitoni”), sono un buon esempio di stati diabatici; in
conseguenza di variazioni nelle coordinate nucleari, le loro energie possono incrociarsi
oppure no; ma certamente sono stati che mantengono la stessa natura fisica, pur
mescolandosi con coefficienti variabili a seconda della geometria molecolare, per
produrre gli stati adiabatici, come vediamo qui di seguito. Indichiamo con H la
matrice che rappresenta l’hamiltoniano elettronico nella base degli stati η:
D E
ˆ (7.147)
Hij (R) = ηi H
el ηj
206
Qui H11 e H22 sono le energie di eccitazione dei due cromofori e H12 `e l’interazione.
Se altri stati non interferiscono, essendo abbastanza lontani in energia, avremo due
ˆ el rappresentabili come
autostati di H
ψ1 = cosθ η1 + sinθ η2
(7.148)
ψ2 = −sinθ η1 + cosθ η2
(7.149)
con autovalori:
E1,2
H11 + H22
=
±
2
s
H11 − H22
2
2
2
+ H12
(7.150)
(il segno + vale per E2 e il − per E1 ). Il “mescolamento” dipende dal parametro
p
2
H22 − H11 − (H22 − H11 )2 + 4H12
(7.151)
tg(θ) =
2H12
Indicando con ψ0 lo stato fondamentale, i dipoli di transizione che determinano
l’assorbimento di radiazione sono
µ0,1 = hψ0 |~
~
µ| ψ1 i = cosθ µ
~ X + sinθ ~µY
(7.152)
µ0,2 = hψ0 |~
~
µ| ψ2 i = −sinθ ~µX + cosθ µ
~Y
(7.153)
(qui per semplicit`a indichiamo con ~µX e ~µY i dipoli di transizione dei due cromofori). La somma delle due forze dell’oscillatore, trascurando la differenza di energia
E2 − E1 che abbiamo supposto piccola, `e proporzionale a µ20,1 + µ20,2 = µ2X + µ2Y e
quindi non dipende dal mescolamento degli stati. Invece, le intensit`
a delle due bande individuali degli stati ψ1 e ψ2 dipendono dal parametro θ (e quindi dal rapporto
(H22 − H11 )/H12 ), nonch´e dall’angolo formato dai due vettori ~µX e ~µY .
Per esempio, nel caso di due cromofori identici, avremo tg(θ) = ±1 ed energie di
transizione H11 ± H12 ; se l’orientazione di X e Y `e tale che ~µX e µ
~ Y sono perpendicolari (p. es., α = 90◦ e α0 = 0◦ in figura 7.11), le due bande avranno uguale intensit`
a,
proporzionale a µ2X = µ2Y . Se invece ~µX e ~µY sono paralleli (p. es., α = α0 = 90◦
e φ = 0◦ ), una transizione avr`a intensit`
a nulla (stato “dark”) e l’altra sar`
a proporzionale a 2µ2X (stato “bright”). In quest’ultimo caso, l’unica banda di assorbimento
sar`
a spostata di H12 rispetto a quella del cromoforo X (o Y) isolato; cio`e, avremo
uno shift batocromo (verso lunghezze d’onda pi`
u lunghe) se H12 < 0 e uno shift
ipsocromo se H12 > 0.
Se l’orientazione dei cromofori `e simmetrica, come nel caso dei bifenili o dei binaftili
(α = α0 in figura 7.11), abbiamo due forze dell’oscillatore proporzionali a
1
(~
µX
2
+~
µY )2 = sin2 α (1 + cosφ) µ2X
207
(7.154)
Figura 7.11: Accoppiamento eccitonico.
µX
~
µ
~Y
φ
~
R
α
α0
e
1
(~
µX
2
−~
µY )2 = [2cos2 α + sin2 α(1 − cosφ)] µ2X
(7.155)
L’associazione di queste espressioni con le due bande, quella con shift ipsocromo e
quella con shift batocromo, dipende dal segno di H12 , come gi`a visto sopra. Dallo
spettro si hanno quindi informazioni pi`
u o meno dirette sulla forza dell’interazione
e sulla geometria molecolare (altre informazioni si possono ottenere dallo spettro di
dicroismo circolare). Notare che il segno di un singolo elemento di matrice normalmente `e arbitrario, in quanto lo sono i segni delle funzioni d’onda; perci`
o, tale segno
non dovrebbe essere collegato a nessun effetto fisico. Tuttavia, qui abbiamo fissato i
segni di ~
µX e di ~
µY , come mostrato in figura 7.11, e di conseguenza i segni di η1 ed
η2 , relativamente a ψ0 ; cos´ı, l’arbitrariet`a sui segni delle funzioni d’onda `e sostanzialmente rimossa. Guardando la questione da un altro punto di vista, un cambiamento
di segno di uno dei tre stati coinvolti (ψ0 , η1 oppure η2 ), fa cambiare due delle tre
grandezze coinvolte (~
µX , ~
µY e H12 ); ma questo, come si verifica facilmente, non ha
nessun effetto sullo spettro.
Consideriamo ora il caso in cui l’interazione tra i due stati eccitati sia di tipo dipolodipolo; l’espressione (7.143) della sezione precedente si pu`
o riscrivere come:
H12 = R−3 µX µY (sinα sinα0 cosφ + 2cosα cosα0 )
(7.156)
Noti i moduli dei due dipoli di transizione (dalle forze dell’oscillatore per le bande di
X e Y separati), le sole variabili che determinano intensit`
a e posizione delle bande
sono quattro grandezze geometriche: la distanza R tra i cromofori, gli angoli tra i
~ e il diedro definito dagli stessi tre vettori (ossia l’angolo φ).
dipoli ~
µX e µ
~ Y ed R,
Per esempio, se la struttura molecolare implica che entrambi i dipoli siano paralleli
~ con versi opposti (α = α0 = 0), avremo:
a R,
H12 = 2R−3 µX µY
(7.157)
208
Nel caso simmetrico (X≡Y) lo stato “bright” avr`a necessariamente uno shift
batocromo.
~ con un angolo φ tra di loro:
Invece, se i dipoli sono perpendicolari a R,
H12 = R−3 µX µY cosφ
(7.158)
Nel caso simmetrico si applicano le formule (7.154-7.155) con α = α0 = 90◦ , e si
trova che la banda spostata a lunghezze d’onda pi`
u corte `e sempre la pi`
u intensa,
qualunque sia l’angolo φ.
Un caso tipico di accoppiamento dipolo-dipolo `e quello della bilirubina IX, costituita da due cromofori dipirrinonici che differiscono leggermente per la posizione del
sostituente vinilico, endo o exo (vedi figura 7.12). Le due bande di assorbimento,
molto larghe e sovrapposte, sono centrate a circa 420 e 460 nm, con una differenza
di energia tra i due stati eccitati stimata in circa 2000 cm−1 . La banda a 460 nm `e
pi`
u intensa (f =0.65) e quella a 420 nm appare come una spalla (f =0.55).
I dati riportati nel seguito sono in parte sperimentali e in parte desunti da calcoli. Presi come unit`
a separate, i due cromofori mostrano una differenza di energia
pi`
u piccola (circa 1200 cm−1 ) e intensit`
a di assorbimento invertite: quello “endo” assorbe a λ pi`
u lunghe con f =0.5, corrispondente ad un dipolo di transizione
µendo =2.6 a.u., mentre quello “exo” assorbe a λ pi`
u corte con f =0.8, ossia µexo =3.2
a.u. Entrambi gli effetti, spostamento delle bande e inversione delle intensit`
a, sono
prevedibili in base al modello di accoppiamento eccitonico.
Nella bilirubina IX, le due met`
a della molecola giacciono approssimativamente su
due piani perpendicolari tra loro, a causa dei legami a idrogeno presenti tra gruppi
ammidici e gruppi carbossilici in catena laterale. I dipoli di transizione sono grosso
modo allineati lungo ciascun sistema di doppi legami coniugati C=C. Stabilite con
una certa arbitrariet`a le posizioni dei baricentri dei due cromofori, risulta che la loro
distanza `e 6.3 ˚
A e gli angoli formati dai dipoli di transizione sono: α(endo) = 54◦ ,
α0 (exo) = 58◦ , φ = 84◦ . La formula (7.156) con questi dati fornisce H12 = 750
cm−1 , in buon accordo con l’aumento della differenza tra le energie di transizione
(da 1200 a 2000 cm−1 ), che si osserva passando dai cromofori isolati alla bilirubina.
Poich´e H12 `e positivo e i due dipoli di transizione formano tra loro un angolo ottuso
(104◦ ), si capisce anche perch´e la banda a lunghezze d’onda pi`
u corte `e indebolita e
l’altra `e intensificata.
Entrambi i sottosistemi endo ed exo possono isomerizzare intorno al loro doppio legame centrale; le rese quantiche sono basse (dell’ordine di 0.01 per l’exo, molto inferiore
per l’endo), probabilmente a causa dei legami a idrogeno che irrigidiscono la struttura. L’accoppiamento eccitonico rende possibile la fotoisomerizzazione del frammento
exo anche a λ > 460 nm, dove tale frammento, se isolato, non assorbirebbe.
209
Figura 7.12: Bilirubina IX. La linea tratteggiata indica l’intersezione tra i piani dei
due cromofori “endo” e “exo”, e le croci (+) i baricentri approssimativi dei cromofori
stessi. Le frecce mostrano i due corrispondenti dipoli di transizione.
H
H
C
H
H
C
C
C
C
C
H
N
C
H
+ ENDO
C
endo
C
H
55.1o
H
H
C
O
H
O
C
H
N
C
C
H
C
C
C
C
H
H
C
H
210
H
+ EXO
H
N
H
C
C
H
O
H
C
H
H
C
H
C
H
C
C
C
C
C
O
O
H
N
H
exo
H
C
H
63.6o
O