שיטת ANN-Z היברידית למידול משטחים חכמים הניתנים להגדרה מחדש המבוססים על ננו-צינוריות פחמן עבור היגוי אלומה בתדר טרה הרץ

Mayameen S. Kadhim; Humam Hussein; Taha A. Elwi

doi:10.3791/70498

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

מאמר שיטה

שיטת ANN-Z היברידית למידול משטחים חכמים הניתנים להגדרה מחדש המבוססים על ננו-צינוריות פחמן עבור היגוי אלומה בתדר טרה הרץ

53 צפיות

⸱

DOI:

10.3791/70498

⸱

28 באוגוסט 2026

Mayameen S. Kadhim¹ , Humam Hussein²^,³ , Taha A. Elwi⁴^,⁵

¹Technical Engineering College, Al-Bayan University, ²College of Computer Science and Mathematics, Tikrit University, ³Avicenna E-Learning Center, Tikrit University, ⁴Department of Automation and Artificial Intelligence Engineering, College of Information Engineering at Al-Nahrain University, ⁵International Applied and Theoretical Research Center (IATRC)

במאמר זה

סיכום

פרוטוקול זה מציג שיטה היברידית המשלבת רשת עצבית מלאכותית והתמרת Z למודל אלקטרומגנטי מדויק של משטחים תבוניים אופטיים ניתנים לתכנון המבוססים על ננו-צינוריות פחמן בעלות דופן אחת, הפועלים בפס התדרים של טרה-הרץ (0.5-30 THz) עבור יישומים של תקשורת אלחוטית בדור שישי (6G), תוך השגת האצה חישובית של 180× עם יכולת כיוונון של מופע ההחזרה העוברת את 310° וטווח היגוי של קרן של ±45°.

תקציר

משטחים אינטליגנטיים בעלי תכונות אופטיות הניתנות להגדרה מחדש המבוססים על ננו-צינוריות פחמן דופן אחת מציעים פתרונות מבטיחים להכוונת אלומות טרה-הרץ ומניפולציה של גלי פוטוניקה במערכות אלחוטיות עתידיות של דור 6 (6G). עם זאת, מידול מדויק של מבנים אלו נותר מאתגר בשל אפקטים של הובלה קוונטית, השראות קינטית והתנהגות אקסיטונית רב-תהודה בטווחים רחבים של תדרים. פרוטוקול זה מתאר מסגרת חישובית היברידית המשלבת מידול מוליכות קוונטית באמצעות פורמליזם קובו (Kubo formalism), החלקת נתונים מבוססת רגרסיה פולינומית וניתוח פונקציית תמסורת במישור Z לאפיון מדויק של תאי יחידה של משטחים אינטליגנטיים בעלי תכונות אופטיות הניתנות להגדרה מחדש מננו-צינוריות פחמן דופן אחת. השיטה מתחילה בתכנון של רזונטורי ננו-פסים מצטלבים של ננו-צינוריות פחמן דופן אחת על מצע קוורץ בעלי כיריליות (10,5) (קוטר 0.60 nm, פער אנרגיה 1.762 eV), ולאחריה סימולציה אלקטרומגנטית של גל מלא ב-CST Microwave Studio בטווח תדרים של 0.5-30 THz. מודל רגרסיה פולינומית מסדר 8 מעבד את פרמטרי S שהופקו כדי להסיר תנודות מספריות ולנבא תגובות אלקטרומגנטיות מוחלקות. פונקציית תמסורת בדידה H(z) עם מונה מסדר 6 ומכנה מסדר 7 מותאמת לאחר מכן באמצעות אופטימיזציית ריבועי פחותים עם פירוק QR, מה שמאפשר ניתוח יציבות של קטבים ואפסים ואימות פסיביות. הפרוטוקול משלב בנוסף כוונון מוליכות קוונטית באמצעות מודולציה של פוטנציאל כימי לאופטימיזציה של הכוונת אלומות. תוצאות מייצגות מדגימות יכולת כוונון של מופע השתקפות העולה על 310°, הגברת בליעה עד 92.3%, טווח הכוונת אלומות של ±45° עם רמות סייד-לוב נמוכות מ-12- dB, והאצה חישובית של 180× בהשוואה לשיטות אופטימיזציה קונבנציונליות של גל מלא. הרגרסיה הפולינומית השיגה שגיאת Root Mean Square במבחן של 0.0688 עם מקדם R² של 0.994, בעוד שהתאמת H(z) השיגה שגיאת Root Mean Square של 0.89 dB. ניתוח היציבות אישר כי כל הקטבים נמצאים בתוך מעגל היחידה. פרוטוקול זה מספק נתיב יעיל וניתן לשחזור לתכנון מטא-משטחים פוטוניים הניתנים לתכנות ומערכות תקשורת טרה-הרץ אינטליגנטיות עבור 6G ומעבר לכך.

