تحليل بيانات حيود الأشعة السينية
Automated X-ray diffraction analysis by Python
يُعد جهاز حيود الأشعة السينية (XRD) بمثابة “قارئ البصمة” للمواد في عالم الأبحاث العلمية. إن أهميته لا تنبع فقط من قدرته على استكشاف المادة، بل في كونه الأداة المعيارية الأولى لفهم “كيفية” ترتيب الذرات داخل أي مادة صلبة، وهو الترتيب الذي يحدد في النهاية جميع خصائصها الفيزيائية والكيميائية والميكانيكية.
حلول برمجية متقدمة بلغة بايثون
1-المعايرة واستخراج معاملات كاليوتي (Caglioti Calibration)
🔬 لماذا نستخدم عينة السيليكون القياسية لمعايرة جهاز الـ XRD؟
عندما نقيس عرض القمة (FWHM) في طيف حيود الأشعة السينية، فإن هذا الاتساع لا يمثل فقط خصائص مادتك (مثل صغر الحجم البلوري أو الإجهاد)، بل إن الجهاز نفسه يضيف اتساعاً إضافياً وهمياً بسبب هندسة الشقوق (Slits) ودقة الكاشف. يُعرف هذا بـ الاتساع الآلي (Instrumental Broadening).
1️⃣ معادلة كاليوتي (Caglioti Equation)
تصف هذه المعادلة كيف يتغير الاتساع الآلي (βinst) باختلاف زاوية الحيود (θ):
المعاملات (U, V, W) هي بصمة جهازك الخاصة. لاستخراجها، نستخدم عينة قياسية مثل السيليكون (NIST 640)؛ لأن حجم بلوراتها ضخم جداً وخالية من الإجهاد، مما يضمن أن أي اتساع نراه في قممها هو ناتج بنسبة 100% عن الجهاز وليس المادة.
2️⃣ تأثير الدومينو على النتائج (لماذا المعايرة حتمية؟)
استخدام قيم (U, V, W) عشوائية أو غير دقيقة سيؤدي إلى انهيار دقة حساباتك بالكامل:
-
خطأ في الاتساع الصافي: لاستخراج اتساع مادتك الفعلي (βphys)، يجب أن نطرح تأثير الجهاز:
βphys ≈ βmeasured – βinst -
حجم بلوري وهمي (معادلة شيرير): الحجم البلوري (D) يتناسب عكسياً مع الاتساع الصافي:
D = (K · λ) / (βphys · cosθ)إذا طرحت قيماً خاطئة لمعاملات كاليوتي وجعلت βphys صغيرة جداً، ستُظهر لك المعادلة حجماً بلورياً ضخماً غير حقيقي! والعكس صحيح.
* سيتم فتح الأداة في نافذة جديدة
2- تحليل بيانات حيود الأشعة السينية (XRD)
الهدف من الأداة: معالجة أطياف حيود الأشعة السينية بطريقة احترافية ترتقي لمستوى النشر العلمي، فصل تداخل أشعة (K-alpha 1 & 2)، حساب نسبة التبلور، استخراج حجم البلورات والإجهاد (Microstrain) مع حساب نسبة الخطأ، وتوليد تقرير PDF شامل تلقائياً.
خطوات العمل:
- تشغيل بيئة العمل: افتح الرابط أعلاه واضغط على زر التشغيل بجوار خلية الكود، أو استخدم الاختصار
Shift + Enter. - رفع البيانات المعملية: عند ظهور زر الرفع، اختر ملف القياس الخاص بك (تدعم الأداة صيغ
.xrdmlالاحترافية بالإضافة إلى.txtأو.csv). - المعالجة التلقائية: سيقوم الكود تلقائياً بفصل الخلفية وحساب نسبة التبلور (Crystallinity).
- النمذجة الرياضية: سيقوم الكود بدراسة كل قمة على حدة باستخدام نموذج Double Pseudo-Voigt وتطبيق تصحيحات Caglioti الآلية.
- التقرير النهائي: سيتم عرض الرسوم البيانية التفاعلية، وتحميل تقرير
PDFجاهز يحتوي على كافة الحسابات بشكل آلي.
ماذا ستقدم لك الأداة؟
- تقرير علمي بصيغة PDF: ملف جاهز للطباعة أو الرفع للأبحاث العلمية يحتوي على المنهجية، نسبة التبلور، الطور المتوقع، وجدول مفصل للقمم.
- جدول بيانات دقيق:
يستخرج لكل قمة:
- زاوية 2θ والمسافة البينية d-spacing (مع نسبة الخطأ ±).
- العرض عند منتصف الارتفاع FWHM (الصافي فيزيائياً بعد خصم اتساع الجهاز).
- حجم البلورة Crystallite Size باستخدام معادلة شيرر مع حساب عدم اليقين (Uncertainty).
- لوحة تفاعلية (Plotly):
شاشة منقسمة تعرض طيف الـ XRD مع المنحنيات المفصلة (Fitted Curves) على جهة، ورسم Williamson-Hall على الجهة الأخرى لتقييم الإجهاد البلوري.
شرح تحليلي للرسوم البيانية (التفاعلية)
تستخدم الأداة مكتبة Plotly لإنشاء لوحة تفاعلية (Interactive Dashboard). يمكنك تحديد أي منطقة لتكبيرها (Zoom)، أو تمرير مؤشر الماوس فوق المنحنيات والنقاط لقراءة القيم بدقة دون الحاجة للرجوع للجدول.
1. طيف الأشعة السينية (الرسم الأيسر – XRD Pattern)
- البيانات الخام (الخط الرمادي الفاتح): هو الطيف الأصلي المتضمن للخلفية غير المتبلورة (Amorphous). يُترك كمرجعية بصرية.
- النمط المعالج (الخط الأزرق المتصل): هو الطيف “النقي” بعد طرح الخلفية، وقاعدته تبدأ من الصفر، مما يجعل حساب المساحات ونسبة التبلور دقيقاً للغاية.
- منحنيات المطابقة (الخطوط الحمراء): تظهر فوق القمم فقط. تمثل دوال (Double Pseudo-Voigt) التي تحاكي شكل القمة بدقة متناهية لحساب العرض (FWHM) والموقع متجاهلة أي تشويش (Noise).
2. رسم Williamson-Hall (الرسم الأيمن)
يُستخدم هذا الرسم لفصل تأثير حجم البلورة عن تأثير الإجهاد الشبكي (Microstrain) على اتساع القمم.
- النقاط (Purple): تمثل القمم المعتمدة. المحور السيني يمثل $4\sin(\theta)$ والمحور الصادي يمثل $\beta\cos(\theta)$.
- خط المطابقة (Orange Fit): يمثل تطبيق نموذج (UDM).
- تقاطع الخط مع المحور الصادي: يُعطي متوسط حجم البلورة (Crystallite Size).
- ميل الخط (Slope): يعبر عن قيمة الإجهاد (Microstrain) داخل الشبكة البلورية الناتجة عن العيوب.
💡 استنتاج التباين الخواص (Anisotropy):
إذا ظهرت لك رسالة (Rejected: R² < 0.75 Data is anisotropic) واختفى الخط البرتقالي، فهذا ليس خطأً برمجياً! بل هو استنتاج فيزيائي دقيق يعني أن نقاط القمم متناثرة ولا تشكل خطاً مستقيماً. هذا يحدث عندما لا تنمو البلورات بشكل متماثل في كل الاتجاهات (مثل الأنابيب أو الصفائح النانوية)، وهنا يرفض الكود بشفافية تطبيق النموذج الخطي الموحد لمنع استخراج نتائج مضللة.
3- مطابقة الأطوار البلورية
أداة تحليل ومطابقة حيود الأشعة السينية (XRD)
الهدف من الأداة: معالجة أطياف حيود الأشعة السينية (XRD) الخام، إزالة تأثيرات الخلفية والضوضاء، ومطابقة القمم المعملية مع الأنماط النظرية المستخرجة من ملفات (CIF) لتحديد الأطوار المتبلورة بدقة.
خطوات العمل:
- تشغيل بيئة العمل: افتح الرابط أعلاه واضغط على زر التنشيط (Run) بجوار خلية الكود، أو استخدم الاختصار
Shift + Enter. - رفع البيانات المعملية: عند ظهور رسالة الطلب، اضغط على "Choose Files" وارفع ملف القياس الخاص بك (بصيغة
XRDMLأوTXT / CSV). - رفع بيانات التركيب البلوري (CIF): بعد ذلك، سيُطلب منك رفع ملفات الـ CIF التي تتوقع وجودها في العينة لإجراء عملية المطابقة.
- الضبط المتقدم (اختياري): إذا لاحظت انزياحاً في القمم المعملية، يمكنك تصحيحه بتعديل قيمة المتغير
ZERO_SHIFTداخل الكود وإعادة التشغيل.
لوحة رسومية تفاعلية تقدم الآتي:
- تصور دقيق لمعالجة الإشارة: يوضح مراحل تصفية البيانات من خلال رسم منحنى البيانات الخام، وخلفية التشتت، والطيف النهائي المعالج الجاهز للتحليل.
- استخراج خصائص القمم (Peak Properties):
عند تمرير مؤشر الماوس فوق أي قمة مستخرجة، ستظهر لك تفاصيلها مباشرة:
- الزاوية الدقيقة 2θ والشدة.
- العرض عند منتصف الارتفاع FWHM (وهو الرقم الأساسي لاستخدامه في معادلة شيرر لحساب حجم الحبيبات).
- المطابقة البصرية الفعّالة للأطوار:
تظهر القمم النظرية كخطوط عمودية (Stick Patterns) ملونة لتسهيل المقارنة وتأكيد الطور البلوري مع القمم المعملية المستخرجة.
شرح تحليلي للرسم البياني
يقدم المخرج الرسومي للأداة تمثيلاً بصرياً ممتازاً لعملية معالجة الإشارة. إليك دلالات المنحنيات والألوان:
1. معالجة الإشارة (Signal Processing)
- البيانات الخام (الخط الرمادي الفاتح المتعرج): هو الطيف الأصلي كما التقطه الجهاز. الارتفاع والتعرجات ناتجة عن تشتت الأشعة العشوائي (Amorphous) والضوضاء الآلية، ويُترك في الرسم كمرجعية بصرية للتأكد من سلامة عملية المعالجة.
- الخلفية المحسوبة (الخط البرتقالي المتقطع): تمثل مسار خوارزمية عزل الخلفية (SNIP) التي تتتبع قاع المنحنى الرمادي بذكاء لتعزل الضوضاء دون تشويه القمم.
- النمط المعالج (الخط الأزرق المتصل): هو الطيف "النقي" الناتج عن طرح الخط البرتقالي من الرمادي. تصبح قاعدته قريبة من الصفر مما يبرز القمم الحقيقية بوضوح.
2. استخراج القمم والمطابقة (Peak Matching)
- القمم المعملية (مثلثات زرقاء داكنة): هي مواقع القمم الحقيقية لعشيرتك التي تمكن الكود من التعرف عليها بعد تصفية الطيف.
- النمط النظري (خطوط ملونة متقطعة): تمثل المواقع الزاوية التي يجب أن تظهر فيها قمم المادة النظرية بناءً على ملف التركيب البلوري (CIF) الذي تم رفعه.
💡 استنتاج المطابقة (Mismatch):
إذا لاحظت أن القمم المعملية (الزرقاء) لا تتطابق مع مواقع القمم النظرية المتوقعة، فهذا يعني فيزيائياً أن العينة لا تحتوي على المادة المرفقة، أو أن ملف الـ CIF لا يخص الطور البلوري الصحيح. في هذه الحالة، يجب البحث عن ملفات CIF لمواد أخرى متوقعة في العينة وتجربتها مجدداً.
