Contrast transfer function

Last updated October 30, 2024

The contrast transfer function (CTF) mathematically describes how aberrations in a transmission electron microscope (TEM) modify the image of a sample.^[1]^[2]^[3]^[4] This contrast transfer function (CTF) sets the resolution of high-resolution transmission electron microscopy (HRTEM), also known as phase contrast TEM.

By considering the recorded image as a CTF-degraded true object, describing the CTF allows the true object to be reverse-engineered. This is typically denoted CTF-correction, and is vital to obtain high resolution structures in three-dimensional electron microscopy, especially electron cryo-microscopy. Its equivalent in light-based optics is the optical transfer function.

Phase contrast in HRTEM

The contrast in HRTEM comes from interference in the image plane between the phases of scattered electron waves with the phase of the transmitted electron wave. Complex interactions occur when an electron wave passes through a sample in the TEM. Above the sample, the electron wave can be approximated as a plane wave. As the electron wave, or wavefunction, passes through the sample, both the phase and the amplitude of the electron beam is altered. The resultant scattered and transmitted electron beam is then focused by an objective lens, and imaged by a detector in the image plane.

Detectors are only able to measure the amplitude, not the phase directly. However, with the correct microscope parameters, the phase interference can be indirectly measured via the intensity in the image plane. Electrons interact very strongly with crystalline solids. As a result, the phase changes due to very small features, down to the atomic scale, can be recorded via HRTEM.

Contrast transfer theory

Contrast transfer theory provides a quantitative method to translate the exit wavefunction to a final image. Part of the analysis is based on Fourier transforms of the electron beam wavefunction. When an electron wavefunction passes through a lens, the wavefunction goes through a Fourier transform. This is a concept from Fourier optics.

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

Take the Fourier transform of the exit wave to obtain the wave amplitude in back focal plane of objective lens
Modify the wavefunction in reciprocal space by a phase factor, also known as the Phase Contrast Transfer Function, to account for aberrations
Inverse Fourier transform the modified wavefunction to obtain the wavefunction in the image plane
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:

\tau (r,z)=\tau _{o}\exp[-i\pi \lambda \int dz'U(r,z')]

\tau _{o}=\tau (r,0)

U(r,z)=2mV(r,z)/h^{2}

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

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

\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:

\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(\theta )=\delta (\theta )+\Phi K(\theta )

$\theta$ = the scattering angle between the transmitted electron wave and the scattered electron wave

$\delta$ = a delta function representing the non-scattered, transmitted, electron wave

$\Phi$ = the Fourier transform of the wavefunction's phase

$K(\theta )$ = the phase shift incurred by the microscope's aberrations, also known as the Contrast Transfer Function:

K(\theta )=\sin[(2\pi /\lambda )W(\theta )]

W(\theta )=-z\theta ^{2}/2+C_{s}\theta ^{4}/4

$\lambda$ = the relativistic wavelength of the electron wave, $C_{s}$ = The spherical aberration of the objective lens

The contrast transfer function can also be given in terms of spatial frequencies, or reciprocal space. With the relationship ${\textstyle \theta =\lambda k}$ , the phase contrast transfer function becomes:

K(k)=\sin[(2\pi )W(k)]

W(k)=-z\lambda k^{2}/2+C_{s}\lambda ^{3}k^{4}/4

$z$ = the defocus of the objective lens (using the convention that underfocus is positive and overfocus is negative), $\lambda$ = the relativistic wavelength of the electron wave, $C_{s}$ = The spherical aberration of the objective lens, $k$ = the spatial frequency (units of m⁻¹)

Spherical aberration

Spherical aberration is a blurring effect arising when a lens is not able to converge incoming rays at higher angles of incidence to the focus point, but rather focuses them to a point closer to the lens. This will have the effect of spreading an imaged point (which is ideally imaged as a single point in the gaussian image plane) out over a finite size disc in the image plane. Giving the measure of aberration in a plane normal to the optical axis is called a transversal aberration. The size (radius) of the aberration disc in this plane can be shown to be proportional to the cube of the incident angle (θ) under the small-angle approximation, and that the explicit form in this case is

