Teoria perturbativa (meccanica quantistica)

studio del comportamento dell'hamiltoniana di un sistema quantistico sotto l'effetto di un potenziale perturbativo

In meccanica quantistica, la teoria perturbativa (o teoria delle perturbazioni) è un insieme di schemi di approssimazione legati all'omonima teoria matematica usati per descrivere un sistema quantistico complicato in termini di uno più semplice.

L'idea è cominciare con un sistema semplice, per il quale è nota una soluzione matematica, e aggiungere all'operatore hamiltoniano un termine "perturbativo", che rappresenti un disturbo lieve del sistema. Se la perturbazione non è troppo grande, le varie grandezze fisiche associate al sistema perturbato (ad esempio, i livelli energetici e gli autostati) possono essere espresse come "correzioni" di quelle del sistema semplice. Queste correzioni, poiché piccole rispetto alle grandezze stesse, possono essere calcolate usando metodi approssimati come lo sviluppo asintotico. Il sistema complicato può essere quindi studiato sulla base della conoscenza di quello semplice. Nei fatti, si descrive un sistema complicato senza soluzione per mezzo di un sistema semplice risolvibile.

Introduzione

modifica

La teoria perturbativa è un metodo di calcolo estremamente importante nella fisica moderna in quanto consente di descrivere i sistemi fisici quantistici reali, la cui quasi totalità è descritta da equazioni differenziali altrimenti difficilmente risolvibili in maniera esatta. Il metodo si basa sull'introduzione, nell'hamiltoniana, di una perturbazione, ovvero un potenziale così piccolo da giustificare uno sviluppo in serie di potenze. Ad esempio, aggiungendo un potenziale elettrico perturbativo all'hamiltoniana di un atomo di idrogeno - per la quale è stata trovata una soluzione esatta - si ottengono delle piccolissime variazioni nelle linee spettrali dell'idrogeno causate proprio dal potenziale perturbativo: questo effetto va sotto il nome di effetto Stark lineare.

Tuttavia, le soluzioni prodotte dalla teoria perturbativa non sono esatte, anche se sono estremamente accurate. Tipicamente i risultati sono espressi in termini di serie di potenze infinite che convergono rapidamente alla soluzione esatta man mano che ci si ferma nello sviluppo ad un ordine sempre più alto. Nella QED, dove l'interazione tra elettrone e fotone è trattata perturbativamente, il calcolo del momento magnetico dell'elettrone è stato determinato in accordo con il dato sperimentale fino all'undicesima cifra decimale. In QED e in altre teorie di campo quantistiche, speciali tecniche di calcolo note come diagrammi di Feynman sono utilizzate per sommare i termini delle serie di potenze.

In certe condizioni, la teoria perturbativa non può essere utilizzata; questo perché il sistema che si vuole descrivere non può essere descritto con l'introduzione di una perturbazione in una situazione ideale libera. In cromodinamica quantistica (QCD), ad esempio, l'interazione tra i quark ed il campo gluonico non può essere trattata perturbativamente a bassa energia a causa del fatto che essa diventa troppo grande. La teoria perturbativa, inoltre, non va bene per descrivere stati che non sono generati con continuità, incluse le condizioni al contorno e i fenomeni collettivi noti come solitoni.

Tra i sistemi che possono essere trattati con la teoria perturbativa vi sono inoltre la struttura fine dell'atomo di idrogeno e degli idrogenoidi, l'effetto Zeeman ed il limite di Paschen-Back. Inoltre, con le tecniche di simulazione moderne, si è in grado di applicare la teoria perturbativa a molti sistemi sempre più complicati, ottenendo delle buone soluzioni numeriche.

A fianco della teoria perturbativa indipendente dal tempo c'è anche la teoria perturbativa dipendente dal tempo, nella quale si considerano sia potenziale sia, soprattutto, soluzioni dipendenti dal tempo. Esistono, infine, altri metodi per ottenere soluzioni approssimate del problema agli autovalori per una data hamiltoniana tra i quali i più importanti sono il metodo variazionale e l'approssimazione WKB.

