دالة نقل التباين

طيف القدرة (تحويل فورييه) لصورة مجهرية إلكترونية نموذجية. يمكن ملاحظة تأثير دالة نقل التباين في الحلقات المتناوبة بين الضوء والظلام (حلقات ثون)، والتي توضح العلاقة بين التباين والتردد المكاني.

تصف دالة نقل التباين ( CTF ) رياضياً كيف تُعدّل الانحرافات في المجهر الإلكتروني النافذ (TEM) صورة العينة. [ 1 ] [ 2 ] [ 3 ] [ 4 ] وتحدد دالة نقل التباين (CTF) هذه دقة المجهر الإلكتروني النافذ عالي الدقة (HRTEM)، المعروف أيضاً باسم المجهر الإلكتروني النافذ ذي التباين الطوري.

باعتبار الصورة المسجلة جسماً حقيقياً متأثراً بتشوه دالة نقل التباين (CTF)، فإن وصف هذه الدالة يسمح بإعادة هندسة الجسم الحقيقي . يُشار إلى هذا عادةً بتصحيح دالة نقل التباين، وهو أمر بالغ الأهمية للحصول على هياكل عالية الدقة في المجهر الإلكتروني ثلاثي الأبعاد، وخاصةً المجهر الإلكتروني المبرد . ويُعادل هذا المفهوم في البصريات الضوئية دالة النقل البصري .

تباين الطور في المجهر الإلكتروني عالي الدقة

ينشأ التباين في المجهر الإلكتروني عالي الدقة (HRTEM) من التداخل في مستوى الصورة بين أطوار موجات الإلكترونات المتناثرة وطور موجة الإلكترونات النافذة. تحدث تفاعلات معقدة عندما تمر موجة إلكترونية عبر عينة في المجهر الإلكتروني النافذ. فوق العينة، يمكن تقريب موجة الإلكترونات كموجة مستوية. أثناء مرور موجة الإلكترونات، أو دالة الموجة ، عبر العينة، يتغير كل من طور وسعة حزمة الإلكترونات. ثم يتم تركيز حزمة الإلكترونات المتناثرة والنافذة الناتجة بواسطة عدسة موضوعية، وتصويرها بواسطة كاشف في مستوى الصورة.

لا تستطيع أجهزة الكشف قياس سوى السعة، وليس الطور مباشرةً. مع ذلك، وباستخدام معايير المجهر الصحيحة، يمكن قياس تداخل الطور بشكل غير مباشر عبر شدة الضوء في مستوى الصورة. تتفاعل الإلكترونات بقوة شديدة مع المواد الصلبة البلورية . ونتيجةً لذلك، يمكن تسجيل تغيرات الطور الناتجة عن تفاصيل دقيقة للغاية، وصولاً إلى المستوى الذري، باستخدام المجهر الإلكتروني عالي الدقة (HRTEM).

نظرية نقل التباين

مخطط الأشعة في المجهر الإلكتروني النافذ مع دالة نقل تباين الطور

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

Contrast transfer theory consists of four main operations:[1]

  1. Take the Fourier transform of the exit wave to obtain the wave amplitude in back focal plane of objective lens
  2. Modify the wavefunction in reciprocal space by a phase factor, also known as the Phase Contrast Transfer Function, to account for aberrations
  3. Inverse Fourier transform the modified wavefunction to obtain the wavefunction in the image plane
  4. Find the square modulus of the wavefunction in the image plane to find the image intensity (this is the signal that is recorded on a detector, and creates an image)

Mathematical form

If we incorporate some assumptions about our sample, then an analytical expression can be found for both phase contrast and the phase contrast transfer function. As discussed earlier, when the electron wave passes through a sample, the electron beam interacts with the sample via scattering, and experiences a phase shift. This is represented by the electron wavefunction exiting from the bottom of the sample. This expression assumes that the scattering causes a phase shift (and no amplitude shift). This is called the Phase Object Approximation.

The exit wavefunction

Following Wade's notation,[1] the exit wavefunction expression is represented by:

τ(r,z)=τoexp[iπλdzU(r,z)]{\displaystyle \tau (r,z)=\tau _{o}\exp[-i\pi \lambda \int dz'U(r,z')]}
τo=τ(r,0){\displaystyle \tau _{o}=\tau (r,0)}
U(r,z)=2mV(r,z)/h2{\displaystyle U(r,z)=2mV(r,z)/h^{2}}

