משוואת המילטון-יעקובי

במכניקה אנליטית, **משוואת המילטון-יעקובי** היא ניסוח חלופי של המכניקה הקלאסית. המשוואה שקולה לניסוחים אחרים כגון חוקי התנועה של ניוטון, המכניקה הלגראנז'ית והמכניקה ההמילטונית. המשוואה נקראת על שמם של ויליאם רואן המילטון וקרל גוסטב יעקובי שפתחו אותה בראשית המאה ה-19.

משוואת המילטון־יעקובי שימושית במיוחד במציאת גדלים נשמרים במערכות בהן מציאת הפתרון המלא של המערכת איננה אפשרית.

## רקע
במהלך סוף המאה ה-17 ותחילת המאה ה-18 נסוב ויכוח בין מדענים שדגלו בכך שהאור הוא חלקיק הנע במסלולים מוגדרים - הקרניים המופיעות באופטיקה גאומטרית, ובין אלו שטענו שהאור הוא גל, וכל נקודה אליה מגיע חזית הגל היא מקור של גל כדורי - עקרון הויגנס. ב-1818 הראה פרנל שניתן להסביר גם את האופטיקה הגאומטרית וגם את תופעות העקיפה באמצעות תיאוריה גלית של האור.

האופטיקה הגאומטרית והתיאור החלקיקי של האור יכול להתבצע בהסתמך על עקרון פרמה. לפי עיקרון זה, קרן אור גאומטרית תטייל בין נקודה A לנקודה B לאורך המסלול בו זמן התנועה של הקרן יהיה מינימלי. מאחר שמהירות האור בתווך היא v=c/n כאשר n הוא מקדם השבירה של החומר הרי שהקרן תנוע לאורך המסלול בו \int_{A}^B n\ ds מינימלי. מצד שני, לפי התיאור הגלי של האור, האור מקיים משוואת גלים מהצורה \nabla^2 \Psi = \frac{n^2}{c^2} \frac{\partial \Psi}{\partial t}, גל אור מונוכרומטי יכול להיכתב בצורה \Psi = a(x,y,z) \exp[2\pi i( f t - \varphi(x,y,z))] כאשר a האמפליטודה של הגל,\varphi הפאזה ו-f התדירות. במקרה כזה, משוואת הגלים הופכת להיות: \nabla ^2\Psi = \frac{4\pi^2 f^2}{c^2} \Psi. אם אורך הגל \lambda= \frac{nf}{c} , קטן אז הפאזה גדולה, a משתנה לאט ביחס להשתנות הפאזה, וניתן להזניח את הנגזרות של a. בנוסף ניתן להזניח את הנגזרת השנייה של הפאזה ביחס לנגזרת של הפאזה בריבוע. מתקבלת המשוואה הבאה עבור הפאזה:

|\nabla \varphi | ^2 = \frac{f^2n^2}{c^2}=\frac{1}{\lambda^2(x,y,z)}

משוואה זו ידועה בתור משוואת האייקונל. קווים עליהם \varphi=C נקראים קווים שווי פאזה, או חזית הגל. האופטיקה הגאומטרית מתקבלת מהאופטיקה הגלית בקירוב זה, וקרן אור גאומטרית תנוע במאונך לקווים שווי פאזה.

ב-1833 הצביע המילטון על העובדה שקיימת הקבלה בין תנועה של מערכת במכניקה קלאסית לתנועה של קרן אור גאומטרית. בשני המקרים קיים פונקציונל שמקבל ערך מינימום לאורך מסלול, במקרה האופטי הפונקציונל הוא זמן התנועה לאורך המסלול ובמקרה המכני הפונקציונל הוא הפעולה. משוואת המילטון־יעקובי מחפשת למצוא את המקבילה המכנית לתיאור הגלי. היא עושה זאת באמצעות חיפוש משטחים שווי "פעולה" שמהווים את חזית הגל המכנית. מרגע שקווים אלו נמצאו, קיים קשר ביניהם לבין המסלול בו המערכת תמשיך להתקדם. ההקבלה הזו מצביעה על מאפיין גלִי של המערכת המכנית, ובכך מהווה בסיס נוח לפיתוח מכניקת הקוונטים.

