În mecanica cuantică , mișcarea într-un câmp central este tipică pentru două particule care interacționează supuse unui potențial dependent de distanța reciprocă a ambelor. Problema poate fi redusă la o singură problemă a particulelor, așa cum se întâmplă în cazul clasic.
Problema celor două corpuri
Să luăm două particule de masă {\ displaystyle m_ {1}} și {\ displaystyle m_ {2}} interacționând printr-un potențial {\ displaystyle V (r)} unde r este distanța relativă a celor două particule. Operatorul hamiltonian este:
- {\ displaystyle H = {\ frac {{\ vec {p}} _ {1} ^ {2}} {2m_ {1}}} + {\ frac {{\ vec {p}} _ {2} ^ { 2}} {2m_ {2}}} + V (| {\ vec {r}} _ {1} - {\ vec {r}} _ {2} |)}
Introducem noile variabile:
- {\ displaystyle {\ vec {R}} = {\ frac {m_ {1} {\ vec {r}} _ {1} + m_ {2} {\ vec {r}} _ {2}} {m_ { 1} + m_ {2}}}}
- {\ displaystyle {\ vec {r}} = {\ vec {r}} _ {1} - {\ vec {r}} _ {2}}
unde este {\ displaystyle {\ vec {R}}} este vectorul de poziție al centrului de masă al celor două particule și {\ displaystyle {\ vec {r}}} este distanța reciprocă a celor două particule. Odată cu introducerea acestor două variabile, trebuie să transformăm operatorii de impuls și astfel Laplacianul energiei cinetice. Raportăm transformările:
- {\ displaystyle {\ frac {\ partial} {\ partial r_ {1}}} = {\ frac {m_ {1}} {m_ {1} + m_ {2}}} {\ frac {\ partial} {\ partial R}} + {\ frac {\ partial} {\ partial r}}}
- {\ displaystyle {\ frac {\ partial} {\ partial r_ {2}}} = {\ frac {m_ {2}} {m_ {1} + m_ {2}}} {\ frac {\ partial} {\ partial R}} - {\ frac {\ partial} {\ partial r}}}
Gradienții (proporționali cu impulsurile) devin:
- {\ displaystyle {\ vec {\ nabla}} _ {1} = {\ frac {m_ {1}} {m_ {1} + m_ {2}}} {\ vec {\ nabla}} _ {R} + {\ vec {\ nabla}} _ {r}}
Și
- {\ displaystyle {\ vec {\ nabla}} _ {2} = {\ frac {m_ {2}} {m_ {1} + m_ {2}}} {\ vec {\ nabla}} _ {R} - {\ vec {\ nabla}} _ {r}}
Luând pătratele gradienților obținem laplacienii:
- {\ displaystyle \ nabla _ {1} ^ {2} = {\ frac {m_ {1} ^ {2}} {(m_ {1} + m_ {2}) ^ {2}}} \ nabla _ {R } ^ {2} + \ nabla _ {r} ^ {2} + {\ frac {2m_ {1}} {m_ {1} + m_ {2}}} {\ vec {\ nabla}} _ {R} \ cdot {\ vec {\ nabla}} _ {r}}
- {\ displaystyle \ nabla _ {2} ^ {2} = {\ frac {m_ {2} ^ {2}} {(m_ {1} + m_ {2}) ^ {2}}} \ nabla _ {R } ^ {2} + \ nabla _ {r} ^ {2} - {\ frac {2m_ {2}} {m_ {1} + m_ {2}}} {\ vec {\ nabla}} _ {R} \ cdot {\ vec {\ nabla}} _ {r}}
apoi Hamiltonianul este transformat:
- {\ displaystyle H = {\ frac {{\ vec {P}} ^ {2}} {2M}} + {\ frac {{\ vec {p}} ^ {2}} {2 \ mu}} + V (r)}
unde este
- {\ displaystyle {\ vec {P}} = - i \ hbar {\ vec {\ nabla}} _ {R}}
Și
- {\ displaystyle {\ vec {p}} = - i \ hbar {\ vec {\ nabla}} _ {r}}
și unde au fost introduse masa totală și masa redusă , respectiv:
- {\ displaystyle M = m_ {1} + m_ {2} \}
- {\ displaystyle \ mu = {\ frac {m_ {1} m_ {2}} {m_ {1} + m_ {2}}}}
Conform teoriei mecanicii cuantice, ecuația Schrödinger dependentă de timp:
- {\ displaystyle H \ Psi ({\ vec {r}}, t) = i \ hbar {\ frac {\ partial \ Psi ({\ vec {r}}, t)} {\ partial t}}}
admite soluții precum:
- {\ displaystyle \ Psi ({\ vec {r}}, t) = \ psi ({\ vec {r}}) \ cdot e ^ {- iEt / \ hbar}}
cu {\ displaystyle e ^ {- iEt / \ hbar}} soluție legată de evoluția temporală a funcției de undă e {\ displaystyle \ psi ({\ vec {r}})} Funcția de undă pentru ecuația Schrödinger independentă de timp:
- {\ displaystyle H \ psi ({\ vec {r}}) = E \ psi ({\ vec {r}})}
Acum, transformând Hamiltonianul, vedem că problema poate fi descompusă în două părți independente, prima parte este legată de soluția ecuației Schrödinger pentru mișcarea centrului de masă al sistemului cu două particule, cealaltă parte la mișcarea particulei {\ displaystyle \ mu} care se deplasează într-un câmp central în raport cu centrul de masă. Putem exprima funcția de undă totală ca:
- {\ displaystyle \ psi ({\ vec {r}} _ {1}, {\ vec {r}} _ {2}) = \ psi _ {R} ({\ vec {R}}) \ cdot \ psi _ {r} ({\ vec {r}})}
Soluția ecuației centrului de masă
Împărțind hamiltonienul în două părți la ecuația pentru centrul de masă obținem din ecuația relativă Schrödinger:
- {\ displaystyle H _ {\ vec {R}} \ psi _ {R} ({\ vec {R}}) = E _ {\ vec {R}} \ psi _ {R} ({\ vec {R} })}
unde este {\ displaystyle H _ {\ vec {R}} = - {\ frac {\ hbar ^ {2}} {2M}} \ nabla _ {\ vec {R}} ^ {2}} . Soluția generală a acestei ecuații este cea a unei particule libere :
- {\ displaystyle \ psi _ {R} ({\ vec {R}}) = Ae ^ {i {\ vec {k}} {\ vec {R}}} + Be ^ {- i {\ vec {k} } {\ vec {R}}}}
adică soluție de undă plană cu energie:
- {\ displaystyle E _ {\ vec {R}} = {\ frac {\ hbar ^ {2} k ^ {2}} {2M}}}
Și {\ displaystyle k} este vectorul de undă: centrul de masă se mișcă ca o particulă liberă, deoarece nu este supus niciunui potențial.
Soluția ecuației radiale pentru particulă {\ displaystyle \ mu} supus câmpului central
Ecuația Schrödinger pentru particula unică, adică pentru masa redusă devine:
- {\ displaystyle \ left (- {\ frac {\ hbar ^ {2}} {2 \ mu}} \ nabla ^ {2} + V (r) \ right) \ psi _ {r} (r) = E \ psi _ {r} (r)}
unde este {\ displaystyle E} reprezintă energia netă a centrului de masă, adică reprezintă energia internă a sistemului. Din moment ce potențialul {\ displaystyle V} este sferic, putem folosi coordonate sferice , noul operator hamiltonian devine:
- {\ displaystyle H = - {\ frac {\ hbar ^ {2}} {2 \ mu}} \ left [{\ frac {1} {r ^ {2}}} {\ frac {\ partial} {\ partial r}} \ left (r ^ {2} {\ frac {\ partial} {\ partial r}} \ right) + {\ frac {1} {r ^ {2} \ sin \ theta}} {\ frac { \ partial} {\ partial \ theta}} \ left (\ sin \ theta {\ frac {\ partial} {\ partial \ theta}} \ right) + {\ frac {1} {r ^ {2} \ sin ^ {2} \ theta}} {\ frac {\ partial ^ {2}} {\ partial \ varphi ^ {2}}} \ right] + V (r)}
Această ecuație poate fi gestionată cu ușurință dacă reconsiderăm operatorul momentului unghiular în coordonate sferice:
- {\ displaystyle {\ mathcal {L}} = - \ hbar ^ {2} \ left [{\ frac {1} {\ sin \ theta}} {\ frac {\ partial} {\ partial \ theta}} \ left (\ sin \ theta {\ frac {\ partial} {\ partial \ theta}} \ right) + {\ frac {1} {\ sin ^ {2} \ theta}} {\ frac {\ partial ^ {2} } {\ partial \ varphi ^ {2}}} \ right]}
Deci putem rescrie ecuația Schrödinger pentru particula unică ca:
- {\ displaystyle \ left (- {\ frac {\ hbar ^ {2}} {2 \ mu}} \ left [{\ frac {1} {r ^ {2}}} {\ frac {\ partial} {\ partial r}} \ left (r ^ {2} {\ frac {\ partial} {\ partial r}} \ right) - {\ frac {\ mathcal {L}} {\ hbar ^ {2} r ^ {2 }}} \ right] + V (r) \ right) \ psi _ {r} (r, \ theta, \ varphi) = E \ psi _ {r} (r, \ theta, \ varphi)}
Soluția acestei ecuații poate fi luată în considerare în continuare:
- {\ displaystyle \ psi _ {r} (r, \ theta, \ varphi) = f (r) \ cdot f (\ theta) \ cdot f (\ varphi)}
de fapt, pe de o parte, soluția părții unghiulare din fizica matematică în termeni de armonici sferice este bine cunoscută:
- {\ displaystyle Y_ {lm} (\ theta, \ varphi) = \ Theta _ {lm} (\ theta) \ cdot \ Phi (\ varphi)}
după cum se știe, armonicele sferice sunt funcții proprii simultane ale operatorului momentului unghiular {\ displaystyle {\ mathcal {L}}} și proiecția sa de-a lungul axei z: {\ displaystyle {\ mathcal {L}} _ {z}} ; indicele {\ displaystyle l} și {\ displaystyle m} în schimb, acestea sunt numerele cuantice unghiulare și magnetice.
Soluția completă pentru particule {\ displaystyle \ mu} care orbitează centrul de masă (indicele {\ displaystyle r} indică funcția sa de undă) este atunci:
- {\ displaystyle \ psi _ {r} (r, \ theta, \ varphi) = R_ {E, l} (r) \ cdot Y_ {lm} (\ theta, \ varphi)}
Trebuie doar să găsim soluția rămasă a părții radiale a ecuației. De fapt, scrierea ecuației radiale de Schrödinger:
- {\ displaystyle \ left (- {\ frac {1} {2 \ mu}} \ left [{\ frac {\ hbar ^ {2}} {r ^ {2}}} {\ frac {d} {dr} } \ left (r ^ {2} {\ frac {d} {dr}} \ right) - {\ frac {l (l + 1) \ hbar ^ {2}} {r ^ {2}}} \ right ] + V (r) \ right) R_ {E, l} (r) = E \ cdot R_ {E, l} (r)}
unde este {\ displaystyle l (l + 1) \ hbar ^ {2}} sunt valorile proprii ale operatorului momentului unghiular {\ displaystyle {\ mathcal {L}}} , noi vedem asta {\ displaystyle R_ {E, l}} depinde și de {\ displaystyle l} dar nu din {\ displaystyle m} , de fapt operatorul nu apare {\ displaystyle {\ mathcal {L}} _ {z}} .
Să rescriem ecuația radială:
- {\ displaystyle \ left [- {\ frac {\ hbar ^ {2}} {2 \ mu}} {\ frac {1} {r ^ {2}}} {\ frac {d} {dr}} r ^ {2} {\ frac {d} {dr}} + V_ {eff} \ right] R_ {E, l} (r) = E \ cdot R_ {E, l} (r)}
unde cu
- {\ displaystyle V_ {eff} = {\ frac {l (l + 1) \ hbar ^ {2}} {2 \ mu r ^ {2}}} + V (r)}
am indicat potențialul efectiv numit energie potențială centrifugă (este respingător), deci vedem că ecuația radială este cea unidimensională a particulei (amintiți-vă că particula în cauză este masa redusă) care se deplasează într-un potențial efectiv. Rescriem în mod explicit ecuația radială ca:
- {\ displaystyle - {\ frac {\ hbar ^ {2}} {2 \ mu}} \ left [{\ frac {1} {r}} {\ frac {d} {dr}} \ left (r {\ frac {d} {dr}} \ right) + {\ frac {1} {r}} {\ frac {d} {dr}} - {\ frac {l (l + 1)} {r ^ {2} }} \ right] R_ {E, l} + \ left (V (r) -E \ right) R_ {E, l} (r) = 0}
Realizăm înlocuirea:
- {\ displaystyle R_ {E, l} (r) = {\ frac {\ chi (r)} {r}}}
asa de:
- {\ displaystyle {\ frac {1} {r}} {\ frac {dR_ {E, l}} {dr}} = {\ frac {1} {r}} \ left ({\ frac {1} {r }} {\ frac {d \ chi} {dr}} - {\ frac {\ chi} {r ^ {2}}} \ right) = {\ frac {1} {r ^ {2}}} \ left ({\ frac {d \ chi} {dr}} - {\ frac {\ chi} {r}} \ right)}
- {\ displaystyle {\ frac {1} {r}} {\ frac {d} {dr}} \ left (r {\ frac {dR_ {E_ {l}}} {dr}} \ right) = {\ frac {1} {r}} {\ frac {d} {dr}} \ left ({\ frac {d \ chi} {dr}} - {\ frac {\ chi} {r}} \ right) = {\ frac {1} {r}} {\ frac {d ^ {2} \ chi} {dr ^ {2}}} - {\ frac {1} {r ^ {2}}} {\ frac {d \ chi } {dr}} + {\ frac {1} {r ^ {3}}} \ chi = {\ frac {1} {r}} \ left ({\ frac {d ^ {2} \ chi} {dr ^ {2}}} - {\ frac {1} {r}} {\ frac {d \ chi} {dr}} + {\ frac {\ chi} {r ^ {2}}} \ right)}
Ecuația Schrödinger pentru {\ displaystyle \ chi} devine
- {\ displaystyle - {\ frac {\ hbar ^ {2}} {2 \ mu}} {\ frac {d ^ {2} \ chi} {dr ^ {2}}} + \ left [{\ frac {\ hbar ^ {2} l (l + 1)} {2 \ mu r ^ {2}}} + V (r) -E \ right] \ chi = 0}
După cum se poate vedea din această ecuație, nu conține {\ displaystyle m} , de fapt este degenerat față de {\ displaystyle {\ mathcal {L}} _ {z}} întrucât într-un câmp central simetric toate direcțiile sunt egale. Funcția radială conține energia și numărul cuantic {\ displaystyle l} , în timp ce partea unghiulară dată de armonicele sferice conține {\ displaystyle l} și {\ displaystyle m} , deci mișcarea este complet definită de cele trei numere {\ displaystyle n, l, m} , care permit calcularea energiei, a momentului unghiular orbital și a componentei sale față de orice axă, care sunt trei operatori care formează un set complet de mărimi fizice pentru mișcare într-un câmp central simetric.
Soluția radială poate fi dată furnizând în mod explicit forma potențialului {\ displaystyle V (r)} cu condiția ca potențialul să nu fie singular la origine:
- {\ displaystyle \ lim _ {r \ to 0} r ^ {2} V (r) = 0} .
Cu toate acestea, putem vedea comportamentul asimptotic al funcției {\ displaystyle \ chi} . Pentru {\ displaystyle r \ to 0} funcția trebuie să nu fie singulară, astfel încât condiția de mai sus trebuie să rămână în continuare:
- {\ displaystyle \ chi (0) = 0}
și ecuația pentru {\ displaystyle \ chi} devine:
- {\ displaystyle {\ frac {d ^ {2} \ chi} {dr ^ {2}}} - {\ frac {l (l + 1)} {r ^ {2}}} \ chi \ simeq 0}
a cărei soluție poate fi găsită pentru:
- {\ displaystyle \ chi (r) \ sim r ^ {s}}
care, înlocuit în precedent, ne oferă ecuația:
- {\ displaystyle s (s-1) -l (l + 1) = 0}
adică o soluție {\ displaystyle s = -l} de la care {\ displaystyle r ^ {- l}} care are o origine neregulată și o soluție regulată {\ displaystyle s = (l + 1)} de la care {\ displaystyle r ^ {l + 1}} care satisface condiția de regularitate {\ displaystyle \ chi (0) = 0} .
Pentru {\ displaystyle r \ to \ infty} ecuația pentru {\ displaystyle \ chi} devine:
- {\ displaystyle {\ frac {d ^ {2} \ chi} {dr ^ {2}}} + {\ frac {2 \ mu E} {\ hbar ^ {2}}} \ chi \ simeq 0}
și folosind condiția de normalizare a funcției de undă:
- {\ displaystyle \ int _ {0} ^ {\ infty} dr \, r ^ {2} | R_ {E, l} | ^ {2} = 1 \, \, \, \ Rightarrow \, \, \ int _ {0} ^ {\ infty} dr \, | \ chi (r) | ^ {2} = 1}
rezultă că funcția de undă trebuie în orice caz să tindă la zero pentru {\ displaystyle r \ to \ infty} . Pentru stări legate, adică pentru {\ displaystyle E <0} avem:
- {\ displaystyle {\ frac {2 \ mu E} {\ hbar ^ {2}}} = - k ^ {2}}
iar comportamentul asimptotic este:
- {\ displaystyle \ chi (r \ to \ infty) \ sim e ^ {- kr}}
În general, putem anticipa doar că soluția ecuației radiale este proporțională cu polinoamele Laguerre :
- {\ displaystyle R_ {E, l} (r) = R_ {nl} (r) \ propto N_ {n, l} L_ {n + l} ^ {2l + 1} (\ cdot)}
unde este {\ displaystyle n} este numărul cuantic principal , {\ displaystyle L_ {n + l} ^ {2l + 1} (\ cdot)} sunt polinoamele Laguerre și {\ displaystyle N_ {nl}} este o constantă de normalizare. Un exemplu tipic este cel al câmpului Coulomb care caracterizează atomul de hidrogen și particula liberă tridimensională.
Considerații suplimentare
După cum putem vedea, valorile proprii ale energiei depind doar de numărul cuantic principal {\ displaystyle n} , prin urmare energia este degenerată în raport cu {\ displaystyle l} și din {\ displaystyle m} : pentru fiecare {\ displaystyle n = 1,2, \ dots, \ infty} , valoarea proprie a momentului unghiular orbital poate lua valorile {\ displaystyle l = 0,1, \ dots, n-1} și, pentru fiecare valoare de {\ displaystyle l} , valoarea proprie a proiecției momentului unghiular poate lua valorile {\ displaystyle m = -l, -l + 1, \ dots, + l} . Degenerarea totală este produsul degenerării în {\ displaystyle l} adică poate lua n valori diferite plus degenerare în {\ displaystyle m} care la rândul său poate lua n valori și, prin urmare, degenerarea totală {\ displaystyle n ^ {2}} . În conformitate cu interpretarea probabilistică a mecanicii cuantice,
- {\ displaystyle \ int dr \, | \ psi _ {r} (r) | ^ {2}}
reprezintă probabilitatea ca particula {\ displaystyle \ mu} orbita centrului de masă al sistemului este undeva în spațiu. Probabilitatea totală a funcției de undă este:
- {\ displaystyle \ int d {\ vec {r}} \, | \ Psi ({\ vec {r}}) | ^ {2}}
Bibliografie
- ( EN ) BH Bransden și Charles J. Joachain, Physics of atoms and molecules , Boston, Addison-Wesley, 2005, ISBN 978-05-82-35692-4 .
Elemente conexe
Portal cuantic : Accesați intrările Wikipedia care se ocupă de cuantică |