Immiscible to miscible quenching instabilities in two-dimensional binary Bose-Einstein condensates

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

تتقصى هذه الدراسة عدم الاستقرار الديناميكي والاضطراب المتولد في مكثف بوز-أينشتاين (BEC) ثنائي الأبعاد (2D) من نظائر الروبيديوم ($^{85}\text{Rb}$ و $^{87}\text{Rb}$) المحصور داخل صندوق دائري منتظم. وبشكل محدد، يركز البحث على انتقالات التبريد من عدم الامتزاج إلى الامتزاج (IMQT).

تتمثل المشكلة الجوهرية في فهم كيف يؤدي الخفض المفاجئ في طول التشتت بين الأنواع ($a_{12}$) إلى دفع النظام من حالة منفصلة طورياً (غير ممتزجة) إلى حالة مختلطة (ممتزجة). يهدف المؤلفون إلى:

• التمييز بين الديناميكيات المدفوعة بتغيرات التفاعل غير الخطي (التبريد/الخمد) وتلك المدفوعة بالقوى الخطية الخارجية.
• تحليل ظهور الاضطراب الكمي، وتوليد الدوامات، وإنتاج موجات الصوت (الفونونات).
• استقصاء قوانين القياس لأطياف الطاقة الحركية (القابلة للانضغاط مقابل غير القابلة للانضغاط) ومقارنتها مع قياس اضطراب كولموغوروف الكلاسيكي ($k^{-5/3}$).
• تحديد تأثير التكوينات المكانية الأولية المختلفة وحجم عملية التبريد على تطور النظام والسلوك التقاربي للنظام.

2. المنهجية

يستخدم المؤلفون نهجاً عددياً يعتمد على صيغة مجال المتوسط لـ "غروس-بيتايفسكي" (GP) المكيفة لنظام ثنائي ثنائي الأبعاد.

• النموذج الفيزيائي:
  • النظام: خليط ثنائي من $^{85}\text{Rb}$ (النوع 1) و $^{87}\text{Rb}$ (النوع 2) بأعداد ذرية متساوية ($N = 3.5 \times 10^4$).
  • الحصر: جهد صندوق دائري ثنائي الأبعاد بنصف قطر $R$ وارتفاع $V_0 \gg \mu$ (الكيميائي).
  • المعادلات: يتم حل معادلات GP ثنائية الأبعاد المترابطة وغير البعدية باستخدام طريقة فورييه للخطوات المنقسمة على شبكة $800 \times 800$.
  • المعلمات: تم تثبيت أطوال التشتت داخل النوع الواحد عند $a_{11} = a_{22} = 90a_0$. طول التشتت بين الأنواع $a_{12}$ هو المعلمة المتحكم بها.
• بروتوكول التبريد (Quenching Protocol):
  • يتم تهيئة النظام في حالة أرضية غير ممتزجة حيث $a_{12} > a_{ii}$.
  • عند $t=0$، يتم خفض $a_{12}$ فجأة (تبريد) إلى قيمة حيث $a_{12} < a_{ii}$، مما يحفز الانتقال إلى الامتزاج.
  • تم اختبار تكوينين مكانين أوليين مختلفين:
    1. مركزي (ثلاثي الأجزاء): ثلاثة نطاقات مكانية غير متداخلة (النوع 2 في المركز، والنوع 1 على الجوانب). $a_{12} = 100a_0 \to 60a_0$ ($\Delta\delta = 4/9$).
    2. محوري (ثنائي الأجزاء): نطاقان منفصلان محورياً. $a_{12} = 110a_0 \to 60a_0$ ($\Delta\delta = 5/9$).
• تقنيات التحليل:
  • تفكيك الطاقة: يتم تفكيك الطاقة الحركية إلى مكونات قابلة للانضغاط (مرتبطة بتقلبات الكثافة/الموجات الصوتية) وغير قابلة للانضغاط (مرتبطة بالدوامية/الحركة الدورانية) باستخدام تفكيك هلمهولتز.
  • التحليل الطيفي: يتم حساب أطياف القدرة للطاقة الحركية في فضاء الزخم ($k$-space) لتحديد قوانين القياس.
  • معلمة الامتزاج ($\Lambda$): تُعرف عبر تكامل التداخل لملفات تعريف الكثافة للنوعين لتتبع الانتقال من الحالة غير الممتزجة ($\Lambda=0$) إلى الحالة الممتزجة ($\Lambda=1$).

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

