End-to-End PDE-Based Quantum Algorithms for Multi-Asset Option Pricing under Local and Stochastic Volatility

بيان المشكلة
تؤدي عملية تسعير الخيارات متعددة الأصول تحت نماذج التقلب المحلي (Black–Scholes) والتقلب العشوائي (Heston) بشكل طبيعي إلى معادلات تفاضلية جزئية (PDEs) ذات أبعاد عالية. تعاني الطرق العددية الكلاسيكية لحل هذه المعادلات، مثل مخططات الفروق المحدودة أو طرق الأس لـلمصفوفات، من "لعنة الأبعاد"، حيث تزداد التكلفة الحسابية بشكل أسي مع عدد الأصول ($d$) ودقة الشبكة ($N$). فبالنسبة لشبكة تحتوي على $N$ من النقاط لكل اتجاه مكاني، تتدرج الطرق الكلاسيكية القائمة على الشبكة كـ $\tilde{O}(N^{d+2})$ لنموذج Black–Scholes و $\tilde{O}(N^{2d+2})$ لنموذج Heston (حيث يتضمن نموذج Heston متغيرات التباين). وبينما تم اقتراح خوارزميات كمومية لتسعير الخيارات، فإن العديد منها يعتمد على محاكاة مسار مونت كارلو أو النهج التبايني (variational approaches) التي تواجه تحديات في إعداد الحالة، والتطبيع (normalization)، والقراءة.

المنهجية
يطور المؤلفون إطار عمل كمومي كامل (end-to-end) للمعادلات التفاضلية الجزئية يأخذ بيانات العقد والنموذج الكلاسيكية كمدخلات ويعيد تقديرات كلاسيكية لقيم الخيارات. يتكون سير العمل من ثلاث مراحل رئيسية:

1. التجزئة والصياغة (Discretization and Formulation): يتم تجزئة معادلات التسعير التفاضلية جزئياً في الفضاء باستخدام فروق محدودة من الدرجة الثانية، مما يحول المعادلة التفاضلية الجزئية المستمرة إلى نظام مكون من معادلات تفاضلية عادية (ODEs) متفرقة وذات أبعاد محدودة. وللتعامل مع الشروط الحدية غير المتجانسة وحدود المصدر، يتم تعزيز نظام الـ ODE الأفيني (affine) إلى نظام خطي متجانس عن طريق إضافة متغير حالة مساعد.
2. التطور الكمومي (Schrödingerisation): بما أن نظام الـ ODE الناتج يكون عموماً غير وحدوي (non-unitary) (بسبب الطبيعة غير الهيرميتية لمؤثر التسعير)، يستخدم إطار العمل تقنية Schrödingerisation. تقوم هذه التقنية بدمج الديناميكيات غير الوحدوية داخل نظام وحدوي أكبر عن طريق إدخال متغير مستمر مساعد ($\xi$) ورسم الديناميكيات كمعادلة شرودنجر. بعد ذلك، يتم تجزئة النظام في بُعد $\xi$ وتحويله عبر طرق فوريه لإنتاج هاميلتوني هيرميتي (Hermitian Hamiltonian).
3. إعداد الحالة والقراءة (State Preparation and Readout):
   • إعداد الحالة: يتم إعداد متجه العائد (مثل Worst-of Call) وملف القطع الناعم المطلوب لعملية الـ Schrödingerisation كحالات كمومية باستخدام ترميز الكتل (block-encodings) لمؤثرات الإحداثيات وتحويل القيمة المفردة الكمومي (QSVT).
   • التطور: يتم محاكاة ديناميكيات التسعير المجزأة عبر التطور الهاميلتوني باستخدام ترميز الكتل وتسلسلات تعديل الطور المتناوبة (Quantum Signal Processing).
   • القراءة: يتم استرداد قيم الخيارات المختارة من الحالة الكمومية النهائية. ولأن السعر المطلوب مشفر في سعة واحدة، يستخدم إطار العمل تقدير السعة الكمومي (QAE) بالاقتران مع اختبار هادامارد. ومن الأهمية بمكان أن عبء القراءة يتدرج مع الجذر التربيعي لحجم الشبكة ($2^{W/2}$)، حيث $W$ هو عدد الكيوبتات في النظام.

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

