Practical Vectorial Mode Solver for Dielectric Waveguides Based on Finite Differences

إليك ملخص تقني مفصل لورقة البحث بعنوان "محلل أنماط متجهي عملي للموجات في الأدلة الموجية العازلة بناءً على الفروق المحدودة" للباحث إرغون سيمشيك.

1. بيان المشكلة

تتناول الورقة البحثية تحدي الحساب الدقيق للأنماط الموجهة (ثوابت الانتشار وملفات الحقل) للأدلة الموجية العازلة غير المغناطيسية والخطية بصرياً.

• الصعوبة الجوهرية: حل معادلة الموجة المتجهية الكاملة للحقل الكهربائي ($\mathbf{E}$) في الهياكل ذات الانقطاعات في السماحية (الواجهات المادية).
• مشكلات محددة:
  • الحلول الزائفة: تولد العديد من الطرق العددية أنماطاً غير فيزيائية.
  • فرض الشروط الحدية: غالباً ما تواجه طرق الفروق المحدودة القياسية صعوبة في فرض شروط الاستمرارية بشكل صحيح عند الواجهات المادية، مما يؤدي إلى عيوب عددية، خاصة عند حساب المكون الطولي للحقل ($E_z$) من المكونات المستعرضة ($E_x, E_y$).
  • التعقيد مقابل الدقة: توجد حلول متجهية كاملة عالية الرتبة، ولكنها قد تكون معقدة حاسوبياً وصعبة التفسير أو التكامل مع النهج القائمة على البيانات.

2. المنهجية

يقترح المؤلف محلاً عددياً قائماً على الفروق المحدودة يصيغ المشكلة كـ مسألة القيمة الذاتية المعممة باستخدام صياغة الحقل الكهربائي.

أ. الصياغة النظرية

• المعادلة الحاكمة: يبدأ المحلل بمعادلة الموجة المتجهية للحقل الكهربائي في منطقة خالية من المصادر:
  $$\nabla^2 \mathbf{E} + \nabla \left( \frac{1}{\varepsilon} \nabla \varepsilon \cdot \mathbf{E} \right) + k_0^2 \varepsilon_r \mathbf{E} = 0$$
• تفكيك الحقل: بافتراض تغيراً متناغماً زمنياً $e^{j(\omega t - \beta z)}$، يتم تفكيك الحقل إلى مكونات مستعرضة ($E_x, E_y$) وطولية ($E_z$).
• النهج المتسق ذاتياً: على عكس الطرق التي تحل $E_x$ و $E_y$ أولاً ثم تشتق $E_z$ عبر علاقة اشتقاق (مما يؤدي إلى أخطاء عند الواجهات)، يقوم هذا النهج بحل جميع المكونات الثلاثة ($E_x, E_y, E_z$) في آن واحد.
• بناء المصفوفة: يتم تحويل المعادلات القياسية المقترنة إلى مسألة قيمة ذاتية مصفوفية من الشكل:
  $$\mathbf{M} \mathbf{E} + \beta \mathbf{L} \mathbf{E} - \beta^2 \mathbf{I} \mathbf{E} = 0$$
  لتحويل هذه المسألة التربيعية إلى مسألة خطية، تم إدخال متغير مساعد $\mathbf{F} = \beta \mathbf{E}$، مما حولها إلى مسألة قيمة ذاتية معممة خطية:
  $$\begin{pmatrix} \mathbf{M} & \mathbf{L} \\ \mathbf{0} & \mathbf{I} \end{pmatrix} \begin{pmatrix} \mathbf{E} \\ \mathbf{F} \end{pmatrix} = \beta \begin{pmatrix} \mathbf{0} & \mathbf{I} \\ \mathbf{I} & \mathbf{0} \end{pmatrix} \begin{pmatrix} \mathbf{E} \\ \mathbf{F} \end{pmatrix}$$
  يضمن هذا التعامل مع الشروط الحدية والمشتقات بشكل متسق عبر متجه الحقل بأكمله.

ب. التنفيذ العددي

• التقطيع (Discretization): تُستخدم شبكة مستطيلة حيث تُقيم مكونات الحقل الكهربائي عند حواف الخلايا والسماحية النسبية ($\varepsilon_r$) عند مراكز الخلايا.
• تقريب الاشتقاق: بينما ذُكرت الفروق المركزية من الدرجة الأولى للوضوح المفاهيمي، يستخدم التنفيذ نماذج فروق مركزية من الدرجة الرابعة للمشتقات المكانية لتحسين الدقة:
  $$\frac{\partial A}{\partial \rho} \approx \frac{-A(\rho+2h) + 8A(\rho+h) - 8A(\rho-h) + A(\rho-2h)}{12h}$$
