בעיית שני הגופים

שני גופים נעים במסלולים אליפטיים סביב [[מרכז מסה משותף]]
במסתם נעים סביב מרכז מסה משותף, הנמצא בתוך הגוף הגדול.]]
בעיית שני הגופים היא בעיה פיזיקלית העוסקת בתיאור תנועתם של שני גופים המפעילים כוח כלשהו זה על גבי זה וללא כוחות חיצוניים. בעיה זו היא אחת הבעיות המרכזיות בפיזיקה, בגלל היישומים החשובים שלה וגם מכיוון שהיא אחת הבעיות היחידות הניתנת לפתרון אנליטי (עבור כוחות מסוימים), במכניקה קלאסית, בתורת הקוונטים, ולפתרון מקורב בתורת היחסות. דוגמאות ליישומים של הבעיה הדו-גופית הם בתיאור תנועת כוכבי לכת מסביב לשמש, בתיאור תנועת האלקטרון סביב הפרוטון באטום המימן, תנועת מולקולה דו אטומית ועוד. 

דוגמאות נוספות לכך כוללות מקרים בהם התנועה אינה חסומה (הגופים מתרחקים למרחק אינסופי זה מזה), כגון אסטרואיד העובר סמוך לשמש, מוסט על ידיה, וממשיך במסלולו מבלי לחזור. מקרים אלו נקראים התנגשויות אלסטיות, והחקר שלהם הוא הבסיס לתורת הפיזורים (scattering theory). להתנגשויות אלסטיות שימושים רבים, החל בתנועת גופים שמימיים, המשך בתיאור תנועת כדורים בביליארד, וכלה בניסויים המנגשים חלקיקים אלמנטריים במאיצים.

יש לציין שבפועל הבעיה הדו-גופית משמשת קירוב למציאות, כיוון שמוזנחים כוחות חיצוניים (לדוגמה השפעת מאדים על תנועת כדור-הארץ סביב השמש), והמבנה הפנימי של שני הגופים (כגון גאות ושפל על פני כדור הארץ, התנגשויות פלסטיות).

==היסטוריה==
הבעיה הדו-גופית במכניקה קלאסית, עבור כח הכבידה, נפתרה לראשונה בצורה אנליטית על ידי אייזק ניוטון, ואפשרה לתת הסבר תאורטי לחוקי קפלר שתיארו את תנועת הכוכבים סביב השמש. הישג זה היווה את גולת הכותרת של המכניקה הניוטונית. 
מודל האטום של בוהר
הבעיה הדו-גופית באלקטרודינמיקה, תנועת האלקטרון סביב הפרוטון באטום המימן, שיחקה תפקיד חשוב בהתפתחות תורת הקוונטים: ניתוח ספקטרום הפליטה של אטום המימן הוביל למודל האטום של בוהר, ומאוחר יותר למשוואת שרדינגר.

תורת הפיזורים שימשה את ארנסט רתרפורד להסיק את מבנה האטום, מתוך פיזורים של קרינת אלפא על אטומי זהב, ובכך הולידה את הפיזיקה הגרעינית. מאז הניסוי של רתרפורד ניסויי התנגשויות הפכו לכלי מרכזי בפיזיקה גרעינית ובפיזיקה של חלקיקים אלמנטריים. פיזורים משחקים תפקיד רב גם בהתפתחות האסטרופיזיקה. בשנת 1977, הצליחה נאס"א למצוא שימוש מקורי בפיזורים: היא הצליחה לשלוח שתי חלליות, וויאג'ר 1 ווויאג'ר 2, אל כוכבי הלכת הרחוקים ביותר של מערכת השמש באמצעות פיזור של החללית מכוכב אחד לכוכב הבא אחריו בתהליך שנקרא מקלעת כבידתית.

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

150px
הפתרון מבוסס על שני שלבים: 
*פירוק הבעיה לשתי בעיות חד-גופיות: תיאור תנועת מרכז המסה של שני הגופים, שהיא למעשה תנועה חופשית; ותיאור התנועה היחסית של שני הגופים. הפירוק הזה מתאפשר בזכות חוק שימור התנע, והפיתוח המתמטי שלו מוצג בסעיף הבא. בעת שאחד הגופים כבד הרבה יותר מהשני, כמו במקרה של תנועת כדור-הארץ סביב השמש, אזי מרכז המסה נמצא סמוך לגוף המסיבי, ולכן הגוף המסיבי כמעט ואיננו נע. במקרה זה הבעיה מצטמצמת לבעיה של תנועת הגוף הקטן סביב הגוף הגדול.
את הפוטנציאל שבין שני גרעינים בתוך מולקולה ניתן לקרב על ידי [[פוטנציאל מורס]]

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