• خط إنتاج متكامل (End-to-End Pipeline): يوفر البحث مساراً كاملاً وصريحاً من البيانات الكلاسيكية المدخلة إلى المخرجات الكلاسيكية، بما في ذلك تحليل مفصل للموارد لإعداد الحالة، ومحاكاة الهاميلتوني، واستعادة التطبيع، والقراءة.
• تحليل التعقيد: يستنتج المؤلفون تعقيد البوابة الرئيسي لاسترداد سعر الخيار عند نقطة واحدة:
  • Black–Scholes متعدد الأصول: $\tilde{O}(d^2 N^{2 + d/2})$.
  • Heston متعدد الأصول: $\tilde{O}(d^2 N^{d + 2})$.
    مقارنة بالأسس الكلاسيكية للشبكات ($\tilde{O}(N^{d+2})$ و $\tilde{O}(N^{2d+2})$ على التوالي)، يقدم الإطار الكمومي تحسينات في عوامل التدرج متعددة الحدود بمقدار $N^{d/2}$ و $N^d$.
• التحقق العددي: تم التحقق من صحة الإطار عددياً مقابل الحلول التحليلية (Black–Scholes)، والمعايير شبه التحليلية (دالة Heston المميزة)، والمحللات العددية الكلاسيكية (الفروق المحدودة والأس لـلمصفوفات).
• استرداد التقلب الضمني: في حالة Heston، نجح الإطار في استرداد أسعار الخيارات عبر مجموعة من مستويات السعر (strikes) وأعاد بناء سطح التقلب الضمني غير المسطح (smile/skew)، مما يثبت قدرته على التقاط ميزات التقلب العشوائي ذات الصلة بالسوق.

النتائج

• الأداء النظري: يؤكد التحليل وجود تفوق كمومي متعدد الحدود في حجم الشبكة $N$ لكل من نموذجي Black–Schches و Heston. يصبح هذا التفوق أكثر وضوحاً مع زيادة عدد الأصول $d$، رغم أن التبعية الأسية لـ $d$ (لعنة الأبعاد) قد تم تخفيفها ولكن لم يتم القضاء عليها.
• الدقة العددية: تظهر المحاكاة أن الإطار الكمومي ينتج أسعار خيارات متوافقة مع المعايير الكلاسيكية وشبه التحليلية. وفي حالة Heston، تظهر أسطح التقلب الضمني المستردة هيكل الميل الصحيح (مثل ميل الأسهم) وتتطابق مع نماذج SSVI شبه التحليلية بفروقات نسبية تقل عن 1%.
• مصادر الخطأ: تُعزى الأخطاء الأساسية في النتائج العددية إلى التجزئة المكانية الخشنة وتقليص النطاق (domain truncation)، وليس إلى المحاكاة الكمومية نفسها. وقد وُجد أن الخطأ الإضافي الناتج عن عملية الـ Schrödingerisation هو خطأ ثانوي.

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

يحافظ المؤلفون على نبرة متواضعة فيما يتعلق بالميزة العملية:

• يوضحون أن المقارنة تتم على مستوى التدرج التقاربي (asymptotic scaling) وليس على مستوى وقت التنفيذ الفعلي، مشيرين إلى أن عدد البوابات المنطقية لا يترجم مباشرة إلى وقت تنفيذ فيزيائي دون مراعاة أعباء تصحيح الخطأ.
• يقرون بأن "لعنة الأبعاد" لم يتم القضاء عليها؛ إذ لا تزال التكلفة تزداد أسياً مع عدد الأصول، ولكن معدل النمو انخفض (من $N^d$ إلى $N^{d/2}$ في نموذج Black–Scholes).
• يشيرون إلى القيود، مثل التركيز الحالي على عوائد الخيارات من النوع الأوروبي مع الشروط الحدية المتوافقة مع الإطار (مثل Worst-of Call)، واستبعاد الميزات المعتمدة على المسار أو ميزات التمرين المبكر التي تتطلب زيادة في الحالة أو طبقات حل تكرارية.

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