מבוא

ההתפתחות המהירה לעבר מערכות אלחוטיות מהדור השישי (6G) האיצה את חקירת פסיי תדרים של טרה-הרץ (THz) ותדרים אופטיים כדי להשיג קצבי נתונים גבוהים במיוחד העולים על 1 Tbps, חישוש חכם, עיצוב אלומה הולוגרפי והנדסת חזית גל אדפטיבית¹^,²^,³. פס התדרים של הטרה-הרץ (0.1-30 THz) מציע רוחב פס רב, אך סובל מניחות חמורה במרחב החופשי (בסביבות 20 dB·km^-1 בתדר של 1 THz), בליעה מולקולרית אטמוספרית על ידי אדי מים ב-557 GHz, 752 GHz, 988 GHz ו-1.13 THz, ורגישות קיצונית לחסימות על ידי חלקיקים אטמוספריים וגשם⁴.

משטחים אינטליגנטיים ניתנים להגדרה מחדש הופיעו כטכנולוגיה מאפשרת המסוגלת להשפיע באופן דינמי על חזיתות גלים אלקטרומגנטיים באמצעות החזרה, בליעה, בקרת פאזה ומיקוד גלים הניתנים לתכנות⁵^,⁶^,⁷. כאשר הם מורחבים לתדרים אופטיים, משטחים אלו מציעים שליטה חסרת תקדים על התקדמות האור, ומאפשרים יישומים כגון LiDAR, תקשורת אופטית במרחב חופשי, תצוגות הולוגרפיות והכוון אלומה אדפטיבי⁸.

מבני משטחים אינטליגנטיים ר可-קונפיגורביליים (reconfigurable intelligent surface) קונבנציונליים המבוססים על רזונטורים מתכתיים נתקלים במגבלות משמעותיות בתדרים של טרה-הרץ ואופטיקה, בשל הגברת הפסדים אומיים הנובעים ממוליכות דרוד (Drude conductivity), שיכוך פלסמוני, מגבלות ייצור עבור גדלי מאפיינים קטנים מ-100 nm, ויכולת כיוונון מוגבלת של תכונות דיאלקטריות קבועות⁹^,¹⁰^,¹¹. מטא-משטחים המבוססים על ננו-חומרים, ובמיוחד אלה המשתמשים בננו-צינורות פחמן בעלי דופן אחת (SWCNTs), מציעים חלופות מבטיחות בשל תכונותיהם הקוונטיות-אלקטרומגנטיות יוצאי הדופן: ממדי ננוסקאלה עם קטרים של 0.7-3 nm, מוליכות קוונטית ניתנת לכיוונון הנשלטת על ידי פוטנציאל כימי, ניידות נשאיים יוצאת דופן העולה על 100,000 cm²·(V·s^-1), מוליכות תרמית גבוהה של כ-3000 W·(m·K^-1), ורזוננסים אקסיטוניים חזקים בתחומי האור הנראה והאינפרא-אדום הקרוב¹²^,¹³^,¹⁴^,¹⁵. SWCNTs מראים אינטראקציות אלקטרומגנטיות חזקות במשטרי טרה-הרץ ופוטוניקה באמצעות צימוד אקסיטון-פוטון, המאפשר יצירת רזונטורים אופטיים קומפקטיים עם תגובות אלקטרומגנטיות ניתנות לתכנות¹⁶. תכונות אלו הופכות את ה-SWCNTs למועמדים אטרקטיביים עבור משטחים אינטליגנטיים ר可-קונפיגורביליים אופטיים, שבהם גודל קומפקטי, תגובה ניתנת לכיוונון ומניפולציה יעילה של גלים הם דרישות תכנון חשובות.