r_{s}=C_{s}\cdot \theta ^{3}\cdot M

where $C_{s}$ is the spherical aberration and $M$ is the magnification, both effectively being constants of the lens settings. One can then go on to note that the difference in refracted angle between an ideal ray and one which suffers from spherical aberration, is

\alpha _{s}=\arctan \left({\frac {b}{R}}\right)-\arctan \left({\frac {b}{R+r_{s}}}\right)

where $b$ is the distance from the lens to the gaussian image plane and $R$ is the radial distance from the optical axis to the point on the lens which the ray passed through. Simplifying this further (without applying any approximations) shows that

\alpha _{s}=\arctan \left({\frac {br_{s}}{R^{2}+Rr_{s}+b^{2}}}\right)

Two approximations can now be applied to proceed further in a straightforward manner. They rely on the assumption that both $r_{s}$ and $R$ are much smaller than $b$ , which is equivalent to stating that we are considering relatively small angles of incidence and consequently also very small spherical aberrations. Under such an assumption, the two leading terms in the denominator are insignificant, and can be approximated as not contributing. By way of these assumptions we have also implicitly stated that the fraction itself can be considered small, and this results in the elimination of the $\arctan()$ function by way of the small-angle approximation;

\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}}

If the image is considered to be approximately in focus, and the angle of incidence $\theta$ is again considered small, then

{\frac {R}{f}}\approx \tan \left(\theta \right)\approx \theta ~~{\text{and}}~~M\approx {\frac {b}{f}}

meaning that an approximate expression for the difference in refracted angle between an ideal ray and one which suffers from spherical aberration, is given by

\alpha _{s}\approx {\frac {C_{s}\cdot R^{3}}{f^{4}}}

Defocus

As opposed to the spherical aberration, we will proceed by estimating the deviation of a defocused ray from the ideal by stating the longitudinal aberration; a measure of how much a ray deviates from the focal point along the optical axis. Denoting this distance $\Delta b$ , it is possible to show that the difference $\alpha _{f}$ in refracted angle between rays originating from a focused and defocused object, can be related to the refracted angle as

