הרמוניות ספריות

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

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

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

 היסטוריה 
הרמוניות ספריות נחקרו לראשונה בהקשר של הפוטנציאל הניוטוני הנובע מחוק הכבידה העולמי בשלושה ממדים. ב-1782, פייר סימון לפלס בספרו "מכניקה שמיימית" עסק בבעיה של קביעת הפוטנציאל הכבידתי  בנקודה כלשהי x, שיוצרות אוסף של מסות נקודתיות mi המוצבות בנקודות xi, שנתון על ידי:

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

כאשר γ היא הזווית בין הווקטורים x ו-x1, והפונקציות  הן פולינומי לז'נדר (ראו גם פיתוח לפלס למולטיפולים אלקטרוסטטיים). שימו לב שאם הנקודה x מוגדרת להיות הראשית אז r=0, כך שקיבלנו את הביטוי הרגיל לפוטנציאל של מסה נקודתית. חשיבות פיתוחו של לז'נדר היא שהוא המפתח לתיאור הפוטנציאל של התפלגויות מסה כלליות, בין אם הן בדידות ובין אם הן רציפות. בעקבות כך, במאמרו מ-1782 לפלס חקר את המקדמים הללו תוך שהוא משתמש בקואורדינטות כדוריות כדי לייצג את הזווית γ בין x1 ל-x. 

ב-1867, ויליאם תומסון ופטר גונטרי טייט הציגו את ההרמוניות הספריות המרחביות בחיבור משותף, שבו גם טבעו את המונח "הרמוניות ספריות" עבור הפונקציות הללו. ההרמוניות המרחביות היו פתרונות פולינומיים הומוגניים  למשוואת לפלס:

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

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

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

כך שהן מייצגות אורביטלים אטומיים שונים.

מבוא
ממוזער|הרמוניות ספריות ממשיות Ylm עבור l מ-0 (שורה עליונה) עד 4 (שורה תחתונה). ההרמוניות השליליות Yl-m מסובבות סביב ציר ה-z בזווית של 90 חלקי m מעלות ביחס להרמוניות החיוביות.

משוואת לפלס בקואורדינטות כדוריות היא

עבור הפרדת משתנים  החלק הזוויתי של משוואת לפלס מקיים:

באמצעות הפרדת משתנים מקבלים ממשוואה דיפרנציאלית חלקית זו שתי משוואות דיפרנציאליות רגילות:

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

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

כאשר  נקראת הרמוניה ספרית ממעלה  וסדר m ו- הוא פונקציית לז'נדר הנלווית (ידועה גם כפולינומי לז'נדר הנלווים). הפקטור N הוא קבוע נרמול והזוויות θ ו-φ מייצגות את זווית הגובה (, היא שווה לאפס בקוטב הצפוני של הכדור, ל  בקו המשווה של הכדור ולפאי בקוטב הדרומי שלו) ואת הזווית המישורית . הקואורדינטות הספריות שמשמשות בערך זה עקביות עם אלה המשמשות פיזיקאים אך שונות מאלה המשמשות מתמטיקאים. 

כאשר פותרים את משוואת לפלס על פני השטח של כדור (על הספירה שלו), תנאי השפה המחזוריים על φ ותנאי רגולריות בקטבים, מבטיחות שהמעלה  והסדר m הם מספרים שלמים המקיימים  ≥ 0 ו-|m| ≤ . תנאי הרגולריות בקטבים הכרחי לצורך הקביעה כי l  מספר שלם, שכן למשוואת לז'נדר נקודות סינגולריות רגולריות ב-, ורק כאשר  הוא מהצורה  עבור l שלם כלשהו מובטח כי פתרון המשוואה בשיטת פרובניוס ייתן טור חזקות סופי של x (כלומר פולינום) המתכנס אל הפתרון עצמו.

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

כאשר  ו-  הם קבועים. האיברים בסכום הראשון שואפים לאפס כאשר r שואף לאינסוף, בעוד שהאיברים בסכום השני שואפים לאפס בראשית.

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

כך שהן אורתונורמליות ביחס למכפלה פנימית:

כאשר הדלתא היא דלתא של קרונקר המקיימת: δaa = 1, δab = 0 אם a ≠ b ו-dΩ = sinθ dφ dθ. 

בחקר הגאודזות והניתוח הספקטרלי משתמשים בנרמול אחר:

שמקיים

בחקר המגנטיות, משתמשים בחצי-נרמול שמידט

שמקיים

גם במכניקת הקוונטים משתמשים לעיתים קרובות בנרמול זה, והוא נקרא על שם הפיזיקאי הישראלי יואל רקח.

