אופטיקה של שכבות דקות

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

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

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

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

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

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

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

על מנת לבחון את אופי החלקיקים יש לפתור את משוואת הגלים במערכות השונות הגורמות להיווצרות החלקיקים.
לשם הפתרון, נניח שהמקדמים הדיאלקטרים (יסומנו \varepsilon ) של החומרים שנדון בהם קבועים, בנוסף נניח תלות זמנית הרמונית של השדה ונקבל את משוואת הלמהולץ. לשם פשטות נגדיר גאומטריית התפשטות חד־ממדית, הגלים מתפשטים בכיוון x בלבד ולכן \varepsilon=\varepsilon(z)
, מישור z=0 מקביל למישור התפשטות הגל ונקבל:

\vec{E}(x,y,z)=\vec{E}(z)e^{ixk_x}

כאשר k_x מסומן גם \beta והוא קבוע ההתפשטות.

ובנוסף:

\frac{\partial^2\vec{E}(z)}{\partial z^2}+(k^2_0\varepsilon -k^2_x )\vec{E}=0

ובאופן דומה עבור השדה המגנטי.

נמצא את המשוואות המצומדות בין השדה החשמלי לבין השדה המגנטי תחת מספר הנחות. האחת, תלות זמנית הרמונית עבורה \frac{\partial}{\partial t}=-i\omega. השנייה, עבור התפשטות בציר x עבורה מתקיים \frac{\partial}{\partial x}=i\beta . שלישית, התפשטות הומוגניות בציר y עבורה \frac{\partial}{\partial y}=0. המשוואות המצומדות שמתקבלות הן:

\frac{\partial E_y}{\partial z}= -i\omega\mu_0H_x

\frac{\partial E_x}{\partial z}-i\beta E_z= i\omega\mu_0H_y

i\beta E_y=i\omega\mu_0H_z

\frac{\partial H_y}{\partial z}= i\omega\varepsilon_0\varepsilon E_x

\frac{\partial H_x}{\partial z}-i\beta H_z= -i\omega\varepsilon_0\varepsilon E_y

i\beta H_y= -i\omega\varepsilon_0\varepsilon E_z

מתקבלים שני סטים של משוואות עם קיטוב שונה של הגל המתפשט, TM ו-TE, כך שעבור כל קיטוב מספר המשואות מצטמצם ונקבל את משוואות הגל בהתאמה:

\frac{\partial^2 H_y}{\partial z^2}+(k^2_0\varepsilon-\beta^2)H_y=0

\frac{\partial^2 E_y}{\partial z^2}+(k^2_0\varepsilon-\beta^2)E_y=0

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

	
### פלסמון-פולריטון SPPs
#### מערכות דו שכבתית
מערכת דו שכבתית היא מערכת של מתכת-חומר דיאלקטרי (IM-isolator metal) או חומר דיאלקטרי-מתכת (MI).	
נסתכל על מערכת דו שכבתית, השקה יחידה, של מוליך בעל מקדם דיאלקטרי שיסומן \varepsilon_1 וחומר דיאלקטרי בעל מקדם דיאלקטרי ממשי חיובי שיסומן \varepsilon_2.

פתרון TM בחומר הדיאלקטרי:


H_y(z)=
A_1e^{ik_xx}e^{k_1z}

E_x(z)=-i
\frac{A_1k_1e^{ik_xx}e^{k_1z}}{\omega\varepsilon_0\varepsilon_1}

E_z(z)=-
\frac{A_1\beta e^{ik_xx}e^{k_1z}}{\omega\varepsilon_0\varepsilon_1}

כאשר A_1 הוא קבוע שערכו לא ידוע עדיין.

פתרון TM במתכת:

H_y(z)=
A_2e^{ik_xx}e^{k_2z}

E_x(z)=i
\frac{A_2k_2e^{ik_xx}e^{k_2z}}{\omega\varepsilon_0\varepsilon_2}

E_z(z)=-
\frac{A_1\beta e^{ik_xx}e^{k_2z}}{\omega\varepsilon_0\varepsilon_2}


יש לדרוש את רציפות השדה החשמלי והשדה המגנטי בנקודת ההשקה z=0 ונקבל:

A_1=A_2

\frac{k_2}{k_1}=-\frac{\varepsilon_2}{\varepsilon_1}

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

יחס הדיספרציה או יחס נפיצה המתקבל:

\beta=k_x=k_0\sqrt{\frac{\varepsilon_1\varepsilon_2}{\varepsilon_1+\varepsilon_2}}