{\sqrt {R^{2}+b^{2}}}\cdot \sin(\alpha _{f})=\Delta b\cdot \sin(\theta '-\alpha _{f})

where $R$ and $b$ are defined in the same way as they were for spherical aberration. Assuming that $\alpha _{f}<<\theta '$ (or equivalently that $|b\cdot \sin(\alpha _{f})|<<|R|$ ), we can show that

\sin(\alpha _{f})\approx {\frac {\Delta b\sin(\theta ')}{\sqrt {R^{2}+b^{2}}}}={\frac {\Delta b\cdot R}{R^{2}+b^{2}}}

Since we required $\alpha _{f}$ to be small, and since $\theta$ being small implies $R<<b$ , we are given an approximation of $\alpha _{f}$ as

\alpha _{f}\approx {\frac {\Delta b\cdot R}{b^{2}}}

From the thin-lens formula it can be shown that $\Delta b/b^{2}\approx \Delta f/f^{2}$ , yielding a final estimation of the difference in refracted angle between in-focus and off-focus rays as

\alpha _{f}\approx {\frac {\Delta f\cdot R}{f^{2}}}

Examples

The contrast transfer function determines how much phase signal gets transmitted to the real space wavefunction in the image plane. As the modulus squared of the real space wavefunction gives the image signal, the contrast transfer function limits how much information can ultimately be translated into an image. The form of the contrast transfer function determines the quality of real space image formation in the TEM.

This is an example contrast transfer function. There are a number of things to note:

The function exists in the spatial frequency domain, or k-space
Whenever the function is equal to zero, that means there is no transmittance, or no phase signal is incorporated into the real space image
The first time the function crosses the x-axis is called the point resolution
To maximize phase signal, it is generally better to use imaging conditions that push the point resolution to higher spatial frequencies
When the function is negative, that represents positive phase contrast, leading to a bright background, with dark atomic features
Every time the CTF crosses the x-axis, there is an inversion in contrast
Accordingly, past the point resolution of the microscope the phase information is not directly interpretable, and must be modeled via computer simulation

Scherzer defocus

The defocus value ( ${\textstyle z}$ ) can be used to counteract the spherical aberration to allow for greater phase contrast. This analysis was developed by Scherzer, and is called the Scherzer defocus.^[7]

$z_{s}=(C_{s}\lambda )^{1/2}$

The variables are the same as from the mathematical treatment section, with $z_{s}$ setting the specific Scherzer defocus, $C_{s}$ as the spherical aberration, and λ as the relativistic wavelength for the electron wave.

The figure in the following section shows the CTF function for a CM300 Microscope at the Scherzer Defocus. Compared to the CTF Function showed above, there is a larger window, also known as a passband, of spatial frequencies with high transmittance. This allows more phase signal to pass through to the image plane.

Envelope function

The envelope function represents the effect of additional aberrations that damp the contrast transfer function, and in turn the phase. The envelope terms comprising the envelope function tend to suppress high spatial frequencies. The exact form of the envelope functions can differ from source to source. Generally, they are applied by multiplying the Contrast Transfer Function by an envelope term Et representing temporal aberrations, and an envelope term Es representing spatial aberrations. This yields a modified, or effective Contrast Transfer Function:

$K_{eff}(k)=E_{t}E_{s}(\sin[(2\pi /\lambda )W(k)]$

Examples of temporal aberrations include chromatic aberrations, energy spread, focal spread, instabilities in the high voltage source, and instabilities in the objective lens current. An example of a spatial aberration includes the finite incident beam convergence.^[8]

As shown in the figure, the most restrictive envelope term will dominate in damping the contrast transfer function. In this particular example, the temporal envelope term is the most restrictive. Because the envelope terms damp more strongly at higher spatial frequencies, there comes a point where no more phase signal can pass through. This is called the Information Limit of the microscope, and is one measure of the resolution.

Modeling the envelope function can give insight into both TEM instrument design, and imaging parameters. By modeling the different aberrations via envelope terms, it is possible to see which aberrations are most limiting the phase signal.

Various software have been developed to model both the Contrast Transfer Function and Envelope Function for particular microscopes, and particular imaging parameters.^[9]^[10]

Linear imaging theory vs. non-linear imaging theory

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]

Related Research Articles

$Diffraction Phenomenon of the motion of waves$

Diffraction is the interference or bending of waves around the corners of an obstacle or through an aperture into the region of geometrical shadow of the obstacle/aperture. The diffracting object or aperture effectively becomes a secondary source of the propagating wave. Italian scientist Francesco Maria Grimaldi coined the word diffraction and was the first to record accurate observations of the phenomenon in 1660.

A hydrogen atom is an atom of the chemical element hydrogen. The electrically neutral hydrogen atom contains a nucleus of a single positively charged proton and a single negatively charged electron bound to the nucleus by the Coulomb force. Atomic hydrogen constitutes about 75% of the baryonic mass of the universe.

$Refractive index Property in optics$

In optics, the refractive index of an optical medium is the ratio of the apparent speed of light in the medium to the speed in air or vacuum. The refractive index determines how much the path of light is bent, or refracted, when entering a material. This is described by Snell's law of refraction, $n 1 sin θ 1 = n 2 sin θ 2$ , where $θ 1$ and $θ 2$ are the angle of incidence and angle of refraction, respectively, of a ray crossing the interface between two media with refractive indices $n 1$ and $n 2$ . The refractive indices also determine the amount of light that is reflected when reaching the interface, as well as the critical angle for total internal reflection, their intensity and Brewster's angle.

Ray transfer matrix analysis is a mathematical form for performing ray tracing calculations in sufficiently simple problems which can be solved considering only paraxial rays. Each optical element is described by a 2 × 2ray transfer matrix which operates on a vector describing an incoming light ray to calculate the outgoing ray. Multiplication of the successive matrices thus yields a concise ray transfer matrix describing the entire optical system. The same mathematics is also used in accelerator physics to track particles through the magnet installations of a particle accelerator, see electron optics.

Angular resolution describes the ability of any image-forming device such as an optical or radio telescope, a microscope, a camera, or an eye, to distinguish small details of an object, thereby making it a major determinant of image resolution. It is used in optics applied to light waves, in antenna theory applied to radio waves, and in acoustics applied to sound waves. The colloquial use of the term "resolution" sometimes causes confusion; when an optical system is said to have a high resolution or high angular resolution, it means that the perceived distance, or actual angular distance, between resolved neighboring objects is small. The value that quantifies this property, θ, which is given by the Rayleigh criterion, is low for a system with a high resolution. The closely related term spatial resolution refers to the precision of a measurement with respect to space, which is directly connected to angular resolution in imaging instruments. The Rayleigh criterion shows that the minimum angular spread that can be resolved by an image-forming system is limited by diffraction to the ratio of the wavelength of the waves to the aperture width. For this reason, high-resolution imaging systems such as astronomical telescopes, long distance telephoto camera lenses and radio telescopes have large apertures.

Fourier optics is the study of classical optics using Fourier transforms (FTs), in which the waveform being considered is regarded as made up of a combination, or superposition, of plane waves. It has some parallels to the Huygens–Fresnel principle, in which the wavefront is regarded as being made up of a combination of spherical wavefronts whose sum is the wavefront being studied. A key difference is that Fourier optics considers the plane waves to be natural modes of the propagation medium, as opposed to Huygens–Fresnel, where the spherical waves originate in the physical medium.

In mathematics, the inverse trigonometric functions are the inverse functions of the trigonometric functions, under suitably restricted domains. Specifically, they are the inverses of the sine, cosine, tangent, cotangent, secant, and cosecant functions, and are used to obtain an angle from any of the angle's trigonometric ratios. Inverse trigonometric functions are widely used in engineering, navigation, physics, and geometry.

$Airy disk Diffraction pattern in optics$

In optics, the Airy disk and Airy pattern are descriptions of the best-focused spot of light that a perfect lens with a circular aperture can make, limited by the diffraction of light. The Airy disk is of importance in physics, optics, and astronomy.

In physics, a wave vector is a vector used in describing a wave, with a typical unit being cycle per metre. It has a magnitude and direction. Its magnitude is the wavenumber of the wave, and its direction is perpendicular to the wavefront. In isotropic media, this is also the direction of wave propagation.

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when plane waves are incident on a diffracting object, and the diffraction pattern is viewed at a sufficiently long distance from the object, and also when it is viewed at the focal plane of an imaging lens. In contrast, the diffraction pattern created near the diffracting object and is given by the Fresnel diffraction equation.

<span class="mw-page-title-main">Envelope (mathematics)</span> Curve external to a family of curves in geometry

In geometry, an envelope of a planar family of curves is a curve that is tangent to each member of the family at some point, and these points of tangency together form the whole envelope. Classically, a point on the envelope can be thought of as the intersection of two "infinitesimally adjacent" curves, meaning the limit of intersections of nearby curves. This idea can be generalized to an envelope of surfaces in space, and so on to higher dimensions.

In quantum physics, the scattering amplitude is the probability amplitude of the outgoing spherical wave relative to the incoming plane wave in a stationary-state scattering process. At large distances from the centrally symmetric scattering center, the plane wave is described by the wavefunction

Optical resolution describes the ability of an imaging system to resolve detail, in the object that is being imaged. An imaging system may have many individual components, including one or more lenses, and/or recording and display components. Each of these contributes to the optical resolution of the system; the environment in which the imaging is done often is a further important factor.

The optical transfer function (OTF) of an optical system such as a camera, microscope, human eye, or projector specifies how different spatial frequencies are captured or transmitted. It is used by optical engineers to describe how the optics project light from the object or scene onto a photographic film, detector array, retina, screen, or simply the next item in the optical transmission chain. A variant, the modulation transfer function (MTF), neglects phase effects, but is equivalent to the OTF in many situations.

<span class="mw-page-title-main">High-resolution transmission electron microscopy</span>

High-resolution transmission electron microscopy is an imaging mode of specialized transmission electron microscopes that allows for direct imaging of the atomic structure of samples. It is a powerful tool to study properties of materials on the atomic scale, such as semiconductors, metals, nanoparticles and sp²-bonded carbon. While this term is often also used to refer to high resolution scanning transmission electron microscopy, mostly in high angle annular dark field mode, this article describes mainly the imaging of an object by recording the two-dimensional spatial wave amplitude distribution in the image plane, similar to a "classic" light microscope. For disambiguation, the technique is also often referred to as phase contrast transmission electron microscopy, although this term is less appropriate. At present, the highest point resolution realised in high resolution transmission electron microscopy is around 0.5 ångströms (0.050 nm). At these small scales, individual atoms of a crystal and defects can be resolved. For 3-dimensional crystals, it is necessary to combine several views, taken from different angles, into a 3D map. This technique is called electron tomography.

In general relativity, a point mass deflects a light ray with impact parameter $by an angle approximately equal to$

<span class="mw-page-title-main">Wrapped exponential distribution</span> Probability distribution

In probability theory and directional statistics, a wrapped exponential distribution is a wrapped probability distribution that results from the "wrapping" of the exponential distribution around the unit circle.

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when the diffraction pattern is viewed at a long distance from the diffracting object, and also when it is viewed at the focal plane of an imaging lens.

In physics and engineering, the radiative heat transfer from one surface to another is the equal to the difference of incoming and outgoing radiation from the first surface. In general, the heat transfer between surfaces is governed by temperature, surface emissivity properties and the geometry of the surfaces. The relation for heat transfer can be written as an integral equation with boundary conditions based upon surface conditions. Kernel functions can be useful in approximating and solving this integral equation.

References

1 2 3 Wade, 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.
↑ Spence, John C. H. (1988 2nd ed) Experimental high-resolution electron microscopy (Oxford U. Press, NY) ISBN 0195054059.
↑ Ludwig Reimer (1997 4th ed) Transmission electron microscopy: Physics of image formation and microanalysis (Springer, Berlin) preview.
↑ Earl J. Kirkland (1998) Advanced computing in electron microscopy (Plenum Press, NY).
↑ "DeBroglie Wavelength". HyperPhysics. Georgia State University. Retrieved 27 April 2017.
↑ "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.
↑ 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.
↑ "Envelope Functions". www.maxsidorov.com. Retrieved 2015-06-12.
↑ "CTF Simulation". Wen Jiang Group. Retrieved 27 April 2017.
↑ Sidorov, Max. "Home of the ctfExplorer" . Retrieved 27 April 2017.
↑ 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.
↑ This page was prepared in part for Northwestern University class MSE 465, taught by Professor Laurie Marks.
↑ Notes prepared by Professor Laurie Marks at Northwestern University.

External links

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[:0-1] 1 2 3 Wade, 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.

[Spence1982-2] Spence, John C. H. (1988 2nd ed) Experimental high-resolution electron microscopy (Oxford U. Press, NY) ISBN 0195054059.

[Reimer97-3] Ludwig Reimer (1997 4th ed) Transmission electron microscopy: Physics of image formation and microanalysis (Springer, Berlin) preview.

[Kirkland1998-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". Wen Jiang Group. Retrieved 27 April 2017.

[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.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]