تصوير الحالات الكمومية لتجمعات اللف المغزلي: كيف نحوّل الطيف إلى حالة كمومية

الأصل: "Computational framework for quantum state tomography of spin ensembles"

arXiv:2609.11868v1 · 2026-09-10 · CC BY 4.0 · 3 دقيقة

فيزياء الكم مادة مكثفة أخرى الفيزياء الحيوية الفيزياء الحاسوبية الأجهزة والكاشفات

حزمة Tomography-NMR المفتوحة تحوّل أطياف الرنين المغناطيسي النووي الخام إلى مصفوفات كثافة بدقة تتجاوز 99%.

السياق

يقف هذا العمل على هامش المشهد: فأقرب الجيران في الأرشيف يشتغلون على الخوارزميات الكمومية وتصحيح الأخطاء والمعلومات الكمومية، دون أن يقتربوا من إعادة بناء الحالات انطلاقاً من الأطياف. وتوجد حوله بالفعل أدوات لمحاكاة ديناميكا اللف المغزلي (PULSEE) ومكتبات عامة (QuTiP وnmrglue)، لكن الحلقة الرابطة بين بيانات المطياف الخام ومصفوفة الكثافة المُتحقَّق منها كانت غائبة. وحزمة Tomography-NMR تملأ هذه الثغرة، واضعةً نصب عينيها معالجة شفافة وقابلة لإعادة الإنتاج لبيانات الرنين المغناطيسي النووي لأغراض التصوير الكمومي.

المنهجية

تنفّذ الحزمة خط تحليل بلغة Python 3.11+، حيث تخضع الإشارات الخام في المجال الزمني لتحويل فورييه مع تصحيح طور تلقائي أو يدوي، ثم تُستخرج شدة قمم الثنائيات بثلاث طرق. الأولى: القياس المباشر لارتفاع القمة عند نقطة واحدة؛ والثانية: التكامل العددي بطريقة سيمبسون داخل نافذة ثابتة؛ والثالثة: مسح منهجي للمعاملات (إزاحات التردد، وعرض النافذة، وتصحيح الانفصال) لتعظيم الدقة قياساً إلى حالة نظرية معلومة. ويُسند مجموع القمتين اليسرى واليمنى للثنائي وفرقهما إلى معاملات مؤثرات الضرب عبر جدول بحث، ثم تُستعاد مصفوفة الكثافة مع فرض قيود فيزيائية (الهرميتية، وتطبيع الأثر، وشبه التعريف الموجب).

النتائج

على 20 حالة تجريبية مُهيّأة على نواتين مرتبطتين من الفوسفور-31 في ADP، بلغ متوسط دقة إعادة البناء وفق مقياس فورتوناتو ما بين 0.975 و0.995 بحسب فئة الحالة. فقد وصلت حالات الأساس إلى 0.995، ومخرجات بوابة هادامارد إلى 0.992، وCNOT إلى 0.983، وحالات بيل إلى 0.980، أما التركيب H+T فبلغ 0.975. ويعطي التكامل المُحسَّن دقة تتجاوز 99% لحالات بعينها (تصل إلى 0.998 للحالة |00⟩)، بينما تظل الطرق التي لا تعرف الحالة الهدف عند حدود 98%. وتُقدَّر أخطاء عناصر مصفوفة الكثافة عند مستوى 10⁻³–10⁻⁴.

الأهمية

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

التطور المستقبلي

سيتطلب التوسّع إلى أنظمة بثلاثة كيوبتات أو أكثر حلّ مشكلة تراكب المتعددات وقابلية تحسين المسح للتوسّع. وقد أُدرج بالفعل نموذج أولي لوحدة خاصة بثلاث لفّات مغزلية، لكن الدعم الكامل يبقى مسألة مفتوحة. كما يُخطَّط لاستبدال البحث الشبكي بتحسين بايزي، وللتحليل في الزمن الحقيقي من أجل البروتوكولات التكيفية، وللمقاييس تحقق جديدة مثل شهود التشابك.

التأثير

سيؤثر العمل في المعلوماتية الكمومية التجريبية، ولا سيما الحوسبة الكمومية بالرنين المغناطيسي النووي، إذ يوفّر أداة موحّدة لتوصيف الحالات. كما سيسرّع تدريب المتخصصين الجدد ويسهّل مقارنة النتائج بين المجموعات البحثية.

الخطوات التالية

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

المسائل المفتوحة الرئيسية

لا تزال المعلوماتية الكمومية تواجه مسألة غير محلولة هي تصوير الحالات القابل للتوسّع: النمو الأُسّي لعدد القياسات مع عدد الكيوبتات. ولا يحلّ هذا العمل تلك المسألة، لكنه يضع معياراً للدقة والشفافية في حالة الكيوبتين، وهو أمر ضروري للتحقق من الطرق المستقبلية.

سُمّيت الحزمة Tomography-NMR وتُوزَّع تحت رخصة MIT، ما يسمح باستخدامها بحرية في المشاريع الأكاديمية والتجارية. وقد أُجري التحقق على جزيء ADP، حيث ترتبط نواتا الفوسفور-31 عبر تفاعل J بنحو 22 هرتز.

\hat{\rho} = \sum_{i,j \in \{o,x,y,z\}} c_{ij} \hat{B}_{ij}

تُمثَّل مصفوفة الكثافة كتركيب خطي من 16 مؤثراً أساسياً بمعاملات c_ij.

\mathrm{Re}(L+R) = \mathrm{Tr}(\hat{\rho}\tilde{\hat{I}}^1_x), \quad \mathrm{Re}(L-R) = \mathrm{Tr}(\hat{\rho}2\tilde{\hat{I}}^1_x \tilde{\hat{I}}^2_z)

تشفّر شدة قمم الطيف مباشرةً القيم المتوقعة للمؤثرات، ما يتيح استخراج معاملات التفكيك.

F_{\mathrm{proj}}(\hat{\rho}_m, \hat{\rho}_{th}) = \frac{\mathrm{Tr}(\hat{\rho}_m \hat{\rho}_{th})}{\sqrt{\mathrm{Tr}(\hat{\rho}_m^2) \mathrm{Tr}(\hat{\rho}_{th}^2)}}

جداء داخلي مُطبَّع بين مصفوفتي الكثافة المقيسة والنظرية، يُستخدم لتقييم جودة إعادة البناء.

أرقام رئيسية

دقة إعادة البناء (المُحسَّنة): حتى 0.998
دقة إعادة البناء (دون معرفة الهدف): نحو 0.980
عدد الحالات المُتحقَّق منها: 20
عدد المعاملات غير الواحدية لكيوبتين: 15
زمن التحسين لحالة واحدة: 2–5 دقائق

أنا المؤلف