Perturbazioni stazionarie

modifica

Si considerano perturbazioni stazionarie, analizzando i casi in cui lo spettro sia degenere e non degenere.

Spettro non degenere

modifica

Questo grafico non è disponibile a causa di un problema tecnico.
Si prega di non rimuoverlo.

Si consideri un sistema libero, sia cioè:

 

dove |n0〉 è un sistema di autostati, espresso con la notazione bra-ket, ortonormale e completo, ossia tale per cui valgono le identità:

  (completezza)
  (ortonormalità)

Supponiamo di inserire un potenziale nel sistema, e si consideri il caso in cui lo spettro dell'hamiltoniana sia non degenere, ossia tale che per ogni autovalore vi sia uno e un solo autostato. Il potenziale rappresenta una perturbazione di tipo additivo rispetto allo stato libero:

 

con λ reale positivo compreso tra 0 e 1.

Il problema degli autovalori diventa quindi:

 

introducendo la quantità Δn = En - En(0), si ottiene:

 

che opportunamente riscritta diventa:

 

A questo punto si pone il problema dell'invertibilità dell'operatore

 

il suo inverso avrà una singolarità sull'autostato |n0〉, mentre

 

non presenta componenti lungo |n0〉 in quanto appartenente ad uno spazio ortogonale a quello di |n0〉. Formalmente questo fatto può essere espresso introducendo l'operatore Φn, che è un proiettore sullo spazio ortogonale ad |n0〉:

 

Con l'introduzione del proiettore, l'autostato diventa quindi:

 

dove è stato aggiunto il ket |n0〉 perché, se λ tende a zero, l'autostato dell'hamiltoniana perturbata deve tendere all'autostato libero. Ovviamente |n〉 andrà opportunamente normalizzato.

A questo punto è necessario calcolare la distanza (scostamento) tra lo stato imperturbato e quello perturbato. Innanzitutto si vede che:

 

poiché stiamo introducendo una perturbazione, il fattore λ sarà piccolo (prossimo allo zero) e quindi Δn può essere espresso attraverso una serie di potenze di λ, della quale interessano solo pochi termini

 

allo stesso modo per l'auto-ket:

 

con il che si può scrivere finalmente:

 

confrontando gli n ordini simili si ottengono altrettante relazioni, una per ciascun ordine:

 

Posto

 

le variazioni agli ordini successivi consentono di scrivere sia lo spostamento in energia tra i livelli, sia l'auto-ket

 
 

A questo punto, essendo Vii nullo e Vij non nullo, e ordinando i livelli in modo tale che Ei(0) > Ej(0), si ottiene che lo scostamento i-esimo è positivo e quello j-esimo negativo: quindi i livelli tendono ad allontanarsi.

 
 

Infine, nel caso dello stato fondamentale n, si può notare che la sua energia si abbassa sempre:

 

Dimostrazione alternativa[1]

modifica

Supponiamo al solito che l'hamiltoniano sia del tipo:

 

dove   è l'hamiltoniano imperturbato, cioè tale che:

 

in cui   sono un insieme completo di autofunzioni dell'operatore imperturbato, e   una perturbazione. Vogliamo trovare la soluzione dell'equazione di Schrödinger:

 

Supponiamo che lo spettro degli autovalori sia non degenere e supponiamo di avere sviluppato la nostra funzione d'onda:

 

Sostituiamo la (4) nella (3):

 

moltiplicando per   otteniamo formalmente i coefficienti  :

 

Finora non abbiamo eseguito approssimazioni. Ora sviluppiamo i valori dell'energia e dei coefficienti in serie:

 
 

dove l'apice indica l'ordine di grandezza, (0) indica l'ordine di  , (1) indica l'ordine di   e così via. Imponiamo le due condizioni:

 .

