Transport-of-intensity equation

Last updated

The transport-of-intensity equation (TIE) is a computational approach to reconstruct the phase of a complex wave in optical and electron microscopy. [1] It describes the internal relationship between the intensity and phase distribution of a wave. [2]

The TIE was first proposed in 1983 by Michael Reed Teague. [3] Teague suggested to use the law of conservation of energy to write a differential equation for the transport of energy by an optical field. This equation, he stated, could be used as an approach to phase recovery. [4]

Teague approximated the amplitude of the wave propagating nominally in the z-direction by a parabolic equation and then expressed it in terms of irradiance and phase:

where is the wavelength, is the irradiance at point , and is the phase of the wave. If the intensity distribution of the wave and its spatial derivative can be measured experimentally, the equation becomes a linear equation that can be solved to obtain the phase distribution . [5]

For a phase sample with a constant intensity, the TIE simplifies to

It allows measuring the phase distribution of the sample by acquiring a defocused image, i.e. .

TIE-based approaches are applied in biomedical and technical applications, such as quantitative monitoring of cell growth in culture, [6] investigation of cellular dynamics and characterization of optical elements. [7] The TIE method  is also applied for phase retrieval in transmission electron microscopy. [8]

Related Research Articles

<span class="mw-page-title-main">Diffraction</span> Phenomenon of the motion of waves

Diffraction is defined as 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.

<span class="mw-page-title-main">Refractive index</span> Ratio of the speed of light in vacuum to that in the medium

In optics, the refractive index of an optical medium is a dimensionless number that gives the indication of the light bending ability of that medium.

<span class="mw-page-title-main">Gaussian beam</span> Monochrome light beam whose amplitude envelope is a Gaussian function

In optics, a Gaussian beam is a beam of electromagnetic radiation with high monochromaticity whose amplitude envelope in the transverse plane is given by a Gaussian function; this also implies a Gaussian intensity (irradiance) profile. This fundamental (or TEM00) transverse Gaussian mode describes the intended output of most (but not all) lasers, as such a beam can be focused into the most concentrated spot. When such a beam is refocused by a lens, the transverse phase dependence is altered; this results in a different Gaussian beam. The electric and magnetic field amplitude profiles along any such circular Gaussian beam (for a given wavelength and polarization) are determined by a single parameter: the so-called waist w0. At any position z relative to the waist (focus) along a beam having a specified w0, the field amplitudes and phases are thereby determined as detailed below.

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×2 ray 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.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

The primitive equations are a set of nonlinear partial differential equations that are used to approximate global atmospheric flow and are used in most atmospheric models. They consist of three main sets of balance equations:

  1. A continuity equation: Representing the conservation of mass.
  2. Conservation of momentum: Consisting of a form of the Navier–Stokes equations that describe hydrodynamical flow on the surface of a sphere under the assumption that vertical motion is much smaller than horizontal motion (hydrostasis) and that the fluid layer depth is small compared to the radius of the sphere
  3. A thermal energy equation: Relating the overall temperature of the system to heat sources and sinks

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.

Geometrical optics, or ray optics, is a model of optics that describes light propagation in terms of rays. The ray in geometrical optics is an abstraction useful for approximating the paths along which light propagates under certain circumstances.

In optics, an ultrashort pulse, also known as an ultrafast event, is an electromagnetic pulse whose time duration is of the order of a picosecond or less. Such pulses have a broadband optical spectrum, and can be created by mode-locked oscillators. Amplification of ultrashort pulses almost always requires the technique of chirped pulse amplification, in order to avoid damage to the gain medium of the amplifier.

<span class="mw-page-title-main">Nonlinear Schrödinger equation</span> Nonlinear form of the Schrödinger equation

In theoretical physics, the (one-dimensional) nonlinear Schrödinger equation (NLSE) is a nonlinear variation of the Schrödinger equation. It is a classical field equation whose principal applications are to the propagation of light in nonlinear optical fibers and planar waveguides and to Bose–Einstein condensates confined to highly anisotropic, cigar-shaped traps, in the mean-field regime. Additionally, the equation appears in the studies of small-amplitude gravity waves on the surface of deep inviscid (zero-viscosity) water; the Langmuir waves in hot plasmas; the propagation of plane-diffracted wave beams in the focusing regions of the ionosphere; the propagation of Davydov's alpha-helix solitons, which are responsible for energy transport along molecular chains; and many others. More generally, the NLSE appears as one of universal equations that describe the evolution of slowly varying packets of quasi-monochromatic waves in weakly nonlinear media that have dispersion. Unlike the linear Schrödinger equation, the NLSE never describes the time evolution of a quantum state. The 1D NLSE is an example of an integrable model.

Self-phase modulation (SPM) is a nonlinear optical effect of light–matter interaction. An ultrashort pulse of light, when travelling in a medium, will induce a varying refractive index of the medium due to the optical Kerr effect. This variation in refractive index will produce a phase shift in the pulse, leading to a change of the pulse's frequency spectrum.

A Van Hove singularity is a singularity in the density of states (DOS) of a crystalline solid. The wavevectors at which Van Hove singularities occur are often referred to as critical points of the Brillouin zone. For three-dimensional crystals, they take the form of kinks. The most common application of the Van Hove singularity concept comes in the analysis of optical absorption spectra. The occurrence of such singularities was first analyzed by the Belgian physicist Léon Van Hove in 1953 for the case of phonon densities of states.

