Mathematical Fundamentals of Modern Linear Optics

There is a fruitful analogy between mechanics and optics. For example, geometrical optics and classical mechanics can be expressed in the same mathematical language of Hamilton’s theory [1]. Moreover, this analogy played a substantial role in the development of wave and statistical optics as a theory that has the same relations to geometric optics and radiometry as quantum mechanics has to classical mechanics. The advent of quantum mechanics radically changed the classical way of physical thinking [2]. Physics that incorporates elements of quantum mechanics is said to be modern physics [3]. Thus, classical optics based on mathematical fundamentals of quantum and statistical mechanics can be called “modern linear optics.” The adjective “linear” emphasizes its difference from modern optics that embraces quantum optical phenomena. Although the traditional wave optics is based on Maxwell’s electrodynamics [4], in the literature there are four approaches to wave optics which are based on the analogy with quantum mechanics.


Introduction
There is a fruitful analogy between mechanics and optics.For example, geometrical optics and classical mechanics can be expressed in the same mathematical language of Hamilton's theory [1].Moreover, this analogy played a substantial role in the development of wave and statistical optics as a theory that has the same relations to geometric optics and radiometry as quantum mechanics has to classical mechanics.The advent of quantum mechanics radically changed the classical way of physical thinking [2].Physics that incorporates elements of quantum mechanics is said to be modern physics [3].Thus, classical optics based on mathematical fundamentals of quantum and statistical mechanics can be called "modern linear optics."The adjective "linear" emphasizes its difference from modern optics that embraces quantum optical phenomena.
Although the traditional wave optics is based on Maxwell's electrodynamics [4], in the literature there are four approaches to wave optics which are based on the analogy with quantum mechanics.
(1) Generalizing the Huygens-Fresnel principle [5], Feynman has created a new mathematical tool: the path integral [6] which can be used as a fundamental principle of wave and statistical optics [7].In the stationary-phase approximation Feynman's integral over all paths transforms into an integral over one path which is satisfied by Fermat's principle [7], that is, the fundamental principle of Hamilton's optics.
(2) It is known that the Legendre transformations are a framework of Hamilton's optics [1,4,8], and the Fourier transformations are a framework of wave optics [9][10][11][12][13][14].Using the stationary phase approximation Walther proved [12][13][14] that these transformations which before were considered independent, are connected to each other: if the complex wave functions are connected by the Fourier transformations, then the phases of these functions are connected by the Legendre transformation.
(3) It is known that geometrical optics is only a shortwave approximation of wave optics and cannot describe optical phenomena as completely as wave optics does, because all equations of geometrical optics can be derived from the wave equation by neglecting certain terms.But Gloge and Marcuse [15,16] proved that the wave equation can indeed be recovered from the Hamilton equation of geometrical optics by quantization, that is, by formal replacing of all of the variables by corresponding operators.
(4) In 1968 Walther proved [17] that the two divisions of optics, which were considered before independently-the theory of the partial spatial coherence (statistical optics) and radiometry of non-Lambertian sources-have in fact the general subject for study, International Journal of Antennas and Propagation namely, the radiation energy transfer.The basic laws of radiometry can be derived from the theory of partial coherence by using the Wigner definition function (WDF) [17,18].
Note that according to Bohr's correspondence principle [19] modern optics in short-wave approximation must reproduce Hamilton's optics and radiometry.
The purpose of the present paper is to unite these four approaches to wave and statistical optics in a single theory, modern optics.In the short-wave approximation the Huygens-Fresnel principle (in form of Feynman's path integral), the Fourier transformations, and concepts (e.g., WDF) of this modern optics must reproduce Fermat's principle, Legendre transformations, and concepts (e.g., radiance) of Hamilton's optics [8] and radiometry [20] in a one-to-one fashion.

Feynman's Path Integral Reproduces Fermat's Principle
It is known that all results of geometrical optics can be derived from only one principle-the Fermat principle [1,4,8], and all results of wave optics (in the scalar quasi-monochromatic approximation) can be derived from the another principle-the Huygens-Fresnel principle, which can be mathematically expressed in the form of Feynman's path integral [5].Let us show that the Fermat principle can be deduced from the Feynman path integral by the method of stationary phase [5,7].
In a scalar quasi-monochromatic approach, that is, neglecting polarisation and spectral structure of radiation, the distribution of a wave field in the plane z can be described by the complex amplitude U(x, y; z).Consider the propagation of a scalar quasi-monochromatic field with wavelength λ o in vacuum along the z-direction through a layer of an inhomogeneous optical medium with the refractive index distribution n(x, y, z) restricted by the input {(x, y)} = {(x, y, z) : z = 0} and the output {(x , y )} = {(x, y, z) : z = Z} planes (see Figure 1).This layer can be described as a linear system which transforms the input complex amplitude U(x, y) into the output complex amplitude U (x , y ) by the superposition integral The kernel K S (x, y; x , y ) of this input-output relation is called an optical propagator.Following Feynman [6,7], the propagator can be written as the integral over all possible virtual paths connecting points (x, y) and (x , y ): where the symbol D(path) is the "path differential measure" signifying integration over all paths, S path (x, y; x , y ) is called the optical length of the virtual path connecting points (x, y) and (x , y ), k o ≡ 2π/λ o is the wave number of the scalar quasi-monochromatic wave in vacuum, and A is a normalization constant.
Let us consider the simplest case.The collection of all virtual paths connecting points (x, y) and (x , y ) can be labelled with a parameter u [8].All these paths can be projected onto the z axis so that they can be described by functions x(z) and y(z).Then the entire assembly of paths running from the input point (x, y) to the output point (x , y ) is given by two functions of two parameters: γ u (z) = {x(z, u), y(z, u)} (see Figure 1).Following Luneburg [1], we use the z-axis as an independent variable, similar to the time axis in analytical mechanics.So the direction in space can be characterized by "velocities": in terms of the angles (θ, ϕ) of the spherical coordinate system (see Figure 2).In this case the optical length of the path takes the following form S u (x, y; x , y ): where is the optical Lagrangian [1,8].In this case the optical propagator (2) can be rewritten as Note that for sufficiently large k o (the short-wavelength approximation), the exponential functions exp[ik o S u (x, y; x , y )] fluctuate so rapidly that they compensate each other, and the final amplitude becomes close to zero.Compensation does not occur in the neighbourhood of the stationary point u (see Figure 1): where the phase function S u (x, y; x , y ) varies relatively slowly, the exponential function constructively interferes, and the final amplitude, the optical propagator equation (6), is rather large.Since outside the neighbourhood of the stationary point u the resulting amplitude nearly equals zero the function S u (x, y; x , y ) can be approximated by a parabola with the vertex at point (u, S u ) (see Figure 1): where Light ray γ u (z) (1) A collection of virtual paths connecting an input point (x, y) with an output point (x , y ) can be labelled with a parameter u : γ : (z, u) → (x, y; z) ∈ R 3 , that is, γ u (0) = (x, y; 0) = (x, y) and γ u (Z) = (x, y; Z) = (x , y ).( 2) All virtual paths going through the optical inhomogeneous medium n(x, y, z) from the input point (x, y) to the output point (x , y ) have the optical length S : u → γu n(x, y, z) The phase function S u (x, y; x , y ) in the neighbourhood of the stationary point u, where ∂S u /∂u| u=u = 0 can be approximated by a parabola with the vertex at point (u, S u ).In this case the optical propagator can be approximated by the function where S u (x , y ; x, y) is a point eikonal, that is, the optical length of the path which satisfies Fermat's principle.ζ) can be characterized by angles (θ, ϕ) of the spherical coordinate system (ϕ ∈ (0, 2π), θ ∈ (0, π/2)), or by "velocities" ( ẋ, ẏ).Thus, in (6) we can substitute the phase function S u (x, y; x , y ) by its approximation (8): As the integral has the closed form [10,21,22] Equation ( 8) can be rewritten as K S x, y; x , y = K S x, y; x , y • exp ik o S u x, y; x , y . ( Equation ( 12) is the stationaryphase approximation of the propagator (6) [5].Note that (12) creates the impression that the wave runs along the path γ u (z) whose optical length is stationary with respect to small perturbations of u, Equation (6), that is, whose optical length satisfies Fermat's principle and the phase function S u (x , y ; x, y) is a point eikonal [5].Thus, using the method of stationary phase, we have shown that Feynman's path integral reproduces Fermat's principle and the point eikonal (see sections 2 and 4 in [8]).

k o -Fourier Transformations Reproduce the Legendre Transformations
For true understanding of problems, it is necessary to consider them from various points of view.This possibility is given by one-to-one mappings such as the Legendre transformations in geometrical optics and the Fourier transformations in wave optics.Note that to understand optics as a whole, one needs to know not only these transformations alone, but also the one-to-one mapping between them which is given by the method of stationary phase [5,12,13,21,22].The direct k o -Fourier transformation of the complexvalued function U(x, y) of real variables (x, y) is a new complex-valued function U(p, q) of new real variables (p, q), which is defined as The inverse k o -Fourier transformation of the complexvalued function U(p, q) of real variables (p, q) is a new complex-valued function U(x, y) of new real variables (x, y), which is defined as The composition of the direct k o -Fourier transformation and the inverse k o -Fourier transformation (or the composition of the inverse k o -Fourier transformation and the direct k o -Fourier transformation) is the identity transformation I: If the complex functions are connected by the k o -Fourier transformations, the phases of these functions are connected by the Legendre transformations [21,22].
Note that if the complex functions U(x, y) and U(p, q) are factorized into the amplitude and phase factors: then their direct and inverse k o -Fourier transformations (13a), and (13b) take the forms The phase functions (x, y) − xp − yqΨ and Ψ(p, q)+xp + yq have stationary points (p, q) and (x, y), respectively, that is, These stationary points are nondegenerate: det Using the stationaryphase approximation (see Appendix A), the Fourier transformations can be estimated by the equations The phase function Ψ(p, q) is the direct Legendre transformation of the continuously differentiable real-valued phase function Ψ(x, y), that is, a new real-valued phase function Ψ(p, q) of new real variables (p, q), which is defined as where p ≡ ∂Ψ ∂x , q ≡ ∂Ψ ∂y .

(19a)
The phase function Ψ(x, y) is the direct Legendre transformation of the continuously differentiable real-valued phase function, that is, a new real-valued phase function Ψ(p, q) of new real variables (x, y), which is defined as If the function Ψ(x, y) is positive definite, det (or negative definite, det If the functions Ψ(x, y) and Ψ(p, q) are positive or negative definite functions, then the composition of the direct Legendre transformation and the inverse Legendre transformation (or the composition of the inverse Legendre transformation and the direct Legendre transformation) is the identity transformation I:

Dual Meaning of Spatial Frequencies
In Fourier optics, the spatial frequency has a dual meaning.
On the one hand it is a sinusoidal grating, and on the other hand it is the momentum of the monochromatic plane wave which goes through this grating.To analyze this in detail, let us consider the simplest optical system, that is, a layer of an optically homogeneous medium with refractive index n and thickness Z bounded by parallel input (x, y) = (x, y, 0) and output planes (x z , y z ) = (x, y, Z) (see Figure 3).It is known that a scalar quasi-monochromatic plane wave travelling through this layer in a direction specified by Figure 4: Direction of the "ray vector" can be characterized by angles (θ, ϕ) of the spherical coordinate system (ϕ ∈ (0, 2π), θ ∈ (0, π/2)), or by momenta (p, q). the vector (p, q, n 2 − p 2 − q 2 ) (see Figure 4) can be described by the equation Using the direction cosines p and q as two independent variables, we can describe the complex amplitude on a plane Z, U(x, y; Z), as a linear superposition of plane waves with weight function U in (p, q) [9-13]: Let us consider the layer of a homogeneous optical medium as a linear system (see Figure 3) [9][10][11][12][13] which transforms the input complex amplitude U in (x, y) ≡ U(x, y; 0) to the output complex amplitude U out (x z , y z ) ≡ U(x, y; Z).Then using the inverse k o -Fourier transformation ( 14), (24a) and (24b) can be rewritten twice where the function is known as the optical Hamiltonian [1,5,8].
International Journal of Antennas and Propagation Thus, the weight function U in (p, q) = ko F x → p y → q {U in (x, y)} can be interpreted as an input frequency spectrum, and the variables (p, q) ∈ R 2 -as spatial frequencies.According to (24b), the input frequency spectrum and the output frequency spectrum U out (p, q) = ko F xz → p yz → q {U out (x z , y z )} are related by the equation where ] is the transfer function of the layer of a homogeneous optical medium with refractive index n of thickness Z.Thus, in region p 2 + q 2 ≤ n 2 variables (p, q) can be interpreted not only as spatial frequencies but also as the direction cosines of the wave vector.This manifests the so called dual meaning of spatial frequencies [11].
According to the convolution theorem [9,10], Equation ( 26) can be rewritten as where is the impulse response function of the layer of a homogeneous optical medium with refractive index n of thickness Z, and is the impulse response function of the layer of a homogeneous optical medium with refractive index n of unit thickness.Using the stationaryphase approximation of the inverse Fourier transformation (18b), the impulse response function (28b) can be rewritten as: Here where Note that (31) can be interpreted as Hamilton's equations which describe propagation of a light ray through an inhomogeneous optical medium with the refractive index distribution n(x, y, z).
In the optically homogeneous medium (n = const) the Legendre transformation (30) of the Hamiltonian (25) is the optical Lagrangian L( ẋ, ẏ) = n 1 + ẋ2 + ẏ2 [1,8,15,16].Since the impulse response function (29) takes the form or, according to (28a), Let us compare the approximation (34) to the well-known exact solution of (28a), and (28b), the so-called Rayleigh-Sommerfeld impulse response function given by [9-13] For large values of Z, the second term between the parentheses is smaller than the first term, and in this most important case approximation ( 34) is very good.
Using the method of stationary phase we have shown that the Fourier transformations between the transfer function and the impulse response function of the layer of a homogeneous optical medium with refractive index n of the unit thickness reproduce the Legendre transformations between the Hamiltonian and the Lagrangian (see Section 3 in [8]).

Input-Output Transformations of Wave Optics Reproduce Eikonals Theory
From the possibility of both coordinate and frequency descriptions of the wave field at the input and output of the optical system follows the existence of four equivalent ways of wave description of this system in terms of "input-output transformations."Using the method of stationary phase, Walter showed that each of these transformations of the wave field reproduces the eikonal which describes the transformation in terms of geometrical optics.Consider a linear optical system bounded by parallel input (x, y) and output (x , y ) planes lying in layers of optically homogeneous media with refractive indices n and n , respectively, and a layer of an optically inhomogeneous medium between them (see Figure 5).Variables with primes refer to the output plane.
As the radiation field in the input and output planes of the linear optical system can be specified either in the coordinate representations U(x, y) and U (x , y ) or in the frequency representations U(p, q) and U (p , q ), according to the superposition principle, the linear optical system can be described by any of the four integral equations [5,[12][13][14]: The kernels of these linear input-output transformations can be called the point K S (x, y; x , y ), point-angle K V (x, y; p , q ), angle-point K V (p, q; x , y ), and angle K T (p, q; p , q ) propagators.
The complex amplitudes, U(x, y) and U (x , y ), and frequency spectra, U(p, q) and U (p , q ), are connected by k o -Fourier transformations.As a result, the propagators of the optical system are connected by the Fourier transformations too [5,[12][13][14]: In these expressions, the frequency representations are transformed into the coordinate representations by the direct Fourier transformations in the input plane and by the inverse Fourier transformations in the output plane.
Note that propagators can be factorized into amplitude and phase factors: K V p, q; x , y = K V p, q; x , y •exp ik o V p, q; x , y , (38c) K T p, q; p , q = K T p, q; p , q • exp ik o T p, q; p , q .(38d) Since the propagators of the optical system are connected by the Fourier transformations (37), then according to (18a), and (18b), their phase functions are connected by the Legendre transformations:
It is known that the point-angle eikonal of the perfect optical system has the form and the point-angle eikonal of the isoplanar (or shift invariant) optical system has the form where W(p , q ) is the wave aberration [23,24].Note that in an optical system with telecentric stop in object space (see Figure 6) the wave aberration is determined in the region of moments limited by the telecentric stop with the radius R: p 2 + q 2 ≤ R/ f , where f is the back focal length (distance between front focal F and front nodal points N of the optical system).In the isoplanar optical system, the point eikonal S isoplanar (x, y; x , y ) and the point propagator K isoplanar s (x, y; x , y ) depend only on the difference variables: In this case the functions K(p q ) and K(x − mx, y − my) are called the transfer function and the impulse response function of the isoplanar optical system, respectively.
Using the method of stationary phase, we have shown that the Fourier transformations between propagators of optical system (37) reproduce the Legendre transformations between its eikonals (39) (see Section 4 in [8]).In the same way, Fourier transformations between the transfer function and the impulse response function of an isoplanar optical system (43) reproduce the Legendre transformations between its wave aberration function and its point eikonal (42).

Quantization: Hamiltonian Reproduces Helmholtz Equation
According to Gloge and Marcuse a formal quantization of the Hamiltonian of geometrical optics leads to the Helmholtz equation of wave optics [15,16].Let as show that this formal procedure can be substantiated by using the relations between the Fourier and the Legendre transformations, (18a), and (18b).Note that the Fourier transformations between inputoutput relations (37) and the Legendre transformations between eikonals (39) are valid for any position of the input and output planes.Let the input and output planes be separated by a layer of a medium of small thickness ε.Then the differential equation of wave optics (the Helmholtz equation) can be derived from the integral superposition (2).
Consider the particular case where the input plane ζ and the output plane ζ +ε are separated by a plane-parallel layer of medium with smooth homogeneities (see Figure 7).Let the thickness ε be less than the characteristic size of the optical homogeneities, but substantially greater than the radiation wavelength λ o /n [5].In this case the medium in a neighbourhood of a point (x, y) can be considered as quasi-homogeneous and the integral superposition equation (36a) with the propagator (38a) can be approximated by the equation [5,6] International Journal of Antennas and Propagation Using the operator representation of the Taylor series (Appendix B), the function U(x − ε ẋ, y − ε ẏ; ζ) can be given as Substituting this expression into the left-hand side of (44), we get or, using the definition of the Fourier transformation (13a), If L( ẋ, ẏ, x, y; ζ) is a smoothly varying function, then, according to the asymptotic (18a), the phases of the functions connected by the Fourier transformation are connected by the Legendre transformation: where The approximation of (48) makes it possible to simplify the right-hand side of (47) and to write it in the differential form Comparing (49) and the optical Hamiltonian ( 25) [1,8,15,16], we can see that the operator H((−i/k o )(∂/∂x), (−i/k o )(∂/∂y), x, y; ζ) is obtained from the optical Hamiltonian by replacing coordinates x, y and momenta p and q with the corresponding operators in the position representation [5,25] x that is, This makes it possible to interpret (50) as the symbolic representation of the Helmholtz equation [5,15,16] So we substantiate the "quantum theory of light rays" (see Appendix C) [15,16]: the operator representation of the Helmholtz equation can be obtained from the optical Hamiltonian by replacing the momenta p and q with the corresponding momentum operators in the position representation.

International Journal of Antennas and Propagation
Note that the symbolic representation of the Helmholtz equation (50) has the complex conjugation analog where the asterisk * denotes complex conjugation.It is shown that formal quantization, that is, replacing coordinates x, y and momenta p and q with the corresponding operators, (52), can be substantiated by using the relation between the Fourier and the Legendre transformations.

Position Representation and Momentum Representation of Operators
It is known that a wave field of radiation can be described in the position and the momentum representations which are connected to each other by Fourier transformations [25].
The fundamental symmetry of the operator theory is the duality principle: the k o -Fourier transformation (13a) connects the complex amplitude U(x, y; z) of the wave field with its spectrum U(p, q; z), and also replaces the operators in position representation (51) with the operators in the momentum representation [25]: More precisely, Therefore the k o -Fourier transformation (13a) replaces the Hamilton operator in the position representation (52) with the Hamilton operator in the momentum representation: That is true for the complex conjugation Hamilton operator as well: Thus, the Hamilton operator in the momentum representation (58) can be derived from the optical Hamiltonian ( 25) by replacing coordinates x, y and momenta p and q with the corresponding operators in the momentum representation (55).
If the Helmholtz equation in the position representation has form (50), then in the momentum representation it takes the following form:
The fundamental equations of optics are the equations describing propagation of light through an optically inhomogeneous medium.They are Hamilton's equations in geometrical optics, (31), and the Helmholtz equation in wave optics, (53).Propagation of the WDF through an optically inhomogeneous (and random) media is described by the socalled transport equation [32,33].Let us show that the duality between the position and momentum representations of the Helmholtz equation, ( 50) and ( 59), allows us to prove it.
If the WDF (61a) is differentiated with respect to the longitudinal variable ζ, then we have ∂W x, y; p, q; ζ ∂ζ so (65) takes the following form: ∂W p, q; x, y; ζ ∂ζ (65) Using (57a), (57b) and (65) can be written as ∂W p, q; x, y; z ∂ζ Using the identity where the symbol Im H denotes the imaginary part of the Hamilton operator H, (66) can be rewritten in the symmetric form [32,33]: ∂W p, q; x, y; z ∂ζ We can approximate the Hamiltonian in (68) via a Taylor series expansion by the Appendix B In that case the transport equation (68) takes the form ∂W p, q; x, y; ζ ∂ζ = −2 Im H p, q; x, y; ζ or, in the symbolic notation of Besieris and Tappert [34], ∂W p, q; x, y; ζ ∂ζ In this expression the arrows over the partial derivatives indicate the direction in which the derivatives act, that is, those with a left arrow act on H(p, q; x, y; ζ), and those with a right arrow act on W(p, q; x, y; ζ).
In a weakly inhomogeneous medium the transport equation (71) reduces to ∂W p, q; x, y; ζ ∂ζ In the special case that H(p, q; x, y; ζ) is a real function of x, y, p, and q, Equation (75) can be formulated in geometric optical terms as follows: along a light ray the WDF has a constant value: It is known that the WDF with some natural assumptions can be interpreted as a wave analogue of the radiance [1].So equation ( 76) is the theorem about invariance of radiance along the light ray, that is, the most fundamental theorem of radiometry.Note that this theorem in the form of the Liouville equation (73) can be derived from the point of view of Hamilton's optics, too [20].

Conclusion
We have unite all known quantum-mechanical approaches to wave and statistical optics in a single theory using Feynman's path integral as a fundamental principle.In particular, we have substantiated the quantization which allows one to translate the Hamiltonian of geometrical optics into the Helmholtz equation of wave optics.

A. The Stationary-Phase Approximation
Wave optics often takes account of oscillatory integrals of the type Then, according to the stationary phase method [5, 12-14, 21, 22], the integral (A.1) can be approximated by the equation

B. The Taylor Series in Operator Form
Using the Taylor series, a function U(x − Δx, y − Δy) can be expanded in a power series up to the linear term:

C. Formal Quantum Theory of Light Rays
In the quantum theory of light rays, the coordinates x, y retain their meaning as numbers, but the canonically conjugate variables p, q (the momenta), ( 52

Figure 1 :
Figure1: Illustration of the path integral applied to a layer of an inhomogeneous medium.The path integral contains three elements.(1)A collection of virtual paths connecting an input point (x, y) with an output point (x , y ) can be labelled with a parameter u : γ : (z, u) → (x, y; z) ∈ R 3 , that is, γ u (0) = (x, y; 0) = (x, y) and γ u (Z) = (x, y; Z) = (x , y ).(2) All virtual paths going through the optical inhomogeneous medium n(x, y, z) from the input point (x, y) to the output point (x , y ) have the optical length S : u → γu n(x, y, z) 1 + ẋ2 + ẏ2 dz.(3) The path integral K S (x, y; x , y ) =

Figure 3 :
Figure 3: A layer of an optically homogeneous medium as a linear system.

Figure 6 :
Figure 6: An isoplanar optical system with telecentric stop in object space.

Figure 7 :
Figure 7: Thin layer of an optically quasi-homogeneous medium.

R 2 g 2 )
x, y exp ik o Ψ x, y dx dy (A.1) in which g(x, y) is a slowly varying complex-valued function and Ψ(x, y) is a slowly varying real-valued function of the real variables (x, y).We assume that the phase functionΨ(x,y) has one stationary point (x o , y o ) in the interior of the region of integration: The stationary point (x o , y o ) is nondegenerate det ⎛ ⎝ Ψ xx x o , y o Ψ xy x o , y o Ψ yx x o , y o Ψ yy x o , y o

R 2 gΨ− 1 / 2 •
x, y exp ik o Ψ x, y dx dy ≈ xx x o , y o Ψ xy x o , y o Ψ yx x o , y o Ψ yy x o , y o ⎞ ⎠ g(x o , x o ) exp ik o Ψ x o , y o .(A.4)

2 n 2
), and the Hamiltonian H, (51), become differential operators p is the free-space wavelength.Thus, the operator equation of the quantum theory of light rays can be obtained from the optical Hamiltonian(25) and operator equation (C.1c) by replacing momenta p and q with the momentum operators (C.1a) and (C.1b) − n 2 x, y; z + U(x, y; z) of the quantum theory of light rays is the usual distribution of the scalar wave field over the plane z.Applying the operator equation (C.2) to the wave function U(x, y; z) results in the equation takes into account the symmetry of the optical medium relative to the 0z axis.Equation (C.3) can be interpreted as the symbolic form of the Helmholtz equation: x, y; z U = 0.(C.4) [5,35]; ζ ∼ = U x, y; ζ −Using the Taylor series, we can expand the exponent exp(−(∂U/∂x)Δx, −(∂U/∂y)Δy) in a power series up to the linear term: Comparing (B.2) with (B.3), we see that the right-hand sides of these equations are equal, so we can write the approximate operator equation[5,35]:U x − Δx, y − Δy; ζ = exp −