נסתכל על הקשר בין וקטור גל k לאורך גל \lambda k=\frac{2\pi}{\lambda}, ניתן לראות שקיים יחס הפוך בין וקטור הגל לאורך הגל, כלומר, העדר חסם עליון על וקטור הגל מעיד על העדר חסם תחתון על אורך הגל. המשמעות היא שחלקיק הפלסמון-פולריטון תומך באורכי גל קצרים מאד וניתן באמצעותו לשבור את גבול הדיפרקציה. גבול אשר מהווה חסם תחתון באופטיקה של חומרים תלת־ממדיים. היתרון באורכי גל קצרים הוא הן לטובת מיקוד הגל והן לטובת הגדלת עוצמתו.

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

L=\frac{1}{2Im{\beta}}

#### מערכת תלת שכבתית

מערכת תלת שכבתית היא מערכת של מתכת-חומר דיאלקטרי-מתכת (MIM) או חומר דילאקטרי-מתכת-חומר דיאלקטרי (IMI).	
בכל מישור של ממשק בין שתי שכבות יכול להיווצר SPPs.
נגדיר אורך דעיכה \hat{z}=\frac{1} המתאר את הדעיכה בכיוון z לתוך נפח החומרים.
כאשר המרחק בין ממשקים סמוכים גדול מאורך הדעיכה \hat{z} ̂ המערכת שקולה לשתי מערכות דו שכבתיות נפרדות אשר ניתן לפתור כל אחת מהן בנפרד ללא תלות בשנייה, כפי שמוצג קודם. אולם, כאשר המרחק שווה או קטן מאורך הדעיכה \hat{z} ̂ נוצרת אינטראקציה בין SPPs הנוצרים על כל ממשק ואינטראקציה זו גורמת ליצירת מודים מצומדים.

נציג פתרון עבור IMI כאשר הפתרון עבור MIM דומה.

עבור קיטוב TM:

בתחום בו z>a:


H_y=
Ae^{ik_xx}e^{-k_3z}


E_x=i
\frac{Ak_3e^{ik_xx}e^{-k_3z}}{\omega\varepsilon_0\varepsilon_3}

E_z=-
\frac{A\beta e^{ik_xx}e^{-k_3z}}{\omega\varepsilon_0\varepsilon_3}


בתחום בו z>-a:

H_y=
Be^{ik_xx}e^{-k_2z}


E_x=-i
\frac{Bk_2e^{ik_xx}e^{-k_2z}}{\omega\varepsilon_0\varepsilon_2}


E_z=-
\frac{iB\beta e^{ik_xx}e^{-k_2z}}{\omega\varepsilon_0\varepsilon_2}



בתחום בו -a:

H_y=
Ce^{ik_xx}e^{k_1z}+De^{ik_xx}e^{-k_1z}


E_x=
\frac{-iCk_1e^{ik_xx}e^{k_1z}}{\omega\varepsilon_0\varepsilon_1}+\frac{iDk_1e^{ik_xx}e^{-k_1z}}{\omega\varepsilon_0\varepsilon_1}


E_z=
\frac{C\beta e^{ik_xx}e^{k_1z}}{\omega\varepsilon_0\varepsilon_2}+\frac{D\beta e^{ik_xx}e^{-k_1z}}{\omega\varepsilon_0\varepsilon_2}



על מנת לפתור את מערכת המשוואות הזו יש לדרוש רציפות על רכיב y של השדה המגנטי ועל רכיב x של השדה החשמלי. כמו כן, להיעזר במשוואת וקטור הגל k^2_i=\beta^2-k^2_0\varepsilon_i כאשר i= 1, 2, 3. באמצעות אלה ניתן לקבל 4 משוואות מצומדות אשר פתרונן מוביל לביטוי עבור יחס הדיספרסיה של המערכת:

e^{-4k_1a}=\frac{\frac{k_1}{\varepsilon_1}+\frac{k_2}{\varepsilon_2}}{\frac{k_1}{\varepsilon_1}-\frac{k_2}{\varepsilon_2}}\frac{\frac{k_1}{\varepsilon_1}+\frac{k_3}{\varepsilon_3}}{\frac{k_1}{\varepsilon_1}-\frac{k_3}{\varepsilon_3}}

אפשר לשים לב שעבור a\to\infin מתקבל הביטוי ליחס הדיספרסיה שך מערכת דו שכבתית.

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

\tanh k_1a=-\frac{k_2\varepsilon_1}{k_1\varepsilon_2}

