תהודה פלסמונית מקומית

בפיזיקה של מצב מוצק, תהודה פלסמונית מקומית (באנגלית: Localized Surface Plasmon Resonance או בקיצור LSPR) היא עירור פלסמוני שיכול להתרחש בחלקיקי מתכת, בדרך כלל בתחום הננומטרי. במצב זה, חלקיק מתכת מוקרן בתדר המתאים לתנודות של אלקטרונים חופשיים שבמתכת. לתנודות אלו קוראים תנודות פלסמוניות והן יכולות לבוא לידי ביטוי בכמה תצורות מרחביות. תחת הקירוב הקוואזי-סטטי, חלקיק בתהודה פלסמונית ניתן לתיאור פשוט כדיפול מתנדנד. כל אימת שתדר השדה האלקטרומגנטי החיצוני נמוך מתדר הפלזמה של המתכת, האלקטרונים החופשיים יעקבו אחר תדר זה ויבצעו את התנודות הפלסמוניות. בקירוב ראשון, ניתן לתאר את האינטראקציה של השדה האלקטרומגנטי עם חלקיקי מתכת בעזרת מודל הפלאסמה. על פי מודל זה, ברקע השדה האלקטרומגנטי נעים (ונדים) אלקטרונים חופשיים של המתכת כנגד הכוחות המחזירים מגרעיני היונים המרכיבים את הגביש. במודל זה אין התייחסות ישירה לפוטנציאל ההרמוני של הגביש ולאינטראקציות בין האלקטרונים החופשיים עצמם, וההתייחסות היחידה למבנה הפסים של החומר באה לידי ביטוי במסה האפקטיבית האופטית של האלקטרון. כאשר ממדי המבנה הפלסמוני גדולים יחסית למרחק הריסון (המרחק הממוצע בין התנגשויות של האלקטרון בגרעין), לפי משוואת התנועה החד־ממדית של אלקטרון בודד מתקיים- m e ∗ x → ¨ = − m e ∗ Γ x → ˙ − | q e | E → 0 e − i ω t ⇒ x → = | q e | m e ∗ ( ω 2 + i Γ ω ) E → 0 e − i ω t {\displaystyle m_{e}^{*}{\ddot {\vec {x}}}=-m_{e}^{*}\Gamma {\dot {\vec {x}}}-\left|{{q}_{e}}\right|{{\vec {E}}_{0}}{{e}^{-i\omega t}}\Rightarrow {\vec {x}}={\frac {\left|{{q}_{e}}\right|}{m_{e}^{*}\left({{\omega }^{2}}+i\Gamma \omega \right)}}{{\vec {E}}_{0}}{{e}^{-i\omega t}}} כאשר m e ∗ {\displaystyle m_{e}^{*}} היא המסה האפקטיבית האופטית של האלקטרון, Γ {\displaystyle \Gamma } מאפיינת את תדירות ההתנגשויות של האלקטרון בגרעיני הגביש, q e {\displaystyle {{q}_{e}}} הוא מטען האלקטרון ו- E → 0 {\displaystyle {{\vec {E}}_{0}}} מייצג את אמפליטודת השדה החשמלי השייך לשדה האלקטרומגנטי המוקרן. ממילא, כיוון שמומנט הדיפול הכולל נתון על ידי P = − n | q e | x → {\displaystyle P=-n\left|{{q}_{e}}\right|{\vec {x}}} כאשר n {\displaystyle n} הוא מספר האלקטרונים ליחידת נפח בחלקיק, אזי כאשר מתחשבים בקשר D = ε 0 E + P = ε 0 ε r E {\displaystyle D={{\varepsilon }_{0}}E+P={{\varepsilon }_{0}}{{\varepsilon }_{r}}E} ניתן לחלץ את הפונקציה הדיאלקטרית של המתכת- ε r = 1 + ( ω 2 − ω p 2 ) τ 2 1 + ω 2 τ 2 + i ω p 2 τ ω ( 1 + ω 2 τ 2 ) {\displaystyle {{\varepsilon }_{r}}={\frac {1+\left({{\omega }^{2}}-\omega _{p}^{2}\right){{\tau }^{2}}}{1+{{\omega }^{2}}{{\tau }^{2}}}}+i{\frac {\omega _{p}^{2}\tau }{\omega \left(1+{{\omega }^{2}}{{\tau }^{2}}\right)}}} כאשר ω p = n q e 2 ε 0 m e ∗ {\displaystyle {{\omega }_{p}}={\sqrt {\frac {nq_{e}^{2}}{{{\varepsilon }_{0}}m_{e}^{*}}}}} הוא תדר הפלאסמה ו- τ = 1 Γ {\displaystyle \tau ={\frac {1}{\Gamma }}} הוא הזמן הממוצע בין התנגשויות של האלקטרון בגרעיני היונים המרכיבים את המבנה הפלסמוני. התלות ε r ( ω ) {\displaystyle {{\varepsilon }_{r}}\left(\omega \right)} מהווה את הבסיס לקיומו של מנגנון התהודה הפלסמונית בחומר. התיאור הפיזיקלי של מנגנון התהודה הפלסמונית תלוי בגאומטריה של החלקיק ובממדיו. האינטראקציה של המבנה הפלסמוני עם השדה האלקטרומגנטי באה לידי ביטוי בפיזור האלקטרומגנטי המתקבל מחלקיק נתון. עבור חלקיקים כדוריים שאינם בהכרח קטנים מאורך הגל המוקרן, פתרון מלא ויחסית מורכב של בעיית הפיזור ניתן באמצעות תאוריית מיי. תחת הקירוב הקוואזי-סטטי ניתן לפשט משמעותית את הפתרון המתקבל על פי תאוריית מיי עבור חלקיקים כדוריים. עבור חלקיקים שאינם כדוריים, תיאור הפיזור והבליעה של פוטונים בהשפעת החלקיק נעשה מורכב בהרבה וברוב המקרים במצב זה לא קיים פתרון אנליטי לבעיית הפיזור. אולם, בגאומטריות מסוימות, תחת הקירוב הקוואזי-סטטי ניתן יהיה לקבל ביטויים יחסית פשוטים לחתכי הפעולה לבליעה ולפיזור אשר בקירוב זה מקסימליים בתדרי התהודה הפלסמונית המקומית של החלקיק הנתון. כאשר חלקיק הנתון בתווך דיאלקטרי בעל קבוע ε m {\displaystyle {{\varepsilon }_{m}}} מוקרן על ידי פוטונים באורך גל λ {\displaystyle \lambda } הגדול ביחס לממדי החלקיק, ניתן לתאר את האינטראקציה בין החלקיק לפוטון כאינטראקציה בין מבנה פלסמוני לשדה אלקטרומגנטי חיצוני תחת הקירוב הקוואזי-סטטי. במקרה של חלקיקי מתכת כדוריים ברדיוס a {\displaystyle a} , משמעות הקירוב a ≪ λ {\displaystyle a\ll \lambda } היא קיומו של תדר רזוננס פלסמוני שאיננו תלוי בממדי החלקיק (בפועל, תדר הרזוננס הפלסמוני גדל ככל שהיחס בין רדיוס החלקיק לאורך הגל גדל, כתוצאה ממנגנון שנקרא "Red shift"). תחת קירוב זה ניתן להניח כי הפרשי הפאזה של השדה החשמלי בסביבת החלקיק זניחים וכי השדה החשמלי באזור זה אחיד. במקרה כזה, פתרון של משוואות מקסוול מנבא פוטנציאל אלקטרוסטטי המתנהג כדיפול, עם פולריזביליות: α = 4 π a 3 ε r − ε m ε r + 2 ε m {\displaystyle \alpha =4\pi {{a}^{3}}{\frac {{{\varepsilon }_{r}}-{{\varepsilon }_{m}}}{{{\varepsilon }_{r}}+2{{\varepsilon }_{m}}}}} כאשר הקשר בין מומנט הדיפול p → {\displaystyle {\vec {p}}} לשדה החשמלי החיצוני E → 0 {\displaystyle {{\vec {E}}_{0}}} (אשר בקירוב אחיד) הוא p → = ε 0 ε m α E → 0 {\displaystyle {\vec {p}}={{\varepsilon }_{0}}{{\varepsilon }_{m}}\alpha {{\vec {E}}_{0}}} . ממילא קיים ω {\displaystyle \omega } מסוים עבורו α {\displaystyle \alpha } מקסימלי. מצב זה נקרא הרזוננס הפלסמוני המקומי של החלקיק והתנאי לקיומו במקרה של חלקיקים כדוריים הוא Re ⁡ { ε r ( ω ) } = − 2 ε m {\displaystyle \operatorname {Re} \left\{{{\varepsilon }_{r}}\left(\omega \right)\right\}=-2{{\varepsilon }_{m}}} . תנאי זה נקרא תנאי פרליך ומהתבוננות בביטוי עבור הפולריזביליות ניתן להראות כי: Re ⁡ { ε r ( ω L S P R ) } = − 2 ε m ⇒ ω L S P R = ω p 1 1 + 2 ε m − 1 ω p 2 τ 2 ≈ ω p 1 + 2 ε m {\displaystyle \operatorname {Re} \left\{{{\varepsilon }_{r}}\left({{\omega }_{LSPR}}\right)\right\}=-2{{\varepsilon }_{m}}\Rightarrow {{\omega }_{LSPR}}={{\omega }_{p}}{\sqrt {{\frac {1}{1+2{{\varepsilon }_{m}}}}-{\frac {1}{\omega _{p}^{2}{{\tau }^{2}}}}}}\approx {\frac {{\omega }_{p}}{\sqrt {1+2{{\varepsilon }_{m}}}}}} כאשר הקירוב האחרון נעשה מתוך הנחה כי ω p 2 τ 2 ≫ 1 + 2 ε m {\displaystyle \omega _{p}^{2}{{\tau }^{2}}\gg 1+2{{\varepsilon }_{m}}} . הנחה זו בדרך כלל מתקיימת במתכות הנפוצות. בכל מקרה, פיתוח זה מותנה בכך ש- ω L S P R < ω c {\displaystyle {{\omega }_{LSPR}}<{{\omega }_{c}}} כאשר ℏ ω c {\displaystyle \hbar {{\omega }_{c}}} היא אנרגיית הסף המקסימלית (של הפוטון המוקרן) ממנה יש לקחת בחשבון את השפעתם הישירה של מעברי האלקטרונים (בין רמות אנרגיה בגביש המתכת) על הפונקציה הדיאלקטרית. במובן מסוים, המודל הנ"ל ניתן להכללה עבור חלקיקים כלליים בעלי אורך, רוחב ועובי נתונים. באופן מדויק ניתן לומר כי עבור חלקיק כללי בדמות אליפסואיד בעל אורך a 1 {\displaystyle {{a}_{1}}} , רוחב a 2 {\displaystyle {{a}_{2}}} ועובי a 3 {\displaystyle {{a}_{3}}} כלשהם- קיים פתרון אנליטי מלא למשוואות מקסוול לפיו הפולריזביליות הדיפולית בקירוב הקוואזי-סטטי יכולה לקבל שלושה ערכים שונים (בהנחה שהשדה האלקטרומגנטי מקוטב בכיוון אחד מצירי האליפסואיד): α i = 4 π 3 a 1 a 2 a 3 ε r − ε m ε m + ( ε r − ε m ) L i , L i = a 1 a 2 a 3 2 ∫ 0 ∞ d φ ( φ + a i 2 ) ( φ + a 1 2 ) ( φ + a 2 2 ) ( φ + a 3 2 ) , i = 1 , 2 , 3 {\displaystyle {\begin{matrix}{{\alpha }_{i}}={\frac {4\pi }{3}}{{a}_{1}}{{a}_{2}}{{a}_{3}}{\frac {{{\varepsilon }_{r}}-{{\varepsilon }_{m}}}{{{\varepsilon }_{m}}+\left({{\varepsilon }_{r}}-{{\varepsilon }_{m}}\right){{L}_{i}}}},&{{L}_{i}}={\frac {{{a}_{1}}{{a}_{2}}{{a}_{3}}}{2}}\int \limits _{0}^{\infty }{{\frac {d\varphi }{\left(\varphi +a_{i}^{2}\right){\sqrt {\left(\varphi +a_{1}^{2}\right)\left(\varphi +a_{2}^{2}\right)\left(\varphi +a_{3}^{2}\right)}}}},}&i=1,2,3\\\end{matrix}}} כאשר α 1 {\displaystyle {{\alpha }_{1}}} , למשל, היא הפולריזביליות הדיפולית של החלקיק המתנהג כדיפול בכיוון ציר x {\displaystyle x} בהשפעת שדה אלקטרומגנטי חיצוני המקוטב בהתאם. ממילא קיימים שלושה מצבי תהודה מקומיים, והתנאי המקביל לתנאי פרליך במקרה זה הוא: Re ⁡ { ε r ( ω L S P R ( i ) ) } = − ( 1 L i − 1 ) ε m ⇒ ω L S P R ( i ) = ω p 1 1 + ( 1 L i − 1 ) ε m − 1 ω p 2 τ 2 ≈ ω p 1 + ( 1 L i − 1 ) ε m {\displaystyle {\begin{matrix}\operatorname {Re} \left\{{{\varepsilon }_{r}}\left(\omega _{LSPR}^{\left(i\right)}\right)\right\}=-\left({\frac {1}{{L}_{i}}}-1\right){{\varepsilon }_{m}}&\Rightarrow &\omega _{LSPR}^{\left(i\right)}={{\omega }_{p}}{\sqrt {{\frac {1}{1+\left({\frac {1}{{L}_{i}}}-1\right){{\varepsilon }_{m}}}}-{\frac {1}{\omega _{p}^{2}{{\tau }^{2}}}}}}\approx {\frac {{\omega }_{p}}{\sqrt {1+\left({\frac {1}{{L}_{i}}}-1\right){{\varepsilon }_{m}}}}}\\\end{matrix}}} כאשר ω L S P R ( i ) {\displaystyle \omega _{LSPR}^{\left(i\right)}} מייצג את תדר הרזוננס הפלסמוני המקומי המתקבל בהשפעת שדה חשמלי המתנדנד בכיוון ציר האליפסואיד ה- i {\displaystyle i} . הקירוב הקוואזי-סטטי מאפשר לתאר תהודה פלסמונית של ננו חלקיקים בצורה יחסית מדויקת עבור חלקיקים שממדיהם קטנים מ- 100 n m {\displaystyle {{100}_{nm}}} . אולם, עבור חלקיקים גדולים יותר יש צורך בניתוח מדויק יותר התקף גם ללא הנחת הקוואזי-סטטיות. ניתוח אלקטרודינמי אנליטי לבעיה זו נעשה מורכב בהרבה כאשר הקירוב הקוואזי-סטטי איננו מוצדק, מאחר שלמעט במקרה של חלקיקים כדוריים, לא קיים פתרון אנליטי מלא לבעיית הפיזור מחלקיק שממדיו אינם קטנים ביחס לאורך הגל המוקרן. הפתרון עבור המקרה הכללי של חלקיקים כדוריים הוצג בשלב יחסית מוקדם. בשנת 1908 פרסם הפיזיקאי הגרמני גוסטב מיי את מאמרו "תרומות לאופטיקה של סביבה טורבידית, בפרט עבור תמיסת קולואידים של מתכות" (“Contributions to the optics of turbid media, particularly of colloidal metal solutions”). במאמר זה הציג מיי פתרון אנליטי מלא לבעיית הפיזור של גלים אלקטרומגנטיים על ידי מבנים כדוריים. פיזור זה נקרא "פיזור מיי", והוא שימושי בעיקר במקרים בהם אורך הגל המפוזר הוא בסדר הגודל של ממדי החלקיק המפזר. ניתן להציג את הפולריזביליות הדיפולית המתקבלת בפיזור מיי בקירוב ראשון על ידי מנת פולינומים: α = 40 π a 3 − 4 π ( ε + ε m ) a 3 x 2 + O ( x 4 ) 10 ( ε + 2 ε m ) ε − ε m − ( ε + 10 ε m ) x 2 − i 160 ε m 1.5 3 x 3 + O ( x 4 ) {\displaystyle \alpha ={\frac {40\pi {{a}^{3}}-4\pi \left(\varepsilon +{{\varepsilon }_{m}}\right){{a}^{3}}{{x}^{2}}+O\left({{x}^{4}}\right)}{{\frac {10\left(\varepsilon +2{{\varepsilon }_{m}}\right)}{\varepsilon -{{\varepsilon }_{m}}}}-\left(\varepsilon +10{{\varepsilon }_{m}}\right){{x}^{2}}-i{\frac {160\varepsilon _{m}^{1.5}}{3}}{{x}^{3}}+O\left({{x}^{4}}\right)}}} כאשר x = π a λ 0 {\displaystyle x={\frac {\pi a}{{\lambda }_{0}}}} ו- λ 0 {\displaystyle {{\lambda }_{0}}} הוא אורך הגל המפוזר כאשר הוא מוקרן בריק. הביטוי לעיל מנבא היסט ספקטרלי של תדר הרזוננס הפלסמוני המקומי, ביחס למיקומו בקירוב הקוואזי-סטטי. היסט זה הוא לכיוון תדרי האינפרא-אדום, ולכן נקרא "Red shift" (היסט אדום). גודל ה- "Red shift" נקבע בהתאם לממדי החלקיק. ככל שהחלקיק גדול יותר, כך ה- "Red shift" גדול יותר. משמעות היסט זה היא קביעתו של אורך הגל עם הנטייה הגבוהה ביותר להתפזר בהשפעת החלקיק. ברמה המאקרוסקופית, המשמעות היא שצבעם של חלקיקי המתכת נקבע לפי גודלם. דוגמה לכך היא צבעו של קולואיד זהב הנקבע לפי גודלם של ננו חלקיקי הזהב אותם הוא מכיל. לחלקיקי מתכת כדוריים יש בדרך כלל תדר דומיננטי יחיד של תהודה פלסמונית מקומית. מצב זה מתקיים כאשר האלקטרונים החופשיים בחלקיק מדמים תצורה מרחבית של דיפול. במקרה של חלקיקים שאינם כדוריים, יכולים להתקבל כמה מצבי תהודה פלסמונית עבור תצורות מרחביות שונות של הדיפול. ישנה תופעה פחות ניכרת, לפיה ייתכנו מצבים נוספים של תהודה פלסמונית מקומית- תהודה פלסמונית של מולטיפולים מסדרים גבוהים. עבור חלקיקי מתכת כדוריים, מצבים אלו באים לידי ביטוי מתמטי באיברים מסדר O ( x 4 ) {\displaystyle O\left({{x}^{4}}\right)} ומעלה אשר הוזנחו בביטוי שהוצג לעיל עבור הפולריזביליות המתקבלת בפיזור מיי. ניתן במקרה כזה להראות תחת הקירוב הקוואזי-סטטי כי לכל חלקיקי מתכת אינפיניטסימלי יש אינסוף תדרי רזוננס פלסמוני מקומיים כאשר התדר (הזוויתי) של הרזוננס הפלסמוני מסדר n {\displaystyle n} נתון על ידי: ω L S P R [ n ] = ω p 1 + ( 1 + 1 n ) ε m {\displaystyle {{\omega }_{LSPR}}\left[n\right]={\frac {{\omega }_{p}}{\sqrt {1+\left(1+{\frac {1}{n}}\right){{\varepsilon }_{m}}}}}} כאשר עבור n ≥ 2 {\displaystyle n\geq 2} התופעה פחות ניכרת ולרוב זניחה. ממילא עבור n = 1 {\displaystyle n=1} מתכנסים לביטוי המוכר עבור ω L S P R {\displaystyle {{\omega }_{LSPR}}} של חלקיקים כדוריים בקירוב הקוואזי-סטטי. פלזמון פיזור דיפול ננו-חלקיק ננו חלקיקי זהב פולריזביליות נקודה קוונטית Craig F. Bohren, Donald R. Huffman (1983) Absorption and Scattering of Light by Small Particles, Wiley Stefan Alexander Maier (2007) Plasmonics, Fundamental and Applications, Springer David K. Cheng (1989) Field and Wave Electromagnetics, Addison-Wesley Uwe Kreibig, Michael Vollmer (1995) Optical Properties of Metal Clusters, Springer Tchiya Zar, Racheli Ron, Omer Shavit, Alon Krause, David Gachet, Adi Salomon (2022) Nanoscopy of Aluminum Plasmonic Cavities by Cathodoluminescence and Second Harmonic Generation, Wiley מדיה וקבצים בנושא תהודה פלסמונית מקומית בוויקישיתוף פרהן רנה, פיזור אלקטרומגנטי, ‏2007 (באנגלית) Craig F. Bohren, Donald R. Huffman (1983) Absorption and Scattering of Light by Small Particles, Wiley Stefan Alexander Maier (2007) Plasmonics, Fundamental and Applications, Springer Uwe Kreibig, Michael Vollmer (1995) Optical Properties of Metal Clusters, Springer