Teorema del viriale

Template:F In meccanica classica, il teorema del viriale è una proposizione che lega la media temporale dell'energia cinetica e dell'energia potenziale di un sistema stabile di N particelle, e che ha importanti risvolti in diverse branche della fisica.

La prima formulazione del teorema è dovuta a Rudolf Clausius, nel 1870. Il nome viriale deriva dal latino vis che significa forza o energia.

Il teorema

Il teorema del viriale afferma che in un sistema di N particelle che si muovono in una regione limitata di spazio, la cui energia cinetica totale sia $T$, vale la relazione

$2⟨T⟩=-\sum _{k=1}^N⟨{𝐅}_k\cdot {𝐫}_k⟩$

dove le parentesi indicano la media temporale ed ${𝐅}_k$ rappresenta la forza che agisce sulla k-esima particella, situata nella posizione ${𝐫}_k$.

Se l'energia potenziale del sistema è una funzione omogenea di grado n delle coordinate, ovvero della forma

$U(r)=\alpha r^n$

cioè proporzionale ad una potenza n della distanza media r tra le particelle, allora il teorema assume la forma

$2⟨T⟩=n⟨U⟩$

dove l'energia potenziale totale media $⟨U⟩$ è la somma dell'energia potenziale tra ogni coppia di particelle.

Nel caso particolare di un potenziale gravitazionale, proporzionale al reciproco della distanza, si ha che

$2⟨T⟩=-⟨U⟩$

dove U è l'energia potenziale gravitazionale.

Dimostrazione

Per dimostrare il teorema si consideri un sistema di masse $m_i$ ognuna indicata da un raggio vettore ${𝐫}_i$ riferito ad una certa origine. Sia ${𝐅}_i$ la forza agente sulla massa i-esima. Indicando con ${𝐩}_i$ la quantità di moto della massa i-esima, allora

$\sum _i{𝐩}_i\cdot {𝐫}_i=\sum _i{m}_i{𝐯}_i\cdot {𝐫}_i=\frac{1}{2}\frac{d}{dt}\sum _i{m}_i{𝐫}_i^2$

L'ultima somma, che si denota con $I$, è pari a metà della traccia del tensore d'inerzia, che corrisponde al momento d'inerzia per un problema bidimensionale, rispetto all'origine del sistema di masse. Derivando questa espressione si ottiene:

$\frac{1}{2}\frac{d^{2}I}{d{t}^2}=\sum _i{\dot{𝐩}}_i\cdot {𝐫}_i+\sum _i{𝐩}_i\cdot {\dot{𝐫}}_i=\sum _i{𝐅}_i\cdot {𝐫}_i+\sum _i{m}_i{𝐯}_i\cdot {𝐯}_i=\sum _i{𝐅}_i\cdot {𝐫}_i+2T$

Dove si è usata la relazione classica ${\dot{𝐩}}_i={𝐅}_i$. Indicata con ${𝐅}_{ij}$ la forza esercitata dalla massa i-esima sulla massa j-esima e tenuto conto della natura gravitazionale delle forze

$\sum _i{𝐅}_i\cdot {𝐫}_i=\sum _i{𝐫}_i\cdot \sum _{j\ne i}{𝐅}_{ij}=\sum _i{𝐫}_i\cdot \sum _{j\ne i}G{m}_i{m}_j\frac{{𝐫}_j-{𝐫}_i}{r_{ij}^3}=\sum _{j>i}\frac{G{m}_i{m}_j}{r_{ij}^3}[{𝐫}_i\cdot ({𝐫}_j-{𝐫}_i)+{𝐫}_j\cdot ({𝐫}_i-{𝐫}_j)]=$

$=\sum _{j>i}\frac{G{m}_i{m}_j}{r_{ij}^3}({𝐫}_j-{𝐫}_i)({𝐫}_i-{𝐫}_j)=-\sum _{j>i}\frac{G{m}_i{m}_j}{r_{ij}}$

L'ultima espressione è quindi semplicemente U, l'energia potenziale gravitazionale totale del sistema di masse.

Siamo quindi giunti alla seguente espressione:

$\frac{1}{2}\frac{d^{2}I}{d{t}^2}=2T+U$

ed il teorema si ottiene quindi mediando entrambi i membri. Vista l'ipotesi di limitatezza dei moti, la media del primo membro è nulla, infatti il valor medio di una qualunque funzione del tempo $f(t)$ è definito come

$\overline{f}=\lim {\mathrm{\backslash limits}}_{T\to +\mathrm{\infty }}\frac{1}{T}\underset{0}{\overset{T}{\int }}f\left(t\right)dt$

Se $f(t)$ è una derivata rispetto al tempo $f\left(t\right)=\frac{dF\left(t\right)}{dt}$ di una funzione limitata $F(t)$ risulta

