עקרון הויגנס

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

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

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

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

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

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

### הניסוח המתמטי של העיקרון
שמאל|250px|ממוזער|המחשה גאומטרית לחישוב של פרנל.

נתייחס למקרה של מקור נקודתי הממוקם בנקודה **P**0, המייצר גלים בתדירות *f*. את ההפרעה ניתן לתאר על ידי המשתנה המרוכב *U*0, הנקרא גם משרעת מרוכבת. הוא מפיק גלים כדוריים עם אורך גל λ ומספר גל *k* = *2π/λ*. המשרעת המרוכבת של הגל הראשי בנקודה **Q** הממוקמת במרחק *r*0 מהנקודה **P**0 ניתנת בביטוי:

U(r_0) = \frac {U_0 e^{ikr_0}}{r_0} 

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

באמצעות התאוריה של הויגנס ועקרון הסופרפוזיציה של גלים, המשרעת המרוכבת בנקודה מרוחקת יותר **P** נמדדת על פי סיכום התרומות מכל נקודה על פני הכדור ברדיוס *r*0. כדי להשיג התאמה עם תצפיות ניסוייות, פרנל מצא שהתרומות מהגלים המשניים שנוצרים מנקודות על פני הכדור חייבים להיות מוכפלים בקבוע -i/\lambda (כאשר i = \sqrt ), ובפקטור הטיה נוסף K(\chi). ההתאמה הראשונה פירושה שהגלים המשניים מאחרים ביחס לגל הראשי ברבע מחזור גל, ושהמשרעת של הגלים המשניים היא ביחס של 1:\lambda  עם משרעת הגל הראשי. הוא הניח גם של-K(\chi) יש ערך מרבי כאשר \chi=0 (כלומר: בכיוון התקדמות הגל הראשי), וערך מאופס כאשר \chi=\pi/2 (כלומר: בניצב לכיוון התקדמות הגל הראשי).

על כן, המשרעת המרוכבת ב-**P** נתונה בנוסחה:

 U(P) = -\frac{i}{\lambda} U(r_0) \int_{S} \frac {e^{iks}}{s} K(\chi)\,dS 

כאשר *S* הוא המשטח של הכדור, ו-*s* הוא הרדיוס-וקטור  בין **Q** ל-**P**.

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

~K(\chi )= \frac{1}{2}(1+\cos \chi)

ל-*K* יש מקסימום ב- \chi=0 כמו בעקרון הויגנס-פרנל; אף על פי כן, *K* לא שווה לאפס כאשר \chi=\pi/2. נוסחת קירכהוף מסבירה מדוע בכלל גלים מתקדמים בכיוון מוגדר מסוים: כאשר \chi=\pi, מקבלים
0 = K(\chi), כלומר הגלים המשניים כלל לא מועברים לאחור.

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

### תיאור מתמטי
המשוואה שמקיים גל ממקור נקודתי היא משוואת הלמהולץ הלא-הומוגנית:
:\nabla^2 U + k^2 U(\vec{r},\vec{r'}) = \delta(\vec{r}-\vec{r'})
כאשר \delta(\vec{r}-\vec{r'}) היא פונקציית דלתא של דיראק התלת-ממדית, המתאפסת בכל המרחב פרט לנקודה \vec{r}=\vec{r'}. הגל המונוכרומטי התלוי בזמן שפותר את משוואת הגלים מתקבל על ידי הכפלת U(\vec{r},\vec{r'}) בפונקציה \ e^{-i \omega t}.
מאחר שפונקציית דלתא תלויה רק במרחק מהראשית ולא בזווית המרחבית, כי אז לאחר החלפת המשתנים: \vec{R} = \vec{r} - \vec{r'} ניתן להשתמש רק בחלק הרדיאלי של אופרטור הלפלסיאן בקואורדינטות כדוריות:
:\nabla ^2U= \frac{1}{R} \frac {\partial ^2}{\partial R^2} (R U) 
כאשר R=|\vec{R}|.
:המשוואה הופכת להיות:
:\frac {\partial ^2}{\partial R^2} (R U) + k^2 R U = R \delta(\vec{R})
מאחר שהפונקציה R מתאפסת בנקודה היחידה שבה פונקציית דלתא איננה מתאפסת, אגף ימין של המשוואה מתאפס באופן זהותי. פתרון המשוואה בקואורדינטות המקוריות הוא:
:U(\vec{r},\vec{r'}) \propto \frac{1} e^{ik|\vec{r}-\vec{r'}|}
זוהי משוואה המתארת גל כדורי, הבוקע מהנקודה \vec{r'}. את קבוע הפרופורציה ניתן למצוא באמצעות תנאי שפה. בנוסף, הפונקציה שהתקבלה היא לפי ההגדרה פונקציית גרין של המשוואה. לפי משפט גרין, משוואת הגל שנוצר ממקור דו-ממדי שאינו נקודתי ניתנת לחישוב באמצעות פונקציה זו על ידי האינטגרל המשטחי:
:U(\vec{r}) = \frac{k}{2 \pi i} \int\!\!\!\int U_{0}(\vec{r'}) \frac{e^{ik|\vec{r} - \vec{r'}|}} \cos \theta d^2r'
כאשר \ U_0 הוא מקור הגל, \theta היא הזווית בין וקטור השטח האינפיניטסימלי לבין הווקטור \vec{r}-\vec{r'}, והאינטגרל הוא על המשטח המחולל את הגל. אינטגרל זה מבטא את הגל הנצפה כסופרפוזיציה של גלים כדוריים, הבוקעים מנקודות שונות על פני המשטח.
:בפרט, אם במישור z=0 ניצב מחסום שלאחריו מתקבל הגל: \ U(x,y,z=0), אז על מסך במרחק z מתקבלת תבנית העקיפה:
:U(x,y,z) = \frac{k}{2 \pi i} \int\!\!\!\int U(x',y',z=0) \frac{e^{ik|\vec{r} - \vec{r'}|}} \cos \theta dx' dy'

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

## לקריאה נוספת
 
*Joseph W. Goodman, **Introduction to Fourier Optics**, 2nd edition, McGraw-Hill, 1996, Chapter 3

* ריצ'רד פיינמן (1988), *אור וחומר*, הוצאת הקיבוץ המאוחד

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

* 

קטגוריה:גלים