• المحلل (Solver): يتم حل نظام المصفوفة المتفرقة الكبيرة الناتجم باستخدام دالة eigs في MATLAB.
• المعالجة اللاحقة:
  • يتم حساب الحقول المغناطيسية ($H_x, H_y, H_z$) مباشرة من مشتقات الحقل الكهربائي باستخدام معادلات "كيرل" لماكسويل (والتي لا تتضمن مشتقات $\varepsilon$، مما يتجنب أخطاء الواجهات).
  • يتم تصفية الأنماط للاحتفاظ فقط بالحلول ذات المعنى الفيزيائي ($n_{eff} < n_{cladding}$).
  • يتم تطبيع الحقول لكثافة قدرة تبلغ $1 \text{ W/m}^2$.

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

1. صياغة متجهية موحدة: المساهمة الأساسية هي تضمين جميع مكونات الحقل الكهربائي الثلاثة في صياغة قيمة ذاتية واحدة ومتسقة ذاتياً. هذا يقلل من عدم الدقة العددية المرتبطة بحساب $E_z$ لاحقاً عند انقطاعات السماحية.
2. العملية والبساطة: يتجنب النهج التعقيد الشديد للمخططات عالية الرتبة مع الحفاظ على دقة كافية للتحليل الاستكشافي والتكامل مع سير العمل القائم على البيانات. إنه يقدم "نقطة مثالية" بين الكفاءة الحسابية والدقة الفيزيائية.
3. القوة عند الواجهات: من خلال معالجة المتجه الكامل للحقل في وقت واحد، يقلل المحلل بفعالية من الحلول الزائفة والعيوب العددية عند حدود المواد.

4. النتائج العددية

تم التحقق من صحة المحلل مقابل Tidy3D (وهو محلل FDTD تجاري) باستخدام هندستين مختلفتين للدليل الموجي عند طول موجي 1550 نانومتر.

• الحالة 1: دليل موجي مستطيل من $\text{Si}_3\text{N}_4$

  • الهندسة: قلب بمقاس $1.6 \mu\text{m} \times 0.7 \mu\text{m}$، محاط بـ $\text{SiO}_2$.
  • الشبكة: 20 نقطة لكل طول موجي (PPW).
  • النتيجة: أظهرت معاملات الانكسار الفعالة ($n_{eff}$) لأول أربعة أنماط فرقاً يقل عن 0.23% مقارنة بـ Tidy3D (الذي استُخدم كحقيقة أرضية بـ 80 PPW).
  • ملفات الحقل: تطابقت توزيعات المقدار لـ $|E_x|$, $|E_y|$, $|E_z|$ والحقول المغناطيسية مع السلوك الفيزيائي المتوقع، مع تركيز الطاقة في القلب.

• الحالة 2: دليل موجي شبه منحرف على غشاء رقيق

  • الهندسة: دليل موجي شبه منحرف (عرض يتراوح من $1.2 \mu\text{m}$ إلى $1.0 \mu\text{m}$) يحاكي نيوبات الليثيوم ذات الغشاء الرقيق، فوق ركيزة.
  • الشبكة: 30 PPW لالتقاط العرض المتغير.
  • النتيجة: كان الفرق في $n_{eff}$ مقارنة بـ Tidy3D أقل من 0.88% لجميع الأنماط الأربعة.
  • الملاحظة: نجح المحلل في التقاط تغيرات الحقل بالقرب من حواف الدليل الموجي والملف المتماثل للنمط الأساسي.

5. الأهمية

• أداة سهلة الوصول: توفر الورقة البحثية أداة "عملية" أسهل في التنفيذ والتفسير من المحللات عالية الرتبة الحديثة، مما يجعلها مناسبة للنماذج الأولية السريعة والأغراض التعليمية.
• التكامل مع البيانات: الطبيعة المبسطة والقابلة للتفسير للكود تجعله مناسباً جداً للتكامل مع حلقات تعلم الآلة أو التحسين حيث تُعطى الأولوية للسرعة الحسابية وشفافية الكود على دقة الدقة القصوى.
• التحقق من النهج: تثبت النتائج أن مخطط الفروق المحدودة المصمم بعناية، والذي يحترم الطبيعة المتجهية للحقل عند الواجهات، يمكنه تحقيق دقة تقل عن 1% دون الأعباء الإضافية للطرق المعقدة القائمة على المعادلات التكاملية أو طرق العناصر المحدودة عالية الرتبة.

في الختن، يقدم هذا العمل نهجاً موثوقاً وفعالاً وسليماً رياضياً للتحليل النمطي للأدلة الموجية العازلة، محققاً توازناً فعالاً بين التعقيد العددي والدقة الفيزيائية.