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

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

חומר דו ממדי והאינטראקציה עם אור

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

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

על מנת לבחון את אופי החלקיקים יש לפתור את משוואת הגלים במערכות השונות הגורמות להיווצרות החלקיקים. לשם הפתרון, נניח שהמקדמים הדיאלקטרים (יסומנו $\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_{0}^{2}\varepsilon -k_{x}^{2}){\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 _{0}H_{x}$

${\frac {\partial E_{x}}{\partial z}}-i\beta E_{z}=i\omega \mu _{0}H_{y}$

$i\beta E_{y}=i\omega \mu _{0}H_{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_{0}^{2}\varepsilon -\beta ^{2})H_{y}=0$

${\frac {\partial ^{2}E_{y}}{\partial z^{2}}}+(k_{0}^{2}\varepsilon -\beta ^{2})E_{y}=0$

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

פלסמון-פולריטון SPPs

מערכות דו שכבתית

מערכת דו שכבתית היא מערכת של מתכת-חומר דיאלקטרי (IM-isolator metal) או חומר דיאלקטרי-מתכת (MI). נסתכל על מערכת דו שכבתית, השקה יחידה, של מוליך בעל מקדם דיאלקטרי שיסומן $\varepsilon _{1}$ וחומר דיאלקטרי בעל מקדם דיאלקטרי ממשי חיובי שיסומן $\varepsilon _{2}$.

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

$H_{y}(z)=A_{1}e^{ik_{x}x}e^{k_{1}z}$

$E_{x}(z)=-i{\frac {A_{1}k_{1}e^{ik_{x}x}e^{k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{1}}}$

$E_{z}(z)=-{\frac {A_{1}\beta e^{ik_{x}x}e^{k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{1}}}$

כאשר $A_{1}$ הוא קבוע שערכו לא ידוע עדיין.

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

$H_{y}(z)=A_{2}e^{ik_{x}x}e^{k_{2}z}$

$E_{x}(z)=i{\frac {A_{2}k_{2}e^{ik_{x}x}e^{k_{2}z}}{\omega \varepsilon _{0}\varepsilon _{2}}}$

$E_{z}(z)=-{\frac {A_{1}\beta e^{ik_{x}x}e^{k_{2}z}}{\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}{|k_{z}|}}$ המתאר את הדעיכה בכיוון z לתוך נפח החומרים. כאשר המרחק בין ממשקים סמוכים גדול מאורך הדעיכה ${\hat {z}}$ ̂ המערכת שקולה לשתי מערכות דו שכבתיות נפרדות אשר ניתן לפתור כל אחת מהן בנפרד ללא תלות בשנייה, כפי שמוצג קודם. אולם, כאשר המרחק שווה או קטן מאורך הדעיכה ${\hat {z}}$ ̂ נוצרת אינטראקציה בין SPPs הנוצרים על כל ממשק ואינטראקציה זו גורמת ליצירת מודים מצומדים.

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

עבור קיטוב TM:

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

$H_{y}=Ae^{ik_{x}x}e^{-k_{3}z}$

$E_{x}=i{\frac {Ak_{3}e^{ik_{x}x}e^{-k_{3}z}}{\omega \varepsilon _{0}\varepsilon _{3}}}$

$E_{z}=-{\frac {A\beta e^{ik_{x}x}e^{-k_{3}z}}{\omega \varepsilon _{0}\varepsilon _{3}}}$

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

$H_{y}=Be^{ik_{x}x}e^{-k_{2}z}$

$E_{x}=-i{\frac {Bk_{2}e^{ik_{x}x}e^{-k_{2}z}}{\omega \varepsilon _{0}\varepsilon _{2}}}$

$<math>E_{z}=-{\frac {iB\beta e^{ik_{x}x}e^{-k_{2}z}}{\omega \varepsilon _{0}\varepsilon _{2}}}$ </math>

בתחום בו $-a<z<a$:

$H_{y}=Ce^{ik_{x}x}e^{k_{1}z}+De^{ik_{x}x}e^{-k_{1}z}$

$E_{x}={\frac {-iCk_{1}e^{ik_{x}x}e^{k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{1}}}+{\frac {iDk_{1}e^{ik_{x}x}e^{-k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{1}}}$

$E_{z}={\frac {C\beta e^{ik_{x}x}e^{k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{2}}}+{\frac {D\beta e^{ik_{x}x}e^{-k_{1}z}}{\omega \varepsilon _{0}\varepsilon _{2}}}$

על מנת לפתור את מערכת המשוואות הזו יש לדרוש רציפות על רכיב y של השדה המגנטי ועל רכיב x של השדה החשמלי. כמו כן, להיעזר במשוואת וקטור הגל $k_{i}^{2}=\beta ^{2}-k_{0}^{2}\varepsilon _{i}$ כאשר $i=1,2,3$. באמצעות אלה ניתן לקבל 4 משוואות מצומדות אשר פתרונן מוביל לביטוי עבור יחס הדיספרסיה של המערכת:

$e^{-4k_{1}a}={\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 \infty$ מתקבל הביטוי ליחס הדיספרסיה שך מערכת דו שכבתית.

משוואה זו לא ניתנת לפתרון אנליטי, לכן על מנת לפשט את המשוואה נניח ששתי השכבות הדיאלקטריות זהות כך ש $\varepsilon _{2}=\varepsilon _{3}$ ו- $k_{2}=k_{3}$. כעת ניתן לפצל את יחס הדיספרציה שהתקבל לזוג משוואות, אחת עבור הפתרונות הזוגיים והשנייה עבור הפתרונות האי-זוגיים, בהתאמה:

$\tanh k_{1}a=-{\frac {k_{2}\varepsilon _{1}}{k_{1}\varepsilon _{2}}}$

$\tanh k_{1}a=-{\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 _{0}d}}\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_{2DEP}^{2}-{\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

הערך באדיבות ויקיפדיה העברית, קרדיט,
רשימת התורמים
רישיון cc-by-sa 3.0

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