Exponentially cheaper coherent phase estimation via uncontrolled unitaries

إليك ملخص تقني مفصل لورقة "تقدير الطور المتماسك الأرخص أسياً عبر الوحدات غير المتحكم بها" (Exponentially cheaper coherent phase estimation via uncontrolled unitaries) للباحث ميركو أميكو (Mirko Amico).

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

يُعد تقدير الطور الكمي (QPE) إجراءً فرعياً أساسياً في العديد من الخوارزميات الكمية (مثل خوارزمية شور، وخوارزمية HHL، وتقدير طاقة الحالة الأرضية). يعتمد بروتوكول QPE القياسي على ارتداد الطور (phase kickback)، حيث تقوم عملية "وحدة متحكم بها" ($C\text{-}U$) بطبع طور الحالة الذاتية على كيوبت مساعد (ancilla).

التحدي الجوهري:

• تكلفة الموارد العالية: تنفيذ بوابة $C\text{-}U$ هو أمر أغلى أسياً من تطبيق $U$ غير المتحكم به. إذا تم تفكيك $U$ إلى $n_1$ من البوابات أحادية الكيوبت و $n_2$ من بوابات الكيوبت الثنائية، فإن جعلها "متحكم بها" يتطلب عادةً تحويل كل بوابة أحادية الكيوبت إلى نسخة متحكم بها (إضافة CNOTs)، وكل بوابة ثنائية الكيوبت إلى هيكل يشبه بوابة Toffoli (إضافة حوالي 6 بوابات CNOT).
• عمق الدائرة: للحصول على دقة قدرها $m$ بت، يتطلب QPE القياسي تطبيق $C\text{-}U^{2^k}$ من أجل $k=0$ إلى $m-1$. إجمالي عدد بوابات الكيوبت الثنائية يتدرج كـ $O(2^m \cdot n_2)$، حيث يضرب العامل الأسي تكلفة التحكم في الوحدة بأكملها.
• قيود أجهزة NISQ: في الأجهزة الكمية ذات المقياس المتوسط المشوبة بالضجيج (NISQ)، تعاني الدوائر العميقة التي تحتوي على العديد من البوابات المتحكم بها من فك الترابط والضجيج، مما يجعل QPE القياسي غير قابل للتنفيذ غالباً. البدائل الموجودة (التقدير الإحصائي للطور) تضحي بالتماسك، مما يفقد القدرة على استخدام QPE كإجراء فرعي في خوارزميات كمية أكبر.

2. المنهجية: ارتداد الطور غير المتحكم به

تقترح الورقة إجراءً جديداً يسمى "ارتداد الطور غير المتحكم به" (Uncontrolled Phase Kickback). يحقق هذا الأسلوب نفس تراكم الطور المتماسك لـ QPE القياسي، ولكنه يستبدل عمليات $C\text{-}U$ المكلفة بتسلسل من تحضيرات الحالة المتحكم بها والوحدات غير المتحكم بها.

الافتراضات الرئيسية:

1. الحالة المرجعية: توجد حالة ذاتية مرجعية $|\phi\rangle$ للوحدة $U$ ذات قيمة ذاتية معروفة $e^{i2\pi\phi}$.
2. دائرة التحضير: توجد وحدة $W$ تقوم بنقل الحالة المرجعية إلى الحالة الذاتية المستهدفة: $W|\phi\rangle = |\psi\rangle$، حيث $|\psi\rangle$ هي الحالة الذاتية ذات الطور المجهول $\theta$ الذي نرغب في تقديره.
3. عدم التبادل: $W$ لا تبدل مع $U$.

أداة الدائرة (الكيوبت المساعد الواحد):
بدلاً من تطبيق $C\text{-}U$، تنفذ الخوارزمية التسلسل التالي على كيوبت مساعد (في البداية في الحالة $|+\rangle$) وعلى النظام (في البداية في الحالة $|\phi\rangle$):