$\overline{f}=\lim {\mathrm{\backslash limits}}_{T\to +\mathrm{\infty }}\frac{1}{T}\underset{0}{\overset{T}{\int }}\frac{dF\left(t\right)}{dt}dt=\lim {\mathrm{\backslash limits}}_{T\to +\mathrm{\infty }}\frac{F\left(T\right)-F\left(0\right)}{T}=0$

Dimostrazione per energia dipendente dal grado delle coordinate

Poiché l'energia cinetica $T$ è una funzione quadratica delle velocità, si ha, per il Teorema di Eulero sulle funzioni omogenee

$\sum _i\frac{\partial T}{\partial {𝐯}_i}\cdot {𝐯}_i=2T$

se ora introduciamo gli impulsi

$\frac{\partial T}{\partial {𝐯}_i}={𝐩}_i$

e le rispettive derivate rispetto al tempo conformemente alle equazioni di Newton

${\dot{𝐩}}_i=-\frac{\partial U}{\partial {𝐫}_i}$

si ottiene

$2T=\sum _i{𝐩}_i\cdot {𝐯}_i=\frac{d}{dt}\left(\sum _i{𝐩}_i\cdot {𝐫}_i\right)-\sum _i{𝐫}_i\cdot {\dot{𝐩}}_i=\frac{d}{dt}\left(\sum _i{𝐩}_i\cdot {𝐫}_i\right)+\sum _i{𝐫}_i\cdot \frac{\partial U}{\partial {𝐫}_i}$

in virtù del Teorema di Eulero sulle funzioni omogenee risulta

$nU=\sum _i{𝐫}_i\cdot \frac{\partial U}{\partial {𝐫}_i}$

mentre per l'ipotesi di limitatezza dei moti il valor medio rispetto al tempo del termine

$\frac{d}{dt}\left(\sum _i{𝐩}_i\cdot {𝐫}_i\right)$

è nullo. Da ciò segue l'asserto

$2⟨T⟩=n⟨U⟩$

che nel caso gravitazionale, in cui $n=-1$, si riduce all'enunciato particolare.

Teorema del viriale in Meccanica Quantistica

Anche in Meccanica Quantistica si ha una variante del teorema del viriale classico.

Denominando con $|E⟩$ un autostato relativo all'autovalore $E$ dell'hamiltoniana

$H=T(𝐩)+U(𝐪)$

dove l'energia cinetica $T(𝐩)$ è sempre una funzione dei quadrati degli impulsi e l'energia potenziale $U(𝐪)$ è ancora una funzione omogenea di grado $n$ delle coordinate $𝐪$, si ha:

$2⟨E|T|E⟩=n⟨E|U|E⟩$

Dimostrazione

In questa dimostrazione, per comodità di scrittura, utilizzeremo la convenzione di Einstein secondo la quale, quando ci sono due indici ripetuti, si sottintende una sommatoria sugli indici stessi, ad esempio:

$q_i{p}_i\equiv \sum _i{q}_i{p}_i$

Per la dimostrazione è utile dimostrare preliminarmente la seguente uguaglianza:

$⟨E|[q_i{p}_i,H]|E⟩=0$.

Infatti, ricordando che $⟨E|H=E⟨E|$, vale:

$⟨E|[q_i{p}_i,H]|E⟩=⟨E|(q_i{p}_{i}H-H{q}_i{p}_i)|E⟩=E⟨E|q_i{p}_i|E⟩-E⟨E|q_i{p}_i|E⟩=0$

Possiamo ora dimostrare la versione quantistica del teorema del viriale:

$0=⟨E|[q_i{p}_i,H]|E⟩=⟨E|q_i[p_i,H]|E⟩+⟨E|[q_i,H]p_i|E⟩=⟨E|q_i[p_i,U]|E⟩+⟨E|[q_i,T]p_i|E⟩$

dove l'ultima uguaglianza segue dal fatto che

$\left[q_i,U\left(𝐪\right)\right]=\left[p_i,T\left(𝐩\right)\right]=0$

Dalle proprietà del commutatore posizione-momento segue che

$\left[q_i,T\left(𝐩\right)\right]=i\hslash \frac{\partial T}{\partial p_i}$

$\left[p_i,U\left(𝐪\right)\right]=-i\hslash \frac{\partial U}{\partial q_i}$

e di nuovo dal teorema di Eulero sulle funzioni omogenee segue

$\frac{\partial T}{\partial p_i}{p}_i=2T$

$\frac{\partial U}{\partial q_i}{q}_i=nU$

Mettendo tutto insieme si ottiene

$-i\hslash ⟨E|nU|E⟩+i\hslash ⟨E|2T|E⟩=0$

da cui l'enunciato

$2⟨E|T|E⟩=n⟨E|U|E⟩$