We compute the eigenvalues of the matrix T {\displaystyle \mathbf {T} } that satisfy eigenequation [ T − λ I ] v = [ A − λ B C D − λ ] v = 0 , {\displaystyle [{\boldsymbol {T}}-\lambda I]\mathbf {v} ={\begin{bmatrix}A-\lambda &B\\C&D-\lambda \end{bmatrix}}\mathbf {v} =0,} by calculating the determinant | A − λ B C D − λ | = λ 2 − ( A + D ) λ + 1 = 0. {\displaystyle {\begin{vmatrix}A-\lambda &B\\C&D-\lambda \end{vmatrix}}=\lambda ^{2}-(A+D)\lambda +1=0.}

As a result, if the input and output planes are located within the same medium, or within two different media which happen to have identical indices of refraction, then the determinant of M is simply equal to 1.

We proceed to calculate the eigenvalues of the transfer matrix: det [ M − λ I ] = 0 , {\displaystyle \det \left[\mathbf {M} -\lambda \mathbf {I} \right]=0,} leading to the characteristic equation λ 2 − tr ⁡ ( M ) λ + det ( M ) = 0 , {\displaystyle \lambda ^{2}-\operatorname {tr} (\mathbf {M} )\lambda +\det(\mathbf {M} )=0,} where tr ⁡ ( M ) = A + D = 2 − d f {\displaystyle \operatorname {tr} (\mathbf {M} )=A+D=2-{\frac {d}{f}}} is the trace of the RTM, and det ( M ) = A D − B C = 1 {\displaystyle \det(\mathbf {M} )=AD-BC=1} is the determinant of the RTM. After one common substitution we have: λ 2 − 2 g λ + 1 = 0 , {\displaystyle \lambda ^{2}-2g\lambda +1=0,} where g = d e f tr ⁡ ( M ) 2 = 1 − d 2 f {\displaystyle g{\overset {\mathrm {def} }{{}={}}}{\frac {\operatorname {tr} (\mathbf {M} )}{2}}=1-{\frac {d}{2f}}} is the stability parameter. The eigenvalues are the solutions of the characteristic equation. From the quadratic formula we find λ ± = g ± g 2 − 1 . {\displaystyle \lambda _{\pm }=g\pm {\sqrt {g^{2}-1}}.}

Paraxial raysin physics

RTM analysis can now be used to determine the stability of the waveguide (and equivalently, the resonator). That is, it can be determined under what conditions light traveling down the waveguide will be periodically refocused and stay within the waveguide. To do so, we can find all the "eigenrays" of the system: the input ray vector at each of the mentioned sections of the waveguide times a real or complex factor λ is equal to the output one. This gives: M [ x 1 θ 1 ] = [ x 2 θ 2 ] = λ [ x 1 θ 1 ] . {\displaystyle \mathbf {M} {\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}}={\begin{bmatrix}x_{2}\\\theta _{2}\end{bmatrix}}=\lambda {\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}}.} which is an eigenvalue equation: [ M − λ I ] [ x 1 θ 1 ] = 0 , {\displaystyle \left[\mathbf {M} -\lambda \mathbf {I} \right]{\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}}=0,} where I = [ 1 0 0 1 ] {\textstyle \mathbf {I} =\left[{\begin{smallmatrix}1&0\\0&1\end{smallmatrix}}\right]} is the 2 × 2 identity matrix.

Paraxialray tracing

For g 2 < 1 {\displaystyle g^{2}<1} let r + {\displaystyle r_{+}} and r − {\displaystyle r_{-}} be the eigenvectors with respect to the eigenvalues λ + {\displaystyle \lambda _{+}} and λ − {\displaystyle \lambda _{-}} respectively, which span all the vector space because they are orthogonal, the latter due to λ + ≠ λ − {\displaystyle \lambda _{+}\neq \lambda _{-}} . The input vector can therefore be written as c + r + + c − r − , {\displaystyle c_{+}r_{+}+c_{-}r_{-},} for some constants c + {\displaystyle c_{+}} and c − {\displaystyle c_{-}} .