<span class="mw-page-title-main">Directivity</span> Measure of how much of an antennas signal is transmitted in one direction

In electromagnetics, directivity is a parameter of an antenna or optical system which measures the degree to which the radiation emitted is concentrated in a single direction. It is the ratio of the radiation intensity in a given direction from the antenna to the radiation intensity averaged over all directions. Therefore, the directivity of a hypothetical isotropic radiator is 1, or 0 dBi.

Ellipsoidal coordinates are a three-dimensional orthogonal coordinate system that generalizes the two-dimensional elliptic coordinate system. Unlike most three-dimensional orthogonal coordinate systems that feature quadratic coordinate surfaces, the ellipsoidal coordinate system is based on confocal quadrics.

Fluorescence interference contrast (FLIC) microscopy is a microscopic technique developed to achieve z-resolution on the nanometer scale.

<span class="mw-page-title-main">Stokes wave</span> Nonlinear and periodic surface wave on an inviscid fluid layer of constant mean depth

In fluid dynamics, a Stokes wave is a nonlinear and periodic surface wave on an inviscid fluid layer of constant mean depth. This type of modelling has its origins in the mid 19th century when Sir George Stokes – using a perturbation series approach, now known as the Stokes expansion – obtained approximate solutions for nonlinear wave motion.

<span class="mw-page-title-main">Mild-slope equation</span> Physics phenomenon and formula

In fluid dynamics, the mild-slope equation describes the combined effects of diffraction and refraction for water waves propagating over bathymetry and due to lateral boundaries—like breakwaters and coastlines. It is an approximate model, deriving its name from being originally developed for wave propagation over mild slopes of the sea floor. The mild-slope equation is often used in coastal engineering to compute the wave-field changes near harbours and coasts.

<span class="mw-page-title-main">Multiple-prism dispersion theory</span> Theory in optics

The first description of multiple-prism arrays, and multiple-prism dispersion, was given by Newton in his book Opticks. Prism pair expanders were introduced by Brewster in 1813. A modern mathematical description of the single-prism dispersion was given by Born and Wolf in 1959. The generalized multiple-prism dispersion theory was introduced by Duarte and Piper in 1982.

Quantum mechanics was first applied to optics, and interference in particular, by Paul Dirac. Richard Feynman, in his Lectures on Physics, uses Dirac's notation to describe thought experiments on double-slit interference of electrons. Feynman's approach was extended to N-slit interferometers for either single-photon illumination, or narrow-linewidth laser illumination, that is, illumination by indistinguishable photons, by Frank Duarte. The N-slit interferometer was first applied in the generation and measurement of complex interference patterns.

In mathematics, the method of steepest descent or saddle-point method is an extension of Laplace's method for approximating an integral, where one deforms a contour integral in the complex plane to pass near a stationary point, in roughly the direction of steepest descent or stationary phase. The saddle-point approximation is used with integrals in the complex plane, whereas Laplace’s method is used with real integrals.

References

  1. Bostan, E. (2014). "Phase Retrieval by Using Transport-of-Intensity Equation and Differential Interference Contrast Microscopy" (PDF). IEEE International Conference on Image Processing (ICIP): 3939–3943. doi:10.1109/ICIP.2014.7025800. ISBN   978-1-4799-5751-4. S2CID   10310598.
  2. Cheng, H. (2009). "Phase Retrieval Using the Transport-of-Intensity Equation". IEEE Fifth International Conference on Image and Graphics: 417–421. doi:10.1109/ICIG.2009.32. ISBN   978-1-4244-5237-8. S2CID   15772496.
  3. Teague, Michael R. (1983). "Deterministic phase retrieval: a Green's function solution". Journal of the Optical Society of America. 73 (11): 1434–1441. doi:10.1364/JOSA.73.001434.
  4. Nugent, Keith (2010). "Coherent methods in the X-ray sciences". Advances in Physics. 59 (1): 1–99. arXiv: 0908.3064 . Bibcode:2010AdPhy..59....1N. doi:10.1080/00018730903270926. S2CID   118519311.
  5. Gureyev, T. E.; Roberts, A.; Nugent, K. A. (1995). "Partially coherent fields, the transport-of-intensity equation, and phase uniqueness". JOSA A. 12 (9): 1942–1946. Bibcode:1995JOSAA..12.1942G. doi:10.1364/JOSAA.12.001942.
  6. Curl, C.L. (2004). "Quantitative phase microscopy: a new tool for measurement of cell culture growth and confluency in situ". Pflügers Archiv: European Journal of Physiology. 448 (4): 462–468. doi:10.1007/s00424-004-1248-7. PMID   14985984. S2CID   7640406.
  7. Dorrer, C. (2007). "Optical testing using the transport-of-intensity equation". Opt. Express. 15 (12): 7165–7175. Bibcode:2007OExpr..15.7165D. doi: 10.1364/oe.15.007165 . PMID   19547035.
  8. Belaggia, M. (2004). "On the transport of intensity technique for phase retrieval". Ultramicroscopy. 102 (1): 37–49. doi:10.1016/j.ultramic.2004.08.004. PMID   15556699.