Time-resolved digital quantum simulation of cosmological particle creation in a de Sitter-radiation transition

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

بيان المشكلة
تتناول الورقة البحثية تحدي محاكاة خلق الجسيمات الكونية في الخلفيات المعتمدة على الزمن، وتحديداً الانتقال من مرحلة "دي سيتر" (de Sitter) إلى مرحلة يهيمن عليها الإشعاع في كون مسطح مكانياً من نوع "فليت-كورت-روبرتسون-ووكر" (FLRW). بينما ركزت المحاكاة الكمومية القائمة على البوابات (gate-based) السابقة للظواهر في الفضاء المنحني (مثل Maceda و Sabín [13]) على تنفيذات "الدفعة الواحدة" (one-shot)—أي تجميع تحويل "بوغوليوبوف" التحليلي النهائي في دائرة واحدة—إلا أن هذا النهج يحجب التاريخ الديناميكي لعملية إنتاج الجسيمات. يهدف المؤلفون إلى تطوير محاكاة كمومية رقمية معتمدة على الزمن تقوم بتجزئة التطور المعتمد على الزمن الكونفي (conformal time)، مما يسمح بمراقبة تراكم إشغال الأزواج في القاعدة الثابتة أثناء الانتقال غير الأديباتي (non-adiabatic). وتعد هذه الطريقة ذات أهمية خاصة في الخلفيات التي لا تتوفر فيها معاملات "بوغوليوبوف" بصيغة مغلقة.

المنهجية
يصيغ المؤلفون المشكلة باستخدام حقل قياسي حقيقي، عديم الكتلة، ومتصل أدنى (minimally coupled)، في نموذج انتقال مفاجئ. يكون عامل القياس $a(\eta)$ مستمراً في مشتقته الأولى، ولكنه غير مستمر في مشتقته الثانية عند زمن كونفي $\eta_e$، مما يقرب عملية انتقال سريع من التضخم إلى هيمنة الإشعاع.

1. صياغة الهاملتوني (Hamiltonian Formulation):

   • يتم رسم ديناميكيات الحقل باستخدام مؤثرات السلم (ladder operators) ثابتة زمنياً مُعرفة بالنسبة لتردد مرجعي $\omega_0 = k$، وذلك ضمن فضاء فوك (Fock space) ثابت.
   • يتم ترميز تطور الخلفية عبر معاملات هاملتوني معتمدة على الزمن $c_Z(\eta)$ و $c_A(\eta)$، والتي تقود التفاعلات الحافظة للعدد وتفاعلات خلق الأزواج، على التوالي.
   • يأخذ الهاملتوني لزوج من الزخم $(+k, -k)$ الشكل $H_k(\eta) = c_Z(\eta)Z + c_A(\eta)A$، حيث $Z$ يعد الكمات الإجمالية و $A$ يولد ضغط الحالة ثنائية النمط (two-mode squeezing).

2. التجزئة المعتمدة على الزمن (Time-Resolved Discretization):

   • يتم تجزئة تطور الزمن الكونفي إلى $N$ من الخطوات المنتظمة باستخدام متغير عديم الأبعاد $y = k\eta$.
   • يتم تقريب مؤثر التطور المرتب زمنياً عبر تقسيم "سترانج" (Strang splitting) من الدرجة الثانية (Trotterization) لكل شريحة زمنية:
     $$U^{(n)}_k \approx e^{-i\tilde{c}_Z Z \Delta y/2} e^{-i\tilde{c}_A A \Delta y} e^{-i\tilde{c}_Z Z \Delta y/2}$$
   • بخلاف مناهج الدفعة الواحدة، تقوم هذه الطريقة بتطبيق كتل دوائر قصيرة المدى بشكل تكراري باستخدام زوايا دوران محددة بواسطة معاملات الهاملتوني اللحظية، بدلاً من المعاملات النهائية المحسوبة مسبقاً.

3. ترميز الكيوبت (Qubit Encoding):

   • تم استخدام ترميز أربعة كيوبتات للإثارة الواحدة لزوج الزخم $(+k, -k)$.
   • يتم قطع كل نمط عند فضاء جزئي يحتوي على إثارة واحدة كحد أقصى، مما يربط حالات القاعدة الفيزيائية $|0_+, 0_-\rangle$، $|1_+, 0_-\rangle$، $|0_+, 1_-\rangle$، و $|1_+, 1_-\rangle$ بحالات حسابية مكونة من أربعة كيوبتات.
   • يتم تفكيك المولدات $Z$ و $A$ إلى سلاسل "بالي" (Pauli strings) $Z_q$ و $A_q$ تعمل على سجل الكيوبتات الأربعة.

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

