Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures

ملخص تقني: إطار عمل التوازن التوافقي كامل الموجة للهياكل الفوتونية المشتتة والمتغيرة زمنياً

بيان المشكلة
تقدم الهياكل الفوتونية المتغيرة زمنياً قدرات فريدة مثل تحويل التردد، وعدم التبادلية، والتحكم الديناميكي في جبهة الموجة من خلال كسر تماثل الإزاحة الزمنية. ومع ذلك، لا يزال النمذجة الدقيقة لهذه الأنظمة تشكل تحدياً كبيراً عندما تكون المواد المعنية مشتتة بقوة (مثل الأكاسيد الموصلة الشفافة، أو الجرافين، أو مواد "إبسيلون القريبة من الصفر" ENZ). في هذه المنصات، لا يمكن معاملة التعديل الزمني كمجرد اضطراب لحظي لمعلمة مستقلة عن التردد؛ بل تقتضي السببية أن الاستجابة المادية تعتمد على تاريخ الاستثارة، مما يعني أن كل توافق فلوكي (Floquet harmonic) مُولَّد يخضع لاستجابة مادية متميزة مع بقائه مرتبطاً بجميع التوافقات الأخرى من خلال التعديل. تواجه الطرق الحسابية الحالية مقايضات: فالطرق في المجال الزمني (مثل FDTD) غالباً ما تتطلب عمليات محاكاة طويلة للوص \text{إلى} الحالة المستقرة وتواجه صعوبة في التحديثات المشتتة السببية، بينما تفتقر الطرق التحليلية أو طرق الموجة المقترنة المتخصصة إلى المرونة الهندسية للتعامل مع الهياكل ثلاثية الأبعاد المعقدة، أو الأوساط غير المتجانسة، أو أطياف التعديل التعسفية.

المنهجية
يقترح المؤلفون إطار عمل للتوازن التوافقي كامل الموجة يحل الاستجابة الكهرومغناطيسية للحالة المستقرة مباشرة في مجال التردد، متجنباً التقدم الزمني. يتضمن جوهر المنهجية إعادة صياغة المشكلة المتغيرة زمنياً كطقم من المشكلات الفرعية المرتبطة في مجال التردد ضمن فضاء فلوكي موسع.

تشمل العناصر الرئيسية للصياغة ما يلي:

• تمثيل المصدر الثانوي: لا يتم تمثيل المكونات المتغيرة زمنياً كشروط حدودية تعتمد على الزمن، بل كمصادر ثانوية مستحثة. تُستخدم كثافات الاستقطاب الحجمية ($P$) للأوساط الحجمية المعدلة، وكثافات التيارات السطحية ($J_s$) للألواح الموصلة المعدلة.
• توسع فلوكي (Floquet Expansion): يتم توسيع جميع المجالات والمصادر المستحثة في سلسلة فلوكي مبتورة ($F(r, t) = \sum F_n(r)e^{j\omega_n t}$)، حيث $\omega_n = \omega_0 + n\Omega$.
• صياغة النظام المقترن: يتم حل المشكلة كنظام خطي عالمي $(L - SC)X = Y_{ext}$.
  • يمثل $L$ مؤثر ماكسويل في مجال التردد للخلفية الساكنة، ويتم تقييمه عند التردد الفيزيائي لكل توافق ($\omega_n$) لمراعاة التشتت المادي.
  • يمثل $S$ مؤثرات حقن المصدر.
  • يمثل $C$ مؤثرات توليد المصدر المشتقة من التوسع فورييه الزمني للمعلمات المادية. تُدخل هذه المؤثرات اقتراناً غير قطري في الفضاء التوافقي، مما يمثل آلية تحويل التردد.
• التنفيذ: تم تنفيذ إطار العمل باستخدام طريقة العناصر المحدودة (FEM). يتم رفع متجه المجهول القياسي لـ FEM إلى متجه في المجال التوافقي، مما يسمح لكل درجة حرية مكانية بحمل معاملات توافقية متعددة. هذا يحافظ على المرونة الهندسية للنمذجة كاملة الموجة مع القدرة على التقاط الاقتران بين التوافقات.

المساهمات الرئيسية

