Cos'è ElectroYou | Login Iscriviti

ElectroYou - la comunità dei professionisti del mondo elettrico

0
voti

L'oscillatore armonico quantistico.

Una particella che vibra oscillando, comandata da una forza di richiamo puo' per esempio essere un elettrone che si avvicina e si allontana dal nucleo atomico e nel suo movimento puo' assumere diversi stati di energia. Una trattazione di tipo quantistico puo' darci evidenza matematica di questo fenomeno. Ricordiamo un elettrone che vibra e' responsabile dei fenomeni d'irradiazione del campo elettromagnetico a lunghezze d'onda dello stesso ordine di grandezza dell' l'oggetto irradiante. L'articolo dopo una prima stesura e' stato abbastanza modificato per renderlo di piu' semplice comprensione e lettura ed e' scaturito da una delle varie lezioni universitarie, che si sono tenute all'universita' di Roma nell'a.a. 2006-2007. Per i polinomi di Hermite si consulti la letteratura matematica, che ne e' nutritissima.


Indice

Trattazione quantistica

Una particella nel suo moto oscillatorio possedere energia potenziale e cinetica. Poiché siamo in meccanica quantistica, andiamo a vedere il suo Hamiltoniano corrispondente.

\hat H = \frac{\hat p}{2m} + \frac{1}{2} K \hat x^2

in cui \hat H, \hat p, \hat x sono gli operatori associati agli osservabili energia totale, quantita' di moto e posizione, mentre K e' la costante elastica dell'oscillatore. Poiche' \hat p = ih\frac{\rm{d}}{\rm{d}x} ed \hat x = x, si puo' costruire l'equazione di Schoendinger , indipendente dal tempo e stazionaria \hat H\psi(x)=E\psi(x). Sostituendo si ottiene:

 \frac{\rm{d}^2 u(x)}{\rm{d}x^2} + \frac{2 m}{\hbar}(E-\frac{K x^2}{2}) u = 0 ,

equazione monodimensionale, supponendo che la vibrazione sia nella direzione delle x.

In base a delle relazioni, che derivano dalle trasformate di Fourier, si ha:

\frac{d^2u}{dx^2} = -4 \pi^2 \int_{-\infty}^{\infty} \eta^2 \tilde{u}(\eta)exp( 2 \pi i \eta x) \rm{d}\eta
x^2 u(x) = -\frac{1}{4 \pi^2} \int_{-\infty}^{\infty} \frac{d^2 \tilde{u}}{d \eta^2} exp( 2 \pi i \eta x ) \rm{d} \eta

Sostituendo nell'equazione differenziale si ottiene:

- 4 \pi^2 \eta^2 \tilde{u} + \frac{2 m E}{\hbar^2}\tilde{u} + \frac{Km}{4 \pi^2 \hbar^2} \frac{d^2 \tilde{u}}{d \eta^2} = 0

relazione che ha la stessa forma di quella di provenienza, prima di fare le sostituzioni.

Le soluzioni sono quelle che si autoriproducono sotto trasformazione di Fourier e tra queste funzioni ci sono i polinomi Hermite - Gauss, di cui la piu' semplice e' la gaussiana. Per l'appunto prendiamo una forma della forma u0(x) = Aexp( − αx2) con A, α > 0.

sostituendo in questa ultima equazione differenziale, si ottiene:

A(-2 \alpha + 4 \alpha^2 x^2 + \frac{2 m E}{\hbar^2} - \frac{mK}{\hbar^2} x^2)exp(-\alpha x^2)

Affinche' valga questa uguaglianza deve essere:

\alpha=\frac{m E}{\hbar^2}
4 \alpha^2= \frac{m K}{\hbar^2}

Si ricava dunque:

\alpha=\frac{\sqrt{m K}}{2 \hbar}

L'energia relativa a questo stato fondamentale, descritto dall'autofunzione u0(x) e' pari a

E=\frac{1}{2}\hbar\sqrt{\frac{K}{m}}=\frac{1}{2}\hbar\omega

dove ω e' la pulsazione dell'oscillatore, relazione che si ottiene dai moti armonici in fisica classica.

Il valore di A si ottiene normalizzando l'integrale della densita' di probabilita' per la posizione, ossia:

\int_{-\infty}^{\infty}|u(x)|^2 \rm{d}x = A^2 \int_{-\infty}^{\infty}exp(-2 \alpha x^2) \rm{d}x=A^2\sqrt{\frac{\pi}{2 \alpha}}=1

risultando alla fine

A=\sqrt{\frac{2 \alpha}{\pi}}


Questo procedimento e' stato fatto per l'autofunzione di ordine zero, per gli ordini superiori l'equazione differenziale sofddisfa i polinimi di Hermite-Gauss di qualunque ordine e l'autofunzione prende la forma di:

u_n = A_n H_n(x \sqrt{2 \alpha} )\cdot  exp(-\alpha x^2)

per n=1,2...ecc.

Imponendo la normalizzazione, si trova An e quindi si ha:

u_n (x)= \sqrt[4]{\frac{2 \alpha}{\pi}} \frac{1}{\sqrt{2^n n!} } H_n(x \sqrt{2 \alpha}) exp(- \alpha x^2)

con i livelli di energia pari a:

E_n=(n+\frac{1}{2}) \hbar \omega