La prima approssimazione è data da  , ponendo  :

 
 

per normalizzare   occorre porre   così:

 

dove nella sommatoria non va considerata la somma per  . Questa prima approssimazione fornisce anche la condizione di approssimazione, cioè deve essere:

 

Determiniamo la seconda approssimazione:

 
 

allora:

 

Spettro degenere

modifica

Nel caso in cui lo spettro sia degenere si ha che ad un valore dell'energia corrispondono più autovettori, ovvero un autospazio di dimensione superiore a 1:

 

con   il sistema di autovettori con autovalore ED(0).

In questo caso la procedura poc'anzi descritta perde la sua validità e non riesce a descrivere correttamente il sistema.

Si supponga che Vm m' = 0, con

 

A questo punto si passa dalla base   a quella  , tale che:

 
 

Quindi, affinché la prima condizione sia valida, deve verificarsi che:

 

Inserendo la completezza nel termine di destra si ottiene alla fine:

 

La perturbazione sarà rappresentabile da una matrice diagonale, quindi in caso di degenerazione ci si può limitare a studiare il comportamento per un multipletto degenere: si costruisce la matrice Vm m', la si diagonalizza e si trovano gli autostati corrispondenti. Questi saranno una nuova base per rappresentare V, che sarà così una matrice diagonale con gli autovalori λ1, ..., λg come elementi della diagonale, i quali sono anche il rango minimo per la Δl.

A questo punto si potrà procedere come nel caso non degenere, utilizzando l'accortezza di sommare sui k che non appartengono alla degenerazione ED(0).

Teoria delle perturbazioni dipendenti dal tempo

modifica
  Lo stesso argomento in dettaglio: Serie di Dyson.

Perturbazioni dipendenti dal tempo possono essere trattate con la tecnica della serie di Dyson. Partendo dall'Equazione di Schrödinger:

 

questa ammette una soluzione nella forma

 

essendo   l'operatore di ordinamento temporale tale che

 

se   e

 

se   cosicché l'esponenziale rappresenti la serie di Dyson

 .

Consideriamo ora il seguente problema perturbativo

 

e assumiamo che il parametro   sia piccolo e di sapere risolvere il problema  .
Operiamo la trasformazione unitaria passando alla cosiddetta rappresentazione di interazione o rappresentazione di Dirac:

 

cosicché l'Equazione di Schrödinger diviene

 

che può essere risolta con la sopracitata serie di Dyson

 
 

che una è serie perturbativa a   piccolo.
Utilizzando le soluzioni del problema non perturbato   e   (possiamo assumere, per esemplificare, uno spettro puramente discreto), avremo fermandoci al primo ordine

 .

Perciò, se il sistema si trovava inizialmente nello stato del problema imperturbato  , l'ampiezza di probabilità che per effetto della perturbazione si vada a trovare nello stato del sistema imperturbato   al primo ordine sarà data da

 

e la probabilità di transizione per unità di tempo sarà data dalla Regola aurea di Fermi.

La teoria delle perturbazioni indipendenti dal tempo può essere derivata dalla teoria delle perturbazioni dipendenti dal tempo.
Per questo scopo scriviamo l'operatore di evoluzione cronologica ottenuto con la serie di Dyson come

 
 

e assumiamo che la perturbazione   sia indipendente dal tempo. Utilizziamo l'identità

 

dove   assumendo per semplicità che lo spettro sia puramente discreto, è possibile riscrivere l'operatore di evoluzione temporale come

 
 

da cui si vede che, al secondo ordine, occorre sommare su tutti gli stati intermedi. Assumiamo ora   e il limite asintotico a tempi grandi. Questo comporta che, per ogni contributo della serie perturbativa, dovremo aggiungere un fattore moltiplicativo   nell'integrando cosicché, nel limite  , noi possiamo ritrovare lo stato finale del sistema eliminando tutti i termini oscillanti ma conservando quelli secolari.   deve essere scelto arbitrariamente piccolo. In questo modo possiamo eseguire gli integrali e, separando i termini diagonali dagli altri, otteniamo

 
 
 