Note that, since the multiplication of matrices is non-commutative, this is not the same RTM as that for a lens followed by free space: S L = [ 1 d 0 1 ] [ 1 0 − 1 f 1 ] = [ 1 − d f d − 1 f 1 ] . {\displaystyle \mathbf {SL} ={\begin{bmatrix}1&d\\0&1\end{bmatrix}}{\begin{bmatrix}1&0\\-{\frac {1}{f}}&1\end{bmatrix}}={\begin{bmatrix}1-{\frac {d}{f}}&d\\-{\frac {1}{f}}&1\end{bmatrix}}.}

Paraxial raysmirror

According to the values of λ 1 {\displaystyle \lambda _{1}} and λ 2 {\displaystyle \lambda _{2}} , there are several possible cases. For example:

Nonparaxial rays

It is often convenient to express this last equation in reciprocal form: 1 q 2 = C + D / q 1 A + B / q 1 . {\displaystyle {\frac {1}{q_{2}}}={\frac {C+D/q_{1}}{A+B/q_{1}}}.}

Paraxialapproximation

Methods using transfer matrices of higher dimensionality, that is 3 × 3, 4 × 4, and 6 × 6, are also used in optical analysis.[9] In particular, 4 × 4 propagation matrices are used in the design and analysis of prism sequences for pulse compression in femtosecond lasers.[5]

The use of transfer matrices in this manner parallels the 2 × 2 matrices describing electronic two-port networks, particularly various so-called ABCD matrices which can similarly be multiplied to solve for cascaded systems.

The ABCD matrix representing a component or system relates the output ray to the input according to [ x 2 θ 2 ] = [ A B C D ] [ x 1 θ 1 ] , {\displaystyle {\begin{bmatrix}x_{2}\\\theta _{2}\end{bmatrix}}={\begin{bmatrix}A&B\\C&D\end{bmatrix}}{\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}},} where the values of the 4 matrix elements are thus given by A = x 2 x 1 | θ 1 = 0 B = x 2 θ 1 | x 1 = 0 , {\displaystyle A=\left.{\frac {x_{2}}{x_{1}}}\right|_{\theta _{1}=0}\qquad B=\left.{\frac {x_{2}}{\theta _{1}}}\right|_{x_{1}=0},} and C = θ 2 x 1 | θ 1 = 0 D = θ 2 θ 1 | x 1 = 0 . {\displaystyle C=\left.{\frac {\theta _{2}}{x_{1}}}\right|_{\theta _{1}=0}\qquad D=\left.{\frac {\theta _{2}}{\theta _{1}}}\right|_{x_{1}=0}.}

Paraxial rayspdf

Consider a beam traveling a distance d through free space, the ray transfer matrix is [ A B C D ] = [ 1 d 0 1 ] . {\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}={\begin{bmatrix}1&d\\0&1\end{bmatrix}}.} and so q 2 = A q 1 + B C q 1 + D = q 1 + d 1 = q 1 + d {\displaystyle q_{2}={\frac {Aq_{1}+B}{Cq_{1}+D}}={\frac {q_{1}+d}{1}}=q_{1}+d} consistent with the expression above for ordinary Gaussian beam propagation, i.e. q = ( z − z 0 ) + i z R {\displaystyle q=(z-z_{0})+iz_{R}} . As the beam propagates, both the radius and waist change.

A different convention for the ray vectors can be employed. Instead of using θ ≈ sin θ, the second element of the ray vector is n sin θ,[2] which is proportional not to the ray angle per se but to the transverse component of the wave vector. This alters the ABCD matrices given in the table below where refraction at an interface is involved.

Yes, there are video tutorials available that demonstrate how to use the measurement tools in WebViewer. One useful resource is the YouTube video titled How to Use Measurement Tools in WebViewer, which provides a comprehensive overview.