הרדוקציות הללו הן כלליות: פירוק לשתי בעיות חד-גופיות, ניתן לביצוע תמיד וצמצום לבעיה חד־ממדית ניתן לביצוע עבור כוח מרכזי. 

כוח גרוויטציוני, וכן כוח קולון (הכוח החשמלי האלקטרוסטטי) בין שני גופים הם מהצורה <math>\frac{C}{r^2} </math>, ובמקרה זה הבעיה נקראת בעיית קפלר. הפתרונות של בעיה זאת הם חתכים קוניים: המסלולים החסומים הם בעלי צורה של אליפסה (כך לדוגמה נע כדור-הארץ סביב השמש) או מעגל, והמסלולים הלא-חסומים הם בעלי צורה של היפרבולה (כך לדוגמה נעים אסטרואידים החולפים בקרבת השמש) או פרבולה. 

בתיאור תנועתם של גרעינים בתוך מולקולות דו אטומיות, ניתן להשתמש בעקרונות של הבעיה הדו-גופית, אלא שהכוחות בין הגרעינים הם מורכבים יותר ונגרמים מתנועתם של האלקטרונים הקושרים אותם. במקרים אלו ניתן לקרב את הפוטנציאל בין שני הגרעינים על ידי פוטנציאל מורס (Morse Potential).

==מבעיית שני הגופים לבעיה של גוף יחיד בשדה כוח==
נסמן את מסת הגופים ב-<math>\ m_1,m_2 </math> ואת מיקומם ב-<math> \vec r_1,\vec r_2 </math>. 
משוואות התנועה הכלליות עבור שני גופים הן:
<div style="text-align: center;">
<math> m_1 \ddot\vec r_1 = \vec F_{21} + \vec F_1^{ext} </math>

<math> m_2 \ddot\vec r_2 = \vec F_{12} + \vec F_2^{ext} </math>
</div>

כאשר <math> \vec F_{ij} </math> הוא הכוח שגוף i מפעיל על גוף j, ו-<math> F_i^{ext} </math> הוא כוח חיצוני שפועל על גוף i. על פי החוק השלישי של ניוטון <math> \vec F_{12} = - \vec F_{21} = \vec F</math> ועבור בעיה דו-גופית מניחים שאין כוחות חיצוניים, כלומר עבור בעיה דו-גופית ניתן לכתוב:
<div style="text-align: center;">
<math> m_1 \ddot\vec r_1 = \vec F </math>

<math> m_2 \ddot\vec r_2 = -\vec F </math>
</div>
על ידי חיבור וחיסור המשוואות הנ"ל ניתן לקבל משוואות המתארות את תנועת מרכז המסה ואת התנועה היחסית בין הגופים:

===תנועת מרכז המסה===
על ידי חיבור המשוואות מקבלים:
<math display="block"> m_1 \ddot \vec r_1 + m_2 \ddot \vec r_2 = 0 </math>
משוואה זו היא משוואת תנועה עבור מרכז המסה, אשר הקואורדינטות שלו נתונות על ידי <math>\vec r_{c.m} = \frac{1}{M} [m_1 \vec r_1 + m_2 \vec r_2] </math> (כאשר <math> M = m_1 +m_2 </math> היא המסה הכוללת). כלומר את המשוואה האחרונה ניתן לכתוב כ: <math>\ M \ddot \vec R = 0</math>. על מרכז המסה לא פועלים כוחות (כיוון שאין כוחות חיצוניים) ולכן על פי החוק הראשון של ניוטון היא נעה במהירות קבועה.

===תנועה יחסית===
משוואות התנועה של שני הגופים <math>\ m_1, m_2 </math>, ניתנות לכתיבה כ:
<div style="text-align: center;">
<math> \ddot \vec r_1 = \frac{\vec F(\vec r_1 - \vec r_2)}{m_1} </math>
;
<math> \ddot \vec r_2 = -\frac{\vec F(\vec r_1 - \vec r_2)}{m_2} </math>
</div>
כאשר כאן נכתבה באופן מפורש העובדה שבשל אינווריאנטיות להזזה הכוח <math>\vec F </math> הפועל בין שני הגופים הוא פונקציה של <math> \vec r_1 - \vec r_2 </math> בלבד ולא של <math> \vec r_1 </math> או <math> \vec r_2 </math> בנפרד.

נגדיר קואורדינטה יחסית <math> \vec r = \vec r_1 - \vec r_2 </math>, המודדת את מיקום הגוף הראשון יחסית למיקום הגוף השני. 
על ידי חיסור המשוואות עבור <math> \vec r_1,\vec r_2 </math> מקבלים משוואה עבור הקואורדינטה היחסית:
<math display="block"> \ddot \vec r = \left(\frac{1}{m_1}+\frac{1}{m_2}\right) \vec F(\vec r) = \frac{m_1+m_2}{m_1 m_2}\vec F(\vec r)</math>

