The solution of linear systems of equations is the basis of many other quantum algorithms, and recent results provided an algorithm with optimal scaling in both the condition number $κ$ and the allowable error $ε$ [PRX Quantum \textbf{3}, 040303 (2022)]. That work was based on the discrete adiabatic theorem, and worked out an explicit constant factor for an upper bound on the complexity. Here we show via numerical testing on random matrices that the constant factor is in practice about 1,200 times smaller than the upper bound found numerically in the previous results. That means that this approach is far more efficient than might naively be expected from the upper bound. In particular, it is about an order of magnitude more efficient than using a randomised approach from [arXiv:2305.11352] that claimed to be more efficient.
- معرّف الورقة: 2312.07690
- العنوان: حل النظام الخطي الكمي الثابت المنفصل له عوامل ثابتة أقل من حل الثابت العشوائي
- المؤلفون: Pedro C. S. Costa, Dong An, Ryan Babbush, Dominic W. Berry
- التصنيف: quant-ph (الفيزياء الكمية)
- وقت النشر: ديسمبر 2023، أحدث نسخة 11 أكتوبر 2025
- رابط الورقة: https://arxiv.org/abs/2312.07690
يعتبر حل الأنظمة الخطية أساساً لعديد من الخوارزميات الكمية. قدمت الأبحاث الحديثة خوارزميات ذات تعقيد أمثل من حيث رقم الشرط κ والخطأ المسموح به ε PRX Quantum 3, 040303 (2022). يستند هذا العمل على النظرية الثابتة المنفصلة ويقدم عوامل ثابتة صريحة لحدود التعقيد. تُظهر هذه الورقة من خلال الاختبارات الرقمية على مصفوفات عشوائية أن العوامل الثابتة العملية أصغر بحوالي 1200 مرة من الحدود العليا المحسوبة رقمياً في النتائج السابقة. هذا يعني أن الطريقة أكثر كفاءة بكثير مما كان متوقعاً بشكل ساذج من الحدود العليا. وعلى وجه الخصوص، فهي أكثر كفاءة بحوالي رتبة حجم واحدة من الطريقة العشوائية التي تدّعي كفاءة أعلى arXiv:2305.11352.
تعتبر مشكلة النظام الخطي الكمي (QLSP) من المشاكل الأساسية في الحوسبة الكمية، وتهدف إلى إنتاج حالة كمية |x⟩ متناسبة مع حل النظام الخطي Ax = b بكفاءة. يشكل حل هذه المشكلة أساساً لعديد من الخوارزميات الكمية الأخرى، بما في ذلك الخوارزميات في مجالات التعلم الآلي الكمي والتحسين الكمي.
- خوارزمية HHL: اقترح Harrow و Hassidim و Lloyd أول خوارزمية للنظام الخطي الكمي بتعقيد O(poly(n)κ²/ε)
- طرق الحوسبة الكمية الثابتة: قدمت الأبحاث اللاحقة تحسينات بناءً على الحوسبة الكمية الثابتة (AQC)، وخاصة Costa وآخرون في 6 الذين حققوا تعقيداً أمثل O(κlog(1/ε)) بناءً على النظرية الثابتة المنفصلة
- الطرق العشوائية: تستخدم طريقة أخرى تطوراً زمنياً عشوائياً لمحاكاة تأثير Zeno الكمي
على الرغم من أن الطريقة الثابتة المنفصلة تحقق تعقيداً تقاربياً أمثل من الناحية النظرية، فإن الحد الأعلى الصارم يعطي عاملاً ثابتاً α = 2305 كبيراً جداً، مما يثير تساؤلات حول كفاءتها العملية. في الوقت نفسه، تدّعي الطريقة العشوائية أنها أكثر كفاءة عملياً بسبب حد أعلى أكثر إحكاماً. تهدف هذه الورقة إلى التحقق من الأداء العملي لهذه الطرق من خلال التجارب الرقمية.
- الكشف عن العامل الثابت العملي للطريقة الثابتة المنفصلة: من خلال تجارب رقمية واسعة النطاق، تم اكتشاف أن العامل الثابت العملي α = 1.84، وهو أصغر بحوالي 1200 مرة من الحد النظري
- إثبات التفوق العملي للطريقة الثابتة المنفصلة: تُظهر الاختبارات الرقمية أن طريقة المشي الكمي أكثر كفاءة بحوالي 7 مرات من الطريقة العشوائية في المتوسط
- توفير مقارنة شاملة للأداء: تشمل حالات المصفوفات المحددة الموجبة والمصفوفات غير الهرميتية العامة، مع اختبارات لأبعاد وأرقام شرط مختلفة
- الأخذ في الاعتبار نفقات الخوارزمية الكاملة: تحليل التعقيد الإجمالي بما في ذلك خطوة التصفية، مما يثبت أن طريقة الثابت المنفصل لا تزال تحقق تحسناً بما لا يقل عن 3 مرات حتى مع الأخذ في الاعتبار جميع النفقات
بالنظر إلى مصفوفة قابلة للعكس A ∈ C^(N×N)، حيث ∥A∥ = 1، وناقل معياري |b⟩ ∈ C^N، الهدف هو تحضير حالة معيارية |x̃⟩ بحيث تقترب من |x⟩ = A^(-1)|b⟩/∥A^(-1)|b⟩∥ بدقة ∥|x̃⟩ - |x⟩∥ ≤ ε.
بالنسبة للمصفوفة الهرميتية المحددة الموجبة A، يتم بناء هاميلتونيان:
- H₀ := (0 Qb; Qb 0)
- H₁ := (0 AQb; QbA 0)
حيث Qb = I_N - |b⟩⟨b| هو مؤثر الإسقاط.
يُعطى هاميلتونيان المعتمد على الزمن بـ:
H(s) = (1 - f(s))H₀ + f(s)H₁
يتم تصميم دالة الجدولة f(s) وفقاً لشرط الفجوة الطاقية:
f(s) = κ/(κ-1)1 - (1 + s(κ^(p-1) - 1))^(1/(1-p))
يتم التحويل إلى شكل هرميتي بمضاعفة بُعد المصفوفة:
A := (0 A; A† 0)
هاميلتونيان المقابل هو:
H(s) = (0 A(f(s))Qb; QbA(f(s)) 0)
حيث A(f) := (1-f)σz⊗I_N + fA.
تطبق الطريقة العشوائية العملية التالية:
e^(-it_q H(v_q)) ⋯ e^(-it₂H(v₂))e^(-it₁H(v₁))
حيث:
- vⱼ = vₐ + j(vb - vₐ)/q هي نقاط زمنية منفصلة
- tⱼ متغيرات عشوائية يتم اختيارها وفقاً لتوزيع احتمالي محدد
دالة كثافة الاحتمال هي:
p_j(t) ∝ (J_p(Δ(vⱼ)|t|/2)/(Δ^(p-1)(vⱼ)|t|^p))²
حيث J_p هي دالة بيسل من النوع الأول، p = 1.165.
- الأبعاد: مصفوفات عشوائية 4×4، 8×8، 16×16
- أرقام الشرط: κ = 10, 20, 30, 40, 50
- أنواع المصفوفات: مصفوفات هرميتية محددة موجبة ومصفوفات غير هرميتية عامة
- عدد العينات: 100 مصفوفة عشوائية مستقلة لكل رقم شرط
- طريقة المشي الكمي: عدد خطوات المشي المطلوبة للوصول إلى خطأ الهدف Δ = 0.4
- الطريقة العشوائية: إجمالي وقت التطور المطلوب ∑|tⱼ| للوصول إلى نفس الخطأ
- تعريف الخطأ: خطأ معيار 2 ∥|x̃⟩ - |x⟩∥
- معامل دالة الجدولة p = 1.4
- استخدام المتوسط الهندسي لحساب التعقيد
- تضمين نطاق الربع والنطاق الكامل لأشرطة الخطأ
- 200 تكرار لكل مثيل من الطريقة العشوائية
بالنسبة لحالة κ = 50:
- الحد النظري: α = 2305
- القياس الفعلي: α = 1.84
- مضاعف التحسن: حوالي 1250 مرة
| رقم الشرط κ | خطوات QW | خطأ QW | وقت RM | خطأ RM | مضاعف التحسن |
|---|
| 10 | 36 | 0.348 | 281 | 0.395 | 7.81 |
| 20 | 76 | 0.381 | 604 | 0.397 | 7.95 |
| 30 | 120 | 0.400 | 963 | 0.398 | 8.03 |
| 40 | 176 | 0.399 | 1330 | 0.397 | 7.56 |
| 50 | 232 | 0.397 | 1722 | 0.399 | 7.42 |
استخدام مصفوفات فعلية من مجموعة SuiteSparse:
- مشكلة الرسم البياني الموجه (ID=168, κ=4.041×10²): QW أسرع من RM بـ 9.58 مرات
- مشكلة محاكاة الدوائر (ID=1199, κ=6.302×10⁵): QW أسرع من RM بـ 457 مرة
تُظهر التجارب أن اعتماد التعقيد على بُعد المصفوفة ضئيل جداً، وهو يتوافق مع التوقعات النظرية، لأن التعقيد يعتمد بشكل أساسي على رقم الشرط وليس على البُعد.
مع الأخذ في الاعتبار التعقيد الإجمالي بعد خطوة التصفية:
- بالنسبة لأخطاء الهدف النموذجية ε > 0.0004، يهيمن الجزء الثابت
- تحتفظ طريقة QW بميزة كبيرة على طريقة RM حتى مع تضمين نفقات التصفية
- حد الخطأ الأمثل Δ حوالي 0.4، وهو يتوافق مع إعداد التجربة
- خوارزمية HHL: عمل رائد، لكن التعقيد O(κ²/ε)
- تضخيم السعة الزمنية المتغيرة: تحسين الاعتماد على الدقة
- طرق الحوسبة الكمية الثابتة: توفير تحجيم تعقيد أفضل
- التصفية متعددة الحدود الأمثل: تحسين إضافي للتنفيذ
أثبت Harrow و Kothari أن الحد الأدنى لمشكلة النظام الخطي الكمي هو Ω(κlog(1/ε))، وهو الحد الذي تم تحقيقه لأول مرة في عمل Costa وآخرون.
اقترح Subaşı وآخرون طريقة عشوائية بتعقيد O(κlog(κ/ε))، وعلى الرغم من وجود عامل log κ إضافي، إلا أنهم يدّعون أنها أكثر كفاءة عملياً بسبب عوامل ثابتة أصغر.
- الفجوة الضخمة بين النظرية والممارسة: العامل الثابت العملي للطريقة الثابتة المنفصلة أصغر بـ 1200 مرة من الحد النظري
- تفوق الطريقة: طريقة المشي الكمي تتفوق على الطريقة العشوائية في جميع حالات الاختبار
- القيمة العملية: الطريقة ليست فقط أمثل من الناحية النظرية بل أيضاً عالية الكفاءة عملياً
- حد الخطأ: استخدام خطأ هدف نسبي كبير Δ = 0.4، قد يؤدي إلى حالات شاذة معينة
- نوع المصفوفة: الاختبارات تركز بشكل أساسي على مصفوفات عشوائية، قد تظهر مصفوفات منظمة في التطبيقات الفعلية أداء مختلفة
- عدالة المقارنة: مقارنة RM تتعلق بوقت التطور وليس بعدد بوابات كمية محددة
- تحليل أكثر دقة للعوامل الثابتة: تطوير حدود نظرية أكثر إحكاماً
- المصفوفات المنظمة: دراسة أداء المصفوفات ذات البنية الخاصة
- التنفيذ على الأجهزة: التحقق من هذه النتائج على أجهزة كمية فعلية
- قيمة عملية عالية: حل الفجوة المهمة بين النظرية والممارسة
- تجارب شاملة: اختبارات واسعة النطاق على مصفوفات عشوائية وحالات تطبيقية فعلية
- تحليل شامل: الأخذ في الاعتبار نفقات الخوارزمية الكاملة بما في ذلك خطوة التصفية
- نتائج موثوقة: جميع حالات الاختبار تُظهر ميزة متسقة
- شرح نظري غير كافٍ: عدم وجود تحليل عميق لسبب صغر العامل الثابت العملي جداً
- معايير المقارنة: قد تكون المقارنة مع طريقة RM غير مباشرة كفاية (الوقت مقابل الخطوات)
- حدود الحجم: أبعاد المصفوفات المختبرة نسبياً صغيرة (الحد الأقصى 16×16)
لهذا العمل تأثير مهم على مجتمع الخوارزميات الكمية:
- إعادة تقييم كفاءة الخوارزمية: تذكير الباحثين بأن الحدود النظرية قد تكون متحفظة جداً
- إرشادات اختيار الخوارزمية: توفير توصيات واضحة لاختيار الخوارزمية للتطبيقات العملية
- اتجاهات البحث المستقبلية: إثارة الحاجة إلى تحليل تعقيد أكثر دقة
هذه الطريقة مناسبة بشكل خاص لـ:
- الخوارزميات الكمية التي تتطلب حل خطي عالي الدقة
- المشاكل العملية برقم شرط معتدل
- سيناريوهات التطبيق الحساسة للعوامل الثابتة
تستشهد هذه الورقة بالمراجع الرئيسية في مجال حل النظام الخطي الكمي، بما في ذلك:
- 1 خوارزمية HHL الأصلية
- 6 طريقة Costa وآخرون الثابتة المنفصلة
- 10 تحسينات Jennings وآخرون للطريقة العشوائية
- 14 تحسينات Berry وآخرون لمحاكاة هاميلتونيان