\tanh k_1a=-\frac{k_1\varepsilon_2}{k_2\varepsilon_1}

### אקסיטון-פולריטון 2DEP
אקסיטון הוא חלקיק הנוצר כאשר בהיווצרות אלקטרון-חור (חשמל) נוצר כוח קולון בין השניים כך שיחדיו הם מהווים חלקיק נייטרלי חשמלית הנקרא אקסיטון. זמן החיים של חלקיק זה מוגבל וכאשר מתרחשת רקומבינציה נפלט פוטון.
קוואזי חלקיק זה ניתן לעירור במוליך למחצה דו־ממדי חד שכבתי מסוג TMDs בתחום הספקטרום הנראה. דוגמאות לחומרי TMDs הם MoS2, MoSe2, WS2 ו־WSe2.
כפי שמפורט בפיתוח עבור פלסמון-פולריטון, גם עבור היווצרות אקסיטונים, בממשק בין מל"מ לבין חומר דיאלקטרי, יש לדרוש שהמל"מ יהיה בעל מקדם דיאלקטרי עם רכיב ממשי שלילי (כיוון שלחומר הדיאלקטרי מקדם דיאלקטרי חיובי). על מנת לבחון דרישה זו, נסתכל על הסוספטיביליות חשמלית χ של ה־TMDs ועל הקשר שלה למקדם הדיאלקטרי \varepsilon :

\chi=1+\varepsilon

\chi(\omega)=\chi_{bg}-\frac{\frac{c}{\omega_0d}\gamma_{r0}}{\omega-\omega_0+i(\frac{\gamma_{nr}}{2}+\gamma_d)}

כאשר \chi_{bg} היא סוספטביליות הרקע, c היא מהירות האור, \omega_0 היא אנרגיית הקשר של האקסיטון, d מייצג את עובי השכבה החד אטומית ו- \gamma_d, \gamma_{nr},\gamma_{r0} הם מקדמי דעיכה רדיואקטיבית, לא רדיואקטיבית ודעיכה טהורה.

נהוג להגדיר עבור אקסיטונים: \gamma_T= 2\gamma_d + \gamma_{nr}+\gamma_{r0}, גודל זה מייצג את הרוחב בכולל של הרזוננס של האקסיטון ומשפיע על האמפליטודה של הסוספטביליות ולכן גם על האמפליטודה של המקדמים הדיאלקטרים. עבור \gamma_T קטן מספיק האמפליטודות של המקדמים הדיאלקטרים יורדות מספיק ומאפשרות למקדם הדיאלקטרי של ה־TMD להגיע לערכים שליליים. ניתן לשנות את \gamma_T באמצעות שינוי טמפ' ובכך לקבל \gamma_T צר מספיק כך שבתחום אנרגיות מסוים המקדם הדיאלקטרי של ה־TMD יהיה שלילי ויוכל לתמוך בהיווצרות אקסיטונים.

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

נפתור את יחס הדיספרסיה של האקסיטון עבור קיטוב TM ונקבל:

\frac{\varepsilon_1}{k_1}+\frac{\varepsilon_2}{k_2}+\frac{\sigma(\omega)}{\omega\varepsilon_0}=0

\sigma(\omega)=-i\chi\varepsilon_0\omega

k_{1,2}=\sqrt{q^2_{2DEP}-\frac{\omega^2\varepsilon_{1,2}}{c^2}}

כאשר q_{2DEP} הוא תנע האקסיטון ו- \sigma(\omega) היא המוליכות.

## יתרונות וחידושים של האופטיקה בחומרים הדו־ממדיים

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

## ראו גם
* התאבכות
* התאבכות משכבות דקות
* שיטת העברת מטריצות

## קישורים חיצוניים
* 
* PLASMONICS: Fundamentals and Applications (2006) (pg 3-50), Stefan A. Maier
*Itai Epstein at el, Highly confined in-plane propagating exciton-polaritons on monolayer semiconductors, IOP publishing 2D Mater. 7 (2020) 035031, https://iopscience.iop.org/article/10.1088/2053-1583/ab8dd4
*מאמר: Joshua D. Caldwell at el, Low-loss, infrared and terahertz nanophotonics using surface phonon polaritons, Nanophotonics 2015); 4: 44–68)
*מאמר: J. A. Dionne at el, Planar metal plasmon waveguides: frequency-dependent dispersion, propagation, localization,
(and loss beyond the free electron model, Physical Review B 72, 075405 (2005

קטגוריה:אופטיקה