Adaptive time-domain simulation of optical cavities with arbitrary dynamics

ملخص تقني: محاكاة زمنية تكيفية للتجاويف البصرية ذات الديناميكيات التعسفية

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

المنهجية
يقدم المؤلفون محاكياً زمنياً سريعاً يعتمد على صيغة تكرارية للمجال الكهربائي داخل التجويف. تقوم نواة الخوارزمية بنشر المجال كمجموع عبر الدورات المتكررة، مع الاحتفاظ صراحة بذاكرة التجويف للمدخلات السابقة ومواضع المرايا. يتم حساب المجال الكهربائي $E(t)$ باستخدام المعادلة التكرارية:

$$E(t) = t_a \sum_{n=0}^{N-1} \left[ (r_a r_b)^n e^{-2ikS_n(t)} E_{in}(t - 2nT) \right] + (r_a r_b)^N e^{-2ikS_N(t)} E(t - 2NT)$$

حيث تمثل $t_a, r_a, r_b$ معاملات المرآة، و $k$ هو العدد الموجي، و $T$ هو زمن الانتشار في اتجاه واحد، و $S_n(t)$ يمثل تغيرات المسار البصيكي المتراكمة نتيجة حركة المرآة.

تشمل الميزات المنهجية الرئيسية ما يلي:

• شروط حدودية تعسفية: يمكن تعديل كل من المجالات الكهربائية المدخلة ومواضع المرايا عند كل خطوة محاكاة، مما يسمح بتعديل سعة وطور (phase) تعسفيين، وتغيرات زمنية في طول التجويف.
• تردد أخذ العينات التكيفي: يحدد المستخدم تردد أخذ العينات المرغوب ($f^{desired}_{calc}$)، لكن المحاكي يعدل هذا داخلياً إلى تردد فعال ($f_{calc}$) يظل متسقاً مع بنية دورة التجويف. ويتم تحقيق ذلك من خلال شجرة قرار قائمة على القواعد تتعامل مع ثلاثة أنظمة:
  1. نظام التردد العالي ($f^{desired}_{calc} > f_{2T}$): تقوم الخوارزمية بدمج عدة سجلات فرعية مزاحة طورياً لزيادة أخذ العينات الفعال.
  2. نظام التردد المنخفض ($f^{desired}_{calc} < f_{2T}/N_{max}$): يتم تحديد نافذة التكامل عند أقصى عدد من الدورات المتكررة ($N_{max} = 5N_{eff}$) لمنع انتشار المساهمات التاريخية الضئيلة، مع الحفاظ على معدل أخذ العينات المطلوب عبر مساهمات جزئية للدورة المتكررة.
  3. النظام المتوسط: تستخدم قاعدة تقريب المنحنى العكسي لاختيار عدد صحيح أمثل من الدورات المتكررة ($N$) لتقليل الخطأ بين الترددات الفعلية والمرغوبة.
• الكفاءة الحسابية: يتجنب المحاكي إعادة تقييم تاريخ المجال الكهربائي الكامل عند كل خطوة. بدلاً من ذلك، يستخدم مخزناً مؤقتاً محدوداً للمدخلات السابقة وعوامل اضمحلال هندسية مسبقة الحساب. يستفيد التنفيذ من التجميع في الوقت المناسب (Just-In-Time compilation) عبر (Numba) لتسريع عملية الجمع عبر سجلات المجال المخزنة.

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

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

النتائج والتحقق
تم التحقق من صحة المحاكي مقابل البيانات التجريبية من تجويف الذراع الشمالي لمرصد "فيرجو" (Virgo) أثناء تجربة إثارة ميكانيكية.

• المقارنة التجريبية: أعادت المحاكاة إنتاج قدرة النفاذ المستمر (DC transmission power) وإشارات خطأ "بود-دريفير-هال" (PDH) أثناء عبور الرنين غير الأديباتي (سرعات المرايا تصل إلى $\approx 7.6$ من السرعة الحرجة). أظهرت النت‌های توافقاً جيداً في كل من السعة والبنية الزمنية، مما يؤكد قدرة النموذج على التقاط الاقتران غير الخطي بين المجال والمرآة.
• الاختبار المعياري: أظهرت اختبارات الأداء على 400 تكوين للتجويف أن الخلفية المدعومة بـ Numba تتفوق بشكل كبير على عمليات التنفيذ بلغة بايثون الصرفة. يتوسع عامل التسريع مع عدد الدورات المتكررة ($N$)، ليصل إلى 17.61 للتجاويف عالية الفاصة ($N=300$)، مع أزمنة تنفيذ منخفضة تصل إلى 16 ميكروثانية.
• العمومية: نجح الإطار في إعادة إنتاج أنظمة خمود التجويف النبضي المبلغ عنها في الأدبيات السابقة، مما يثبت قابليته للتطبيق بما يتجاوز تجاويف فابري-بيرو القياسية لتشمل أنظمة مثل الرنانات الحلقية.

الأهمية
يضع هذا البحث المحاكي كأداة متعددة الاستخدامات لدراسات المجال الزمني للمرنانات البصرية حيث تكون الديناميكيات العابرة والآثار غير الأديباتية حاسمة. وتكمن أهميته الأساسية في كفاءته ومرونته، مما يجعله مناسباً لـ:

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

يخلص المؤلفون إلى أن هذا الإطار يسد الفجوة بين النمذجة الفيزيائية المفصلة والمتطلبات الحسابية لاستراتيجيات التحكم الحديثة القائمة على البيانات.