2025-11-15T23:07:18.927821

Green's function expansion for multiple coupled optical resonators with finite retardation using quasinormal modes

Fuchs, Ren, Hughes et al.
The electromagnetic Green's function is a crucial ingredient for the theoretical study of modern photonic quantum devices, but is often difficult or even impossible to calculate directly. We present a numerically efficient framework for calculating the scattered electromagnetic Green's function of a multi-cavity system with spatially separated open cavities (with arbitrary shape, dispersion and loss) and finite retardation times. The framework is based on a Dyson scattering equation that enables the construction of the Green's function from the quasinormal modes of the individual resonators within a few-mode approximation and a finite number of iteration steps without requiring nested integrals. The approach shows excellent agreement with the full numerical Green's function for the example of two coupled dipoles located in the gaps of two metal dimers serving as quasinormal mode cavities, and is easily extended to arbitrarily large separations and multiple cavities.
academic

توسيع دالة جرين للرنانات البصرية المقترنة المتعددة مع التأخير المحدود باستخدام الأنماط شبه الطبيعية

المعلومات الأساسية

  • معرّف الورقة: 2510.12511
  • العنوان: توسيع دالة جرين للرنانات البصرية المقترنة المتعددة مع التأخير المحدود باستخدام الأنماط شبه الطبيعية
  • المؤلفون: روبرت فوكس¹، جوانجوان رين²، ستيفن هيوز²، مارتن ريختر¹
  • المؤسسات: ¹جامعة برلين التقنية، ²جامعة كوين
  • التصنيف: فيزياء المواد المكثفة، البصريات، الفيزياء الكمية
  • تاريخ النشر: 15 أكتوبر 2025
  • رابط الورقة: https://arxiv.org/abs/2510.12511

الملخص

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

خلفية البحث والدافع

خلفية المشكلة

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

المساهمات الأساسية

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

شرح الطريقة

تعريف المهمة

الإدخال: نظام N من التجاويف المفتوحة المفصولة مكانياً، كل تجويف له شكل وتشتت وخسائر عشوائية الإخراج: دالة جرين الكهرومغناطيسية الكاملة للنظام متعدد التجاويف G^(N)(r,r',ω) القيود: ضمن تقريب قليل الأنماط، مع الأخذ في الاعتبار تأثيرات التأخير المحدود، تجنب التكاملات المتداخلة المعقدة

الإطار النظري

1. وصف النظام

يتم تمثيل ثابت العزل الكهربائي الكامل لنظام متعدد التجاويف كما يلي:

ε(r,ω) = ε_back(r,ω) + Σ(i=1 to N) V_i(r,ω)

حيث V_i(r,ω) هو الاضطراب في ثابت العزل الكهربائي الناجم عن التجويف i.

2. الأنماط شبه الطبيعية للتجويف الفردي

تحقق الأنماط شبه الطبيعية f̃_iμ(r) معادلة هيلمهولتز بدون مصادر:

∇×∇×f̃_iμ(r) - (ω̃²_iμ/c²)ε_i(r,ω̃_iμ)f̃_iμ(r) = 0

مع ترددات مميزة معقدة ω̃_iμ = ω_iμ - iγ_iμ، حيث γ_iμ > 0.

3. حقول QNM المنتظمة

لحل مشكلة تباعد QNM في المجال البعيد، يتم إدخال حقول QNM المنتظمة:

F̃_iμ(r,ω) = ∫_V_i d³r' Δε(r',ω)G_back(r,r',ω)·f̃_iμ(r')

4. البناء التكراري لدالة جرين متعددة التجاويف

الابتكار الأساسي هو البناء التكراري من خلال سلسلة من معادلات التناثر:

G^(n)(r,r',ω) = G^(n-1)(r,r',ω) + ∫_V_n d³s Δε(s,ω)G^(n-1)(r,s,ω)·G^(n)(s,r',ω)

شروط الإنهاء:

G^(k)(r,r',ω)|_(r,r'∈V_i) ≈ G_i(r,r',ω)

5. تمثيل التكامل السطحي

لتحسين الكفاءة العددية، يتم تحويل التكاملات الحجمية إلى تكاملات سطحية:

ω²/c² G^(n)(r,r')|_(r∉V_n) = χ_V_n(r') ω²/c² G^(n-1)(r,r') + [حدود التكامل السطحي]

نقاط الابتكار التقني

  1. استراتيجية التحليل: يتم تحليل التناثر متعدد التجاويف بشكل طبيعي إلى منتج عمليات التناثر ثنائي التجاويف، مما يتجنب التكاملات المتداخلة
  2. معالجة تأثيرات التأخير: يتم تضمين أوقات التأخير المحدودة بشكل صحيح من خلال عامل الطور e^(iωR_ij/c)
  3. تقريب القطب: يتم تطبيق تقريب القطب على دوال الغلاف البطيئة التغير ω→ω_i، مما يحسن الكفاءة الحسابية
  4. تحسين التكامل السطحي: استخدام التكاملات السطحية لتجنب مشاكل عدم استمرارية المجال الكهربائي عند حدود الوسط

الإعداد التجريبي

تكوين النظام

  • النظام المختبر: ثنائي معادن كرنان يعمل كتجويف QNM
  • نموذج المادة: نموذج درود، ε_j(ω) = 1 - ω_p²/(ω² + iωγ_pj)
  • إعدادات المعاملات:
    • تردد البلازما: ℏω_p = 8.2934 eV
    • معاملات التخميد: ℏγ_p1 = 0.0928 eV, ℏγ_p2 = 0.3ℏγ_p1
  • ترددات QNM المميزة:
    • ثنائي المعادن 1: ℏω̃_1 = (1.6904 - i0.0652) eV (λ_1 ≈ 733.46 nm)
    • ثنائي المعادن 2: ℏω̃_2 = (1.6482 - i0.0388) eV (λ_2 ≈ 752.24 nm)

تكوين الاختبار

  • إعداد ثنائي القطب: ثنائي قطب متطابق مستقطب في اتجاه z: d_a = d_b = dê_z
  • الموضع: يقع في فجوات ثنائي المعادن على التوالي (r_a ∈ V_1, r_b ∈ V_2)
  • مسافة الفصل: R_12 = 2020 nm و 760 nm في تكوينين

مؤشرات التقييم

قوة الاقتران ثنائي القطب:

g_ba(ω) = (1/ε_0) d_b · G(r_b, r_a, ω) · d_a

طرق المقارنة

  1. دالة جرين العددية الكاملة (G_full): حل صارم لمعادلات ماكسويل ثلاثية الأبعاد
  2. توسيع QNM (G_QNM): طريقة هذه الورقة، الصيغ (12)+(14)
  3. تقريب القطب (G_QNM^pole): الصيغ (12)+(15)، يتم تضمين الطور في تقريب القطب
  4. التوسيع أحادي التجويف (G_1): الطريقة التقليدية التي تأخذ في الاعتبار فقط QNM للتجويف الفردي

النتائج التجريبية

النتائج الرئيسية

1. حالة الفصل الكبير (R_12 = 2020 nm)

  • تُظهر طريقة توسيع QNM اتفاقاً ممتازاً مع الحل العددي الكامل
  • يعتبر المعالجة الصحيحة لتأثيرات التأخير (e^(iωR_12/c)) حاسمة للدقة
  • تنحرف طريقة تقريب القطب (G_QNM^pole) عن الحل العددي لعدم تضمينها تأثيرات التأخير بشكل كامل
  • يفشل التوسيع أحادي التجويف في تقدير قوة الاقتران بشكل كبير، ولا يمكنه التقاط تأثيرات اقتران QNM بين التجاويف

2. حالة الفصل المتوسط (R_12 = 760 nm)

  • تستمر طريقة توسيع QNM في المطابقة المثالية مع الحل العددي
  • تُظهر طريقة تقريب القطب أداءً أفضل بسبب ضعف تأثيرات التأخير
  • يستمر التوسيع أحادي التجويف في الفشل من الناحية الكمية والنوعية

3. النتائج العددية الرئيسية

بالنسبة إلى R_12 = 2020 nm: N_21 = (2.0694 - i0.1357) × 10^(-7) nm^(-2) بالنسبة إلى R_12 = 760 nm: N_21 = (5.3494 + i0.2984) × 10^(-7) nm^(-2)

التحقق من الطريقة

  • بدون معاملات ملائمة: تستند جميع الحسابات إلى المبادئ الأولى، بدون الحاجة إلى ضبط المعاملات
  • الاعتماد على التردد: تم التقاط الخصائص المعتمدة على التردد لقوة الاقتران بشكل صحيح
  • تأثيرات الطور: يعتبر طور التأخير e^(iωR_12/c) حاسماً للدقة في مسافات الفصل الكبيرة

الكفاءة الحسابية

  • تجنب التكاملات المتداخلة: يتم تحليل التناثر متعدد التجاويف إلى عمليات ثنائية التجاويف
  • خطوات تكرارية محدودة: يتطلب نظام بـ N تجويف N خطوة تكرارية
  • تحجيم تربيعي: ضمن تقريب أحادي النمط، يتناسب عدد تكاملات التناثر بشكل تربيعي مع عدد التجاويف

الأعمال ذات الصلة

الاتجاهات البحثية الرئيسية

  1. نظرية الأنماط شبه الطبيعية: نظرية الانبعاث التلقائي للتجاويف المفتوحة لـ Sauvan وآخرين
  2. تكميم المجال: طرق تكميم المجال الكهرومغناطيسي بناءً على QNM
  3. نظرية التجاويف المقترنة: نظرية الأنماط شبه الطبيعية المقترنة (CQT) وقيودها
  4. طرق دالة جرين: توسيع QNM للتجويف الفردي وتقنيات التنظيم

مزايا هذه الورقة

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

الخلاصة والنقاش

الاستنتاجات الرئيسية

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

القيود

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

الاتجاهات المستقبلية

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

التقييم المتعمق

المزايا

  1. الابتكار النظري قوي: التطبيق التكراري لمعادلة التناثر الديسون جديد وذكي
  2. القيمة العملية عالية: يحل احتياجات الحساب العملية لأنظمة البصريات الكمية متعددة التجاويف
  3. التحقق كافٍ: تم التحقق من دقة الطريقة من خلال المقارنة مع الحل العددي الصارم
  4. الكتابة واضحة: الاشتقاق النظري دقيق والصورة الفيزيائية واضحة

أوجه القصور

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

التأثير

  1. المساهمة الأكاديمية: توفير أداة مهمة للنظرية البصرية الكمية متعددة التجاويف
  2. التطبيق التقني: لها قيمة تطبيقية محتملة في مجالات الشبكات الكمية والحوسبة الكمية
  3. الأهمية المنهجية: توضح إمكانيات تطبيق نظرية QNM في الأنظمة المعقدة
  4. قابلية إعادة الإنتاج: الإطار النظري واضح، مما يسهل على الباحثين الآخرين التطبيق والتوسع

السيناريوهات المناسبة

  1. تصميم أجهزة البصريات الكمية: التحليل النظري لأنظمة التجاويف المقترنة
  2. دراسة البلازمونيات: حساب الاستجابة الكهرومغناطيسية للهياكل المعدنية النانوية
  3. الشبكات الكمية: تحليل الاقتران بين العقد الكمية المفصولة مكانياً
  4. البحث في بلورات الفوتون: دراسة أنظمة التجاويف المعيبة المقترنة

المراجع

تستشهد الورقة بـ 47 مرجعاً مهماً، تغطي نظرية الأنماط شبه الطبيعية والبصريات الكمية وتكميم المجال الكهرومغناطيسي وغيرها من المجالات ذات الصلة، مما يوفر أساساً نظرياً متيناً لهذا البحث.

التقييم الشامل: هذه ورقة بحثية عالية الجودة في الفيزياء النظرية، تقترح طريقة مبتكرة لحل المشكلة المهمة المتمثلة في حساب دالة جرين لأنظمة بصرية متعددة التجاويف. الاشتقاق النظري دقيق والتحقق العددي كافٍ، مما يجعلها ذات أهمية كبيرة لمجالات البصريات الكمية والفوتونيات. تجعل الشمولية والقابلية للتوسع للطريقة لها آفاق تطبيقية واسعة.