There exist infinite ways to decompose a ray transfer matrix T = [ A B C D ] {\displaystyle \mathbf {T} ={\begin{bmatrix}A&B\\C&D\end{bmatrix}}} into a concatenation of multiple transfer matrices. For example in the special case when n 1 = n 2 {\displaystyle n_{1}=n_{2}} :

Consider a beam traveling through a thin lens with focal length f. The ray transfer matrix is [ A B C D ] = [ 1 0 − 1 / f 1 ] . {\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}={\begin{bmatrix}1&0\\-1/f&1\end{bmatrix}}.} and so q 2 = A q 1 + B C q 1 + D = q 1 − q 1 f + 1 {\displaystyle q_{2}={\frac {Aq_{1}+B}{Cq_{1}+D}}={\frac {q_{1}}{-{\frac {q_{1}}{f}}+1}}} 1 q 2 = − q 1 f + 1 q 1 = 1 q 1 − 1 f . {\displaystyle {\frac {1}{q_{2}}}={\frac {-{\frac {q_{1}}{f}}+1}{q_{1}}}={\frac {1}{q_{1}}}-{\frac {1}{f}}.} Only the real part of 1/q is affected: the wavefront curvature 1/R is reduced by the power of the lens 1/f, while the lateral beam size w remains unchanged upon exiting the thin lens.

This beam can be propagated through an optical system with a given ray transfer matrix by using the equation[further explanation needed]: [ q 2 1 ] = k [ A B C D ] [ q 1 1 ] , {\displaystyle {\begin{bmatrix}q_{2}\\1\end{bmatrix}}=k{\begin{bmatrix}A&B\\C&D\end{bmatrix}}{\begin{bmatrix}q_{1}\\1\end{bmatrix}},} where k is a normalization constant chosen to keep the second component of the ray vector equal to 1. Using matrix multiplication, this equation expands as q 2 = k ( A q 1 + B ) 1 = k ( C q 1 + D ) . {\displaystyle {\begin{aligned}q_{2}&=k(Aq_{1}+B)\\1&=k(Cq_{1}+D)\,.\end{aligned}}}

(R, w, and q are functions of position.) If the beam axis is in the z direction, with waist at z0 and Rayleigh range zR, this can be equivalently written as[8] q = ( z − z 0 ) + i z R . {\displaystyle q=(z-z_{0})+iz_{R}.}

This technique, as described below, is derived using the paraxial approximation, which requires that all ray directions (directions normal to the wavefronts) are at small angles θ relative to the optical axis of the system, such that the approximation sin θ ≈ θ remains valid. A small θ further implies that the transverse extent of the ray bundles (x and y) is small compared to the length of the optical system (thus "paraxial"). Since a decent imaging system where this is not the case for all rays must still focus the paraxial rays correctly, this matrix method will properly describe the positions of focal planes and magnifications, however aberrations still need to be evaluated using full ray-tracing techniques.[1]

The same matrices can also be used to calculate the evolution of Gaussian beams[7] propagating through optical components described by the same transmission matrices. If we have a Gaussian beam of wavelength λ 0 {\displaystyle \lambda _{0}} , radius of curvature R (positive for diverging, negative for converging), beam spot size w and refractive index n, it is possible to define a complex beam parameter q by:[8] 1 q = 1 R − i λ 0 π n w 2 . {\displaystyle {\frac {1}{q}}={\frac {1}{R}}-{\frac {i\lambda _{0}}{\pi nw^{2}}}.}

Creating measurement annotations in a web application is straightforward with Apryse. You can add various types of measurement annotations such as distance, perimeter, and area. For a detailed guide, check out the Measurement Annotations Documentation.