Esistoni quindi infiniti autostati con energia che variano a partire da un livello fondamentale per n=0, detta energia minima di punto zero. Il fatto che l'energia minima non sia nulla per lo stato fondamentale, porta alla concordanza con il principio d'indeterminazione, perche' ad ogni livello n non si sa perfettamente quanto valgano precisamente la posizione x e la quantita' di moto p. SE fossero determinate, esse dovrebbero essere entrambe zero. Si puo' dimostrare che il prodotto delle indeterminazioni Δx e Δp ha come valore minimo \frac{\hbar}{2}

Se scriviamo la densita' di probabilita' allo stato zero per la posizione x, come:

|u_0 (x)|^2 = A^2 exp[-\frac{x^2}{2(\Delta x)^2}], si vede che normalizzando risulta:

\Delta x=\frac{1}{2 \sqrt{\alpha}}

La trasformata di Fourier per la d.d.p. dell'impulso si ottiene dividendo per \hbar il modulo quadro della trasformata di Fourier di u0(x). Si ricorda che la trasformata di exp( − αx2) e' proporzionale ad exp(-\frac{\pi^2 \eta^2}{\alpha}) , per cui la d.d.p dell'impulso e' proporzionale a exp(-\frac{2 \pi^2 p^2}{\alpha \hbar^2})

ed eguagliando

exp (-\frac{p^2}{2 \Delta p^2}) = exp (-\frac{2 \pi^2 p^2}{\alpha \hbar^2})

Si ricava che

\Delta p = \hbar\sqrt{\alpha}

e

\Delta x \cdot \Delta p = \frac{\hbar}{2}

Le d.d.p.calcolate in modulo, piu' n e' elevato e piu' approssimano l'andamento parabolico deterministico.

Si osservino, per l'appunto le 3 figure sottostanti (Rif.: Università' degli studi di Padova, anno 2011, Corso di Chimica Fisica II, Prof.ssa Maria Brustolon, 7. L'oscillatore quantistico).

Fig.1. Funzioni d

Fig.1. Funzioni d'onda e densita' di probabilita'


Fig.2. Livelli energetici di densita

Fig.2. Livelli energetici di densita' di probabilita'


Fig.3. Rappresentazione delle d. d. p. All

Fig.3. Rappresentazione delle d. d. p. All'aumentare di n, approssimanti la parabola deterministica


Appendice

I polinimi di Hermite

Si distinguono due differenti tipi di polinomi di Hermite che sono entrambi equivalenti e ricavabili l'uno dall'altro con formule di carattere generale. Riportamo solo quelli che arrivano a n = 10. Si distinguono in polinomi probabilistici:

HE0(x) = 1
HE1(x) = x
HE2(x) = 4x2 − 1
HE3(x) = x3 − 3x
HE4(x) = x4 − 6x2 + 3
HE5(x) = x5 − 10x3 + 15x
HE6(x) = x6 − 15x4 + 45x2 − 15
HE7(x) = x7 − 21x5 − 105x3 − 105x
HE8(x) = x8 − 28x6 + 210x4 − 420x2 + 105
HE9x) = x9 − 36x7 + 378x5 − 1260x3 + 945x
HE10(x) = x10 − 45x8 + 630x6 − 3150x4 + 4726x2 − 946

ed in polinomi fisici:

HE0(x) = 1
HE1(x) = 2x
HE2(x) = 4x2 − 2
HE3(x) = 8x3 − 12x
HE4(x) = 16x4 − 48x2 + 120x
HE5(x) = 32x5 − 160x3 + 15x
HE6(x) = 64x6 − 480x4 + 720x2 − 120
HE7(x) = 128x7 − 1344x5 − 3360x3 − 1680x
HE8(x) = 2568 − 3584x6 + 13440x4 − 13440x2 + 1680
HE9x) = 512x9 − 9216x7 + 48384x5 − 30240x3 + 1680x
HE10(x) = 1024x10 − 23040x8 + 161280x6 − 403200x4 − 302400

L'oscillatore armonico classico

Un particella che vibra attorno ad un punto centrale, che e' l'origine in posizione x, puo' assumersi avente una legge di spostamento di tipo dinusoidale del tipo:

x = Asint + φ)

dove A e' l'elongazione, φ la fase iniziale per t=0 ed ω la pulsazione legata alla frequenza, con la relazione

ω = 2πf

Una vibrazione rispetto ad un punto centrale e' possibile perche esiste in qualche modo una forza di richiamo, che tende a riportare la particella verso l'origine. Per questo l'equazione dinamico - fisica di questo moto e' data da:

-Kx = m\frac{d^2 x}{dt^2}

dove K e' la costante elastica di moto armonico data da

K = mω2

Poiche la forza -Kx e' il gradiente del potenziale (assumiamo di stare in campo conservativo) ne risulta che l'energia potenziale relativa a questo moto oscollatorio vale

E_p = \frac{1}{2} K x^2

mentre l'energia cinetica e'

E_c = \frac{p^2}{2 m} = \frac{1}{2} m v^2 = \frac{1}{2}  m A^2\omega^2 cos(\omega t + \phi).

Conoscendo l'espressione di x e dopo una valutazione di una proprieta' trigonometrica, si ottiene che l'energia totale di un oscollatore armonico classico vale

E_t =\frac{1}{2}  m \omega^2 A^2 = \frac{1}{2} K A^2

Ossia, l'energia di una particella che si muove di moto armonico, a causa di una forza di richiamo e' proporzionale al quadrato dell'ampiezza massima.

0

Commenti e note

Inserisci un commento

Inserisci un commento

Per inserire commenti è necessario iscriversi ad ElectroYou. Se sei già iscritto, effettua il login.