נהוג להגדיר מסה מצומצמת <math> \mu = \frac{m_1 m_2}{m_1+m_2} </math> בעזרתה ניתן לכתוב:
<math display="block"> \mu \ddot \vec r = \vec F(\vec r) </math>
כלומר זו פשוט משוואת תנועה של גוף בעל מסה <math>\ \mu </math> עליו פועל כוח <math> \vec F </math>.

===סיכום===
בבעיית שני הגופים שקולה מתמטית לשתי בעיות בלתי תלויות פשוטות יותר:
* תנועת מרכז המסה - כיוון שעל מרכז המסה לא פועל כל כוח, זוהי בעיה טריוויאלית.
* תנועה יחסית סביב מרכז המסה - זו בעיה השקולה לבעיה של גוף יחיד הנע בהשפעת כוח נתון. במקרים רבים הכוח הנ"ל הוא כוח מרכזי. במקרים אלו הבעיה פשוטה יחסית.

לאחר פתרון שתי הבעיות הנ"ל, ניתן לשלבן באופן הבא על מנת לקבל את תנועת שני הגופים:
<div style="text-align: center;">
<math>\vec r_1(t) = \vec R(t) + \frac{m_2}{m_1+m_2}\vec r(t) </math>

<math>\vec r_2(t) = \vec R(t) - \frac{m_1}{m_1+m_2}\vec r(t) </math>
</div>

ניתן לראות שכאשר אחת המסות גדולה בהרבה מהמסה השנייה, כמו במקרה של כדור-הארץ והשמש, המסה הכבדה יותר מבצעת תנועה קטנה מאוד, בקירוב טוב היא נשארת נייחת, וכן ניתן לראות שהמסה המצומצמת כמעט שווה למסה הקטנה יותר. 

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

==מבעיה תלת־ממדית לבעיה דו־ממדית==

בהינתן גוף הנע בהשפעת שדה כוחות מרכזי במרחב (זאת אומרת שדה כוח שמפעיל בכל נקודה כוח על הגוף בכיון ראשית הצירים), ניתן להראות, באמצעות חוק שימור התנע הזוויתי, שהגוף נע במישור קבוע.

התנע הזוויתי (של גוף נקודתי ביחס לראשית) נתון על ידי הנוסחה:
<math>\vec L=\vec r \times m \vec v</math>,
כאשר:
* <math>\vec r</math> - הווקטור המחבר מהראשית לגוף
* <math>\vec v</math> - מהירות הגוף 
* <math>m</math> - מסת הגוף
* <math>\times</math> - מכפלה וקטורית
מהחוק השני של ניטון נובע כי קצב השינוי של התנע הזוויתי הוא מומנט הכוח הפועל על הגוף:
<math display="block">\dot{\vec L}=\frac{d}{dt}(\vec r\times m \vec v)= \dot \vec r\times m \vec v+ \vec r\times m \dot \vec v=0+ \vec r\times \vec F=\vec r\times \vec F</math>
לכן, אם הכוח הוא מרכזי (ז"א מקביל ל- <math>\vec r</math>) אז התנע הזוויתי לא משתנה. 
מאידך הווקטור <math>\vec r</math> תמיד מאונך לתנע הזוויתי. לכן הגוף נמצא תמיד במשור המאונך ל- <math>\vec L</math>.

==מבעיה דו־ממדית לבעיה חד־ממדית==

אם הכוח הפועל על הגוף תלוי רק במרחק בין הגוף לראשית, ניתן לפשט את הבעיה על ידי מעבר לקואורדינטות קוטביות במישור. ניתן להביע את גודל התנע הזוויתי על ידי:
<math display="block">L=\mu r^2 \dot \theta </math>
ואת האנרגיה הכוללת על ידי:
<math display="block"> E = \frac{1}{2} \mu \dot r^2 +\frac{L^2}{2\mu r^2}+V(r)</math>
בביטוי האחרון שני האיברים הראשונים הם האנרגיה הקינטית (המתאימים לתנועה הרדיאלית והסיבובית) והשלישי הוא האנרגיה הפוטנציאלית.
הפוטנציאל האפקטיבי בבעיה הדו-גופית]]
כיוון שהתנע הזוויתי <math>L</math> קבוע ניתן להתייחס אל <math> \frac{L^2}{2\mu r^2} </math> כאל חלק מהאנרגיה הפוטנציאלית, ולהתייחס לגוף כנע במימד אחד (הקואורדינטה r) בהשפעת פוטנציאל אפקטיבי <math> V_{eff}(r) = V(r) +\frac{L^2}{2\mu r^2} </math>. ניתן לראות שהפוטנציאל הנובע מהתנע הזוויתי הולך לאין-סוף כאשר r הולך ל-0 ולכן הגופים המושכים זה את זה לא יכולים להתקרב עד למרחק 0 אחד לשני.
משוואת התנועה עבור הקואורדינטה r תהיה:
<math display="block">\mu \frac{d^2r}{dt^2}=-\frac{dV_{eff}}{dr} </math>

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

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