1. إطار عمل موحد: توفر الطريقة بنية جبرية واحدة لمعالجة كل من الأوساط الحجمية المعدلة والألواح الموصلة المعدلة، مما يدعم أطياف التعديل الزمني التعسفية (أحادية النغمة، أو متعددة النغمات، أو غير الجيبية).
2. معالجة التشتت: على عكس الطرق التي تفترض تشتتاً ضعيفاً أو تقريبات الغلاف، يقوم هذا الإطار بتقييم التشتت المادي ومؤثرات الإشعاع عند التردد المحدد لكل توافق مُوَلَّد.
3. الملاحظات المحللة توافقياً: يعطي الحل مباشرة سعة الحالة المستقرة، والطور، وملفات تعريف المجال لجميع توافقات فلوكي المحتفظ بها، مما يتيح حساب التشتت، والامتصاص، وتدفق القدرة، وتوزيعات المجال القريب الخاص بكل توافق دون الحاجة لمعالجة بيانات المجال الزمني.

النتائج
تم التحقق من صحة الإطار وتجربته من خلال مثالين ممثلين:

1. أسطوانة جرافين ملفوفة ذات ضخ بارامتري:

   • الإعداد: أسطوانة عازلة مغلفة بطبقة جرافين مع مستوى فيرمي معدل زمنياً ($\Omega = 2\omega_0$)، تعمل في نظام التضخيم البارامتري المتدهور.
   • التحقق: أظهرت النتائج العددية للتشتت والامتصاص المحلل توافقياً توافقاً ممتازاً مع حل فلوكي التحليلي المستمد من الهندسة الأسطوانية.
   • النتائج: نجحت المحاكاة في رصد إعادة توزيع الطاقة من الموجة الحاملة إلى التوافق المخفض ($\omega_{-1}$)، مما أكد قدرة الإطار على نمذجة التشتت القوي والتعديل السريع في أنظمة البلازمون السطحي الرنينية.

2. سطح ميتا-سطحي بصري (يعتمد على ITO):

   • الإعداد: سطح ميتا-سطحي عاكس يتكون من محزز $TiO_2$ ثنائي مع طبقة (ITO) ذات "إبسيلون قريبة من الصفر" (ENZ) معدلة زمنياً ومدعومة بالذهب. تم تعديل تردد البلازما لـ ITO لدفع تغييرات فائقة السرعة في السماحية.
   • التطبيق: استُخدم الإطار لتصميم وتحسين السطح الميتا-سطحي من أجل انحراف الحزمة المحول للتردد. كان الهدف هو تحويل الموجة الضوئية الساقطة إلى نطاقات جانبية وإعادة توجيهها إلى رتبة حيود مكانية محددة ($m=1$).
   • النتائج: نجح التصميم الأمثل في إعادة توجيه النطاقات الجانية المولدة ($\omega_{\pm 1}$) إلى رتبة الحيود المكانية الأولى. أظهرت النتائج التفاعل بين تشتت ENZ، واقتران درود (Drude) المعتمد على الزمن، والحيود المكاني، مما يؤكد فائدة الطريقة لتصميم أجهزة بصرية واقعية ومعقدة.

الأهمية
يزعم البحث أن هذا الإطار يضع منصة عامة لنمذجة وتصميم الهياكل الفوتونية النشطة المشتتة. ومن خلال الجمع بين الشفافية الفيزيائية لمنهجيات التوازن التوافقي والمرونة الهندسية لـ FEM كاملة الموجة، يتغلب هذا الأسلوب على قيود الحلول الحالية للأنظمة المتغيرة زمنياً. فهو يسمح للباحثين بقياس كيفية إعادة توزيع التعديل للطاقة الكهرومغناطيسية بين القنوات المزاحة ترددياً بشكل مباشر، مما يوفر ملاحظات أساسية (القدرة، الامتصاص، والحيود المحلل توافقياً) لتصميم الأنظمة الفوتونية الزمانية-المكانية التالية القائمة على الجرافين، والأكاسيد الموصلة الشفافة، والأسطح الميتا-سطحية الرنينية. يضع هذا العمل الإطار ليس فقط كأداة للتحليل، بل أيضاً للتصميم العكسي للأجهزة التي تتطلب تحويل التردد، والتضخيم البارامتري، والتحكم الديناميكي في جبهة الموجة.