تقدم هذه الورقة البحثية محاكاة هاميلتونيان عشوائية مسارية (PRHS)، وهي خوارزمية مبتكرة تمدد طريقة qDrift إلى رتب تعسفية عبر توظيف شرائح زمنية مترابطة وتوزيعات احتمالية شبه مقدارية، مما يحقق تعقيد استعلام دون حدود متعددة الحدود بالنسبة للدقة المستهدفة دون الحاجة إلى كيوبتات مساعدة، وتتفوق بشكل كبير على طريقة qDrift القياسية في المحاكاة العددية.
في مجال الحوسبة الكمومية، تعد إحدى المهام الأكثر جوهرية هي محاكاة كيفية تغير نظام فيزيائي بمرور الوقت. تخيل جزيئاً، وهو تجمع معقد من الذرات والإلكترونات، يغير شكله أو يتفاعل مع بيئة جديدة. للتنبؤ بهذا السلوك، يجب على العلماء حساب كيفية دفع تطوره بواسطة طاقته، والتي توصف بكائن رياضي يسمى "الهاميلتوني" (Hamiltonian). هذا ليس مجرد تمرين نظري؛ بل هو المحرك وراء العديد من الخوارزميات الكمومية القوية المستخدمة في الكيمياء، والتحسين، وحل المعادلات المعقدة. وتكمن الصعوبة في أن أوصاف الطاقة هذه غالباً ما تتكون من أجزاء مختلفة ومتفاعلة لا تكتفي بمجرد الجمع البسيط. وعندما لا تتعاون هذه الأجزاء، يصبح حساب التغيير الإجمالي صعباً للغاية، مما يتطلب تقريباً خطوة بخطوة يمكن أن يصبح سريعاً جداً أو عرضة للأخطاء بالنسبة للأجهزة الحالية.
لسنوات، اعتمد الباحثون على طريقة تسمى تفكيك "تروتر-سوزوكي" (Trotter–Suzuki decomposition) لتقسيم هذه المشكلة. يقوم هذا النهج بتقطيع الوقت إلى فترات زمنية ضئيلة وتطبيق قواعد الطاقة لكل جزء بشكل متتالٍ. ورغم فعاليته، إلا أن لهذه الطة عيباً كبيراً: فكلما زاد عدد أجزاء النظام، زادت الخطوات المطلوبة، وكلما وجب أن تكون النتيجة أكثر دقة، استهلكت الطريقة موارد حاسوبية أكبر. ولتجاوز ذلك، تم تقديم تقنية أحدث تسمى "qDrift". فبدلاً من اتباع تسلسل مرتب وصارم، تختار "qDrft" جزءاً واحداً من الطاقة عشوائياً لتطبيقه في كل مرة. وهذا يلغي العقوبة المرتبطة بوجود أجزاء كثيرة، مما يجعل العملية أكثر كفاءة. ومع ذلك، فإن لـ "qDrift" حدها الخاص: فهي بطبيعتها طريقة من الدرجة الأولى، مما يعني أن دقتها تتحسن خطياً فقط مع إضافة المزيد من الخطوات. وللحصول على نتيجة عالية الدقة، لا تزال بحاجة إلى عدد هائل من الخطوات، مما يستنزف المكاسب في الكفاءة.
لقد قدم باحث الآن نهجاً جديداً يسمى "محاكاة الهاميلتوني المساري العشوائي" (Pathwise Random Hamiltonian Simulation أو PRHS)، والذي يتغلب على هذا القصور دون الحاجة إلى أجهزة إضافية. الفكرة الجوهرية هي أخذ الطبيعة العشوائية لـ "qDrift" وإضافة طبقة من التنسيق الذكي. فبدلاً من اختيار قاعدة طاقة واحدة لكل شريحة زمنية بشكل مستقل، تقوم الطريقة الجديدة بتقسيم كل خطوة زمنية إلى سلسلة من الشرائح الأصغر المترابطة. وضمن الشريحة الواحدة، يكون اختيار أي قاعدة طاقة سيتم تطبيقها في الشريحة الأولى مرتبطاً رياضياً بالاختيارات في الشرائح اللاحقة. هذا الارتباط ليس عشوائياً بالمعنى الفوضوي؛ بل تحكمه توزيعات رياضية فريدة ومحددة صاغها المؤلف وأثبت أنها الوحيدة من نوعها. ومن خلال نسج هذه الاختيارات معاً بعناية، تقوم الطريقة بإلغاء الأخطاء التي قد تتراكم لولا ذلك، مما يسمح للمحاكاة بالوصول إلى درجة أعلى بكثير من الدقة.
وقد أظهر الباحث أن هذه الطريقة يمكن ضبطها للعمل في نظامين متميزين. في الحالات التي تعمل فيها المحاكاة لفترة طويلة، تعود الاستراتيجية المثلى طبيعياً إلى بروتوكول "qDrift" القياسي، مما يؤكد أن الطريقة الجديدة هي تعميم حقيقي للطريقة القديمة. ومع ذلك، في النظام الذي تكون فيه الدقة العالية هي الهدف، تتصرف الط {الطريقة} بشكل مختلف. فمن خلال زيادة عدد الشرائح المترابطة داخل كل خطوة، تنمو تكلفة تحقيق دقة معينة ببطء أكبر بكثير من أي طريقة سابقة. وفي الواقع، أظهر الباحث أن عدد الخطوات المطلوبة للوصول إلى دقة مستهدفة ينمو بشكل أبطأ من أي قوة لقلب تلك الدقة. وهذا يعني أنه بالنسبة للدقة العالية جداً، تصبح الطريقة الجديدة أكثر كفاءة بشكل ملحوظ من الأساليب القياسية، وكل ذلك مع بقائها خالية من "الكيوبتات المساعدة" (ancilla qubits) الإضافية التي تتطلبها التقنيات عالية الدقة الأخرى.
للتحقق من هذه التنبؤات النظرية، أجرى الباحث محاكاة عددية على خمسة أنظمة جزيئية مختلفة، بما في ذلك الهيدروجين، وهيدريد الهيليوم، والماء. وقارن أداء طريقته الجديدة ببروتوكول "qDrift" القياسي، مع ضمان استخدام نفس القدر من الجهد الحسابي، مقاساً بعدد العمليات الأساسية المنجزة. كانت النتائج مذهلة؛ ففي هذه المحاكاة، حققت الط {الطريقة} الجديدة دقة أفضل بمقدار اثنين إلى أربعة مراتب عشرية من "qDrift" لنفس التكلفة. وبينما تفرض الطريقة عبئاً إحصائياً لأنها تعتمد على توزيع يتضمن قيماً سالبة - وهي ميزة تتطلب إعادة قياس النتائج النهائية - إلا أن المكسب في الدقة كان كبيراً لدرجة أنه فاق التكلفة الإضافية بكثير. تؤكد الدراسة أنه من خلال إدخال ارتباط محدد ومحسوب بين الاختيارات العشوائية، يمكن دفع حدود ما يمكن أن تحققه عمليات المحاكاة الكمومية العشوائية، مما يوفر أداة قوية للكيمياء الكمومية عالية الدقة دون الحاجة إلى موارد أجهزة إضافية.
ملخص تقني: محاكاة هاميلتونيان عشوائية عبر المسارات (PRHS)
بيان المشكلة تعد محاكاة الهاميلتونيان (Hamiltonian simulation) عملية أولية أساسية في الحوسبة الكمومية، وهي ضرورية لخوارزميات مثل تقدير الطور، والتحسين، وحل الأنظمة الخطية. يعتمد النهج القياسي على تقنيات تقسيم المؤثرات (operator-splitting)، مثل تفكيكات "تروتر-سوزوكي" (Trotter–Suzuki)، والتي تقرب مؤثر التطور e−iHt عبر تركيب عمليات تطور ناتجة عن حدود هاميلتونيان أبسط. وبينما تعد هذه الطرق بسيطة مفاهيميًا ولا تتطلب كيو بتات مساعدة (ancillary qubits)، إلا أن صيغ "تروتر-سوزوكي" الحتمية تظهر اعتمادًا حدوديًا على عدد حدود الهاميلتونيان (Γ) وعلى مقلوب الدقة المستهدفة (1/ϵ).
ولمعالجة الاعتماد على Γ، تم تقديم بروتوكول qDrift كبديل عشوائي. يقوم qDrift بأخذ عينات من حدود الهاميلتونيان من توزيع احتمالي كلاسيكي، مما يلغي الاعتماد الصريح على Γ. ومع ذلك، فإن qDrift محدود جوهريًا بالدقة من الدرجة الأولى في زمن التطور، مما يعني أن تعقيد الاستعلام الخاص به يظل خطيًا في 1/ϵ. وقد عانت المحاولات السابقة لتوسيع الطرق العشوائية إلى درجات أعلى من مقايضات كبيرة:
يحقق qSWIFT دقة من درجات تعسفية ولكنه يتطلب كيو بت مساعدًا (ancilla qubit).
يظل Composite qDrift خاليًا من الكيو بتات المساعدة، ولكنه يعيد إدخال الاعتماد على Γ.
يحقق qSHIFT درجات أعلى، ولكنه يتطلب حل نظام من المعادلات الكلاسيكية ينمو بسرعة مع الدرجة المستهدفة.
وبناءً على ذلك، ظل بناء تعميم لـ qDrift من درجة تعسفية يحافظ على بساطة الخلو من الكيو بتات المساعدة، والاستقلال عن Γ، وقاعدة أخذ عينات قابلة للحساب بكفاءة، مشكلة مفتوحة.
المنهجية يقدم البحث محاكاة هاميلتونيان عشوائية عبر المسارات (PRHS)، وهي توسيع منهجي لـ qDrift إلى درجة تعسفية M في معامل زمن التطور Λt (حيث Λ=∑∣λγ∣).
أخذ العينات المترابط (Correlated Sampling): على عكس qDrift القياسي، الذي يأخذ عينات مستقلة لكل خطوة زمنية، يقوم PRHS بتقسيم كل خطوة زمنية t إلى M من الشرائح المترابطة. ويُحكم التطور بواسطة تسلسل من M من حدود الهاميلتونيان، α=(α1,…,αM)، المأخوذة من توزيع احتمالي شبه حقيقي (quasi-probability distribution)qM(α).
بناء التوزيع الفريد: يثبت المؤلف وجود توزيع احتمالي شبه حقيقي فريد qM(α) يحقق دقة من الدرجة M. وقد اشتق تعبيرًا مغلق الشكل لهذا التوزيع. يعتمد البناء على تفكيك أي متجه مؤشر α إلى أقصى سلاسل من المدخلات المتساوية المتتالية. ويُعرف التوزيع باستخدام أعداد ستيرلينج الأولى الموقعة، s(c,b)، وأوزان qDrift القياسية μγ=∣λγ∣/Λ: qM(α)=b1,…,bL≥1∑n!M!l=1∏Lbl!μglblcl!s(cl,bl) حيث n=∑bl، وتمثل (gl,cl) قيمة وعدد السلسلة رقم l في المسار.
الانحياز وحدود الخطأ: القناة الكمومية الناتلة تقرب التطور الحقيقي بانحياز يتلاشى بشكل عامي (factorially) في M. يتم تحديد خطأ مسافة الماس (diamond-distance error) لكتلة واحدة بـ: δM≤(M+1)!(2Λt)M+1e2Λt بالنسبة لإجمالي زمن التطور T=Nt المكون من N من الكتل المستقلة، يتدرج الخطأ الإجمالي كـ: δ(N,M)≤N(M+1)!(2ΛT/N)M+1e2ΛT/N.
التعامل مع الاحتمالات شبه الحقيقية: نظرًا لأن qM(α) يمكن أن تأخذ قيمًا سالبة، يستخدم البروتوكول أخذ العينات بالأهمية (importance sampling). يتم أخذ العينات من التوزيع المطلق المعياري ∣qM(α)∣/∥qM∥1، وتتم إعادة تحجيم نتائج القياس بواسطة ∥qM∥1. وهذا يؤدي إلى عامل تضخيم لعدد الطلقات (shot-number amplification factor) قدره ∥qM∥12N.
المساهمات الرئيسية والتحسين يؤكد البحث أن عيوب أخذ العينات شبه الحقيقي (زيادة عدد الطلقات) يمكن التخفيف منها من خلال التحسين المشترك لحجم الكتلة M وعدد التكرارات N.
حد الزمن الطويل: عندما ΛT→∞، تتقارب الاستراتيجية المثلى نحو M∗=1 و N∗∼2Λ2T2/ϵ، مما يستعيد بروتوكول qDrift القياسي مع تدرج فرعي غير مثالي تربيعي في الزمن.
حد الدقة العالية: عندما ϵ→0 (عند ثبات ΛT)، تنتقل الاستراتيجية المثلى إلى N∗=1 و M∗∼ln(1/ϵ)/lnln(1/ϵ). في هذا النظام، يتدرج إجمالي تكلفة الاستعلام كـ: C(N∗,M∗)∼(ϵ1)lnln(1/ϵ)lnη حيث η هو معدل خاص بالمشكلة. هذا التدرج أبطأ من أي قوة لـ 1/ϵ، مما يمثل تحسنًا كبيرًا عن التدرج الخطي لـ qDrift وصيغ "تروتر-سوزوكي" ذات الدرجة الثابتة.
النتائج تحقق المؤلف من صحة PRHS من خلال محاكاة عددية لخمسة هاميلتونيانات جزيئية (H2، HeH+، LiH، BeH2، و H2O) تم إسقاطها على الكيو بتات عبر تحويل جوردان-وينجر (Jordan–Wigner transformation).
الدقة: بالنسبة لتكلفة استعلام ثابتة (تطابق عدد قنوات التطور الأولية المطبقة)، حقق PRHS دقة أعلى بمقدار درجتين إلى أربع درجات عشرية من qDrift عبر النطاق المختبر لـ ΛT.
عبء أخذ العينات: بينما ينمو عبء أخذ العينات ∥qM∥1 أسيًا مع M، تشير النتائج العددية إلى أن هذا النمو متواضع مقارنة بالربح الأسي في الدقة. ظل ∥qM∥1 المرصود أقل بكثير من الحد العالمي النظري المشتق في البحث.
الأهمية يزعم البحث أنه يحل المشكلة المفتوحة المتمثلة في بناء بروتوكول محاكاة هاميلتونيان عشوائي، ذي درجة تعسفية، وخالٍ من الكيو بتات المساعدة، ومستقل عن Γ. يوفر PRHS مسارًا منهجيًا لمحاكاة عالية الدقة دون العبء المترتب على الكيو بتات مساعدة أو ترميزات كتل معقدة. ومن خلال تحسين عمق الارتباط (M) مقابل عدد التكرارات (N)، يتنقل PRHS بين نظام qDrift القياسي ونظام الدقة العالية حيث ينمو تكلفة الاستعلام دون كثير من القوى (sub-polynomially) بالنسبة لمقلوب الدقة. يشير المؤلف إلى أنه بينما يقتصر تحليلهم على المقارنة مع تفكيكات نوع "تروتر-سوزوكي" (باستثناء الطرق التي تقدم تدرجًا تقاربيًا أفضل مثل تحويل القيمة المفردة الكمومي QSVT والتي تأتي على حساب الكيو بتات المساعدة)، فإن PRHS يوفر بديلًا عمليًا وفعالًا من حيث الموارد للأجهزة الكمومية القريبة من الحقبة الحالية. ويُقترح في العمل المستقبلي دمج إعادة الوزن الواعي بالبنية التحتية (architecture-aware reweighting) لتقليل التكاليف الحسابية الفعلية على مستوى الأجهزة.