A ray transfer matrix can be regarded as a linear canonical transformation. According to the eigenvalues of the optical system, the system can be classified into several classes.[3] Assume the ABCD matrix representing a system relates the output ray to the input according to [ x 2 θ 2 ] = [ A B C D ] [ x 1 θ 1 ] = T v . {\displaystyle {\begin{bmatrix}x_{2}\\\theta _{2}\end{bmatrix}}={\begin{bmatrix}A&B\\C&D\end{bmatrix}}{\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}}=\mathbf {T} \mathbf {v} .}

This relates the ray vectors at the input and output planes by the ray transfer matrix (RTM) M, which represents the optical component or system present between the two reference planes. A thermodynamics argument based on the blackbody radiation [citation needed] can be used to show that the determinant of a RTM is the ratio of the indices of refraction: det ( M ) = A D − B C = n 1 n 2 . {\displaystyle \det(\mathbf {M} )=AD-BC={\frac {n_{1}}{n_{2}}}.}

If the waveguide is stable, no ray should stray arbitrarily far from the main axis, that is, λN must not grow without limit. Suppose g 2 > 1 {\displaystyle g^{2}>1} . Then both eigenvalues are real. Since λ + λ − = 1 {\displaystyle \lambda _{+}\lambda _{-}=1} , one of them has to be bigger than 1 (in absolute value), which implies that the ray which corresponds to this eigenvector would not converge. Therefore, in a stable waveguide, g 2 ≤ 1 {\displaystyle g^{2}\leq 1} , and the eigenvalues can be represented by complex numbers: λ ± = g ± i 1 − g 2 = cos ⁡ ( ϕ ) ± i sin ⁡ ( ϕ ) = e ± i ϕ , {\displaystyle \lambda _{\pm }=g\pm i{\sqrt {1-g^{2}}}=\cos(\phi )\pm i\sin(\phi )=e^{\pm i\phi },} with the substitution g = cos(ϕ).

Paraxial raysexamples

Now, consider a ray after N passes through the system: [ x N θ N ] = λ N [ x 1 θ 1 ] . {\displaystyle {\begin{bmatrix}x_{N}\\\theta _{N}\end{bmatrix}}=\lambda ^{N}{\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}}.}

Another simple example is that of a thin lens. Its RTM is given by: L = [ 1 0 − 1 f 1 ] , {\displaystyle \mathbf {L} ={\begin{bmatrix}1&0\\-{\frac {1}{f}}&1\end{bmatrix}},} where f is the focal length of the lens. To describe combinations of optical components, ray transfer matrices may be multiplied together to obtain an overall RTM for the compound optical system. For the example of free space of length d followed by a lens of focal length f: L S = [ 1 0 − 1 f 1 ] [ 1 d 0 1 ] = [ 1 d − 1 f 1 − d f ] . {\displaystyle \mathbf {L} \mathbf {S} ={\begin{bmatrix}1&0\\-{\frac {1}{f}}&1\end{bmatrix}}{\begin{bmatrix}1&d\\0&1\end{bmatrix}}={\begin{bmatrix}1&d\\-{\frac {1}{f}}&1-{\frac {d}{f}}\end{bmatrix}}.}

R e = R / cos ⁡ θ {\displaystyle R_{e}=R/\cos \theta } effective radius of curvature in the sagittal plane (vertical direction) R = radius of curvature, R > 0 for concave, valid in the paraxial approximation θ is the mirror angle of incidence in the horizontal plane.

Dividing the first equation by the second eliminates the normalization constant: q 2 = A q 1 + B C q 1 + D , {\displaystyle q_{2}={\frac {Aq_{1}+B}{Cq_{1}+D}},}

You can explore a live demo of the measurement tools and see how they work in real-time by visiting the Measurement Tools Showcase.

n2 = refractive index of the lens itself (inside the lens). R1 = Radius of curvature of First surface. R2 = Radius of curvature of Second surface. t = center thickness of lens.

Paraxial raysformula

