În mecanica clasică , teorema virială este o propoziție care leagă media în timp a energiei cinetice și a energiei potențiale a unui sistem stabil de particule N și care are implicații importante în diferite ramuri ale fizicii .
Prima formulare a teoremei se datorează lui Rudolf Clausius , în 1870 . Numele virial derivă din latina vis care înseamnă putere sau energie .
Teorema
Teorema virială afirmă că într-un sistem de N particule care se mișcă într-o regiune limitată a spațiului, a cărei energie cinetică totală este {\ displaystyle T} , raportul merită
- {\ displaystyle 2 \ left \ langle T \ right \ rangle = - \ sum _ {k = 1} ^ {N} \ left \ langle \ mathbf {F} _ {k} \ cdot \ mathbf {r} _ {k } \ right \ rangle}
unde parantezele indică media în timp și {\ displaystyle \ mathbf {F} _ {k}} reprezintă forța care acționează asupra k-a particulă, situată în poziție {\ displaystyle \ mathbf {r} _ {k}} .
Dacă energia potențială a sistemului este o funcție omogenă a gradului n al coordonatelor, adică a formei
- {\ displaystyle U (r) = \ alpha r ^ {n} \}
care este proporțională cu o putere n a distanței medii r dintre particule, atunci teorema ia forma
- {\ displaystyle 2 \ langle T \ rangle = n \ langle U \ rangle}
unde energia potențială totală medie {\ displaystyle \ left \ langle U \ right \ rangle} este suma energiei potențiale dintre fiecare pereche de particule.
În cazul particular al unui potențial gravitațional , proporțional cu reciprocitatea distanței, avem acest lucru
- {\ displaystyle 2 \ langle T \ rangle = - \ langle U \ rangle}
unde U este energia potențială gravitațională .
Demonstrație
Pentru a demonstra teorema, luați în considerare un sistem de mase {\ displaystyle m_ {i}} fiecare indicat printr-o rază vectorială {\ displaystyle \ mathbf {r} _ {i}} referindu-se la o anumită origine. Este {\ displaystyle \ mathbf {F} _ {i}} forța care acționează asupra primei mase. Indicând cu {\ displaystyle \ mathbf {p} _ {i}} impulsul masei i, atunci
- {\ displaystyle \ sum _ {i} \ mathbf {p} _ {i} \ cdot \ mathbf {r} _ {i} = \ sum _ {i} m_ {i} \ mathbf {v} _ {i} \ cdot \ mathbf {r} _ {i} = {\ frac {1} {2}} {\ frac {d} {dt}} \ sum _ {i} m_ {i} \ mathbf {r} _ {i} ^ {2}}
Ultima sumă, care este notată cu {\ displaystyle I} , este egal cu jumătate din urma tensorului de inerție , care corespunde momentului de inerție pentru o problemă bidimensională, în raport cu originea sistemului de masă. Derivând această expresie obținem:
- {\ displaystyle {\ frac {1} {2}} {\ frac {d ^ {2} I} {dt ^ {2}}} = \ sum _ {i} {\ dot {\ mathbf {p}}} _ {i} \ cdot \ mathbf {r} _ {i} + \ sum _ {i} \ mathbf {p} _ {i} \ cdot {\ dot {\ mathbf {r}}} _ {i} = \ sum _ {i} \ mathbf {F} _ {i} \ cdot \ mathbf {r} _ {i} + \ sum _ {i} m_ {i} \ mathbf {v} _ {i} \ cdot \ mathbf { v} _ {i} = \ sum _ {i} \ mathbf {F} _ {i} \ cdot \ mathbf {r} _ {i} + 2T}
Unde a fost folosită relația clasică {\ displaystyle \ mathbf {\ dot {p}} _ {i} = \ mathbf {F} _ {i}} . Indicat cu {\ displaystyle \ mathbf {F} _ {ij}} forța exercitată de masa a i-a pe masa a j-a și luând în considerare natura gravitațională a forțelor
- {\ displaystyle \ sum _ {i} \ mathbf {F} _ {i} \ cdot \ mathbf {r} _ {i} = \ sum _ {i} \ mathbf {r} _ {i} \ cdot \ sum _ {j \ neq i} \ mathbf {F} _ {ij} = \ sum _ {i} \ mathbf {r} _ {i} \ cdot \ sum _ {j \ neq i} Gm_ {i} m_ {j} {\ frac {\ mathbf {r} _ {j} - \ mathbf {r} _ {i}} {r_ {ij} ^ {3}}} = \ sum _ {j> i} {\ frac {Gm_ { i} m_ {j}} {r_ {ij} ^ {3}}} [\ mathbf {r} _ {i} \ cdot (\ mathbf {r} _ {j} - \ mathbf {r} _ {i} ) + \ mathbf {r} _ {j} \ cdot (\ mathbf {r} _ {i} - \ mathbf {r} _ {j})] =}
- {\ displaystyle = \ sum _ {j> i} {\ frac {Gm_ {i} m_ {j}} {r_ {ij} ^ {3}}} (\ mathbf {r} _ {j} - \ mathbf { r} _ {i}) (\ mathbf {r} _ {i} - \ mathbf {r} _ {j}) = - \ sum _ {j> i} {\ frac {Gm_ {i} m_ {j} } {r_ {ij}}}}
Ultima expresie este deci pur și simplu U, energia potențială gravitațională totală a sistemului de masă.
Prin urmare, am ajuns la următoarea expresie:
- {\ displaystyle {\ frac {1} {2}} {\ frac {d ^ {2} I} {dt ^ {2}}} = 2T + U}
iar teorema se obține prin urmare prin medierea ambilor membri. Având în vedere ipoteza limitării mișcărilor, media primului membru este zero, de fapt valoarea medie a oricărei funcții a timpului {\ displaystyle f (t)} este definit ca
- {\ displaystyle {\ bar {f}} = \ lim \ limits _ {T \ to + \ infty} {1 \ over T} \ int _ {0} ^ {T} {f \ left (t \ right)} dt}
De sine {\ displaystyle f (t)} este un derivat în raport cu timpul {\ displaystyle f \ left (t \ right) = {{dF \ left (t \ right)} \ over {dt}}} a unei funcții limitate {\ displaystyle F (t)} se dovedește
- {\ displaystyle {\ bar {f}} = \ lim \ limits _ {T \ to + \ infty} {1 \ over T} \ int _ {0} ^ {T} {{dF \ left (t \ right) } \ over {dt}} dt = \ lim \ limits _ {T \ to + \ infty} {{F \ left (T \ right) -F \ left (0 \ right)} \ over T} = 0}
Dovadă pentru energie dependentă de gradul coordonatelor
Deoarece energia cinetică {\ displaystyle T} este o funcție pătratică a vitezei, pe care o avem, prin teorema lui Euler asupra funcțiilor omogene
- {\ displaystyle \ sum \ limits _ {i} {\ frac {\ partial T} {\ partial \ mathbf {v} _ {i}}} \ cdot \ mathbf {v} _ {i} = 2T}
dacă acum introducem impulsurile
- {\ displaystyle {\ frac {\ partial T} {\ partial \ mathbf {v} _ {i}}} = \ mathbf {p} _ {i}}
și derivatele respective cu privire la timp conform ecuațiilor lui Newton
- {\ displaystyle \ mathbf {\ dot {p}} _ {i} = - {\ frac {\ partial U} {\ partial \ mathbf {r} _ {i}}}}
primesti
- {\ displaystyle 2T = \ sum \ limits _ {i} {\ mathbf {p} _ {i} \ cdot} \ mathbf {v} _ {i} = {\ frac {d} {dt}} \ left ({ \ sum \ limits _ {i} {\ mathbf {p} _ {i} \ cdot \ mathbf {r} _ {i}}} \ right) - \ sum \ limits _ {i} {\ mathbf {r} _ {i}} \ cdot \ mathbf {\ dot {p}} _ {i} = {\ frac {d} {dt}} \ left ({\ sum \ limits _ {i} {\ mathbf {p} _ { i} \ cdot \ mathbf {r} _ {i}}} \ right) + \ sum \ limits _ {i} {\ mathbf {r} _ {i}} \ cdot {\ frac {\ partial U} {\ parțial \ mathbf {r} _ {i}}}}
în virtutea teoremei lui Euler privind funcțiile omogene rezultă
- {\ displaystyle nU = \ sum \ limits _ {i} {\ mathbf {r} _ {i}} \ cdot {\ frac {\ partial U} {\ partial \ mathbf {r} _ {i}}}}
în timp ce pentru ipoteza limitării mișcărilor valoarea medie în raport cu timpul termenului
- {\ displaystyle {\ frac {d} {dt}} \ left ({\ sum \ limits _ {i} {\ mathbf {p} _ {i} \ cdot \ mathbf {r} _ {i}}} \ right )}}
este nul. Din aceasta rezultă afirmația
- {\ displaystyle 2 \ langle T \ rangle = n \ langle U \ rangle}
decât în cazul gravitațional, în care {\ displaystyle n = -1} , se reduce la declarația specială.
Teorema virială în mecanica cuantică
Tot în Mecanica cuantică există o variantă a teoremei viriale clasice.
Sunând cu {\ displaystyle | E \ rangle} o eigenstate referitoare la valoare proprie {\ displaystyle E} a hamiltonianului
- {\ displaystyle H = T (\ mathbf {p}) + U (\ mathbf {q})}
unde energia cinetică {\ displaystyle T (\ mathbf {p})} este întotdeauna o funcție a pătratelor impulsurilor și a energiei potențiale {\ displaystyle U (\ mathbf {q})} este încă o funcție omogenă a gradului {\ displaystyle n} coordonate {\ displaystyle \ mathbf {q}} , avem:
- {\ displaystyle 2 \ langle E | T | E \ rangle = n \ langle E | U | E \ rangle}
Demonstrație
În această demonstrație, pentru comoditatea scrierii, vom folosi convenția conform căreia, atunci când există doi indici repetați, este implicată o însumare pe indici înșiși, de exemplu:
- {\ displaystyle q_ {i} p_ {i} \ equiv \ sum \ limits _ {i} {q_ {i} p_ {i}}}
Pentru dovadă este util să se demonstreze preliminar următoarea egalitate:
- {\ displaystyle \ langle E | [{q} _ {i} {p} _ {i}, H] | E \ rangle = 0} .
Într-adevăr, amintindu-mi asta {\ displaystyle \ langle E | H = E \ langle E |} , se aplică următoarele:
- {\ displaystyle \ langle E | [{q} _ {i} {p} _ {i}, H] | E \ rangle = \ langle E | ({q} _ {i} {p} _ {i} HH {q} _ {i} {p} _ {i}) | E \ rangle = E \ langle E | {q} _ {i} {p} _ {i} | E \ rangle -E \ langle E | { q} _ {i} {p} _ {i} | E \ rangle = 0}
Acum putem demonstra versiunea cuantică a teoremei viriale:
- {\ displaystyle 0 = \ langle E | [{q} _ {i} {p} _ {i}, H] | E \ rangle = \ langle E | {q} _ {i} [{p} _ {i }, H] | E \ rangle + \ langle E | [{q} _ {i}, H] {p} _ {i} | E \ rangle = \ langle E | {q} _ {i} [{p } _ {i}, U] | E \ rangle + \ langle E | [{q} _ {i}, T] {p} _ {i} | E \ rangle}
unde ultima egalitate rezultă din faptul că
- {\ displaystyle \ left [{q_ {i}, U \ left ({\ mathbf {q}} \ right)} \ right] = \ left [{p_ {i}, T \ left ({\ mathbf {p} } \ right)} \ right] = 0}
Din proprietățile comutatorului poziție-moment rezultă că
- {\ displaystyle \ left [{q_ {i}, T \ left ({\ mathbf {p}} \ right)} \ right] = i \ hbar {\ frac {\ partial T} {\ partial p_ {i}} }}
- {\ displaystyle \ left [{p_ {i}, U \ left ({\ mathbf {q}} \ right)} \ right] = - i \ hbar {\ frac {\ partial U} {\ partial q_ {i} }}}
și din nou din teorema lui Euler asupra funcțiilor omogene, urmează
- {\ displaystyle {\ frac {\ partial T} {\ partial p_ {i}}} p_ {i} = 2T}
- {\ displaystyle {\ frac {\ partial U} {\ partial q_ {i}}} q_ {i} = nU}
Punând totul împreună, îl obții
- {\ displaystyle -i \ hbar \ langle E | nU | E \ rangle + i \ hbar \ langle E | 2T | E \ rangle = 0}
de aici și afirmația
- {\ displaystyle 2 \ langle E | T | E \ rangle = n \ langle E | U | E \ rangle}