1. المتحكم بـ $W$ (Controlled-$W$): إذا كان الكيوبت المساعد هو $|1\rangle$، يتم تطبيق $W$ على النظام ($|\phi\rangle \to |\psi\rangle$). إذا كان $|0\rangle$، لا يتم فعل شيء.
   • الحالة: $\frac{1}{\sqrt{2}}(|0\rangle|\phi\rangle + |1\rangle|\psi\rangle)$.
2. غير المتحكم بـ $U$ (Uncontrolled-$U$): تطبيق $U$ على النظام.
   • الحالة: $\frac{1}{\sqrt{2}}(e^{i2\pi\phi}|0\rangle|\phi\rangle + e^{i2\pi\theta}|1\rangle|\psi\rangle)$.
3. المتحكم بـ $W^\dagger$ المفتوح (Controlled-$W^\dagger$ - Open-Controlled): إذا كان الكيوبت المساعد $|0\rangle$، يتم تطبيق $W$ (أو إذا كان $|1\rangle$ يتم تطبيق $W^\dagger$ اعتماداً على تصميم الكتلة المحدد) لإعادة مسارات النظام إلى حالة مشتركة.
   • الآلية: المسار $|0\rangle$ (الحالي $|\phi\rangle$) يتم نقله إلى $|\psi\rangle$ عبر $W$. المسار $|1\rangle$ (الحالي $|\psi\rangle$) يبقى $|\psi\rangle$.
   • الحالة الناتجة: $\frac{1}{\sqrt{2}}(e^{i2\pi\phi}|0\rangle + e^{i2\pi\theta}|1\rangle) \otimes |\psi\rangle$.
4. فك الترابط (Disentanglement): النظام الآن في الحالة الذاتية $|\psi\rangle$ بغض النظر عن حالة الكيوبت المساعد. الطور النسبي $(\theta - \phi)$ يقع بالكامل على الكيوبت المساعد.

التوسع لتقدير $m$ بت:
للدقة متعددة البتات، تُكرر هذه الأداة $m$ من المرات.

• الكيوبتات المساعدة الوسيطة: بعد خطوة $U^{2^k}$ غير المتحكم بها، يتم تطبيق 1-controlled-$W^\dagger$. هذا يفك ترابط الكيوبت المساعد الحالي ويعيد النظام إلى الحالة المرجعية $|\phi\rangle$ من أجل التكرار التالي.
• الكيوبت المساعد النهائي: يمكن استخدام $W$ مفتوح التحكم لترك النظام في الحالة المستهدفة $|\psi\rangle$ إذا لزم الأمر لخطوات الخوارزمية اللاحقة.

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

• تعريف الإجراء: صياغة "ارتداد الطور غير المتحكم به" كبديل جاهز للاستخدام لبوابات $C\text{-}U$ في الخوارزميات التي تتطلب تقدير طور متماسك.
• تحليل الموارد: إثبات أن هذا النهج يحقق تخفيضاً أسياً في عدد بوابات الكيوبت الثنائية (CNOTs) مقارنة بـ QPE القياسي في الحد ذي الصلة.
• الحفاظ على التماسك: على عكس الطرق الإحصائية، يحافظ هذا الأسلوب على التماسك الكمي الكامل، مما يسمح باستخدام تحويل فورييه الكمي العكسي (IQFT) والاندماج في خوارزميات كمية أكبر.
• القابلية للتعميم: إثبات أن التقنية تنطبق على أي خوارزمية تستخدم ارتداد الطور، بشرط وجود حالة ذاتية مرجعية ودائرة تحضير $W$ فعالة.

4. النتيات وتقديرات الموارد

تقدم الورقة مقارنة كمية لعدد بوابات الكيوبت الثنائية ($GateCount_{2q}$) لتقدير الطور بـ $m$ بت.

• تكلفة QPE القياسية:
  $$GateCount_{2q}^{std} \approx (2^m - 1)(2n_1 + 6n_2)$$
  حيث $n_1$ و $n_2$ هما بوابات أحادية وثنائية الكيوبت في $U$. العامل الأسي $2^m$ يضرب تكلفة التحكم في كل بوابة في $U$.