בתורת הקוונטים המקובלת היום, לפי הניסוח שלה על ידי משוואת שרדינגר (או משוואת הייזנברג), ניתן לפתור את הבעיה הדו-גופית באופן דומה לאופן שבה היא נפתרת במכניקה הקלאסית: פירוק הבעיה לשתי בעיות חד-גופיות, צמצום הבעיה למימד אחד, ולבסוף פתרון משוואת הגל במימד הנותר. 

===צמצום הבעיה לשתי בעיות חד-גופיות===
ההמילטוניאן של בעיה דו-גופית (עם פוטנציאל מרכזי וללא התחשבות בספין) יהיה מן הצורה:
<math display="block"> \mathcal{H} = \frac{\vec p_1^2}{2m_1}+\frac{\vec p_2^2}{2m_2} + V(|\vec r_1 - \vec r_2|) </math>
נגדיר אופרטורים חדשים המתאימים למיקום ולתנע של מרכז המסה:
<math display="block"> \vec R = \frac{m_1 \vec r_1 +m_2 \vec r_2}{m_1+m_2}\, \ \vec P = \vec p_1 + \vec p_2 </math>
ולמיקום והתנע היחסיים:
<math display="block"> \vec r = \vec r_1 - \vec r_2 \, \ \vec p = \frac{m_2 \vec p_1 - m_1 \vec p_2}{m_1+m_2}</math>
האופרטורים הנ"ל מקיימים יחסי חילוף קנוניים ובעזרתם ניתן לכתוב את ההמילטוניאן כ:
<math> \mathcal{H} = H_{cm}+H_{rel}</math>,
כאשר:
*<math> H_{cm} = \frac{\vec P^2}{2M} </math> הוא החלק המתאר את מרכז המסה (המתנהג כחלקיק חופשי).
*<math> H_{rel} = \frac{\vec p^2}{2m} + V(r)</math> הוא החלק המתאר את התנועה היחסית, הנע בהשפעת הפוטנציאל המרכזי <math>\ V(r) </math>.

האופרטורים <math>\ H_{cm} </math> ו-<math>\ H_{rel} </math> חילופיים ולכן ניתן לטפל בהם בנפרד. בדומה למקרה הקלאסי, החלק של מרכז המסה הוא טריביאלי.

===הקוארדינטות הזוויתיות והקוארדינטה הרדיאלית===
השלב השני (המוסבר בפירוט רב יותר כאן), רדוקציה לבעיה חד־ממדית, שוב מתאפשר הודות לשימור תנע זוויתי, או ליתר דיוק מכך שהבעיה סימטרית לסיבובים. משום כך ניתן להפריד את משוואת שרדינגר לשתי משוואות, אחת עבור המרחק בין הגרעינים, r, ואחת עבור החלק הזוויתי. הרדוקציה מתבצעת על ידי הפרדת משתנים: כותבים את פונקציית הגל באופן <math>\ \psi(\vec{r}) = R(r) Y_{l,m}(\theta,\phi)</math> ואזי משוואת שרדינגר נפרדת לשתי משוואות:

* משוואה עבור החלק הזוויתי:
**: <math>\ L^2 = \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 \phi^2} </math>
*: פתרון משוואה זאת הן ההרמוניות הספריות <math>\ Y_{l,m}</math>.
* משוואה עבור החלק הרדיאלי:
**: <math>\ \frac{1}{r} \frac{\partial ^2 }{\partial r^2} \left( r R(r) \right) + \frac{\hbar^2 l (l+1)}{2 m r^2} R(r) + V(r)R(r) = E \ R(r)</math>
*: פתרון המשוואה הרדיאלית לגבי <math>\ u(r) = rR(r)</math> מערב בדרך כלל פונקציות בסל.

==ראו גם==
* בעיית קפלר - בעיה דו-גופית עם כוח ריבועי הפוך (כגון כוח הכבידה)
* בעיית שלושת הגופים
* משפט ברטראן
* כוח מרכזי
* שיטת החתכים הטלואים

==קישורים חיצוניים==

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

קטגוריה:ערכים שבהם תבנית בריטניקה אינה מתאימה
קטגוריה:מכניקה
קטגוריה:מכניקת הקוונטים