لا يمكن توصيف المادة توصيفاً علمياً متكاملاً بالاعتماد على تقنية واحدة؛ فكل تقنية تقيس استجابة فيزيائية مختلفة ضمن حجم تفاعل وعمق تحليل وحساسية ونموذج رياضي محدد له حدود صلاحية.
القيمة العلمية لا تكمن في كمية البيانات، بل في القدرة على تحويل قياسات مستقلة إلى استنتاج فيزيائي قابل للدفاع عنه، مرفوقاً بهامش عدم يقين كمي:
$$\text{Measurement} \neq \text{Interpretation}$$
$$\boxed{\text{Measurement} + \text{Valid Model} + \text{Calibration} + \text{Uncertainty} + \text{Cross-validation} \rightarrow \text{Defensible Physical Interpretation}}$$
1. حيود الأشعة السينية (XRD)
1.1 الأساس والأطوار البلورية
$$n\lambda = 2d\sin\theta$$
مطابقة مواقع القمم وحدها لا تكفي لإثبات طور معين عند تداخل الأطوار أو ضعف القمم أو وجود مكوّن غير متبلور. حد الكشف النموذجي لطور بلوري في جهاز مخبري قياسي يتراوح بين 1% و5% وزناً (انظر القسم 12.8): غياب قمة ليس دليلاً على غياب الطور، بل يعني فقط أن تركيزه دون هذا الحد.
1.2 حجم المجال المتماسك (معادلة Scherrer)
$$D = \frac{K\lambda}{\beta\cos\theta}, \qquad K \approx 0.89\text{–}0.94$$
يمثل \(D\) حجم المجال المتماسك حيودياً، لا حجم الحبيبة المجهرية؛ فالحبيبة الواحدة قد تحوي عدة مجالات تفصل بينها حدود أو عيوب بلورية.
تعتمد قيمة \(K\) على شكل البلورية المفترض (كروية تقريباً 0.9 عند استخدام عرض القمة عند نصف الارتفاع FWHM، مقابل قيمة تقارب 1.0 عند استخدام عرض القمة المتكامل)، وعلى افتراض توزيع أحادي الحجم. توزيعات lognormal للحجوم تتطلب معالجة بالالتفاف الرياضي (convolution) لا مجرد تعديل قيمة \(K\). مثال محسوب كامل لعدم اليقين في هذه المعادلة يرد في القسم 12.6.
1.3 التصحيح الآلي لاتساع القمة
$$\beta_{\text{sample}} = \sqrt{\beta_{\text{obs}}^2 – \beta_{\text{inst}}^2} \quad \text{(تقريب Gaussian)}$$
$$\beta_{\text{sample}} = \beta_{\text{obs}} – \beta_{\text{inst}} \quad \text{(تقريب Lorentzian بحت)}$$
القمم الضعيفة ذات الأجنحة الواسعة في المواد النانوية تقترب من الشكل Lorentzian البحت، فيصبح الطرح المباشر أدق من الجذر التربيعي في هذه الحالة. أما الحالة العامة (ملف Voigt أو pseudo-Voigt) فتتطلب تحليل ملف مركّب ضمن Rietveld refinement، دون معادلة مبسطة واحدة تُعمَّم على كل الحالات.
1.4 تحليل Williamson-Hall
$$\beta\cos\theta = \frac{K\lambda}{D} + 4\epsilon\sin\theta \qquad \text{(نموذج UDM)}$$
هذه المعادلة تمثل نموذج الانفعال المنتظم (UDM). عند تباين المعامل المرن مع الاتجاه البلوري، تُستبدل \(\epsilon\) بـ \(\epsilon/Y_h\) ضمن نماذج USDM أو UDEDM، أو يُستخدم مخطط Size-Strain Plot بدلاً منها. الخطية الظاهرة في مخطط Williamson-Hall لا تثبت تلقائياً صحة النموذج؛ فعند وجود إجهاد غير متجانس أو عدم تناظر اتجاهي أو تداخل قمم، يصبح تحليل Warren-Averbach (تحليل فورييه لاتساع القمم) أو Rietveld refinement أكثر ملاءمة. الطريقتان تكامليتان لا بديلتان: يحلل Rietveld النمط الكامل والمعلمات البنيوية دفعة واحدة، بينما يفصل Warren-Averbach تحديداً مساهمتي الحجم والانفعال عبر معاملات فورييه.
2. المجهر الإلكتروني الماسح (SEM)
تتميز الإلكترونات الثانوية (SE) بحساسية عالية للتضاريس (شكل الجسيمات، الخشونة، التشققات، المسامية). أما الإلكترونات المرتدة (BSE) فتنتج تبايناً يعتمد على العدد الذري المتوسط، لكنه تباين تركيبي نوعي لا تحديد كمي مباشر؛ يحتاج دعماً من EDS للتحقق منه.
حجم تفاعل الإلكترونات (معادلة Kanaya-Okayama)
$$\boxed{R_{KO}\,[\mu\text{m}] = \frac{0.0276 \times A\,[\text{g/mol}] \times E_0^{1.67}\,[\text{keV}]}{Z^{0.89} \times \rho\,[\text{g/cm}^3]}}$$
اتفاقية الوحدات إلزامية: \(A\) بوحدة g/mol، \(E_0\) بوحدة keV، \(\rho\) بوحدة g/cm³، والناتج بوحدة μm. صيغت العلاقة أصلاً لمجال طاقة يقارب 1 إلى 30 keV، وتُستخدم كتقدير أولي فقط خارج هذا المجال.
| المادة | \(A\) (g/mol) | \(Z\) | \(\rho\) (g/cm³) | \(R_{KO}\) (μm) |
|---|---|---|---|---|
| C (جرافيت) | 12.01 | 6 | 2.26 | ≈ 4.4 |
| Al | 26.98 | 13 | 2.70 | ≈ 4.2 |
| Si | 28.09 | 14 | 2.33 | ≈ 4.7 |
| Au | 196.97 | 79 | 19.32 | ≈ 0.86 |
هذا تقدير لمدى الإلكترونات الكلي، وليس الدقة المكانية الفعلية لتحليل EDS. عمق توليد الأشعة السينية المميزة أصغر من \(R_{KO}\) فعلياً، ويعتمد على مقطع التأين وعلى علاقات من نوع Anderson-Hasler، كما يعتمد الحجم الفعلي على زاوية الميل وهندسة الكاشف وامتصاص الأشعة داخل المادة. في العينات المعقدة أو متعددة الطبقات، تُستخدم محاكاة Monte Carlo (مثل برنامج CASINO) لتقدير واقعي لتوزيع التفاعل. النتيجة العملية: دقة EDS المكانية في العينات السميكة من رتبة الميكرومترات، لا حجم بقعة الإلكترونات النانوية.
3. المجهر الإلكتروني النافذ (TEM)
3.1 معادلة الكاميرا (SAED)
$$R \cdot d = \lambda L \quad \Longrightarrow \quad d = \frac{\lambda L}{R}$$
حيث \(R\) نصف قطر الحلقة الحيودية، ويجب عدم الخلط بينه وبين قطرها عند الحساب.
3.2 الطول الموجي النسبوي للإلكترون
$$\lambda = \frac{h}{\sqrt{2m_0eV\left(1+\dfrac{eV}{2m_0c^2}\right)}}$$
عند الجهود العالية يجب أخذ التأثيرات النسبوية في الاعتبار. عند \(V=200\) kV تنتج المعادلة \(\lambda \approx 2.51\) pm، وهي القيمة المعيارية المستخدمة في معايرة طول الكاميرا.
3.3 تحضير العينة بتقنية FIB
يُحدث استخدام أيونات Ga⁺ عالية الطاقة زرع أيونات وعيوباً إشعاعية وإعادة ترسيب وطبقة سطحية متضررة. لا توجد قيمة ثابتة عامة لعمق هذا الضرر؛ فهو يعتمد على نوع المادة وطاقة الأيونات والجرعة وزاوية السقوط. أي تعميم رقمي (مثل افتراض ضرر ثابت بين 5 و20 نانومتراً دائماً) غير دقيق علمياً. يمكن تقليل الضرر بتنظيف نهائي عند جهد منخفض بين 2 و5 kV، وهذا مهم خصوصاً عند إجراء EDS أو EELS قرب السطح المتضرر.
4. مطيافية الإلكترونات الضوئية بالأشعة السينية (XPS)
4.1 الأساس الفيزيائي
$$E_k = h\nu – E_b – \phi$$
4.2 عمق التحليل: من IMFP إلى EAL
$$\boxed{d_{95\%} \approx 3\,\lambda_{EAL}\cos\theta}$$
العامل 3 يعني أن نحو 95% من الإلكترونات الضوئية المقاسة (بما أن \(1-e^{-3}\approx0.95\)) تنشأ من هذا العمق. يجب التمييز بين كميتين تُخلَطان خطأً في كثير من الأدبيات:
- المسار الحر الوسطي غير المرن IMFP: متوسط المسافة بين حدثين غير مرنين متتاليين، وهو خاصية جوهرية للمادة نفسها.
- الطول الفعّال للإضعاف EAL: يصف الشدة المقاسة فعلياً بعد إدخال تأثير التشتت المرن؛ وعملياً \(\lambda_{EAL}\) يقارب 0.7 إلى 0.9 من قيمة \(\lambda_{IMFP}\) لأغلب المواد، وفق تعريفات معيار ISO 18115-1 وقاعدة بيانات NIST EAL (بحسب أعمال Powell وJablonski).
استخدام IMFP بدل EAL يُبالغ في تقدير عمق التحليل بنسبة تتراوح بين 10% و30%، وهو خطأ منهجي يتجاهل سماكات الأغشية الرقيقة والبروفايلات العمقية. للنمذجة الكمية الدقيقة في الطبقات المتعددة، تُستخدم أداة SESSA الصادرة عن NIST.
صيغة \(\cos\theta\) تفترض سطحاً مستوياً ناعماً، وتتفكك عند زوايا \(\theta\) أكبر من 60 إلى 70 درجة، حيث تسيطر الخشونة والتشتت المرن وتأثيرات الانكسار. كما يجب التمييز بين عمق اختراق الأشعة السينية، وهو أكبر بكثير ويُقاس بالميكرومترات، وعمق التحليل الإلكتروني الذي لا يتجاوز بضعة نانومترات.
4.3 التحليل الزاوي مقابل الحفر الأيوني
يعطي تغيير زاوية القياس في تقنية Angle-Resolved XPS (ARXPS) بروفايلاً عمقياً غير متلف نسبياً. أما الحفر بالرش الأيوني (sputtering) فقد يُحدث اختزال أكاسيد وإزالة تفضيلية لعناصر معينة وخلط طبقات وعيوباً مصطنعة؛ فلا يُفسَّر البروفايل الناتج عنه كصورة مثالية للتوزيع الأصلي للعناصر.
4.4 حساب النسب الذرية
$$\frac{n_A}{n_B} = \frac{I_A/S_A}{I_B/S_B}$$
تعتمد الدقة على تصحيح الخلفية بطريقة Shirley أو Tougaard، وعلى اختيار مناطق القمم، وعلى شحن العينة. قمة C 1s القريبة من 284.8 eV مرجع شائع، لكن أعمال Greczynski وHultman أظهرت أنها قد تُدخل انحيازات منهجية كبيرة، لأن الكربون السطحي قد لا يمثل تركيب العينة نفسها. يُستحسن التحقق بمراجع بديلة، مثل حافة فيرمي للمعادن أو قمة Au 4f⁷/₂ عند 84.0 eV، أو مناقشة الحد المرجعي منهجياً في كل تقرير.
5. مطيافية الأشعة السينية المشتتة للطاقة (EDS/EDX)
5.1 قانون موزلي الكامل
$$\nu = \frac{3}{4}R_\infty c\,(Z-\sigma)^2 \quad \Longrightarrow \quad \sqrt{\nu}\propto(Z-\sigma)$$
حيث \(\sigma\) ثابت الحجب الفعّال، ويقارب 1 لسلسلة \(K_\alpha\)، وأكبر من ذلك لسلاسل \(L\). هذه الصياغة تفسيرية تاريخياً، وليست الأساس الكامل لخوارزميات التحليل الكمي الحديثة في EDS.
5.2 التصحيح الكمي
$$C_i = K_i\,Z_i\,A_i\,F_i \qquad \text{(نموذج ZAF، للعينات السميكة الماصة)}$$
في العينات الرقيقة أو متعددة الطبقات، تصبح نماذج \(\phi(\rho z)\) أكثر ضرورة، لأن افتراض العينة المتجانسة الكافية السمك يفشل عندها. أما في TEM-EDS للعينات الرقيقة فالإطار المعياري هو نسبة Cliff-Lorimer:
$$\frac{C_A}{C_B} = k_{AB}\frac{I_A}{I_B}$$
حيث تُعايَر معاملات \(k_{AB}\) تجريبياً أو حسابياً حسب الجهاز والمادة.
5.3 حدود كشف العناصر الخفيفة (\(Z<5\))
ضعف الكشف ناتج عن عوامل متعددة مجتمعة: انخفاض طاقة الأشعة المميزة، وامتصاصها داخل العينة والنافذة، وضعف كفاءة الكاشف عند الطاقات المنخفضة، والمنافسة مع عمليات Auger. لا يصح اختزال السبب في عامل واحد. عملياً، حد الكشف يتراوح بين 0.1% و1% وزناً للعناصر المتوسطة والثقيلة، ويرتفع إلى 1% إلى 5% للعناصر الخفيفة، ويصبح غير عملي للعناصر ذات العدد الذري أقل من 5 (انظر القسم 12.8).
6. مطيافية تشتت رامان
6.1 الإزاحة الأساسية
$$\Delta\nu = \nu_0 – \nu_s$$
6.2 الإجهاد والانفعال: الصيغة الموترية العامة
$$\Delta\omega = -\omega_0\,\gamma_{ij}\,\epsilon_{ij}$$
حيث \(\gamma_{ij}\) موتر Grüneisen للنمط الاهتزازي، و\(\epsilon_{ij}\) موتر الانفعال. هذه الصيغة تمثل الحد الخطي منخفض الانفعال، وتختزل إلى الحالة القياسية البسيطة فقط عند انفعال هيدروستاتيكي متماثل في بلورة مكعبية، حيث \(\gamma_{ij}=\gamma\,\delta_{ij}\):
$$\Delta\omega = -\omega_0\,\gamma\,\text{Tr}(\epsilon_{ij}) = -\omega_0\,\gamma\,(\epsilon_{xx}+\epsilon_{yy}+\epsilon_{zz})$$
خارج هذه الحالة، كما في الجرافين أو MoS₂ أو عند انفعال غير متماثل، تظهر مساهمات قصية منفصلة؛ فمثلاً للنمط \(E_{2g}\):
$$\Delta\omega \sim \gamma\,\omega_0(\epsilon_{xx}+\epsilon_{yy}) + \beta\,\omega_0(\epsilon_{xx}-\epsilon_{yy})\cos 2\theta$$
وتصبح تقنية Polarized Raman بزوايا استقطاب متعددة ضرورية لفصل هذه المساهمات. عند إجهادات كبيرة خارج الحد الخطي، تُوصف الإزاحات بمعاملات تشوه الفونون \(p\)، \(q\)، \(r\) ضمن إطار Anastassakis مع حل معادلة secular، وهو المستوى الأدق للبلورات المكعبية كالسيليكون.
6.3 علاقة Tuinstra-Koenig وحدود صلاحيتها
$$\frac{I_D}{I_G} = \frac{C(\lambda)}{L_a} \qquad (L_a \gtrsim 2\ \text{nm}, \text{ عيوب متفرقة})$$
عند اضطراب عالٍ، حين يصبح \(L_a\) صغيراً جداً، ينعكس الاتجاه وفق نموذج Ferrari-Robertson:
$$\frac{I_D}{I_G} \propto L_a^2 \qquad (L_a \lesssim 2\ \text{nm}, \text{ عيوب متراكبة})$$
تطبيق العلاقة الأولى خارج نطاقها ينتج استنتاجاً معكوساً تماماً بشأن حجم المجال البنيوي.
7. المجهر بالقوة الذرية (AFM)
7.1 قوى التفاعل والخشونة
$$U(r) = 4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}-\left(\frac{\sigma}{r}\right)^6\right] \qquad \text{(Lennard-Jones، نموذج مبسط جزئي)}$$
$$F = -k\Delta z, \qquad f_0 = \frac{1}{2\pi}\sqrt{\frac{k}{m^*}}$$
$$R_a = \frac{1}{N}\sum_{i=1}^{N}|z_i-\bar z|, \qquad R_q = \sqrt{\frac{1}{N}\sum_{i=1}^{N}(z_i-\bar z)^2}$$
تشمل القوى الفعلية أيضاً القوى الكهروستاتيكية والشعرية وقوى الالتصاق، ويمر تحليل منحنيات القوة الكمي عبر نماذج Hertz أو DMT أو JKR، فليس Lennard-Jones نموذجاً شاملاً لكل تفاعلات AFM.
تنبيه حاسم حول الخشونة: \(R_a\) و\(R_q\) ليستا خاصيتين جوهريتين مستقلتين عن القياس؛ فكلتاهما تعتمدان على مساحة المسح وكثافة البكسل، ومقارنة عينتين بمساحتي مسح مختلفتين بلا معنى فيزيائي. الخشونة خاصية مرتبطة بالمقياس (scale-dependent)، وتُبلَّغ دائماً مع مساحة المسح وحجم البكسل وظروف القياس، وفق معيار ISO 25178.
7.2 تأثير شكل رأس المجس
$$I(x,y) = S(x,y)\oplus T(x,y) \qquad \text{(عملية Dilation المورفولوجية)}$$
استرجاع الشكل الحقيقي \(S(x,y)\) ليس مجرد عملية Erosion عكسية بسيطة؛ إنها مسألة عكسية غير مثالية (ill-posed). خوارزمية Villarrubia تُقدِّر الشكل الحقيقي، لكنها تعتمد على معرفة أو تقدير مسبق لهندسة الرأس، ولا تضمن استرجاعاً كاملاً خصوصاً عند معالم أصغر من نصف قطر انحناء الرأس. تتطلب المعالجة الموثوقة معايرة الرأس، واستخدام رؤوس مرجعية معروفة الهندسة، ومقارنة النتائج مع TEM أو SEM عند الحاجة.
8. مطيافية الامتصاص UV-Vis-NIR
8.1 قانون Beer-Lambert والتصحيح الانعكاسي
$$A = \log_{10}\!\left(\frac{I_0}{I}\right) = \varepsilon\,c\,t, \qquad \alpha = \frac{2.303\,A}{t}$$
تفترض هذه الصيغة انعكاساً مهملاً (\(R\approx0\))، وهو افتراض يفشل غالباً في الأغشية والصفائح. عند قياس النفاذية \(T\) والانعكاس \(R\) معاً، تُستخرج \(\alpha\) من العلاقة غير المتداخلة الكاملة:
$$T = \frac{(1-R)^2\,e^{-\alpha t}}{1 – R^2 e^{-2\alpha t}}$$
$$\boxed{\alpha = \frac{1}{t}\ln\!\left[\frac{\sqrt{(1-R)^4 + 4R^2T^2} + (1-R)^2}{2T}\right]}$$
وتختزل إلى الصيغة الشائعة عند تحقق الشرط \(R^2e^{-2\alpha t}\ll1\)، أي في حالة امتصاص معتدل أو انعكاسية منخفضة:
$$\alpha \approx \frac{1}{t}\ln\!\left[\frac{(1-R)^2}{T}\right]$$
يجب تصحيح \(T\) بقياس عينة الأساس (substrate) فارغة أولاً. في الأغشية الرقيقة المتماسكة التي تظهر فيها أهداب التداخل (interference fringes)، تُستخدم طريقة المغلف لـ Swanepoel لاستخراج معامل الانكسار \(n(\lambda)\) ومعامل الامتصاص \(\alpha\) وسماكة الغشاء في آن واحد. وفي القياس بالمجال المشتت، يجب فصل المكوّن المرتد المنتشر باستخدام كرة تكامل (integrating sphere).
8.2 العينات المسحوقة: تحويل Kubelka-Munk
$$F(R_\infty) = \frac{(1-R_\infty)^2}{2R_\infty} = \frac{K}{S}$$
الشرط الحاسم الذي يُغفل كثيراً: يساوي \(F(R_\infty)\) النسبة \(K/S\)، وبما أن \(K\) يتناسب مع \(\alpha\)، فإن \(F(R_\infty)\) يتناسب مع \(\alpha\) فقط إذا كان معامل التشتت \(S\) مستقلاً عن الطول الموجي. هذا الشرط يُخترق عادة قرب حافة الامتصاص، حيث يتغير \(S(\lambda)\) مع البنية الدقيقة، مما يُحيِّز قيمة \(E_g\) المستخرجة من مخططات \([F(R_\infty)h\nu]^{1/n}\) بشكل منهجي.
من توصيات التخفيف: استخدام مرجع عاكس موحد مثل BaSO₄ أو MgO، مع طحن وكبس متطابقين بين جميع العينات؛ واعتبار المقارنات بين العينات نسبية تحت شروط تحضير متطابقة لا مطلقة؛ وعند الحاجة إلى دقة مطلقة، اللجوء إلى نماذج التشتت المنتشر المتقدمة بدل K-M المبسط. تبقى صلاحية K-M نفسها مشروطة بتشتت متعدد متماثل الاتجاه تقريباً، وعمق شبه لانهائي للعينة، ونسبة امتصاص إلى تشتت ضمن مجال معقول.
8.3 مخطط Tauc لفجوة الطاقة البصرية
$$[\alpha h\nu]^{1/n} = B(h\nu-E_g) \qquad \text{أو} \qquad [F(R_\infty)h\nu]^{1/n} = B(h\nu-E_g)$$
يُحدَّد الأس \(n=\tfrac12\) للانتقال المباشر المسموح، و\(n=2\) للانتقال غير المباشر المسموح، من طبيعة بنية النطاقات الإلكترونية للمادة نفسها، لا من أفضل خطية رسومية فحسب؛ ويتطلب دعماً مستقلاً من حسابات DFT أو الأدبيات المرجعية. معلمة الانحناء \(B\) ليست ثابتة تافهة: قيمتها النسبية بين العينات عند نفس \(n\) تعكس الاضطراب البنيوي وتفاوت المدى القصير، وينبغي أن تُبلَّغ مع \(E_g\) لا أن تُهمَل. وفي مواد ثنائية الأبعاد مثل MoS₂، تُعقّد رنانات الإكسيتون استخراج قيمة \(E_g\) أحادية بسيطة.
8.4 طاقة أورباخ
$$\alpha = \alpha_0\exp\left(\frac{h\nu-E_0}{E_U}\right)$$
حيث \(E_0\) طاقة مرجعية تُختار عادة قريبة من حافة الامتصاص. اختيار \(E_0\) تعسفياً يؤثر على قيمة \(\alpha_0\) المستخرجة، لكنه لا يؤثر على \(E_U\) نفسها، لأنها تُستخرج من الميل فقط:
$$\ln\alpha = \ln\alpha_0 + \frac{h\nu – E_0}{E_U}$$
لا تمثل \(E_U\) مقياساً معزولاً لنقص التبلور، بل تتأثر بأنواع متعددة من الاضطراب البنيوي والعيوب والاهتزازات الحرارية معاً.
9. مصفوفة التكامل بين التقنيات
| التقنية | المعلومة الأساسية | النطاق المكاني أو الفيزيائي | تحقق متبادل نموذجي |
|---|---|---|---|
| XRD | الأطوار، البنية، حجم المجالات، الإجهاد | بلوري إلى متوسط | TEM، Raman |
| SEM | المورفولوجيا السطحية | نانومتر إلى ميكرومتر | TEM، AFM |
| TEM | البنية الداخلية النانوية، العيوب | دون-نانومتري | XRD، EDS |
| EDS | التركيب العنصري الكمي | حجم التفاعل الإلكتروني، ميكرومترات في العينات السميكة | XPS |
| XPS | الكيمياء السطحية، حالات التأكسد | سطحي، بعمق \(\approx 3\lambda_{EAL}\cos\theta\) | EDS، Raman |
| Raman | الاهتزازات، الإجهاد، العيوب | محلي، بحجم بقعة الليزر | XRD، XPS |
| AFM | الطبوغرافيا، الخشونة | سطحي ثلاثي الأبعاد، مرتبط بمقياس المسح | SEM، TEM |
| UV-Vis-NIR | الخواص البصرية، فجوة الطاقة | إلكتروني وبصري | XRD، Raman |
10. من البيانات إلى الاستنتاج: لماذا لا يكفي توافق النتائج
توافق نتيجتين ظاهرياً لا يثبت العلاقة السببية بينهما.
المثال الأول: انخفاض عرض قمة XRD لا يعني تلقائياً زيادة حجم الحبيبة؛ فالاتساع يتأثر أيضاً بالإجهاد المجهري والعيوب والاتساع الآلي وطبيعة ملف القمة نفسها. عزل مساهمة الحجم وحدها يتطلب تحليلاً كاملاً بطريقة Williamson-Hall أو Warren-Averbach، لا قراءة اتساع القمة مباشرة.
المثال الثاني: انخفاض \(E_g\) لا يعني تلقائياً زيادة العيوب. فقد ينتج هذا الانخفاض عن تغيّر الطور، أو تغيّر التركيب الكيميائي، أو تأثير الحصر الكمي (Quantum Confinement) في الأبعاد النانوية، أو الإجهاد، أو تغيّر الترتيب البنيوي. ومن التفسيرات البديلة أيضاً تأثير Burstein-Moss: في المواد شديدة التطعيم يمتلئ أسفل نطاق التوصيل، فيتسع \(E_g\) الظاهري بصورة مصطنعة؛ واختلاف مستوى التطعيم بين عينتين يُغيّر مقدار هذا الاتساع الظاهري، فقد يُقرأ خطأً كتغيّر في كثافة العيوب بينما هو في الحقيقة فرق في انزياح Burstein-Moss نفسه. كل هذه التفسيرات البديلة يجب استبعادها أو مناقشتها قبل الجزم بالعيوب كسبب وحيد.
$$\boxed{\text{Measurement} \rightarrow \text{Model} \rightarrow \text{Validation} \rightarrow \text{Uncertainty} \rightarrow \text{Physical Interpretation}}$$
11. سلسلة الربط الفيزيائي
$$\boxed{\text{Processing} \rightarrow \text{Structure} \rightarrow \text{Chemistry} \rightarrow \text{Morphology} \rightarrow \text{Defects} \rightarrow \text{Properties}}$$
مثال تطبيقي: معالجة حرارية معينة تُغيّر التبلور وتركيب الأطوار، فتظهر آثارها عبر XRD في الأطوار واتساع القمم، وعبر TEM في المجالات المتماسكة، وعبر SEM في مورفولوجيا الجسيمات، وعبر Raman في الأنماط الاهتزازية والعيوب، وعبر XPS في حالات التأكسد، وعبر UV-Vis في \(E_g\) أو \(E_U\). ولتحويل هذا الترابط إلى علاقة سببية فعلية، يلزم مبدأ تصميم التجارب: تغيير متغير معالجة واحد مع تثبيت بقية المتغيرات. عندها تصبح التقنيات أدوات مستقلة لاختبار الفرضية الفيزيائية نفسها، لا قياسات معزولة.
12. تصنيف مصادر الخطأ وانتشار عدم اليقين الكمي
12.1 أخطاء تحضير العينة
ضرر FIB، والانفعال الميكانيكي الناتج عن طحن المساحيق، والتلوث السطحي، وتغيّر التركيب أثناء التحضير.
12.2 أخطاء حجم التفاعل وأداة القياس
حجم تفاعل SEM-EDS، وتشوه رأس AFM، وعمق تحليل XPS الذي يُقدَّر بـ \(3\lambda_{EAL}\cos\theta\)، والامتصاص والتشتت في القياسات البصرية.
12.3 أخطاء النموذج الرياضي
استخدام Scherrer كقياس مباشر لحجم الحبيبة، وتطبيق Williamson-Hall على مادة شديدة عدم التناظر دون تبرير، واستخدام Tuinstra-Koenig خارج نطاق صلاحيته، وتطبيق Tauc دون أساس فيزيائي لاختيار \(n\)، وتفسير طاقة أورباخ كمقياس وحيد للعيوب، واستخدام Kubelka-Munk دون التحقق من افتراضاته، والإبلاغ عن نتائج بلا هامش عدم يقين وبلا فترة ملاءمة معلنة.
12.4 أخطاء المعايرة والمرجع
عدم تصحيح الاتساع الآلي في XRD، وخطأ معايرة طول الكاميرا في TEM، وانحراف طاقة XPS بسبب الشحن، وعدم دقة عوامل الحساسية، ومعايرة غير كافية لكاشف EDS، وجهل هندسة رأس AFM.
12.5 القانون العام لانتشار عدم اليقين
لكمية \(f(x_1,\dots,x_n)\)، وعندما تكون المتغيرات غير مترابطة:
$$\boxed{\delta f = \sqrt{\sum_{i=1}^{n}\left(\frac{\partial f}{\partial x_i}\,\delta x_i\right)^2}}$$
وعند وجود ارتباطات، تُضاف الحدود التبادلية \(2\sum_{i المدخلات: قمة عند \(2\theta=40°\) بأشعة Cu Kα (\(\lambda=1.5406\) Å)، \(\beta_{obs}=0.30°\)، \(\beta_{inst}=0.08°\)، عدم يقين الملاءمة \(\delta\beta_{obs}=\delta\beta_{inst}=\pm0.01°\)، وعدم يقين موضع القمة \(\delta(2\theta)=\pm0.02°\). الخطوة الأولى، التصحيح الآلي عبر \(\beta_s=\sqrt{\beta_{obs}^2-\beta_{inst}^2}\): $$\frac{\partial\beta_s}{\partial\beta_{obs}}=\frac{\beta_{obs}}{\beta_s}=1.038, \qquad \frac{\partial\beta_s}{\partial\beta_{inst}}=-\frac{\beta_{inst}}{\beta_s}=-0.277$$ $$\delta\beta_s = \sqrt{(1.038\times0.01)^2+(0.277\times0.01)^2} = \pm0.011° \quad\Rightarrow\quad \frac{\delta\beta_s}{\beta_s}\approx3.7\%$$ الخطوة الثانية، القيمة المركزية عند \(\beta_s = 0.289° = 5.05\times10^{-3}\) راديان: $$D = \frac{0.9\times1.5406}{5.05\times10^{-3}\times\cos20°} \approx 29.2\ \text{nm}$$ الخطوة الثالثة، ميزانية عدم اليقين الكاملة من العلاقة \(\delta D/D=\sqrt{(\delta K/K)^2+(\delta\beta/\beta)^2+(\tan\theta\,\delta\theta)^2+(\delta\lambda/\lambda)^2}\): $$\boxed{D = (29 \pm 2)\ \text{nm}}$$ يهيمن على عدم اليقين هنا ملاءمة الملف وعدم يقين \(K\) المنهجي، بينما يبقى حد الزاوية مهملاً عند الزوايا المتوسطة، وينمو مع \(\tan\theta\) عند الزوايا العالية. القاعدة الذهبية: رقم واحد بلا هامش عدم يقين ليس قياساً، بل ادعاء. $$R_{wp} = \sqrt{\frac{\sum_i w_i\,(y_i^{obs}-y_i^{calc})^2}{\sum_i w_i\,(y_i^{obs})^2}}, \qquad R_{exp} = \sqrt{\frac{N-P}{\sum_i w_i\,(y_i^{obs})^2}}, \qquad \chi^2 = \left(\frac{R_{wp}}{R_{exp}}\right)^2$$ لا يجوز مقارنة الأرقام التالية باعتبارها «حدود كشف» من النوع نفسه؛ فبعضها يمثل حد كشف تركيبي، وبعضها حساسية سطحية، وبعضها دقة طبوغرافية أو فوتومترية. وهي جميعاً رتب تقريبية تعتمد على الجهاز والعينة والتحضير وزمن القياس والخلفية ونموذج التحليل. القاعدة المنهجية: لا يمثل «حد الكشف» خاصية ثابتة للتقنية وحدها، بل نتيجة لتفاعل التقنية مع العينة والجهاز وظروف القياس ونموذج التحليل. كما أن غياب الإشارة تحت حد الكشف لا يثبت غياب المكوّن، بل يضع قيداً على مقدار ما يمكن اكتشافه ضمن شروط القياس المحددة. تجيب XRD عن سؤال البنية البلورية، وتجيب SEM وTEM عن سؤال المورفولوجيا والبنية الداخلية، وتجيب EDS وXPS عن سؤال التركيب والكيمياء، وتجيب Raman عن سؤال الروابط والإجهاد والاضطراب، وتجيب AFM عن سؤال الطبوغرافيا، وتجيب UV-Vis-NIR عن سؤال الخواص البصرية والإلكترونية. التوصيف المتقدم للمواد ليس مجموعة أجهزة أو معادلات منفصلة، بل منظومة استدلال فيزيائي متعددة الأدلة. لا تُقاس قوتها الحقيقية بعدد التقنيات المستخدمة، بل بالالتزام بحدود صلاحية كل نموذج، وبمعرفة وحدات واتفاقيات كل صيغة كوحدات Kanaya-Okayama، وبالتمييز بين الكميات المتشابهة والمختلفة مثل IMFP مقابل EAL، أو \(F(R_\infty)\) مقابل \(\alpha\). والأهم من كل ذلك: رفض القفز من توافق النتائج الظاهري إلى استنتاج سببي، ورفض تقديم أي رقم بلا هامش عدم يقين، لأن الرقم بلا هامش ليس قياساً، بل ادعاء. تنويه: استُعين بأدوات الذكاء الاصطناعي في صياغة وتنظيم محتوى هذا المقال، مع مراجعة المعادلات والمفاهيم وحدود الصلاحية الفيزيائية قبل النشر. ونرحب بأي ملاحظات أو تصحيحات علمية تسهم في إثراء المحتوى وتحسين دقته.12.6 مثال محسوب بالكامل: عدم اليقين في معادلة Scherrer
المصدر القيمة المساهمة في \(\delta D/D\) عرض القمة المصحح \(\delta\beta_s\) ±0.011° 3.7% موضع القمة \(\delta(2\theta)=\pm0.02°\) \(\tan20°\times1.75\times10^{-4}\) 0.006% (مهمل) الطول الموجي \(\delta\lambda/\lambda\) نحو \(10^{-5}\) مهمل المعامل \(K\) (منهجي: شكل البلورية، التوزيع الحجمي) ±5% 5% الإجمالي التربيعي — ≈ 6.2%
12.7 مؤشرات جودة التحليل في Rietveld وأمثاله
12.8 الحدود النموذجية للكشف والحساسية والدقة
التقنية المقياس أو النطاق النموذجي المعنى الفيزيائي أهم القيود XRD نحو 1%–5% وزناً لطور بلوري في القياس المخبري التقليدي حد كشف تقريبي لطور بلوري ثانوي يعتمد بشدة على التبلور، وشدة القمم، والخلفية، وتداخل القمم، وطبيعة الطور. يمكن الوصول إلى حدود أقل في ظروف محسنة، لذلك لا يمثل الرقم حدًا فيزيائيًا ثابتًا. EDS-SEM نحو 0.1%–1% وزناً في كثير من التطبيقات العملية، مع إمكانية الوصول إلى حدود أقل لبعض العناصر في ظروف محسنة حد كشف عنصري في منطقة التفاعل يعتمد على العنصر، وطاقة الحزمة، وزمن الجمع، ومعدل العد، والخلفية، والكاشف، وتصحيحات ZAF أو نماذج φ(ρz). تكون العناصر الخفيفة أكثر صعوبة. EDS-TEM لا يوجد حد تركيبي عام ثابت صالح لجميع الأنظمة تحليل عنصري موضعي على مقياس نانوي يتأثر بسماكة العينة، والكاشف، والعنصر، وزمن القياس، والخلفية، وتداخل القمم، وشبكة الحامل. الميزة الأساسية هي الدقة المكانية النانوية وليس رقم حد كشف موحد. XPS نحو 0.1%–1% ذرياً في كثير من التطبيقات حد كشف تقريبي للتركيب الكيميائي السطحي يعتمد على المقطع العرضي، وشدة القمة، والخلفية، وتداخل القمم، والشحن، وهندسة القياس. عمق التحليل يجب تقديره باستخدام EAL المناسبة وليس IMFP وحده. Raman لا يوجد حد كشف تركيزي عام لـ Raman التقليدي حساسية طيفية تعتمد على شدة الإشارة ونسبة الإشارة إلى الضوضاء تتأثر بشدة بالفلورية، وطول موجة الليزر، والقدرة، وزمن التكامل، وحجم البقعة، والمادة. يمكن لتقنيات مثل Resonance Raman وSERS تحقيق حساسية أعلى بكثير في تطبيقات متخصصة. AFM حساسية رأسية قد تصل إلى رتبة 0.01–0.1 nm في أنظمة وظروف مناسبة حساسية قياس الارتفاع، وليس حد كشف تركيزي تعتمد الدقة الجانبية على هندسة الرأس، والتفاعل مع السطح، وحجم البكسل، ومساحة المسح، والاستقرار الميكانيكي والحراري. نصف قطر الرأس لا يساوي بالضرورة الدقة الجانبية. UV-Vis-NIR قد تصل الحساسية/الدقة الفوتومترية إلى رتبة 10⁻³–10⁻⁴ Abs في أجهزة وظروف مناسبة أصغر تغير قابل للقياس في الامتصاص تحت شروط محددة لا يمثل هذا الرقم حد كشف تركيزي. يعتمد حد الكشف الفعلي على ε وطول المسار والضوضاء والاستقرار الضوئي وstray light وطريقة المعايرة.
13. الإطار المنهجي الكامل: خمس مراحل
الخلاصة