למרות יתרונות אלו, מידול מדויק של משטחים אינטליגנטיים ניתנים להגדרה מחדש (RIS) המבוססים על SWCNT נותר מאתגר עקב השפעות של תובעת קוונטית הדורשות גישות של פונקציית גרין לא-שיווי משקל, השראות קינטית השולטת בתדרים מעל 100 GHz, התנהגות אקסיטונית מרובת-תהודה עם אנרגיות קשירה של עד 0.4 eV, תופעות דיספרסיה לא-ליניאריות הנובעות מצימוד אלקטרון-פונון, תכונות אופטיות התלויות בכיראליות הדורשות הגדרה של אינדקסים כיראליים (n, m), ומוליכות התלויה בטמפרטורה הדורשת מידול תרמי עקבי¹⁷^,¹⁸^,¹⁹^,²⁰. גישות מסורתיות של מעגלים שווי ערך ושיטות התאמה אלקטרומגנטית קונבנציונליות, כגון קירוב רציונלי (rational approximation) והתאמה וקטורית (vector fitting), נכשלות לעיתים קרובות בייצוג מדויק של אינטראקציות מורכבות אלו על פני טווחים רחבים של תדרים המכסים 0.5-30 THz, רוחב פס של 60×²¹^,²²^,²³. מגבלות אלו עלולות להפחית את דיוק המידול ולהגביר את המורכבות החישובית במהלך תכנון משטחים אינטליגנטיים ניתנים להגדרה מחדש בתדר טרה-הרץ, דבר המדגיש את הצורך בגישות מידול יעילות השומרות על פרשנות פיזיקלית תוך לכידה מדויקת של התנהגות אלקטרומגנטית רחבת פס.

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

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

פרוטוקול

1. תכנון תא יחידה של יחידת RIS אופטית מבוססת SWCNT

בחירת כיראליות של SWCNT
1. בחר כיראליות SWCNT של (10,5) בהתבסס על ניתוח מוליכות קוונטית.
2. חשב את קוטר הננו-צינור באמצעות הנוסחה המצוינת להלן
  
  כאשר a_cc=0.142 nm הוא אורך הקשר פחמן-פחמן. כיראליות של (10,5) מניבה קוטר של 0.60 nm ופער אנרגיה (bandgap) של 1.762 eV, אופטימליים לפעולה בתדר טרה-הרץ.
3. חשב את הזווית הכיראלית באמצעות
הגדרת גיאומטרית תא היחידה
1. תכנן רזונטורי ננו-פסים של SWCNT מצולבים על מצע קוורץ עם הפרמטרים הבאים: מיתנות מצע של 3.8, עובי מצע של 500 nm, ממדי תא יחידה של 1.5 × 1.5 µm², רוחב פס SWCNT של 50 nm, אורך פס SWCNT של 700 nm, ומרווח בין הפסים של 120 nm. ממדים אלו מבטיחים פעולה בתת-אורך גל לאורך פס התדרים 0.5-30 THz.
הגדרת סימולציה ב-CST Microwave Studio
1. הפעל את CST Microwave Studio וצור פרויקט חדש באמצעות ה-Microwave and RF frequency domain solver.
2. בנה את גיאומטריית תא היחידה באמצעות הפרמטרים שהוגדרו. החל תנאי שפה של תא יחידה בכיווני x ו-y עם גבולות מחזוריים. הגדר פורטי Floquet לאורך כיוון z לצורך עירור גל מישורי.
3. הגדר את טווח סריקת התדרים מ-0.5 THz עד 30 THz עם גודל צעד של 0.05 THz, ליצירת 590 נקודות תדר. הגדר את ה-time-domain solver עם עידון רשת אדפטיבי (adaptive mesh refinement) והגדר יעד התכנסות של -40 dB.
4. הרץ את הסימולציה האלקטרומגנטית של הגל המלא. ייצא את מקדם ההחזרה המורכב S₁₁, מקדם המעבר S₂₁, פאזת ההחזרה וספקטרומי הבליעה כקבצי CSV.

2. מידול מוליכות קוונטית באמצעות פורמליזם קובו

אתחול קבועים פיזיקליים
1. אתחול קבועים פיזיקליים: מטען אלמנטרי e=1.602×10⁻¹⁹C, קבוע פלאנק רדוקציוני ħ=1.0546×10⁻³⁴J·s, קבוע בולצמן k_B=1.3806×10⁻²³J·K^-1, טמפרטורה T=300K, מהירות פרמי v_F=8×10⁵m·s^-1, ועכבת חלל ריק η₀=377Ω.
חישוב מוליכות תוך-פסית (Intra-band)
1. חישוב המוליכות התוך-פסית עבור כל נקודת תדר באמצעות:
  
  כאשר τ=0.5ps הוא זמן הרפיית הנושאי מטען ו-E_F היא אנרגיית פרמי (המשתנה בין 0.1e.V ל-0.4eV).
חישוב מוליכות בין-פסית (Inter-band)
1. חישוב המוליכות הבין-פסית באמצעות:
2. חישוב מוליכות פני השטח הכוללת כ- σ_total(ω) = σ_intra(ω) + σ_inter(ω).
חישוב תגובה אופטית
1. חישוב השהיית הפאזה של התשתית עבור קוורץ (ε_r=3.8, עובי=500nm) באמצעות φ_sub=2πf(n_sub)d_sub/c, כאשר n_sub=.
2. חישוב מקדם השתקפות מהמוליכות באמצעות
3. חישוב מקדם השתקפות כ- R(ω)=|S₁₁(ω)|². חישוב בליעה כ- A(ω)=1-R(ω)-T(ω).