## ניסוח מתמטי
בהינתן סט של קואורדינטות קנוניות \{q_i,p_i\}_{i=1}^n והמילטוניאן H(\mathbf{q};\mathbf{p};t) עבור מערכת פיזיקלית כך שמתקיימים משוואות המילטון \dot{q_i}=\frac{\partial H}{\partial p_i},\dot{p_i}=-\frac{\partial H}{\partial q_i} ה**פונקציה המנהלת של המילטון <math>S(\mathbf{q};\alpha_1...,\alpha_n;t)</math>** (Hamilton's principal function, פירוט בהמשך) הפותרת את **משוואת המילטון-יעקובי**:

\frac{\partial S}{\partial t} + H\left(\mathbf{q};\frac{\partial S}{\partial q_1},...,\frac{\partial S}{\partial q_n};t\right)=0

מקיימת p_i = \frac{\partial S}{\partial q_i} ו- \frac{d\beta_i}{dt} =\frac{d}{dt} \frac{\partial S}{\partial \alpha_i} = 0. כלומר הגדלים \beta_i = \frac{\partial S}{\partial \alpha_i} הם גדלים נשמרים. גם הגדלים \alpha_i הם קבועים במהלך התנועה. אם ידועים הערכים הללו בזמן t=t_0 ניתן להפוך את סט המשוואות f_i(\mathbf{q};\alpha_1,...\alpha_n;t) = \frac{\partial S}{\partial \alpha_i} = \beta_i ולקבל q_i = q_i(\boldsymbol\alpha,\boldsymbol\beta,t).

אם ההמילטוניאן לא תלוי ישירות בזמן H(\mathbf{q};\mathbf{p}) ה**פונקציה המאפיינת של המילטון <math>W(\mathbf{q};\gamma_1,...,\gamma_n)</math>** (Hamilton's characteristic function) הפותרת את **משוואת המילטון-יעקובי הלא תלויה בזמן**:

f(\gamma_1, ..., \gamma_n) = H\left(\mathbf{q};\frac{\partial W}{\partial q_1},...,\frac{\partial W}{\partial q_n}\right)

מקיימת p_i = \frac{\partial W}{\partial q_i} וכן \frac{d}{dt} \frac{\partial W}{\partial \gamma_i} = \frac{\partial f}{\partial \gamma_i} = v_i \implies \frac{\partial W}{\partial \gamma_i } = v_i t + \beta_i . הפונקציה f ניתנת לבחירה כדי לפשט את הביטויים בפתרון, ובחירות שונות לא ישנו את פתרון הבעיה. אם נעשה שימוש בפונקציה f(\alpha_1,...,\alpha_n)= \alpha_1, מתקיים \frac{\partial W}{\partial \alpha_i} = \left\{ \begin{matrix} t+\beta_1 & i=1 \\ \beta_i & i\neq 1 \end{matrix} \right.  ומתקיים הקשר הבא בין הפונקציה המנהלת לפונקציה המאפיינת S(\mathbf{q};\alpha_1,...,\alpha_n;t) = W(\mathbf{q};\alpha_2,...,\alpha_n) - \alpha_1 t, אם הלגרנז'יאן המקורי יכול להיכתב בצורה L(\mathbf{q},\mathbf{\dot{q}})= L_0(\mathbf{q}) + \sum_{i,j=1}^n L_2^{ij}(\mathbf{q}) \dot{q_i}\dot{q_j} אז \alpha_1=E היא האנרגיה של המערכת .

## פיתוח באמצעות פונקציה יוצרת
טרנספורמציות קנוניות מאפשרות לעבור ממערכת אחת של קואורדינטות קנוניות (\mathbf{q},\mathbf{p}) המקיימת את משוואת המילטון עבור ההמילטוניאן H(\mathbf{q},\mathbf{p},t), לסט אחר (\mathbf{Q},\mathbf{P}) המקיים את משוואת המילטון עבור המילטוניאן אחר K(\mathbf{Q},\mathbf{P},t). ניתן לפתח טרנספורמציה קנונית באמצעות שימוש בפונקציה יוצרת מסוג שני G(\mathbf{q},\mathbf{P},t) כך שיתקיים p_i = \frac{\partial G}{\partial q_i}; Q_i = \frac{\partial G}{\partial P_i}; K = H + \frac{\partial G}{\partial t}. משוואת המילטון־יעקובי מגיעה מטרנספורמציה קנונית שמטרתה ליצור מצב בו ההמילטוניאן החדש שווה זהותית ל-0 K(\mathbf{Q},\mathbf{P},t)=0. במצב כזה, מתקיים \dot Q_i = \frac{\partial K}{\partial P_i} = 0;\dot P_i = -\frac{\partial K}{\partial Q_i} = 0 כלומר הקואורדינטות החדשות (\mathbf{Q},\mathbf{P}) הן קבועות במהלך התנועה, ותנועת המערכת במרחב הפאזה כפי שהוא מתואר על ידי הקואורדינטות האלו הוא נקודה (בניגוד למסלול חד־ממדי המתאר את התנועה במרחב הפאזה כפי שהוא נפרש על ידי קואורדינטות אחרות). בהתאם למשוואות של הפונקציה היוצרת נקבל  0= K = H(\mathbf{q},\mathbf{p},t) + \frac{\partial G(\mathbf{q}, \mathbf{P},t)}{\partial t} אבל לפי המשוואות הללו מתקיים גםp_i = \frac{\partial G}{\partial q_i} ולכן  0= H\left(\mathbf{q};\frac{\partial G}{\partial q_1}...,\frac{\partial G}{\partial q_n};t\right) + \frac{\partial G(\mathbf{q}, \mathbf{P},t)}{\partial t} שהיא משוואת המילטון־יעקובי עבור הפונקציה היוצרת G. התוצאה לפיה \frac{d}{dt} \frac{\partial G}{\partial P_i}=0 מתקבלת מהמשוואות Q_i = \frac{\partial G}{\partial P_i};\dot{Q}_i = 0.

## פעולה ופונקציות המילטון
הנגזרת המלאה לפי הזמן של הפונקציה המנהלת של המילטון מקיימת:

\frac{d S(\mathbf{q},\boldsymbol{\alpha},t)}{dt} = \frac{\partial S}{\partial t} + \sum_{i=1}^n \frac{\partial S}{\partial q_i} \dot{q}_i + \frac{\partial S}{\partial \alpha_i} \dot{\alpha}_i =\frac{\partial S}{\partial t} + \sum_{i=1}^n \frac{\partial S}{\partial q_i} \dot{q}_i = -H + \sum_{i=1}^n p_i \dot{q}_i = L 

כלומר מתקיים S(\mathbf{q};\boldsymbol{\alpha};t) = \int L\ dt + c כאשר L הוא הלגראנז'יאן. משוואה זו מזכירה את הגדרת הפעולה של המערכת S_{action}[\mathbf{q}(t)] = \int_{t_0}^{t_1} L(\mathbf{q}(t);\dot\mathbf{q}(t);t)\ dt המשמשת בעקרון הפעולה המינימלית כדי לפתח את משוואות אוילר־לגראנז'. על אף שהפונקציה המנהלת והפעולה מזכירות אחת את השנייה, הן אובייקטים מתמטיים שונים - הפונקציה המנהלת היא פונקציה בעלת תחום במרחב הפאזה והזמן, לעומתה הפעולה היא פונקציונל בעלת תחום במרחב המסלולים הפיזיקליים במרחב הקונפיגורציות.

הקשר בין הפונקציה המנהלת לפעולה נובע מכך שאם נסתכל במערכת, שהתחילה בזמן t=t_0 בנקודה במרחב הקונפיגורציות \mathbf{q}=\mathbf{q}_0 והמקיימת \dot\mathbf{q}(t_0)= \mathbf{v}_0 המסלול \mathbf{q}(t)  המייצר נקודת אקסטרמום לפעולה נקבע באופן יחיד מתנאי ההתחלה. בנוסף, אם הפרש הזמנים בין שתי נקודות קטן מספיק, קיים רק מתנאי התחלה אחד (ולפיכך רק מסלול אחד) שיכול להביא את המערכת מהנקודה \mathbf{q}(t_0) = \mathbf{q_0} לנקודה \mathbf{q}(t) = \mathbf{q} כך שהפעולה תקבל ערך אקסטרימום. לפיכך, בהינתן נקודת התחלה במרחב הקונפיגורציות, ניתן להגדיר פונקציה של הקוארדינטות בסיום המסלול \mathbf{q} כך ש-S(\mathbf{q},t;\mathbf{q_0},t_0) = S_{action}[\mathbf{q}(t)] = \int_{t_0}^{t_1} L(\mathbf{q},\mathbf{\dot{q}},t)\ dt, כאשר המסלול הוא המסלול היחיד בו הפעולה מקבלת ערך אקסטרימום. הפונקציה הזו מקיימת את משוואת המילטון-יעקובי, והיא למעשה הפונקציה המנהלת של המילטון. בפיתוח זה, התנעים הקנוניים \boldsymbol{\alpha} שהפונקציה המנהלת תלויה בהן מופיעים כקבועי אינטגרציה של פתרון משוואת המילטון-יעקובי. הקואורדינטות הראשוניות והסופיות \mathbf{q}_0,\mathbf{q} קובעות את המהירויות הראשוניות ובאמצעותן את הערכים של \boldsymbol{\alpha}.

באופן דומה, הנגזרת המלאה לפי הזמן של הפונקציה המאפיינת מקיימת:

\frac{d W(\mathbf{q},\boldsymbol{\alpha})}{dt} = \sum_{i=1}^n \frac{\partial W}{\partial q_i} \dot{q}_i + \frac{\partial W}{\partial \alpha_i} \dot{\alpha}_i =\sum_{i=1}^n \frac{\partial W}{\partial q_i} \dot{q}_i = \sum_{i=1}^n p_i \dot{q}_i 

ולכן מתקיים W = \int \sum_{i=1}^n p_i \dot{q_i}\ dt = \int \sum_{i=1}^n p_i\ dq_i הביטוי הזה נקרא קשור לפעולה המקוצרת W[\mathbf{q}(t)] = \int_{\mathbf{q}_0}^{\mathbf{q}_1} \mathbf{p}\ d\mathbf{q} באופן דומה לאופן בו הפונקציה המנהלת קשור לפעולה.

## דוגמאות
### מתנד הרמוני חד־ממדי
מתנד הרמוני חד־ממדי מתואר על ידי ההמילטוניאן: H(q,p,t) = \frac{1}{2m}\left(p^2 + m^2\omega^2 q^2 \right) לכן משוואת המילטון יעקובי הלא תלויה בזמן עבור המתנד היא:\frac{1}{2m} \left[\left(\frac{\partial W}{\partial q}\right)^2+ m^2\omega^2 q^2 \right]=\alpha לאחר אלגברה נקבל W= \int \sqrt{2m\alpha - m^2 \omega^2 q^2}\ dq ומכאן t+\beta' = \frac{\partial W}{\partial \alpha} = m \int \frac{1}{\sqrt{2m\alpha-m^2\omega^2q^2}}\ dq= \sqrt{\frac{m}{2\alpha}} \int \frac{1}{\sqrt{1-\frac{m\omega^2q^2}{2\alpha} }}\ dq= \frac{1}{\omega} \arcsin\left(q\sqrt{\frac{m\omega^2}{2\alpha}}\right)

ומכאן q = \sqrt{\frac{2\alpha}{m\omega^2}}\sin(\omega t + \beta) וכן p = \frac{\partial W}{\partial q} = \sqrt{2m\alpha-m\omega^2 q^2} = \sqrt{2m\alpha} \cos(\omega t+\beta). \alpha היא האנרגיה הנשמרת של המערכת, ו-\beta היא הפאזה הנשמרת אף היא.

## הקשר למשוואת שרדינגר
נכתוב את משוואת המילטון־יעקובי הלא תלויה בזמן עבור מערכת קואורדינטות קרטזיות:

\frac{1}{2m} |\nabla W(x,y,z)|^2 + V(x,y,z) = E \implies |\nabla W(x,y,z)|^2 = 2m(E-V(x,y,z))

משוואה זו, היא משוואת האייקונל עבור הפונקציה W. לכן, באופן מקביל לתיאוריה האופטית, אם נניח שתנועת החלקיק היא תנועה של גל, והפונקציה \Psi = a\exp\left[\frac{i}{\hbar} (Et-W(x,y,z))\right] מתארת את הגל המונוכרומטי, נקבל שמשוואת הגלים שמקיימת הפונקציה \Psi היא:

|\nabla\Psi|^2 = \frac{2m(E-V)}{\hbar^2}\Psi 

וזוהי משוואת שרדינגר הלא תלויה בזמן. המשוואה התלויה בזמן יכולה להתקבל על ידי ההחלפה של E\Psi= i\hbar \frac{\partial \Psi}{\partial t} .

לחלופין, ניתן להתחיל ממשואת שרדינגר ולהשתמש בהצבה \Psi = a\exp\left[\frac{iS(x,y,z,t)}{\hbar} \right]. שימוש זה יניב:

\frac{|\nabla S|^2}{2m}+ \frac{i\hbar\nabla ^2 S}{2m}+ V = -\frac{\partial S}{\partial t}, משוואת המילטון־יעקובי מתקבלת בגבול הקלאסי \hbar\to0, בו מתקבל \frac{|\nabla S|^2}{2m}+ V = H = -\frac{\partial S}{\partial t} גבול זה מקביל לגבול של האופטיקה הגאומטרית בו אורך הגל קטן והפאזה גדולה.

## הערות שוליים

קטגוריה:מכניקה
קטגוריה:מכניקה אנליטית