Where the exit wavefunction τ is a function of both r{\displaystyle r} in the plane of the sample, and z{\displaystyle z} perpendicular to the plane of the sample. τo{\displaystyle \tau _{o}} represents the wavefunction incident on the top of the sample. λ{\displaystyle \lambda } is the wavelength of the electron beam,[5] which is set by the accelerating voltage. U{\displaystyle U} is the effective potential of the sample, which depends on the atomic potentials within the crystal, represented by V{\displaystyle V}.

Within the exit wavefunction, the phase shift is represented by:

ϕ(r)=πλdzU(r,z){\displaystyle \phi (r)=\pi \lambda \int dz'U(r,z')}

This expression can be further simplified taken into account some more assumptions about the sample. If the sample is considered very thin, and a weak scatterer, so that the phase shift is << 1, then the wave function can be approximated by a linear Taylor polynomial expansion.[6] This approximation is called the Weak Phase Object Approximation.

The exit wavefunction can then be expressed as:

τ(r,z)=τo[1+iϕ(r)]{\displaystyle \tau (r,z)=\tau _{o}[1+i\phi (r)]}

The phase contrast transfer function

Passing through the objective lens incurs a Fourier transform and phase shift. As such, the wavefunction on the back focal plane of the objective lens can be represented by:

I(θ)=δ(θ)+ΦK(θ){\displaystyle I(\theta )=\delta (\theta )+\Phi K(\theta )}

θ{\displaystyle \theta } = the scattering angle between the transmitted electron wave and the scattered electron wave

δ{\displaystyle \delta } = a delta function representing the non-scattered, transmitted, electron wave

Φ{\displaystyle \Phi } = the Fourier transform of the wavefunction's phase

K(θ){\displaystyle K(\theta )}= إزاحة الطور الناتجة عن انحرافات المجهر، والمعروفة أيضًا باسم دالة نقل التباين:

ك(θ)=الخطيئة[(2π/λ)دبليو(θ)]{\displaystyle K(\theta )=\sin[(2\pi /\lambda )W(\theta )]}دبليو(θ)=-zθ2/2+جsθ4/4{\displaystyle W(\theta )=-z\theta ^{2}/2+C_{s}\theta ^{4}/4}

λ{\displaystyle \lambda }= الطول الموجي النسبي لموجة الإلكترون،جs{\displaystyle C_{s}}= الانحراف الكروي للعدسة الشيئية

يمكن أيضًا التعبير عن دالة نقل التباين بدلالة الترددات المكانية، أو الفضاء المتبادل. مع العلاقةθ=λك{\textstyle \theta =\lambda k}تصبح دالة نقل تباين الطور كما يلي:

ك(ك)=الخطيئة[(2π)دبليو(ك)]{\displaystyle K(k)=\sin[(2\pi )W(k)]}دبليو(ك)=-zλك2/2+جsλ3ك4/4{\displaystyle W(k)=-z\lambda k^{2}/2+C_{s}\lambda ^{3}k^{4}/4}

z{\displaystyle z}= عدم وضوح العدسة الشيئية (باستخدام الاصطلاح القائل بأن عدم وضوح الصورة موجب وزيادة وضوحها سالبة)،λ{\displaystyle \lambda }= الطول الموجي النسبي لموجة الإلكترون،جs{\displaystyle C_{s}}= الانحراف الكروي للعدسة الشيئية،ك{\displaystyle k}= التردد المكاني (بوحدات م −1 )

الانحراف الكروي

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

رs=جsθ3م{\displaystyle r_{s}=C_{s}\cdot \theta ^{3}\cdot M}

أينجs{\displaystyle C_{s}}هو الانحراف الكروي وم{\displaystyle M}يمثل التكبير، وكلاهما ثابتان فعليًا لإعدادات العدسة. ويمكن للمرء أن يلاحظ بعد ذلك أن الفرق في زاوية الانكسار بين شعاع مثالي وشعاع يعاني من الزيغ الكروي هو

αs=دالة الظل العكسي(بR)-دالة الظل العكسي(بR+رs){\displaystyle \alpha _{s}=\arctan \left({\frac {b}{R}}\right)-\arctan \left({\frac {b}{R+r_{s}}}\right)}