• تكلفة QPE غير المتحكم بها:
  $$GateCount_{2q}^{unctrl} \approx m(4n_{W,1} + 12n_{W,2}) + (2^m - 1)n_2$$
  هنا، العامل الأسي $2^m$ يضرب فقط تكلفة الكيوبت الثنائي المجرد لـ $U$ ($n_2$)، لأن $U$ يُطبق بشكل غير متحكم به. العمليات المتحكم بها تقتصر على $W$ (و $W^\dagger$) والتي تتدرج خطياً مع $m$.

• الميزة الأسية:
  في الحد الذي تهيمن فيه البوابات أحادية الكيوبت على $U$ (وهو أمر شائع في نماذج Trotter، أي $n_2 \to 0$)، تصبح النسبة بين التكاليف:
  $$\frac{GateCount_{2q}^{std}}{GateCount_{2q}^{unctrl}} \approx \frac{2^m \cdot 2n_1}{m \cdot 4n_{W,1}} \propto \frac{2^m}{m}$$
  يمثل هذا تخفيضاً أسياً في عمق الدائرة وعدد بوابات الكيوبت الثنائية.

5. التطبيقات والأمثلة

توضح الورقة الطريقة من خلال مثالين محددين:

1. طاقة الحالة الأرضية لنموذج هايزنبرغ (Heisenberg Model):

   • السياق: تقدير طاقة الحالة الأرضية لسلسلة سبين.
   • الإعداد: الحالة المرجعية $|\phi\rangle = |0\dots0\rangle$ هي حالة ذاتية ذات طور معروف (يتم تحديده بواسطة حدود $ZZ$). الحالة المستهدفة هي $|\psi\rangle$ (الحالة الأرضية).
   • الفائدة: عملية $C\text{-}U$ (التطور الزمني) مكلفة للغاية بسبب العديد من التفاعلات ثنائية الكيوبت. بينما $C\text{-}W$ (تحضير الحالة الأرضية من $|0\rangle$) غالباً ما يكون تسلسلاً بسيطاً من دورات الكيوبت المنفردة أو بوابات Pauli-X، مما يقلل بشكل كبير من العبء.

2. خوارزمية شور للتحليل إلى عوامل (إيجاد الرتبة):

   • السياق: إيجاد الرتبة $r$ لعدد صحيح $a$ بمقياس $N$.
   • الإعداد: الوحدة هي الضرب النمطي $M_a$. الحالة المرجعية $|\phi\rangle = |0\dots0\rangle$ (أو $|1\rangle$ في ترميزات معينة) لها طور معروف. $W$ هي مجرد بوابة قلب البت (X gate) لتحضير الحالة المدخلة.
   • القيود: في خوارزمية شور، تكون الحالة المدخلة غالباً عبارة عن تراكب من الحالات الذاتية. تشير الورقة إلى أنه بالنسبة للتراكبات، يمكن تطبيق الطريقة غير المتحكم بها فقط على البت الأول من سجل الطور (حيث يمكن إعادة ضبط النظام)، بينما قد تظل البتات اللاحقة تتطلب عمليات متحكم بها قياسية ما لم يتم العثست على مسار تحضير حالة ذاتية محدد.

6. الأهمية والخلاصة

• الجدوى لأجهزة NISQ: من خلال إزالة الحاجة إلى دوائر "وحدة متحكم بها" عميقة، تجعل هذه الطريقة تقدير الطور المتماسك عالي الدقة ممكناً على الأجهزة الحالية ذات أوقات التماسك المحدودة.
• الأثر الخوارزمي: تسمح هذه الطريقة لخوارزميات مثل HHL، وتقدير السعة (amplitude estimation)، ومحاكاة الكيمياء الكمية، بالعمل بدوائر أقل عمقاً بكثير، مما قد يمكّن من حل مشكلات كانت تُعتبر سابقاً أعمق من أن تستوعبها الأجهزة الحالية.
• المقايضات: تتطلب الطريقة معرفة مسبقة بحالة ذاتية مرجعية ودائرة $W$ فعالة لنقلها إلى الحالة المستهدفة. إذا كانت الحالة المستهدفة هي تراكب معقد من حالات ذاتية دون مسار تحضير معروف، فإن إمكانية تطبيق الطريقة تكون محدودة.

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