3. הפקת נתונים ועיבוד מקדים של CST

יצירת נתוני CST מדומים
1. ייצור תוצאות סימולציית גל מלא (full-wave) של CST מדומות על ידי הוספת תהודות ריאליסטיות ל-S₁₁ הקוונטי.
2. הגדרת חמש תהודות: אקסיטון E₁₁ ב-2.8 THz (עוצמה -15.2 dB, מופע 45°), אקסיטון E₂₂ ב-5.6 THz (-22.8 dB, -120°), תהודה פלסמונית ב-12.4 THz (-8.5 dB, 60°), מוד חלל (cavity mode) ב-18.9 THz (-12.1 dB, -30°), ותהודה מסייעת פונונית ב-24.7 THz (-6.8 dB, 15°).
3. הוספת רעש נומרי עם רמת רעש 0.03 (SNR=30.5dB) כדי לדמות תנודות נומריות של CST.
עיבוד מקדמי של נתונים
1. נורמליזציה של ציר התדר לטווח [0,1] באמצעות
  $נוסחת תדר מנורמל $z_{f_{norm}}=\frac{f-f_{min}}{f_{max}-f_{min}}$ המציגה את המשוואה.$
2. הפקה של החלקים הממשיים והמדומים של S₁₁ והעוצמה של S₂₁.
3. פיצול מערך הנתונים לאימון (80%), תיקוף (10%) ובדיקה (10%) באמצעות פרמוטציה אקראית. נורמליזציה של הקלטים לממוצע אפס ושונות יחידה.

4. רגרסיה פולינומיאלית להחלקת נתונים

התאם 8^th התאם פולינום מסדר 8 לחלק הממשי של S₁₁ באמצעות polyfit: P_real = polyfit(X_train, Y_train(:,1), 8).^th סדר פולינום עבור החלק המדומה: P_imag = polyfit(X_train, Y_train(:,2), 8). התאם פולינום מסדר 8^th סדר פולינום עבור גודל S₂₁: P_mag = polyfit(X_train, Y_train(:,3), 8). הפק תחזיות מוחלקות באמצעות polyval על פני טווח התדרים המלא. חשב S₁₁ מוחלק כ- S₁₁_ann = S11_real_ann + i × S11_imag_ann.

5. התאמת פונקציית תמסורת במישור Z

קבע את תדירות הדגימה F_s=60THz (פי 2 מהתדירות המקסימלית לפי קריטריון ניקוויסט). מפה את התדרים למישור Z באמצעות z=e^{i2πf/Fs}.
הגדר את פונקציית התמסיירת הבדידה במישור Z באופן הבא:

כאשר סדר המונה n=6 וסדר המכנה m=7.
עבור אופטימיזציה בריבועים הפחותים, הגדר משוואות ליניאריות A × x= B, כאשר A מכיל את איברי המונה והמכנה ו-B מכיל את S₁₁ שנחזה על ידי רשת נוירונים מלאכותית (ANN). פתור באמצעות פירוק QR ליציבות נומרית:
ניתוח קטבים-אפסים ויציבות
1. חלץ את הקטבים על ידי פתרון פולינום המכנה. הבטח יציבות על ידי השלכת קטבים שערכם |pole| ≥ 1 אל תוך מעגל היחידה באמצעות poles(p) = poles(p) / (|poles(p)| + 0.1).
2. שחזר את המכנה מהקטבים המיוצבים. חלץ את האפסים על ידי פתרון פולינום המונה. העריך את H(z) על גבי רשת תדרים וחשב את שגיאת השורש הממוצע של הריבוע (RMSE).

6. אופטימיזציה של היגוי האלומה

גורם המערך הנורמלי עבור מערך פאזי ליניארי נתון על ידי:

כאשר I_n = 1 עבור עירור אחיד, k = 2π/λ, d = מרווח של λ/2, ו-β_n הוא השיתה הפרוגרסיבית. עבור זוויות היגוי למטרה של ±45°, ±30°, ו-0°, חשב את השיתה בין אלמנטים סמוכים כ-.EQUATION
1. חשב את דפוסי השדה הרחוק והערך את רמות האונות הצדדיות ואת רוחב האלומה בחצי עוצמה.

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

תוצאות

