În matematică, ecuația Hamilton–Jacobi descrie o condiție necesară de extrem geometric în generalizarea problemelor calculului variațional. În fizică, ea este o reformulare a mecanicii clasice și ca atare echivalentă cu alte formulări, precum mecanica newtoniană, mecanica lagrangiană și mecanica hamiltoniană. În particular, ecuația Hamilton–Jacobi este folositoare la identificarea mărimilor care se conservă într-un sistem mecanic, ceea ce este posibil chiar și în cazul în care problema mecanică nu poate fi rezolvată complet.
De asemenea, ecuația Hamilton–Jacobi este singura formulare din mecanică în care mișcarea unui sistem de particule este descrisă într-un formalism asemănător cu propagarea unei unde. În acest sens, a fost atins un obiectiv al fizicii teoretice (datând din secolul 18 de la Johann Bernoulli): găsirea unei analogii între propagarea luminii și mișcarea unei particule. Ecuația de undă pentru sistemele mecanice este similară, dar nu identică, cu ecuația lui Schrödinger; din acest motiv, ecuația Hamilton–Jacobi înlesnește abordarea mecanicii cuantice, pornind de la mecanica clasică.
Ecuația Hamilton–Jacobi este o ecuație cu derivate parțiale neliniare de ordinul întâi pentru o funcție S ( q 1 , … , q N ; t ) {\displaystyle S(q_{1},\dots ,q_{N};t)} , numită funcția principală a lui Hamilton:
Această ecuație derivă din mecanica Hamiltoniană prin tratarea funcției S {\displaystyle S} ca funcție generatoare pentru o tranformare canonică a Hamiltonianului H ( q 1 , … , q N ; p 1 , … , p N ; t ) {\displaystyle H(q_{1},\dots ,q_{N};p_{1},\dots ,p_{N};t)} . Impulsul generalizat corespunzător primei derivate a funcției S {\displaystyle S} , în funcție de coordonatele generalizate este:
Schimarea în acțiune de la o traiectorie la alta apropiată este dată de ecuația:
Deoarece traiectoria mișcării actuale satisface ecuația Euler-Lagrange, variația δ S {\displaystyle \delta S} este zero. În primul termen al ecuației punem δ q k ( t 1 ) = 0 {\displaystyle \delta q_{k}(t_{1})=0} , iar valoarea δ q k ( t 2 ) {\displaystyle \delta q_{k}(t_{2})} o notăm simplu prin δ q k {\displaystyle \delta q_{k}} . Înlocuind ∂ L / ∂ q ˙ k {\displaystyle \partial L/\partial {\dot {q}}_{k}} prin p k {\displaystyle p_{k}} , obținem în final:
Pornind de la această relație urmează că, derivatele parțiale ale acțiunii în funcție de coordonate sunt egale cu impulsurile generalizate corespunzătoare.
Similar, coordonatele generalizate pot fi obținute prin derivarea acțiunii S {\displaystyle S} , în funcție de impulsurile generalizate.
Prin inversarea acestor ecuatii, se poate determina evoluția unui sistem mecanic, adică, putem determina coordonatele generalizate ca funcții de timp. Pozițiile și vitezele inițiale apar drept constante de integrare ale soluției S {\displaystyle S} , corespunzând mărimilor care se conservă în timpul evoluției sistemului mecanic, precum energia, momentul unghiular, sau vectorul Laplace-Runge-Lenz.
Ecuația Hamilton-Jacobi este o ecuație cu derivate parțiale de ordinul întâi a acțiunii S {\displaystyle S} în funcție de N {\displaystyle N} coordonate generalizate q 1 , … , q N {\displaystyle q_{1},\dots ,q_{N}} și timpul t. Impulsurile generalizate nu apar în S {\displaystyle S} , ci numai în derivatele lui S {\displaystyle S} . Remarcabil este faptul că, funcția S {\displaystyle S} este egală cu acțiunea clasică.
Pentru comparație, nici în ecuația de mișcare echivalentă Euler-Lagrange din mecanica Lagrangiană nu apar impulsurile generalizate. Totuși, ele formează un sistem de N ecuații cu derivate de ordinul al doilea în funcție de coordonatele generalizate și timp. O altă comparație este aceea cu ecuația de mișcare a lui Hamilton, care este de fapt un sistem de 2N ecuații de ordinul întâi în funcție de coordonatele generalizate, impulsurile generalizate p 1 , … , p N {\displaystyle p_{1},\dots ,p_{N}} și timp.
Deoarece ecuația Hamilton-Jacobi este o expresie echivalentă a problemelor de miminizare a integralelor, precum principiul lui Hamilton, ea poate fi folositoare și în alte probleme de calcul variațional, sau mai general, în alte ramuri ale matematicii sau fizicii, precum sistemele dinamice, geometria simplectică sau haosului cuantic. De exemplu, ecuația Hamilton-Jacobi este folositoare la determinarea geodezicelor pe o mulțime Riemanniană, care este o problemă importantă variațională din geometria Riemanniană.
Pentru a fi conciși, folosim variabile îngroșate, precum q {\displaystyle \mathbf {q} } , pentru a reprezenta cele N {\displaystyle N} coordonate generalizate:
care nu se transformă neapărat printr-o rotație ca un vector. Produsul scalar este definit aici drept suma produselor componentelor corespunzătoare, adică:
Orice transformare canonică implică o funcție generatoare G 2 ( q , P , t ) {\displaystyle G_{2}(\mathbf {q} ,\mathbf {P} ,t)} , care conduce la relațiile:
Pentru a deriva ecuația Hamilton-Jacobi, alegem o funcție generatoare S ( q , P , t ) {\displaystyle S(\mathbf {q} ,\mathbf {P} ,t)} care face noul Hamiltonian K {\displaystyle K} egal cu zero. Astfel că, toate derivatele sale sunt de asemenea zero, iar Hamiltonianul devine trivial:
adică, noile coordonate și impulsuri generalizate sunt constante. Noul impuls generalizat P {\displaystyle \mathbf {P} } este notat prin α 1 , α 2 , … , α N − 1 , α N {\displaystyle \alpha _{1},\alpha _{2},\ldots ,\alpha _{N-1},\alpha _{N}} , adică, P m = α m {\displaystyle P_{m}=\alpha _{m}} .
Atunci, ecuația Hamilton-Jacobi rezultă din ecuația transformată K {\displaystyle K} :
care este echivalentă cu ecuația:
deoarece p = ∂ S / ∂ q {\displaystyle \mathbf {p} =\partial S/\partial \mathbf {q} } .
Noile coordonate generalizate Q {\displaystyle \mathbf {Q} } sunt de asemenea constante, notate cu β 1 , β 2 , … , β N − 1 , β N {\displaystyle \beta _{1},\beta _{2},\ldots ,\beta _{N-1},\beta _{N}} . Odată ce le-am rezolvat pentru S ( q , α , t ) {\displaystyle S(\mathbf {q} ,{\boldsymbol {\alpha }},t)} , se obțin ecuațiile:
sau, pentru claritate, scrise pentru componentele lui Q {\displaystyle \mathbf {Q} } :
Aceste N {\displaystyle N} ecuații pot fi inversate pentru a găsi coordonatele generalizate originale q {\displaystyle \mathbf {q} } ca funcții de constantele α {\displaystyle {\boldsymbol {\alpha }}} și β {\displaystyle {\boldsymbol {\beta }}} , astfel putând rezolva problema originală.
Ecuația Hamilton-Jacobi este foarte folositoare când poate fi rezolvată via variabilelor separabile aditive, care identifică direct constantele de mișcare. De exemplu, timpul t poate fi separat dacă Hamiltonianul nu depinde explicit de timp. Î acest caz, derivata funcție de timp ∂ S ∂ t {\displaystyle {\frac {\partial S}{\partial t}}} trebuie să fie o constantă, notată cu − E {\displaystyle -E} , dând soluția:
în care, funcția independentă de timp W ( q ) {\displaystyle W(\mathbf {q} )} este numită uneori și funcția caracteristică a lui Hamilton. Atunci, ecuația Hamilton–Jacobi redusă poate fi scrisă astfel:
Pentru a ilustra separabilitatea și pentru alte variabile presupunem că, oricare coordonată generalizată q k {\displaystyle q_{k}} și derivata ei ∂ S ∂ q k {\displaystyle {\frac {\partial S}{\partial q_{k}}}} apar împreună în Hamiltonian ca o singură funcție ψ ( q k , ∂ S ∂ q k ) {\displaystyle \psi \left(q_{k},{\frac {\partial S}{\partial q_{k}}}\right)} , iar H se scrie:
În acest caz, funcția S {\displaystyle S} poate fi despărțită în două funcții, una care depinde numai de q k {\displaystyle q_{k}} și alta care depinde numai de coordonatele generalizate rămase:
Substituția acestor formule în ecuația Hamiton-Jacobi arată că funcția ψ {\displaystyle \psi } trebuie să fie o constantă, aici notată cu Γ k {\displaystyle \Gamma _{k}} , obținând o ecuație diferențială ordinară de ordinul întâi pentru S k ( q k ) {\displaystyle S_{k}(q_{k})} :
În anumite cazuri, funcția S {\displaystyle S} poate fi separată complet în N {\displaystyle N} funcții S m ( q m ) {\displaystyle S_{m}(q_{m})} , obținând:
În astfel de cazuri, problema se rezolvă prin N {\displaystyle N} ecuații diferențiale ordinare.
Separabilitatea funcției S {\displaystyle S} depinde de Hamiltonian și de modul în care sunt alese coordonatele generalizate. Pentru coordonate ortogonale și Hamiltonian care nu depinde de timp și este pătratic pentru impulsurile generalizate, S {\displaystyle S} este complet separabilă dacă energia potențială este separabilă aditiv pentru fiecare coordonată, caz în care, termenul energiei potențiale pentru fiecare coordonată este multiplicat corespunzător printr-un factor dependent de coordonate în termenul impulsului Hamiltonianului (condiția Staeckel). Pentru a ilustra acest lucru, în secțiunea următoare sunt date câteva exemple în coordonate ortogonale.
Hamiltonianul în coordonate sferice poate fi scrie sub forma:
Ecuația Hamilton-Jacobi este complet separabilă în aceste coordonate, demonstrându-se că U {\displaystyle U} are forma analoagă cu:
în care U r ( r ) {\displaystyle U_{r}(r)} , U θ ( θ ) {\displaystyle U_{\theta }(\theta )} și U ϕ ( ϕ ) {\displaystyle U_{\phi }(\phi )} sunt funcții arbitrare. Substituind soluția complet separată S = S r ( r ) + S θ ( θ ) + S ϕ ( ϕ ) − E t {\displaystyle S=S_{r}(r)+S_{\theta }(\theta )+S_{\phi }(\phi )-Et} în ecuația Hamilton-Jacobi, obținem:
Acestă ecuație poate fi rezolvată prin integrări succesive de ecuații diferențiale ordinare, începând cu ecuația ϕ {\displaystyle \phi } :
unde Γ ϕ {\displaystyle \Gamma _{\phi }} este o constantă de mișcare care elimină dependența de ϕ {\displaystyle \phi } din ecuația Hamilton–Jacobi:
Următoarea ecuație diferențială ordinară implică coordonata generalizată θ {\displaystyle \theta } :
în care Γ θ {\displaystyle \Gamma _{\theta }} este o altă constantă de mișcare care elimină dependența de θ {\displaystyle \theta } și reduce ecuația Hamilton-Jacobi la ecuația finală diferențială:
a cărei integrare completează soluția pentru S {\displaystyle S} .
Hamiltonianul în coordonate cilindrice eliptice poate fi scris astfel:
unde focarul elipsei este localizat în ± a {\displaystyle \pm a} , pe axa x {\displaystyle x} . Ecuația Hamilton-Jacobi este complet separabilă în aceste coordonate, demonstrându-se că U {\displaystyle U} are o forma analoagă cu:
în care U μ ( μ ) {\displaystyle U_{\mu }(\mu )} , U ν ( ν ) {\displaystyle U_{\nu }(\nu )} și U z ( z ) {\displaystyle U_{z}(z)} sunt funcții arbitrare. Substituind soluția complet separabilă S = S μ ( μ ) + S ν ( ν ) + S z ( z ) − E t {\displaystyle S=S_{\mu }(\mu )+S_{\nu }(\nu )+S_{z}(z)-Et} în ecuația Hamilton-Jacobi, obținem:
Separând prima ecuație diferențială ordinară (funcție numai de z):
obținem ecuația Hamilton–Jacobi redusă (după ce multiplicăm cu 2m și rearanjăm ecuația):
care poate fi separată în două ecuații diferențiale ordinare:
care rezolvate, conduc la o soluționare completă a lui S {\displaystyle S} .
Hamiltonianul în coordonate cilindrice parabolice poate fi scris sub forma:
Ecuația Hamilton-Jacobi este complet separabilă în aceste coordonate, demonstrându-se că U {\displaystyle U} are o forma analoagă cu:
în care U σ ( σ ) {\displaystyle U_{\sigma }(\sigma )} , U τ ( τ ) {\displaystyle U_{\tau }(\tau )} și U z ( z ) {\displaystyle U_{z}(z)} sunt funcții arbitrare. Substituind soluția complet separabilă S = S σ ( σ ) + S τ ( τ ) + S z ( z ) − E t {\displaystyle S=S_{\sigma }(\sigma )+S_{\tau }(\tau )+S_{z}(z)-Et} în ecuația Hamilton-Jacobi, obținem:
Izosuprafața funcției S ( q ; t ) {\displaystyle S(\mathbf {q} ;t)} poate fi determinată la oricare timp t. Mișcarea unei izosuprafețe S {\displaystyle S} ca funcție de timp este definită prin mișcarea unei particule dintr-un punct q {\displaystyle \mathbf {q} } al izosuprafeței. Mișcarea unei astfel de izosuprafețe poate fi gândită ca o undă care se mișcă prin spațiul q {\displaystyle \mathbf {q} } , cu toate că nu se subordonează exact ecuației de undă. Pentru a arăta acest lucru considerăm că S {\displaystyle S} reprezintă faza unei unde
în care ℏ {\displaystyle \hbar } este o constantă introdusă pentru a face exponențiala adimensională. Schimbarea în amplitudine a undei poate fi reprezentată printr-un număr complex S {\displaystyle S} . Putem rescrie ecuația Hamilton-Jacobi astfel:
care este o variantă neliniară a ecuației Schrödinger.
Invers, plecând de la ecuația lui Schrödinger și de la Ansatzul lui ψ {\displaystyle \psi } , obținem:
Limita clasică ( ℏ → 0 {\displaystyle \hbar \rightarrow 0} ) a ecuației Schrödinger de mai sus devine identică cu următoarea variantă a ecuației Hamilton-Jacobi:
în care g i k {\displaystyle g^{ik}} sunt componentele contravariante ale tensorului metric, m este masa de repaus a particulei, iar c este viteza luminii.