• تعزيز مقدار الـ IMQT: على عكس الدراسات السابقة (المرجع [1]) التي استخدمت مقادير تبريد أصغر، تستخدم هذه الدراسة انتقالات أكبر بكثير ($\Delta\delta = 4/9$ و $5/9$). وهذا يسمح بمراقبة استجابات ديناميكية أسرع وفصل أوضح للمقاييس الزمنية.
• الارتباط الخطي بين التردد والتبريد: يضع المؤلفون علاقة خطية جديدة بين تردد التذبذب التقاربي لمعلمة الامتزاج ($\nu_\Lambda$) ومقدار التبريد الأولي ($\Delta\delta$).
  • الصيغة: $\Delta\delta \approx \alpha \nu_\Lambda + 0.15$.
  • يعتمد الميل $\alpha$ على الهندسة الأولية ($\alpha \approx 0.37$ لـ 3 نطاقات، $\alpha \approx 1.17$ لـ نطاقين)، بينما يمثل الجزء المقطوع ($0.15$) عتبة IMQT مستقلة عن الهندسة.
• تحديد تأثير عنق الزجاجة (Bottleneck Effect): تؤكد الدراسة وجود "تأثير عنق الزجاجة" في أطياف الطاقة قبل الوصول إلى منطقة التبديد فوق البنفسجية، مما يشير إلى انحراف عن سلوك القياس الكلاسيكي في المقاييس المتوسطة.
• هيمنة الأنماط القابلة للانضغاط: وُجد أنه في النظام التقاربي طويل الأمد، تهيمن الطاقة الحركية القابلة للانضغاط (الموجات الصوتية) على الطاقة غير القابلة للانضغاط (الدوامات) بعامل يصل إلى 5 تقريباً، بغض النظر عن التكوين الأولي.

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

• الديناميكيات والمقاييس الزمنية:
  • تؤدي مقادير IMQT الأكبر إلى بدء عدم استقرار وتطور أسرع بشكل ملحوظ. يصل النظام إلى نظام الامتزاج التقاربي في وقت أبكر بكثير (على سبيل المثال، $t \approx 10$) مقارنة بالدراسات السابقة ذات التبريد الأصغر ($t \approx 20$).
  • يولد النظام العديد من ثنائيات أقطاب الدوامات والتدفقات المضطربة فور عملية التبريد.
• القياس الطيفي (كولموغوروف):
  • خلال بدء عدم الاستقرار (الفترات الزمنية القصيرة)، تظهر أطياف الطاقة غير القابلة للانضغاط والقابلة للانضغاط قياس كولموغوروف الكلاسيكي ($k^{-5/3}$)، وهو ما يميز الاضطراب.
  • يُلاحظ هذا القياس في نوافذ زمنية محددة (على سبيل المثال، $0.5 \le t \le 0.8$ للحالة المركزية؛ $0.7 \le t \le 1.2$ للحالة المحورية).
  • مع تطور النظام، تنحرف الأطياف عن هذا القياس، مظهرةً تأثير عنق الزجاجة قبل دخول منطقة التبديد.
• الاستقرار التقاربي:
  • في الحد التقاربي طويل الأمد ($t > 20$)، يستقر إنتاج الدوامات والموجات الصوتية.
  • يظل عدد الدوامات ($N_v$) ثابتاً (غير صفري)، بينما يتذبذب تداخل الكثافة $\Lambda$ بتردد يتناسب طردياً مع قوة التبريد الأولية.
• الاختلافات في التكوين:
  • المركزي (3 نطاقات): يكون بدء عدم الاستقرار أسرع قليلاً، ويكون عنق الزجاجة الطيفي أكثر وضوحاً في النوع الخارجي ($^{85}\text{Rb}$).
  • المحوري (نطاقان): يظهر النظام أهداب تداخل موضعية بالقرب من حدود المصيدة، وهي متميزة عن الحالة المركزية.

5. الأهمية

يعزز هذا العمل فهم الاضطراب الكمي والديناميكيات غير المتوازنة في السوائل الكمية الثنائية.

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

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