• إطار عمل معتمد على الزمن: تقدم الورقة صياغة لتتبع التراكم الديناميكي لإشغال الأزواج عبر الانتقال، مما يميزها عن عمليات المحاكاة الساكنة التي تركز على الحالة النهائية فقط.
• تسلسل التحقق (Validation Hierarchy): تحقق المؤلفون من المحاكاة عبر أربعة مستويات:
  1. تطور "ماتريكس-تروتر" (Matrix-Trotter) في الفضاء الجزئي الفيزيائي (مرجع تحليلي).
  2. محاكاة حالة المتجه (statevector) لـ Qiskit بدون ضجيج (للتحقق من تفكيك بالي).
  3. محاكاة Qiskit Aer ذات الدفعات المحدودة (لتقييم التقلبات الإحصائية).
  4. التنفيذ على أجهزة صلبة ضحلة (shallow hardware execution) على معالجات IBM الكمومية.
• اختبار جدوى الأجهزة: نفذ العمل كتلة دائرة ضحلة ممثلة ($N=1$) على خلفية ibm_fez (IBM Quantum Heron r2) لتوصيف الضجيج، والتسرب (leakage)، وفعالية تخفيف الخطأ في هذا السياق المحدد.

النتائج

• اتساق المحاكي: في البيئات الخالية من الضجيج، تعيد محاكاة "ماتريكس-تروتر" ومحاكاة حالة المتجه باستمرار المرجع التحليلي للانتقال المفاجئ لعدد الجسيمات المتأخر، $n_k = 1/[4(k\eta_e)^4]$، ضمن نظام الإثارة الواحدة المحكوم ($n_k \ll 1$، أو $x = |k\eta_e| \gtrsim 2.2$).
• إحصائيات الدفعات المحدودة: تُظهر محاكاة الدفعات المحدودة على محاكي Aer أن التقلبات الإحصائية تنخفض مع زيادة عدد الدفعات، متبعةً الاتجاه المتوقع، مع حدوث انحرافات أساسية عند قيم $x$ الكبيرة حيث يكون إنتاج الجسيمات منخفضاً.
• أداء الأجهزة:
  • نجحت تجربة الأجهزة لـ $N=1$ في تنفيذ كتلة الدائرة الضحلة.
  • ومع ذلك، أظهرت النتائج خطأً متبقياً في الأجهزة من رتبة $10^{-2}$، تهيمن عليه أخطاء القراءة وضجيج البوابات.
  • تم تطبيق تقنيات تخفيف الخطأ، وتحديداً "تخفيف القياس الخالي من المصفوفة" (M3) و"الاستقراء صفري الضجيج" (ZNE). قلل M3 من التسرب من ~17% إلى ~14% لكن كان تأثيره ضئيلاً على عدد الجسيمات المستخلص. وأظهرت تقديرات ZNE عدم استقرار وعدم يقين كبيرين.
  • ظلت أعداد الجسيمات المقاسة متجمعة حول أرضية الضجيج ($10^{-2}$)، وهي أعلى بكثير من إشارة $N=1$ المثالية ($\sim 2.5 \times 10^{-3}$)، وبعيدة عن المرجع التحليلي الكامل.

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

• التحقق: يكمن النجاح الرئيسي في الاتساق الداخلي للمحاكاة في المحاكيات الكمومية، مما يؤكد أن البناء المعتمد على الزمن يعيد إنتاج المرجع التحليلي بشكل صحيح في النظام المقطوع.
• القيود الهاردويرية: تنص الورقة صراحة على أن إعادة البناء الكمي لطيف الجسيمات على أجهزة NISQ الحالية يظل خارج نطاق الأداء الحالي. لا تُعتبر تجربة IBM كتحقق من خلق الجسيمات الكونية، بل كاختبار جدوى لكتلة الدائرة المجمعة وكأداة تشخيص للضجيج والتسرب.
• الآفاق المستقبلية: يقترح المؤلفون أن نموذج الانتقال المفاجئ يعمل كمرجع قابل للتطبيق، مع وجود خطوات تالية طبيعية تشمل عمليات الاستكمال السلس (smooth interpolations) حيث لا تتوفر المعاملات التحليلية، وزيادة قطع فضاء فو (Fock space cutoff) للتعامل مع أعداد الإشغال الأعلى، بشرتماد توفر الأجهزة المستقبلية لعمق الدوائر المطلوب.

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