בחירת כיראליות של SWCNT
הפרוטוקול המתואר יושם עבור תא יחידה אופטי של RIS מסוג (10,5) SWCNT הפועל בטווח התדרים 0.5-30 THz. תוצאות מייצגות מדגימות את היעילות של גישת מידול ה-hybrid polynomial-Z לאפיון אלקטרומגנטי מדויק ולאופטימיזציה של ניתוב קרן.

ניתוח מוליכות קוונטית
פורמליזם קובו (Kubo) גילה כי ה-SWCNT (10,5) מציג מוליכות שטח מורכבת הנשלטת על ידי רכיבים דמיוניים (אינדוקטיביים) לאורך פס התרהרצ. ב-5 THz עם E_F = 0.2 eV, החלק הממשי של המוליכות היה 1.19 × ...

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

דיון

שלבים קריטיים בפרוטוקול דורשים תשומת לב מדויקת כדי להבטיח יישום מוצלח. ראשית, בחירה מדויקת של האינדקסים הכיראליים של ה-SWCNT חיונית, מאחר שפער הפס (bandgap) והתגובה האופטית תלויים מאוד בכיראליות. הכיראליות (10,5) המפורטת בפרוטוקול זה מספקת פער פס אופטימלי של 1.762 eV עבור פעולה בתדר טרהרץ, אך משתמשים המכוונים לפס תדרים שונה צריכים לחשב את הכיראליות המתאימה באמצעות נוסחת פער הפס E_g = 2ħv_F/d = (2 × 1.0546×10⁻³⁴ × 8×10⁵)/(d) eV. לדוגמה, כיוון לפעולה ב-10 THz (אנרגיית פוטון של 41 meV) דורש ננו-צי...

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

גילויים

המחברים מצהירים כי אין ניגודי עניינים.

תודות

המחברים מבקשים להביע את תודתם הכנה למרכז הבינלאומי למחקר תיאורטי וприкладתי (IATRC), רובע בגדאד, עיראק, על תמיכה מדעית וטכנית רבה. מחקר זה לא קיבל מענק ספציפי מגופי מימון במגזר הציבורי, המסחרי או הבלתי-מסחרי. משאבים חישוביים סופקו על ידי אוניברסיטת Al-Bayan.

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

חומרים

רשימת החומרים שנעשה בהם שימוש במאמר זה
שם	חברה	מספר קטלוג	הערות
CST Microwave Studio	Dassault Systèmesèמדידה של רמות חלבון בנסיוב של תאי ביופסיה של בלוטת הלבלב באמצעות בדיקת ELISA	לא רלוונטי	גרסה 2024, פותר במישור התדר
MATLAB	MathWorks	לא רלוונטי	גרסה R2014a או מאוחרת יותר
מצע קוורץ	פרוסת סיליקון אוניברסיטאית	4526	עובי של 500 ננומטר, ε_r = 3.8
כירליות של SWCNT (10,5)	NanoIntegris	SWCNT-105	בקוטר של 0.60 ננומטר, >90% מוליכים למחצה
מחשב אישי	לא רלוונטי	לא רלוונטי	32 GB RAM, מינימום 8 ליבות CPU

מקורות

Xiao, M., et al. Millimeter wave communications for future mobile networks. IEEE J Sel Areas Commun. 35, 1909-1935 (2017).
Kumar, A., et al. RIS-assisted terahertz communications for 6G networks: A comprehensive overview. IEEE Access. , (2025).
Liaskos, C., et al. A new wireless communication paradigm through software-controlled metasurfaces. IEEE Commun Mag. 56, 162-169 (2018).
Wu, Q., et al. Intelligent reflecting surface-aided wireless communications: A tutorial. IEEE Trans Commun. 69, 3313-3351 (2021).
Rafique, A., et al. Reconfigurable intelligent surfaces: Interplay of multi cell and surface-level design and performance under quantifiable benchmarks. IEEE Open J Commun Soc. 4, 1583-1599 (2023).
Mayaram, K., et al. Computer-aided circuit analysis tools for RFIC simulation: algorithms, features, and limitations. IEEE Trans Circuits Syst II. 47, 274-286 (2000).
Yesilyurt, O., Turhan-Sayan, G. Metasurface lens for ultra-wideband planar antenna. IEEE Trans Antennas Propag. 68, 719-726 (2019).
Gustavsen, B., Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv. 14, 1052-1061 (2002).
Zhang, Y., et al. Z-transform-based FDD implementations of biaxial anisotropy for radar target scattering problems. Remote Sens. 14, 2397(2022).
Hall, S. H., Heck, H. L. Advanced Signal Integrity for High-Speed Digital Designs. , John Wiley & Sons. (2011).

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

הדפסות חוזרות והרשאות

תגיות

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