The ray tracing technique is based on two reference planes, called the input and output planes, each perpendicular to the optical axis of the system. At any point along the optical train an optical axis is defined corresponding to a central ray; that central ray is propagated to define the optical axis further in the optical train which need not be in the same physical direction (such as when bent by a prism or mirror). The transverse directions x and y (below we only consider the x direction) are then defined to be orthogonal to the optical axes applying. A light ray enters a component crossing its input plane at a distance x1 from the optical axis, traveling in a direction that makes an angle θ1 with the optical axis. After propagation to the output plane that ray is found at a distance x2 from the optical axis and at an angle θ2 with respect to it. n1 and n2 are the indices of refraction of the media in the input and output plane, respectively.

The theory of Linear canonical transformation implies the relation between ray transfer matrix (geometrical optics) and wave optics.[6]

The measurement count tool in PDFTron for Android allows you to count and measure objects within a PDF document. This feature is especially useful for tasks such as quantity takeoffs in construction and engineering projects. For more details, you can read the blog post on Measurement Count Tool and More with PDFTron Android 9.2.3

Let m = ( A + D ) 2 {\displaystyle m={\frac {(A+D)}{2}}} , and we have eigenvalues λ 1 , λ 2 = m ± m 2 − 1 {\displaystyle \lambda _{1},\lambda _{2}=m\pm {\sqrt {m^{2}-1}}} .

As one example, if there is free space between the two planes, the ray transfer matrix is given by: S = [ 1 d 0 1 ] , {\displaystyle \mathbf {S} ={\begin{bmatrix}1&d\\0&1\end{bmatrix}},} where d is the separation distance (measured along the optical axis) between the two reference planes. The ray transfer equation thus becomes: [ x 2 θ 2 ] = S [ x 1 θ 1 ] , {\displaystyle {\begin{bmatrix}x_{2}\\\theta _{2}\end{bmatrix}}=\mathbf {S} {\begin{bmatrix}x_{1}\\\theta _{1}\end{bmatrix}},} and this relates the parameters of the two rays as: x 2 = x 1 + d θ 1 θ 2 = x 1 + d θ 1 {\displaystyle {\begin{aligned}x_{2}&=x_{1}+d\theta _{1}\\\theta _{2}&={\hphantom {x_{1}+d}}\theta _{1}\end{aligned}}}

RTM analysis is particularly useful when modeling the behavior of light in optical resonators, such as those used in lasers. At its simplest, an optical resonator consists of two identical facing mirrors of 100% reflectivity and radius of curvature R, separated by some distance d. For the purposes of ray tracing, this is equivalent to a series of identical thin lenses of focal length f = R/2, each separated from the next by length d. This construction is known as a lens equivalent duct or lens equivalent waveguide. The RTM of each section of the waveguide is, as above, M = L S = ( 1 d − 1 f 1 − d f ) . {\displaystyle \mathbf {M} =\mathbf {L} \mathbf {S} ={\begin{pmatrix}1&d\\{\frac {-1}{f}}&1-{\frac {d}{f}}\end{pmatrix}}.}

After N waveguide sectors, the output reads M N ( c + r + + c − r − ) = λ + N c + r + + λ − N c − r − = e i N ϕ c + r + + e − i N ϕ c − r − , {\displaystyle \mathbf {M} ^{N}(c_{+}r_{+}+c_{-}r_{-})=\lambda _{+}^{N}c_{+}r_{+}+\lambda _{-}^{N}c_{-}r_{-}=e^{iN\phi }c_{+}r_{+}+e^{-iN\phi }c_{-}r_{-},} which represents a periodic function.

Thus the matrices must be ordered appropriately, with the last matrix premultiplying the second last, and so on until the first matrix is premultiplied by the second. Other matrices can be constructed to represent interfaces with media of different refractive indices, reflection from mirrors, etc.

Ray transfer matrix analysis (also known as ABCD 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 (surface, interface, mirror, or beam travel) 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.