We present four quantum algorithms for solving a multidimensional drift-diffusion equation. They rely on a quantum linear system solver, a quantum Hamiltonian simulation, a quantum random walk, and the quantum Fourier transform. We compare the complexities of these methods to their classical counterparts, finding that diagonalization via the quantum Fourier transform offers a quantum computational advantage for solving linear partial differential equations at a fixed final time. We employ a multidimensional amplitude estimation process to extract the full probability distribution from the quantum computer.
معرّف الورقة : 2505.21221العنوان : Quantum algorithms for solving a drift-diffusion equation: A complexity analysisالمؤلفون : Ellen Devereux (جامعة وارويك وفوجيتسو المملكة المتحدة)، Animesh Datta (جامعة وارويك)التصنيف : quant-phتاريخ النشر : 16 أكتوبر 2025 (نسخة arXiv المسبقة)رابط الورقة : https://arxiv.org/abs/2505.21221 تقترح هذه الورقة أربع خوارزميات كمية لحل معادلة الانجراف-الانتشار متعددة الأبعاد، تستند إلى محللات الأنظمة الخطية الكمية، محاكاة هاميلتونيان الكمية، المسارات العشوائية الكمية، وتحويل فورييه الكمي. من خلال تحليل التعقيد، تتم مقارنة هذه الطرق مع نظيراتها الكلاسيكية، حيث يُكتشف أن طريقة القطرية القائمة على تحويل فورييه الكمي توفر ميزة حسابية كمية في حل معادلات تفاضلية جزئية خطية في وقت نهائي ثابت. تستخدم الورقة عملية تقدير السعة متعددة الأبعاد لاستخراج التوزيع الاحتمالي الكامل من الحاسوب الكمي.
معادلة الانجراف-الانتشار (DDE) هي فئة مهمة من المعادلات التفاضلية الجزئية، بالصيغة المحددة:
∂ p ( x , t ) ∂ t = ∑ i = 1 d [ a ∂ ∂ x i [ p ( x , t ) ] + D ∂ 2 ∂ x i 2 [ p ( x , t ) ] ] \frac{\partial p(x,t)}{\partial t} = \sum_{i=1}^d \left[a\frac{\partial}{\partial x_i}[p(x,t)] + D\frac{\partial^2}{\partial x_i^2}[p(x,t)]\right] ∂ t ∂ p ( x , t ) = ∑ i = 1 d [ a ∂ x i ∂ [ p ( x , t )] + D ∂ x i 2 ∂ 2 [ p ( x , t )] ]
حيث x = { x 1 , . . . , x d } ∈ R d x = \{x_1, ..., x_d\} \in \mathbb{R}^d x = { x 1 , ... , x d } ∈ R d متجه بـ d d d بُعد، و a a a و D D D هما معاملات الانجراف والانتشار الموجبة على التوالي.
الأهمية النظرية : معادلة الانجراف-الانتشار كمعادلة فوكر-بلانك تصف سرعة الجسيمات، وترتبط ارتباطاً وثيقاً بمعادلات بلاك-شولز وناڤييه-ستوكسالتطبيقات العملية : تُستخدم في نمذجة المخاطر المالية والتنبؤ بقوة الرياح في عدة صناعات لدعم القراراتالتحديات الحسابية : تتطلب الطرق العددية التقليدية تقسيماً لمجالات مشاكل كبيرة ومعقدة، مما يستهلك ذاكرة وموارد حسابية ضخمةتشمل الطرق الكلاسيكية لحل المعادلات التفاضلية الجزئية:
محللات المعادلات الخطية مثل طريقة التدرج المرافق طرق المسارات العشوائية طرق القطرية القائمة على تحويل فورييه السريع تواجه هذه الطرق تحديات في معالجة المشاكل عالية الأبعاد حيث ينمو التعقيد الحسابي بشكل أسي مع البُعد.
اقتراح أربع خوارزميات كمية : تستند إلى محللات الأنظمة الخطية الكمية، التطور الزمني الكمي، المسارات العشوائية الكمية، وتحويل فورييه الكميتحليل نظري للتعقيد : توفير تحليل تفصيلي لتعقيد الوقت، مما يثبت شروط وجود الميزة الكميةطريقة تقدير السعة متعددة الأبعاد : تطبيق أول لتقدير السعة متعدد الأبعاد في حل المعادلات التفاضلية الجزئية، مما يحقق استخراج التوزيع الاحتمالي الكاملالتحقق من الجدوى العملية : التحقق من قيمة التطبيق التجاري للطريقة من خلال أمثلة النمذجة الماليةالبحث عن حل تقريبي p ~ ~ ( x , t ) \tilde{\tilde{p}}(x,t) p ~ ~ ( x , t ) لمعادلة الانجراف-الانتشار بحيث يرضي في اللحظة t = T t=T t = T :
∣ ∣ p ~ ~ ( x , t ) − p ( x , t ) ∣ ∣ ∞ ≤ ϵ ||\tilde{\tilde{p}}(x,t) - p(x,t)||_\infty \leq \epsilon ∣∣ p ~ ~ ( x , t ) − p ( x , t ) ∣ ∣ ∞ ≤ ϵ
حيث ϵ ∈ ( 0 , 1 ) \epsilon \in (0,1) ϵ ∈ ( 0 , 1 ) خطأ معطى، و x ∈ [ − L , L ] d x \in [-L,L]^d x ∈ [ − L , L ] d .
استخدام صيغة الفرق الأمامي الزمني والمركزي المكاني:
خطوة المسافة: Δ x = 2 L / n x \Delta x = 2L/n_x Δ x = 2 L / n x خطوة الزمن: Δ t = T / n t \Delta t = T/n_t Δ t = T / n t شرط الاستقرار: Δ t ≤ Δ x 2 / ( 2 d D ) \Delta t \leq \Delta x^2/(2dD) Δ t ≤ Δ x 2 / ( 2 d D ) الفكرة الأساسية : تحويل معادلة الانجراف-الانتشار المقسمة إلى نظام معادلات خطية A p ~ = b A\tilde{p} = b A p ~ = b تعقيد الوقت : O ~ ( d 5 T 2 ζ ( a L + D ) 2 ϵ q ϵ c ) \tilde{O}\left(\frac{d^5T^2\zeta(aL+D)^2}{\epsilon_q\epsilon_c}\right) O ~ ( ϵ q ϵ c d 5 T 2 ζ ( a L + D ) 2 ) شروط التطبيق : يتطلب عدد شرط محدود (κ L = 5 \kappa_L = 5 κ L = 5 عندما D Δ t / Δ x 2 ≤ 1 / 5 D\Delta t/\Delta x^2 \leq 1/5 D Δ t /Δ x 2 ≤ 1/5 و a / D < 2 10 a/D < 2\sqrt{10} a / D < 2 10 )الفكرة الأساسية : استخدام محاكاة هاميلتونيان والتركيب الخطي للمؤثرات الأحادية لتحقيق التطور الزمنيتعقيد الوقت : O ~ ( d d / 2 + 3 T d / 2 + 2 ζ d / 4 + 1 L 2 ( a L + D ) d / 2 ϵ q ϵ c d / 4 + 2 ) \tilde{O}\left(\frac{d^{d/2+3}T^{d/2+2}\zeta^{d/4+1}L^2(aL+D)^{d/2}}{\epsilon_q\epsilon_c^{d/4+2}}\right) O ~ ( ϵ q ϵ c d /4 + 2 d d /2 + 3 T d /2 + 2 ζ d /4 + 1 L 2 ( a L + D ) d /2 ) الخصائص : محاكاة مباشرة لعملية التطور الزمني الكميالفكرة الأساسية : الاستفادة من الطبيعة العشوائية لمعادلة الانجراف-الانتشار، محاكاة من خلال المسار العشوائي الكميتعقيد الوقت : O ~ ( d ( d + 7 ) / 2 T d / 2 + 1 ζ d / 4 + 1 / 2 ( a L + D ) d / 2 + 1 ϵ q ϵ c d / 4 + 1 / 2 ) \tilde{O}\left(\frac{d^{(d+7)/2}T^{d/2+1}\zeta^{d/4+1/2}(aL+D)^{d/2+1}}{\epsilon_q\epsilon_c^{d/4+1/2}}\right) O ~ ( ϵ q ϵ c d /4 + 1/2 d ( d + 7 ) /2 T d /2 + 1 ζ d /4 + 1/2 ( a L + D ) d /2 + 1 ) الميزة : قابل للتوسع إلى معادلات تفاضلية جزئية عشوائية غير خطيةالفكرة الأساسية : استخدام تحويل فورييه الكمي لقطرية المصفوفات الدورية، حساب L n t L^{n_t} L n t مباشرةتعقيد الوقت : O ~ ( d ( d / 2 + 2 ) T d / 2 ζ d / 4 ( a L + D ) d / 2 ϵ q ϵ c d / 4 ) \tilde{O}\left(\frac{d^{(d/2+2)}T^{d/2}\zeta^{d/4}(aL+D)^{d/2}}{\epsilon_q\epsilon_c^{d/4}}\right) O ~ ( ϵ q ϵ c d /4 d ( d /2 + 2 ) T d /2 ζ d /4 ( a L + D ) d /2 ) الميزة الرئيسية : تطبيق متزامن لجميع القيم الذاتية في حالة التراكب الكميتوسيع مبتكر لتقدير السعة أحادي البُعد إلى الحالة متعددة الأبعاد:
تحديد الإحداثيات ذات الاحتمالية ≥ ϵ q \epsilon_q ϵ q من خلال العينات بناء الدالة f ( x ) = ⟨ x , p ⟩ f(x) = \langle x,p \rangle f ( x ) = ⟨ x , p ⟩ استخدام تقدير الطور لاستخراج التوزيع الاحتمالي التعقيد: O ( log ( n x d / δ ) log ( 1 / ϵ q ) ϵ q ) O\left(\frac{\log(n_x^d/\delta)\log(1/\epsilon_q)}{\epsilon_q}\right) O ( ϵ q l o g ( n x d / δ ) l o g ( 1/ ϵ q ) ) بناءً على عملية أورنشتاين-أولينبيك في النمذجة المالية:
T = 5000 T = 5000 T = 5000 يوم، L = 10 L = 10 L = 10 دولاراتa = 0.2366 a = 0.2366 a = 0.2366 ، D = 0.2455 D = 0.2455 D = 0.2455 d = 3 d = 3 d = 3 أبعاد، ζ = 1 \zeta = 1 ζ = 1 دولار− 4 ^{-4} − 4 مقارنة تعقيد الوقت تحليل تعقيد المسافة شروط الميزة الكمية شروط الميزة الكمية: عندما ϵ q ≥ O ~ ( ϵ c d / 4 d D d d / 2 ζ d / 4 L d ( a L + D ) d / 2 ) \epsilon_q \geq \tilde{O}\left(\frac{\epsilon_c^{d/4}d^D d^{d/2}}{\zeta^{d/4}L^d(aL+D)^{d/2}}\right) ϵ q ≥ O ~ ( ζ d /4 L d ( a L + D ) d /2 ϵ c d /4 d D d d /2 ) ، تتفوق طريقة القطرية الكمية على الطريقة الكلاسيكية.
الطريقة التعقيد الكلاسيكي التعقيد الكمي الميزة الكمية المعادلة الخطية 4.85 × 10 22 / ϵ c 3 4.85 \times 10^{22}/\epsilon_c^3 4.85 × 1 0 22 / ϵ c 3 4.14 × 10 10 / ( ϵ c ϵ q ) 4.14 \times 10^{10}/(\epsilon_c\epsilon_q) 4.14 × 1 0 10 / ( ϵ c ϵ q ) ✓ خطوة زمنية 1.24 × 10 18 / ϵ c 2.5 1.24 \times 10^{18}/\epsilon_c^{2.5} 1.24 × 1 0 18 / ϵ c 2.5 5.23 × 10 17 / ( ϵ c 2.75 ϵ q ) 5.23 \times 10^{17}/(\epsilon_c^{2.75}\epsilon_q) 5.23 × 1 0 17 / ( ϵ c 2.75 ϵ q ) × مسار عشوائي 1.24 × 10 18 / ϵ c 2.5 1.24 \times 10^{18}/\epsilon_c^{2.5} 1.24 × 1 0 18 / ϵ c 2.5 4.73 × 10 12 / ( ϵ c 1.25 ϵ q ) 4.73 \times 10^{12}/(\epsilon_c^{1.25}\epsilon_q) 4.73 × 1 0 12 / ( ϵ c 1.25 ϵ q ) ✓ قطرية 8.07 × 10 11 / ϵ c 1.5 8.07 \times 10^{11}/\epsilon_c^{1.5} 8.07 × 1 0 11 / ϵ c 1.5 6.98 × 10 7 / ( ϵ c 0.75 ϵ q ) 6.98 \times 10^7/(\epsilon_c^{0.75}\epsilon_q) 6.98 × 1 0 7 / ( ϵ c 0.75 ϵ q ) ✓
يبلغ تعقيد المسافة لجميع الطرق الكمية O ~ ( d / ϵ q ) \tilde{O}(d/\epsilon_q) O ~ ( d / ϵ q ) ، يتحدد بشكل أساسي من خلال ترميز الحالة الكمية وبروتوكولات القياس.
طريقة تعيين هاميلتونيان : تعيين المعادلات التفاضلية الجزئية إلى هاميلتونيان، الحل من خلال معادلة شرودنجرطريقة النظام الخطي : بناء نظام معادلات خطية بعد التقسيم للحل الكميالخوارزميات الكمية المتغيرة : طرق متغيرة مناسبة لأجهزة NISQأول مقارنة منهجية لأربع طرق كمية لحل معادلة الانجراف-الانتشار متعددة الأبعاد إدخال تقدير السعة متعدد الأبعاد لاستخراج التوزيع الاحتمالي الكامل توفير تحليل نظري صارم للتعقيد قطرية تحويل فورييه الكمي هي الطريقة الكمية الأكثر فعالية لحل معادلة الانجراف-الانتشار في وقت ثابتوجود ميزة كمية لكن يتطلب استيفاء شروط معاملات محددةتقدير السعة متعدد الأبعاد يوسع بنجاح نطاق تطبيق حل المعادلات التفاضلية الجزئية الكميةالاعتماد على المعاملات : الميزة الكمية تعتمد بشدة على معاملات المشكلة ومتطلبات الخطأنطاق التطبيق : بعض الطرق تنطبق فقط على نطاقات معاملات محددة (مثل طريقة النظام الخطي الكمي)تعقيد التنفيذ : يتطلب ذاكرة الوصول العشوائي الكمية (QRAM) وأجهزة كمية متقدمة أخرىالتوسع إلى معادلات تفاضلية جزئية ذات معاملات متغيرة مكانياً دراسة طرق حل المعادلات التفاضلية الجزئية غير الخطية الكمية تحسين التنفيذ العملي للخوارزميات الكمية الصرامة النظرية : توفير تحليل تعقيد كامل وإثبات رياضيشمولية الطريقة : مقارنة منهجية لأربع طرق كمية مختلفةالقيمة العملية : التحقق من القيمة التجارية للطريقة من خلال التطبيقات الماليةالابتكار التقني : أول تطبيق لتقدير السعة متعدد الأبعاد في حل المعادلات التفاضلية الجزئيةشروط الميزة الكمية صارمة : يتطلب ϵ q \epsilon_q ϵ q استيفاء شروط محددة لتحقيق ميزة كميةمتطلبات الأجهزة عالية : يتطلب حاسوب كمي متسامح مع الأخطاء و QRAMقيود التطبيق : ينطبق بشكل أساسي على المعادلات التفاضلية الجزئية الخطية، التوسع للحالات غير الخطية محدودالمساهمة الأكاديمية : توفير أساس نظري مهم لحل المعادلات التفاضلية الجزئية الكميةالآفاق العملية : تطبيقات محتملة في النمذجة المالية والحسابات العلمية وغيرهاالتقدم التكنولوجي : دفع تطور الخوارزميات الكمية في حل المعادلات التفاضلية الجزئيةحل المعادلات التفاضلية الجزئية الخطية عالية الأبعاد نمذجة المخاطر المالية وتسعير الخيارات محاكاة عمليات الانتشار في الحسابات العلمية التطبيقات التي تتطلب معلومات التوزيع الاحتمالي الكامل تستشهد الورقة بـ 43 مرجعاً ذا صلة، تغطي بشكل أساسي:
الأسس النظرية للخوارزميات الكمية الطرق العددية للمعادلات التفاضلية الجزئية محللات الأنظمة الخطية الكمية المسارات العشوائية الكمية وتحويل فورييه الكمي العمليات العشوائية في النمذجة المالية التقييم الإجمالي : هذه ورقة عالية الجودة في نظرية الخوارزميات الكمية، تقدم مساهمات مهمة في مجال حل المعادلات التفاضلية الجزئية الكمية. على الرغم من أن التطبيقات العملية تواجه حالياً قيوداً على الأجهزة، فإنها تضع أساساً نظرياً متيناً لتطبيقات الحوسبة الكمية المستقبلية في مجال الحسابات العلمية.