أينب{\displaystyle b}هي المسافة من العدسة إلى مستوى الصورة الغاوسية وR{\displaystyle R}هي المسافة الشعاعية من المحور البصري إلى النقطة على العدسة التي مر بها الشعاع. وبتبسيط ذلك أكثر (دون تطبيق أي تقريبات) يتضح أن

αs=دالة الظل العكسي(برsR2+Rرs+ب2){\displaystyle \alpha _{s}=\arctan \left({\frac {br_{s}}{R^{2}+Rr_{s}+b^{2}}}\right)}

يمكن الآن تطبيق تقريبين للمضي قدمًا بطريقة مباشرة. يعتمدان على افتراض أن كليهمارs{\displaystyle r_{s}}وR{\displaystyle R}أصغر بكثير منب{\displaystyle b}وهذا يُعادل القول بأننا نأخذ في الاعتبار زوايا سقوط صغيرة نسبيًا، وبالتالي انحرافات كروية صغيرة جدًا. وبناءً على هذا الافتراض، يكون الحدّان الرئيسيان في المقام ضئيلين، ويمكن تقريبهما على أنهما لا يُساهمان في النتيجة. ومن خلال هذه الافتراضات، ذكرنا ضمنيًا أيضًا أن الكسر نفسه يُمكن اعتباره صغيرًا، وهذا يؤدي إلى حذفه.دالة الظل العكسي(){\displaystyle \arctan()}الدالة عن طريق تقريب الزاوية الصغيرة؛

αsدالة الظل العكسي(برsب2)برsب2=رsب=جsθ3مب{\displaystyle \alpha _{s}\approx \arctan \left({\frac {br_{s}}{b^{2}}}\right)\approx {\frac {br_{s}}{b^{2}}}={\frac {r_{s}}{b}}={\frac {C_{s}\cdot \theta ^{3}\cdot M}{b}}}

إذا اعتُبرت الصورة واضحة تقريبًا، فإن زاوية السقوطθ{\displaystyle \theta }ثم يُعتبر صغيراً مرة أخرى،

Rولون برونزي(θ)θ  و  مبو{\displaystyle {\frac {R}{f}}\approx \tan \left(\theta \right)\approx \theta ~~{\text{and}}~~M\approx {\frac {b}{f}}}

وهذا يعني أن التعبير التقريبي للفرق في زاوية الانكسار بين شعاع مثالي وشعاع يعاني من الزيغ الكروي يُعطى بالصيغة التالية:

αsجsR3و4{\displaystyle \alpha _{s}\approx {\frac {C_{s}\cdot R^{3}}{f^{4}}}}

عدم التركيز

على عكس الزيغ الكروي، سنتابع بتقدير انحراف الشعاع غير المركز عن الوضع المثالي من خلال تحديد الزيغ الطولي؛ وهو مقياس لمدى انحراف الشعاع عن النقطة البؤرية على طول المحور البصري.Δب{\displaystyle \Delta b}، من الممكن إثبات أن الفرقαو{\displaystyle \alpha _{f}}يمكن ربط زاوية الانكسار بين الأشعة الصادرة من جسم مركز وجسم غير مركز بزاوية الانكسار كما يلي:

R2+ب2الخطيئة(αو)=Δبالخطيئة(θ-αو){\displaystyle {\sqrt {R^{2}+b^{2}}}\cdot \sin(\alpha _{f})=\Delta b\cdot \sin(\theta '-\alpha _{f})}

أينR{\displaystyle R}وب{\displaystyle b}يتم تعريفها بنفس الطريقة التي تم تعريفها بها بالنسبة للانحراف الكروي. بافتراض أنαو<<θ{\displaystyle \alpha _{f}<<\theta '}(أو ما يعادل ذلك)|بالخطيئة(αو)|<<|R|{\displaystyle |b\cdot \sin(\alpha _{f})|<<|R|}يمكننا أن نبين أن

الخطيئة(αو)Δبالخطيئة(θ)R2+ب2=ΔبRR2+ب2{\displaystyle \sin(\alpha _{f})\approx {\frac {\Delta b\sin(\theta ')}{\sqrt {R^{2}+b^{2}}}}={\frac {\Delta b\cdot R}{R^{2}+b^{2}}}}

بما أننا كنا بحاجةαو{\displaystyle \alpha _{f}}أن تكون صغيرة، ولأنθ{\displaystyle \theta }كون الشيء صغيرًا يعنيR<<ب{\displaystyle R<<b}، لدينا تقريبًا لـαو{\displaystyle \alpha _{f}}مثل

αوΔبRب2{\displaystyle \alpha _{f}\approx {\frac {\Delta b\cdot R}{b^{2}}}}

يمكن إثبات ذلك من خلال صيغة العدسة الرقيقة.Δب/ب2Δو/و2{\displaystyle \Delta b/b^{2}\approx \Delta f/f^{2}}مما ينتج عنه تقدير نهائي للفرق في زاوية الانكسار بين الأشعة البؤرية والأشعة غير البؤرية كما يلي:

αوΔوRو2{\displaystyle \alpha _{f}\approx {\frac {\Delta f\cdot R}{f^{2}}}}

أمثلة

تحدد دالة نقل التباين مقدار إشارة الطور التي تُنقل إلى دالة الموجة في الفضاء الحقيقي في مستوى الصورة. وبما أن مربع القيمة المطلقة لدالة الموجة في الفضاء الحقيقي يُعطي إشارة الصورة، فإن دالة نقل التباين تحد من كمية المعلومات التي يمكن ترجمتها في النهاية إلى صورة. ويحدد شكل دالة نقل التباين جودة تكوين الصورة في الفضاء الحقيقي في المجهر الإلكتروني النافذ.

تم إعداد وظيفة CTF عبر تطبيق ويب صغير أنشأه جيانغ وشيو، وهو متاح على الرابط التالي: https://ctfsimulation.streamlit.app/

هذا مثال على دالة نقل التباين. هناك عدة أمور يجب ملاحظتها:

  • توجد هذه الدالة في مجال التردد المكاني ، أو فضاء k
  • عندما تكون قيمة الدالة مساوية للصفر، فهذا يعني عدم وجود نفاذية، أو عدم وجود إشارة طور مدمجة في صورة الفضاء الحقيقي.
  • تُسمى المرة الأولى التي تعبر فيها الدالة المحور السيني بنقطة التحليل.
  • لتحقيق أقصى قدر من إشارة الطور، من الأفضل عمومًا استخدام ظروف تصوير تدفع دقة النقطة إلى ترددات مكانية أعلى.
  • عندما تكون الدالة سالبة، فإن ذلك يمثل تباينًا طوريًا موجبًا، مما يؤدي إلى خلفية ساطعة، مع ميزات ذرية داكنة.
  • في كل مرة يتقاطع فيها منحنى CTF مع المحور السيني، يحدث انعكاس في التباين
  • وبناءً على ذلك، فإنه بعد دقة النقطة للمجهر، لا يمكن تفسير معلومات الطور بشكل مباشر، ويجب نمذجتها عبر محاكاة الحاسوب.

عدم وضوح شيرزر

قيمة عدم التركيز (z{\textstyle z}يمكن استخدام هذه التقنية لمعادلة الانحراف الكروي، مما يسمح بتباين طوري أكبر. وقد طوّر شيرزر هذا التحليل، ويُطلق عليه اسم "عدم تركيز شيرزر". [ 7 ]

zs=(جsλ)1/2{\displaystyle z_{s}=(C_{s}\lambda )^{1/2}}

المتغيرات هي نفسها المستخدمة في قسم المعالجة الرياضية، معzs{\displaystyle z_{s}}ضبط خاصية عدم التركيز الخاصة بـ Scherzer،جs{\displaystyle C_{s}}باعتبارها الانحراف الكروي، وλ باعتبارها الطول الموجي النسبي لموجة الإلكترون.

يوضح الشكل في القسم التالي دالة نقل التباين (CTF) لمجهر CM300 عند ضبط بؤرة شيرزر. بالمقارنة مع دالة نقل التباين الموضحة أعلاه، توجد نافذة أكبر، تُعرف أيضًا بنطاق التمرير، للترددات المكانية ذات النفاذية العالية. وهذا يسمح بمرور المزيد من إشارة الطور إلى مستوى الصورة.

دالة المغلف

دالة نقل التباين (CTF) لمجهر CM300 يتم تخميدها بواسطة دوال الغلاف الزمني والمكاني.

تمثل دالة الغلاف تأثير الانحرافات الإضافية التي تُضعف دالة نقل التباين، وبالتالي الطور. تميل حدود الغلاف المكونة لدالة الغلاف إلى كبح الترددات المكانية العالية. قد يختلف الشكل الدقيق لدوال الغلاف من مصدر لآخر. عمومًا، يتم تطبيقها بضرب دالة نقل التباين بحد غلاف Et يمثل الانحرافات الزمنية، وحد غلاف Es يمثل الانحرافات المكانية. ينتج عن ذلك دالة نقل تباين معدلة، أو فعالة.

كهـوو(ك)=هـتهـs(الخطيئة[(2π/λ)دبليو(ك)]{\displaystyle K_{eff}(k)=E_{t}E_{s}(\sin[(2\pi /\lambda )W(k)]}

تشمل أمثلة الانحرافات الزمنية الانحرافات اللونية، وانتشار الطاقة، وانتشار البؤرة، وعدم استقرار مصدر الجهد العالي، وعدم استقرار تيار العدسة الشيئية. ومن أمثلة الانحرافات المكانية تقارب الشعاع الساقط المحدود. [ 8 ]

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

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

تم تطوير برامج متنوعة لنمذجة كل من دالة نقل التباين ودالة الغلاف لأنواع معينة من المجاهر، ومعايير تصوير محددة. [ 9 ] [ 10 ]

نظرية التصوير الخطي مقابل نظرية التصوير غير الخطي

The previous description of the contrast transfer function depends on linear imaging theory. Linear imaging theory assumes that the transmitted beam is dominant, there is only weak phase shift by the sample. In many cases, this precondition is not fulfilled. In order to account for these effects, non-linear imaging theory is required. With strongly scattering samples, diffracted electrons will not only interfere with the transmitted beam, but will also interfere with each other. This will produce second order diffraction intensities. Non-linear imaging theory is required to model these additional interference effects.[11][12]

Contrary to a widespread assumption, the linear/nonlinear imaging theory has nothing to do with kinematical diffraction or dynamical diffraction, respectively.

Linear imaging theory is still used, however, because it has some computational advantages. In Linear imaging theory, the Fourier coefficients for the image plane wavefunction are separable. This greatly reduces computational complexity, allowing for faster computer simulations of HRTEM images.[13]

See also

References

  1. 123Wade, R. H. (October 1992). "A brief look at imaging and contrast transfer". Ultramicroscopy. 46 (1–4): 145–156. doi:10.1016/0304-3991(92)90011-8.
  2. Spence, John C. H. (1988 2nd ed) Experimental high-resolution electron microscopy (Oxford U. Press, NY) ISBN 0195054059.
  3. Ludwig Reimer (1997 4th ed) Transmission electron microscopy: Physics of image formation and microanalysis (Springer, Berlin) preview.
  4. Earl J. Kirkland (1998) Advanced computing in electron microscopy (Plenum Press, NY).
  5. "DeBroglie Wavelength". HyperPhysics. Georgia State University. Retrieved 27 April 2017.
  6. "Weak-phase-objects (WPO) in TEM observations - Practical Electron Microscopy and Database - An Online Book - EELS EDS TEM SEM". www.globalsino.com. Retrieved 2015-06-12.
  7. Scherzer (1949). "The theoretical resolution limit of the electron microscope". Journal of Applied Physics. 20 (1): 20–29. Bibcode:1949JAP....20...20S. doi:10.1063/1.1698233.
  8. "Envelope Functions". www.maxsidorov.com. Retrieved 2015-06-12.
  9. "CTF Simulation Web App". Wen Jiang Group. Retrieved 24 January 2026.
  10. Sidorov, Max. "Home of the ctfExplorer". Retrieved 27 April 2017.
  11. Bonevich, Marks (May 24, 1988). "Contrast Transfer Theory for Non-Linear Imaging". Ultramicroscopy. 26 (3): 313–319. doi:10.1016/0304-3991(88)90230-6.
  12. This page was prepared in part for Northwestern University class MSE 465, taught by Professor Laurie Marks.
  13. Notes prepared by Professor Laurie Marks at Northwestern University.