כאשר משתמשים בזהות הבאה (ראו פונקציות לז'נדר הנילוות)

ניתן להראות שכל ההרמוניות המנורמלות מקיימות

כאשר הכוכבית מסמנת לקיחת צמוד מרוכב. משוואה זו נובעת גם מהיחס בין ההרמוניות הספריות עם מטריצת D של ויגנר (ראו: Wigner D-matrix#Relation with spherical harmonic functions).

פאזת קונדון-שורטלי

מקור אחד לבלבול בהגדרה של הרמוניות ספריות עוסק בפקטור פאזה של m(1-) שבדרך כלל קוראים לו "פאזת קונדון-שורטלי" (Condon-Shortley) בספרות העוסקת במכניקת הקוונטים. בקרב החוקרים העוסקים בפיזיקה קוונטית, ישנה פרקטיקה נפוצה או לכלול את הפאזה בהגדרת פונקציות לז'נדר הנלווים, או לבלוע אותה בהגדרה של ההרמוניות הספריות. אין דרישה להשתמש בפאזת קונדון-שורטלי בהגדרת ההרמוניות הספריות, אך כלילתה יכולה לפשט מספר פעולות קוונטיות, ובפרט העלאה והורדה של אופרטורי תנע זוויתי. חוקרים העוסקים בגאודזות ומגנטיות לא כוללות את פאזת קונדון-שורטלי בהגדרות שלהם עבור הרמוניות ספריות.

פיתוח בהרמוניות ספריות

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

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

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

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

ראו גם רשימת הרמוניות ספריות ממשיות עד ל-l = 5. פונקציות אלה נבדלות בפאזה של (-1)m מאלה הניתנות בערך זה.

אנליזה ספקטרלית

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

כאשר

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

כאשר

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

וזה קירוב טוב בדרך כלל. כאשר β = 0 הספקטרום "לבן" כאשר כל מעלה היא בעלת הספק שווה. כאשר β < 0, אומרים שהספקטרום "אדום" ויש יותר הספק במעלות הנמוכות עם אורכי גל ארוכים מאשר במעלות הגבוהות. כאשר β > 0 אומרים שהספקטרום "כחול".

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

.

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

.

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



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

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

ראו גם: '''הרמוניות ספריות נוספות (עד l=10)

הכללות
ניתן לראות את ההרמוניות הספריות כהצגה של חבורת הסימטריה של הסיבובים סביב נקודה (SO(3)) או החבורה עם האלגברה הזהה לה SU(2). ככאלה, חבורות אלה תופסות את הסימטריה של ספירה דו-ממדית. כל סט של הרמוניות ספריות עם l נתון ממפות את החבורה להצגה אי-פריקה (irreducible representation) של SO(3).

בנוסף לכך, הספירה הדו-ממדית שקולה לספירת רימן. קבוצת הסימטריות השלמה של ספירת רימן מתוארת על ידי חבורת העתקות מביוס , שהיא איזומורפית כחבורת לי ממשית לחבורת לורנץ. האנלוג עבור ההרמוניות ספריות עבור חבורת לורנץ הם טורים היפרגאומטריים; ואכן, ניתן לבטא את ההרמוניות הספריות במונחים של טורים היפרגאומטריים, כשם ש-SO(3) היא תת-חבורה של PSL(2,C).

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

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

לקריאה נוספת
 A.R. Edmonds, Angular Momentum in Quantum Mechanics, (1957) Princeton University Press, .
 E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra, (1970) Cambridge at the University Press, , See chapter 3.
 J.D. Jackson, Classical Electrodynamics, 
 Albert Messiah, Quantum Mechanics, volume II. (2000) Dover. .
 D. A. Varshalovich, A. N. Moskalev, V. K. Khersonskii Quantum Theory of Angular Momentum,(1988) World Scientific Publishing Co., Singapore, 

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

קישורים כלליים

 Spherical harmonics applied to Acoustic Field analysis on Trinnov Audio's research page
 Spherical harmonics on Mathworld
 Spherical Harmonic Models of Planetary Topography
 Spherical harmonics generator in OpenGL
 OpenGL Spherical harmonics demo
 General Solution to LaPlace's Equation in Spherical Harmonics (Spherical Harmonic Analysis). Solid Earth Geophysics.

תוכנות
 SHTOOLS: Fortran 95 software archive
 HEALPIX: Fortran 90 and C++ software archive
 SpherePack: Fortran 77 software archive
 SpharmonicKit: C software archive
 Frederik J Simons: Matlab software archive
 NFFT: C subroutine library (fast spherical Fourier transform for arbitrary nodes)
 Shansyn: spherical harmonics package for GMT/netcdf grd files

הערות שוליים

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