dove la serie secolare nel tempo produce la serie relativa agli autovalori del problema perturbato e la parte restante dà la correzione alle autofunzioni. L'operatore di evoluzione cronologica viene applicato ad uno autostato qualsiasi del problema non perturbato e nel caso in questione la serie perturbativa produce una serie secolare, ossia valida a tempi piccoli.

Teoria delle perturbazioni forti

modifica

In maniera del tutta analoga al caso delle piccole perturbazioni, è possibile sviluppare una teoria delle perturbazioni forti. Si consideri l'Equazione di Schrödinger:

 

La serie adiabatica è una serie di Dyson duale a quella del caso precedente che si applica nel limite in cui una perturbazione diviene infinitamente grande.[2][3] L'approccio è generale e può essere esemplificato nel modo seguente. Si consideri il problema perturbativo:

 

con  . Il nostro scopo è di trovare una soluzione del tipo

 

ma una sostituzione diretta nell'equazione considerata non produce risultati utili. La situazione può essere accomodata effettuando un riscalamento della variabile tempo come   e questo produce la serie di equazioni perturbative

 
 
 

che si risolve conoscendo la soluzione dell'equazione all'ordine principale. Ma già abbiamo visto che per questa possiamo utilizzare l'approssimazione adiabatica. Nel caso particolare in cui   non dipenda dal tempo si ha la serie di Wigner-Kirkwood spesso utilizzata in meccanica statistica. In questo caso infatti, possiamo introdurre una trasformazione unitaria come

 

che definisce una rappresentazione libera, poiché stiamo cercando di eliminare il termine di interazione. A questo punto, in modo duale al caso delle piccole perturbazioni, dobbiamo risolvere l'equazione di Schrödinger

 

in cui vediamo che il parametro di espansione   compare solo nell'esponenziale e dunque, la corrispondente serie di Dyson, una serie di Dyson duale, è significativa per grandi valori di   ed è

 
 

che, dopo il cambiamento di scala nel tempo   scopriamo essere una serie in   giustificando il nome di serie duale di Dyson. Tale serie si è infatti ottenuta semplicemente cambiando la scelta della perturbazione scambiando   con  . Questo principio è detto principio di dualità in teoria delle perturbazioni. La scelta   produce, come detto, una serie di Wigner-Kirkwood, che è una serie di gradiente. La serie di Wigner-Kirkwood è una serie semiclassica con gli autovalori determinati allo stesso modo che per l'approssimazione WKB[4].

  1. ^ Landau Lev D., Lifsits Evgenij M., Fisica Teorica 3 - Meccanica quantistica: teoria non relativistica, in Editori Riuniti.
  2. ^ Ali Mostafazadeh, Quantum adiabatic approximation and the geometric phase, Phys. Rev. A 55, 1653 (1997).
  3. ^ Marco Frasca, Duality in Perturbation Theory and the Quantum Adiabatic Approximation, Phys. Rev. A 58, 3439 (1998).
  4. ^ Marco Frasca, A strongly perturbed quantum system is a semiclassical system, Proc. R. Soc. A 463, 2195 (2007).

Bibliografia

modifica
  • Simmonds, Mann, A First Look at Perturbation Theory, Dover reprint of Krieger, Malabar (FL) 1986.
  • Bender, Orszag, Advanced mathematical methods for scientits and engineers: Asymptotic Methods and Perturbation Theory, Springer 1999.

Voci correlate

modifica

Collegamenti esterni

modifica
Controllo di autoritàLCCN (ENsh85100182 · J9U (ENHE987007538619405171
  Portale Quantistica: accedi alle voci di Wikipedia che trattano di quantistica