## Abstract

A theoretical analysis is made of Gaussian beam propagation in media where the refractive index and/or the gain constant varies quadratically with the distance from the optic axis. Use of a complex beam parameter simplifies the analysis. A differential equation for the complex beam parameter is deduced from the wave equation for the complex beam parameter is deduced from the wave equation, and various of its solutions are discussed. This includes a discussion of light propagation in media with a gain profile, which are found capable of supporting stationary beams of constant diameter.

© 1965 Optical Society of America

## Introduction

This paper considers the propagation of optical modes or Gaussian beams of light of the kind produced by laser oscillators or used in microwave beam waveguides. It discusses the propagation of these beams through optical systems like lenses, gas lenses, and quite general lenslike media including media exhibiting a nonuniform loss or gain profile. An example of a medium with a nonuniform gain profile is the medium of a gas laser amplifier where the gain decreases with the distance from the axis of the laser tube.

The laws governing the propagation of the light beam will be formulated in terms of a complex valued beam parameter. This complex parameter gives information on the width of the beam (spot size), and on the curvature of the phase front in each beam cross section of interest. The advantage of using a complex beam parameter is apparent from the recent literature. Following the initial proposal by Collins various circle diagrams for Gaussian beams were introduced,, which are (or can be regarded as) plotted in the complex plane of a beam parameter of this kind. It was also found that the laws of beam propagation in homogeneous media can be written in a very simple form if a complex notation is used. It will be seen that the use of a complex parameter also simplifies the analysis of beam propagation in lenslike media.

First, the paper recalls the laws of propagation for spherical waves in lossless lenslike media in the geometrical optics approximation. The variation of the phase fronts of these waves along the optic axis is described in terms of the elements of the ray matrix. These matrix elements can be obtained by tracing paraxial rays through the medium. Then Gaussian beams in lossless lenslike media are considered, and it is found that the law for the transformation of the curvature radius of spherical waves is formally the same as the transformation law for the complex beam parameter (ABCD law).

The paper continues by discussing wave propagation in lenslike media where loss or gain varies with the square of the distance from the optic axis. Assuming paraxial waves of small wavelengths, a parabolic equation is obtained from the scalar wave equation for these media. This parabolic equation is a generalization of the parabolic equation for a homogeneous medium which is used in diffraction theory. Light beams with a Gaussian amplitude profile are obtained as a solution of this parabolic equation together with a differential equation of first order. It can be transformed into a linear differential equation of second order which has the form of the paraxial ray equation for the medium. The connection between the beam parameter equation and the differential equation used by Tien et al. is studied and several solutions are discussed. Beam propagation in media with a gain profile is investigated in greater detail.

The paper is concluded by a brief discussion of analog electric circuits which correspond to the various lenslike media. The input and output admittances of these circuits (or impedances in the dual circuits) correspond to the complex parameters of the beam in the input and output planes of the optical system as in the lumped circuit analogs described by Deschamps and Mast.

## Propagation of Spherical Waves in the Geometrical Optics Approximation

Some well-known properties of spherical waves are recalled here for later comparison with the laws obeyed by Gaussian beams. Consider first the curvature of two phase fronts of a spherical wave propagating in the z direction. Call the radii of curvature of the two phase fronts R1 and R2, respectively, and z the distance between the intersections of these phase fronts with the optic axis. It is quite obvious from Fig. 1 that the curvature radii are related by

$R2=R1+z.$
Here and throughout this paper a curvature radius is defined as positive if the phase front is concave as seen from the left (z = −∞).

If a spherical wave passes through a lens of focal length f, it is transformed into another spherical wave. For an incoming wave with a radius of curvature R1 at the lens the wave going out from the lens has a curvature radius R2 given by

$1R2=1R1−1f.$

Consider now paraxial light rays passing through a medium whose refractive index near the optic axis is described by

$n(r,z)=n0(z)−12n2(z)·r2.$
This kind of medium is called a lenslike medium. Examples are inhomogeneous crystals and gas lenses. The propagation of paraxial light rays in a lenslike medium is described by the paraxial ray equation,,
$ddz(n0dxdz)=−n2x,$
where x(z) indicates the x coordinate of the ray position in the plane z = const. A similar equation holds for the y coordinate. In Fig. 2 there is schematically indicated a ray path through a lenslike medium and an input and an output plane. The quantities of interest are also indicated. They are the ray position x2 and the ray slope x2′ in the output plane and the corresponding ray parameters x1 and x1′ in the input plane. The solution of the paraxial ray equation yields a linear relation between the input and output parameters which is conveniently written in matrix form,
$(x2x2′)=(ABCD)(x1x1′).$

The passage of rays through lenses and lens combinations can also be described in this matrix form, and the elements A, B, C, and D of the ray matrix are known for several optical structures. The matrix elements are determined by the over-all focal length of the optical structure of interest and by the position of its principal planes.,

The rays associated with a spherical wave are perpendicular to the wavefront. The position x and slope x′ of a paraxial ray are related to the radius of curvature R of the wavefront by

$R=xx′·$
By tracing a paraxial ray through a given optical structure one can determine the transformation of a wavefront of a spherical wave passing through this structure. For a spherical wave with a radius of curvature R1 at the input plane of the structure, one finds for the curvature radius R2 at the output plane
$R2=AR1+BCR1+D·$

## Propagation of Gaussian Beams, ABCD Law

The properties of Gaussian beams of light in homogeneous media are well known. Near the optic axis the field distribution E(r,z) of a fundamental Gaussian mode is described by

$E(r,z)=w0wexp{−j(kz+φ)−r2(1w2+jk2R)},$
where a constant amplitude factor was dropped. The light is assumed to propagate in z direction, r measures the distance from the z axis (optic axis), k = 2π/λ is the propagation constant in the medium, and φ (z) indicates an additional phase shift due to the geometry of the beam. The radius of the beam (spot size) is w(z), and the radius of curvature of the phase front at z is R(z).

As the light beam propagates in space it expands due to diffraction, but the transverse field distribution remains Gaussian. The law of expansion is

$w2=w02[1+(λzπw02)2],$
where z is measured from the beam waist. At the waist the radius of the beam is w0, and the phase front is plane. The radius of the phase-front curvature changes along the optic axis according to the law
$R=z[1+(πw02λz)2],$
and the additional phase shift is
$φ=tan−1(λzπw02).$
It is the main object of this paper to determine the parameters w, R, and φ of a Gaussian beam which has passed through a lenslike medium of the most general type.

It is convenient to combine the beam parameters w and R in a complex-valued parameter q defined by

$1q=1R−jλπw2.$
For homogeneous media the transformation laws (9) and (10) can then be written as one equation in the simple form
$q2=q0+z,$
where
$q0=jπw02λ.$
Instead of defining z as the distance from the beam waist one can now redefine z as the distance between any two reference planes in the medium. Because of the linearity of (13) the relation between the input parameter q1 and the output parameter q2 can then be written as
$q2=q1+z.$

The parameters of a Gaussian beam are transformed when it passes through an ideal thin lens. The beam radii remain unchanged, i.e., the beam radius w1 immediately to the left of the lens will be equal to the beam radius w2 immediately to the right of it. The curvature radii R1 and R2 of the corresponding phase fronts are transformed in the same way as the curvature radii of spherical waves, i.e., they are related by (2). Therefore, the lens transformation law for the complex beam parameter is ,

$1q2=1q1−1f.$

The laws for the transformation of Gaussian beams are seen to be formally the same as the transformation laws obeyed by the spherical waves of geometrical optics. Equation (15) is formally the same as (1), and (16) is formally the same as (2). The complex beam parameter q plays a role corresponding to the one played by the radius R of the spherical waves and might be called a complex curvature radius.

Like other optical structures a lenslike medium can be regarded as composed of a set of thin lenses. The number of lenses required for a fair representation of the medium may, of course, be very large. But the fact remains that (15) and (16) can be used to trace a Gaussian beam through each lens element, while (1) and (2) are used to trace a spherical wave. As a result of ray tracing one has the law (7) for the transformation of a spherical wave by the optical structure. Because of the formal equivalence discussed above, the transformation law for Gaussian beams must be

$q2=Aq1+BCq1+D.$
This law has been called the ABCD law. It allows one to trace Gaussian beams through any optical structure for which the elements A, B, C, and D of the ray matrix (or the focal length and the principal planes) are known. Its application to some optical structures of interest has been given in [ref. 3].

## Almost-Plane Waves in Lenslike Media

Pierce starts from the wave equation to derive the propagation laws for Gaussian beams in free space and in sequences of lenses. A modification of his results was used by Tien et al. for their study of beam propagation in various lenslike media. The present derivation of the propagation laws for Gaussian beams in quite general lenslike media starts also with the scalar wave equation*

$ΔE+k2E=0,$
where E is a scalar field or potential, and k is defined by
$k2=ω2∊μ(1−jσ∊ω).$
Here ω is the angular frequency of the field, is the dielectric constant of the medium, σ the conductivity, and μ the permeability. For fields oscillating like exp (jωt) the values of σ are positive for lossy media, and negative for media with gain. A lenslike medium of the most general type is characterized by a quadratic radial dependence of k2 near the optic axis, i.e., one can write
$k2(r,z)=k02−k0k2r2.$
The medium constants k0(z) and k2(z) are generally complex valued to allow for loss or gain. They are assumed to be slow functions of z.

The main interest here is in waves that propagate primarily in the z direction, or almost-plane waves. Accordingly, one writes

$E=ψ(x,y,z)·exp(−jΦ),$
with
$Φ′(z)=k0(z),$
and the prime indicating differentiation with respect to z. For k0 = const, the exponential in (21) becomes the usual exp (−jk0z). The function ψ is assumed to vary so slowly with z that its second derivative can be neglected. If one inserts (20), (21), and (22) into the wave equation (18) and neglects ψ″ one obtains
$Δtψ−2jk0ψ′−jk0′ψ−k0k2r2ψ=0,$
where
$Δt=Δ−∂2∂z2=∂2∂x2+∂2∂y2.$
The parabolic equation (23) is a generalized form of the parabolic equation for lossless homogeneous media which was used by Leontovich et al. in their analyses of various diffraction phenomena. This parabolic equation takes into account the transverse diffusion of the field, but it neglects longitudinal diffusion. It yields good results for paraxial waves whose wavelengths are much smaller than their transverse dimensions.

A wave with a Gaussian transverse field distribution is an exact solution of Eq. (23). To determine the parameters of the Gaussian beam one can write

$ψ=exp−j(p+12Qr2).$
using a complex phase parameter P, and a complex beam parameter Q which corresponds to the parameter 1/q used above. Both parameters are functions of z. By inserting (25) into (23) one gets
$Q2r2+2jQ+k0(2P′+Q′r2)+jk0′+k0k2r2=0.$
Now one compares terms with equal powers of r, and obtains an expression for the complex phase P
$P′=−j(Qk0+k0′2k0),$
as well as a differential equation for the complex beam parameter Q
$Q2+k0Q′+k0k2=0.$

Equation (28) is a nonlinear differential equation of the Ricatti type. Its connection to the paraxial ray equation, and several of its solutions will be discussed below.

## Beam Parameter Equation and Paraxial Ray Equation

Solutions of (28) describe the propagation properties of Gaussian beams. These solutions can also be constructed from solutions of the paraxial ray equation. The connection between these two equations is made by

$Q=k0x′x.$
If this is inserted into (28) one arrives at
$(k0x′)′+k2x=0.$
This equation is formally the same as the paraxial ray equation (4) except that one has, in general, complexvalued k0 and k2 instead of the real n0 and n2 which describe lossless media. In media with loss or gain the variable x is also complex-valued, and its physical meaning is not so clear as its meaning as a ray position in lossless media. However, the general solutions of (30) can be combined according to (29) to yield solutions of (28) for the complex beam parameter which is a measure for the width and the phase-front curvature of Gaussian beams. Equation (29) is essentially the ABCD law (17) extended in validity to media with loss or gain.

The example of a medium with k0 = const and k2 = const will help to illustrate the preceding remarks. For k2 ≠ 0 such a medium can support Gaussian beams whose diameter does not change along the optic axis. The Q parameter of this kind of a beam is obtained from (28) by postulating that Q′ = 0. One finds

$Qm=−jk0k2.$
The sign of the square root has been so chosen that a real and positive beam radius w is obtained. The general solution for a beam of varying diameter is easily derived from (30). The general solution of this equation is
$x=aejγz+be−jγz,$
where a and b are constants which depend on the initial conditions at z = 0, and
$γ=k2k0.$
To get a general solution for the complex beam parameter one inserts (32) and its derivative with respect to z into (29). After some rearranging one arrives at
$Q2=QmQ1+Qm+(Q1−Qm)e2jγzQ1+Qm−(Q1−Qm)e2jγz,$
where Q1 is the beam parameter at the input plane at z = 0 and Q2 is the output parameter at z. Use of this general formula will be made later when the effect of a medium with a gain profile on Gaussian beams is studied in greater detail.

## Refractive Indices and Gain Constants

For a more detailed study of beam propagation the complex propagation constant is separated into its real and imaginary part

$k=β+jα,$
where β is the real propagation constant, and α is the gain constant (which is negative for loss). At optical and infrared frequencies the gain per wavelength is usually very small (αβ), and one can write to a good approximation
$k2=β2+2jαβ.$
To find a relation between α and the conductivity one compares (36) with (19) and gets
$α=−σ2μ∊.$
In a medium described by (35) the power of a plane wave grows like
$P2=P1:e2αz.$
This power gain is often measured in decibels (dB). The gain constant α is related to the dB gain per distance z by
$α=dB gain8.686(z)$
as one easily computes.

Near the optic axis the propagation constant of a lenslike medium changes generally as in

$β=β0−12β2r2.$
Closely related to β is the refractive index of the medium
$n=λ2πβ=n0−12n2r2,$
where λ is the wavelength in free space. In the most general type of a lenslike medium one has also a gain constant that varies quadratically with the distance from the optic axis
$α=α0−12α2r2.$

It is not difficult to relate the above gain and propagation constants to the complex constants used in the beam parameter equations. It follows from (20) that, near the axis, the complex propagation constant is

$k=k0−12k2r2.$
Consequently,
$k0=β0+jα0,$
and
$k2=β2+jα2.$

Keeping in mind that the gain per wavelength is small, one can write for the product k0k2

$k0k2=β0β2+j(α0β2+α2β0).$

## Solution for Beams in Free Space

In this section the solution of the beam parameter equation for Gaussian beams in free space is discussed briefly. This will serve to connect the present method to the published literature on the subject, and to demonstrate that the parabolic equation (23) leads to the same results as the theory of Fresnel diffraction.

In free space one has α0 = α2 = β2 = 0 and β0 = 2π/λ. The beam parameter q defined in (12) is related to the parameter Q defined in (25) by

$1q=λ2πQ.$
If this is introduced in (28), the beam parameter equation becomes
$dqdz=1.$
The solution of (48) is, clearly,
$q2=q1+z,$
which is seen to be identical to (15).

The complex phase shift on the axis is computed by means of (27). Together with (47) and (49) this gives

$p2−p1=−j∫0zdzq1+z=jlnq1q1+z.$
If one measures z from the beam waist, one has a q1 = q0 that is purely imaginary and given by (14). As is well known, the real and imaginary parts of the phase shift can then be separated by using the identity
$ln(u+jυ)=Inu2+υ2+jtan−1υu.$
The imaginary part of the phase shift produces the factor w0/w appearing in (8), and the real part predicts the same additional phase shift as (11).

## The Lossless Lenslike Medium

For a lossless medium one has α0 = α2 = 0. If the definitions of the refractive index (41) and of the beam parameter q (47) are introduced, the beam parameter equation (28) can be written in the form

$1q2+n0(1q)′+n0n2=0.$
Because of the recent interest in gas lenses, several studies have been made,, of Gaussian beam propagation in lenslike media. A familiar result for media with constant n0 and n2 is the guided beam of constant diameter 2wm whose parameter qm is obtained from (52) by putting q′ = 0.
$1qm=−jn0n2=−jλπwm2.$
Tien et al. have studied the propagation of Gaussian beams through various media, each characterized by refractive indices n0(z) and n2(z) with a certain z dependence of practical interest. They have based their work on a differential equation for the beam radius w. In order to see the connection between their differential equation and the beam parameter equation (52) one introduces in the latter the definition of the q parameter given in (12). By comparing the imaginary parts of (52) one then obtains
$1R=n0w′w,$
which is a relation that has already been given by Pierce. One then compares the real parts of (52), eliminates R by means of (54), and arrives at
$n0w(n0w′)′=−n0n2+λ2π2w4.$
This is exactly the differential equation (12) of Tien et al. It is interesting to note that the differential equation (52) for the complex beam parameter q has a simpler form than the differential equation (55) for the real beam radius w.

## Medium with a Gain Profile

In this section the propagation of Gaussian beams in a medium with a quadratic gain profile is studied in greater detail. This is a medium where the gain decreases with the distance from the axis according to (42). The plasma medium of various gas laser tubes is of that nature.

For simplicity it is assumed here that the gain constants are independent of z, i.e., α0 = const and α2 = const. The refractive index of the medium is taken to be that of free space, β0 = 2π/λ and β2 = 0. With these assumptions (46) becomes

$k0k2=j2πλα2.$
From this and (31) it follows that a medium with a gain profile is capable of supporting a beam of constant diameter even though there is no refractive index variation. The q parameter of this stationary beam is computed from (31) and (56) as
$1qm=λ2πQm=12(1−j)α2λπ.$
Inspecting this expression one concludes that the stationary beam of constant width has phase fronts whose curvature is finite, positive, and independent of z. It is indicated by (57) that the radius of curvature Rm of these wavefronts is
$Rm=πwm2λ=2πα2λ⋅$
and that the radius of the beam wm is given by
$wm2=2λπα2.$
One notes that the real and imaginary parts of 1/qm are of equal magnitude. This reminds one of a Gaussian beam in free space where the parameter 1/q has the same property at a distance z = πw02/λ (the confocal distance) from the beam waist. In this particular beam cross section the phase-front curvature is stationary and the corresponding radius of curvature goes through its minimum value. In free space the beam diameter expands at the confocal distance, but the gain profile of the present medium counteracts this tendency. If the diffractional and geometrical expansion of the beam balances the size-reducing effect of the gain profile, the beam diameter remains stationary. This leads to a stationary beam with the properties discussed above.

For several practical laser tubes it is reasonable to assume that the gain constant has its largest value at the tube axis, decreases quadratically, and is zero at the tube walls. If one calls the tube radius r0 one can, then, write for the gain constant

$α=α0(1−r2r02).$
Comparing with (42) one has
$α2=2α0/r02,$
while (58) and (59) become
$Rm=r02πλα,wm2=r02λπα0.$

The reported gain of a He–Xe amplifier tube of 4-mm diam (r0 = 2 mm) at a wavelength of 3.5 μ is about 100 dB/m. From this value one computes with (39) a gain constant on the axis of approximately

$α0=11.5m−1.$
By inserting these values in (62) one calculates for Rm and wm the values
$Rm=79cm,wm=0.94mm.$
One notes that, at these high gain values, the stationary beam radius is smaller than the radius of the tube, and the tube can guide its own beam, provided there is no saturation. For small gain values the stationary radius wm is larger than r0, and the tube cannot support a stationary beam. But even then wm can be used as a parameter to characterize the effect of the gain profile on the expansion or contraction of Gaussian beams with radii w that are smaller than r0.

For the stationary beam the formula (27) for the complex phase shift becomes

$p′=−jqm=−(1+j)zRm.$
If this is integrated, one gets
$p2=p1−(1+j)1Rm.$
The imaginary part of this phase shift indicates the reduction in gain due to the gain profile. One concludes that the amplitude of the stationary beam grows along the axis like
$exp(α0−1Rm)z.$

A stationary beam is obtained only when the light is so injected into the medium that the incoming beam matches the stationary beam in diameter and phase-front curvature. If these matching conditions are not satisfied, the beam will generally fluctuate along z in diameter and curvature. To study this situation in more detail one has to evaluate the constant γ given in (33). For the medium with a gain profile one gets

$γ=k2k0=jqm=1+jRm.$
This is substituted in (34) together with (47) to yield the result
$1q2=1qm(q1+qm)−(q1−qm)exp2(j−1)zRm(q1+qm)+(q1−qm)exp2(j−1)zRm.$
The above formula describes the passage of a Gaussian beam through a medium with a gain profile for any input condition. For q1 = qm one gets q2 = qm which is, of course, exactly what has just been said about matching of the input beam. For very large z(z → ∞) one also gets q2 = qm no matter what the input conditions. This means that the parameters of a beam will first fluctuate more or less around the stationary values wm and Rm but will finally settle down to equilibrium. All these remarks assume, of course, that the light beam remains a paraxial wave while fluctuating, and that there is no interference from the tube walls or from other obstacles.

Equation (69) can be so cast as to show directly the fluctuations of both the wavefront curvature and the beam diameter. In order to do that one defines two angles δ and φ by

$q1−qmq1+qm=exp2(δ−jφ).$
The angles δ and φ are determined by the initial conditions of the beam. If one uses, in addition, the abbreviations
$σ=2(zRm−δ),ρ=2(zRm−φ),$
one can write (69) in the form
$1q2=1❘qmeσ−ejρeσ+ejρ=1qmshσ−jsinρchσ+cosρ.$
Comparison of the real parts of (72) leads to an expression for the fluctuation of the curvature radius R2 of the beam phase front
$R2Rm=chσ+cosρshσ−sinρ.$
The imaginary parts of (73) indicate the fluctuations of the beam radius w2
$w22w22=chσ+cosρshσ+sinρ.$

## Electric Circuit Analogs

The equations (15) and (16) for the complex beam parameter have the same form as the well-known formulas of the input impedance of a series circuit of two impedances and of a parallel circuit of two impedances. On this basis Deschamps and Mast have discussed lumped electric circuits which are analogs of systems of lenses. Two dual circuit analogs of a system of lenses are shown in Fig. 3.

The first is an admittance analog. Here the admittance Y of each section is related to the complex beam parameter in the corresponding beam cross section by

$Y=jq=λπw2+jR.$
The admittance Y1 corresponds to the beam parameter q1 in the input plane. The lens spacing d1 is represented by a capacitance C1 = 1/ωd1 connected in series. The first lens of focal length f1 is represented by an inductance L1 = f1/ω connected in parallel, and so on. Finally, the admittance Y2 of the circuit analog is a measure for the beam parameter q2 in the output plane.

Below the admittance analog there is shown the dual impedance analog. In the latter all admittances from above are replaced by impedance

$z=jq,$
series circuits are replaced by parallel circuits, and vice versa.

The analog of a lenslike medium is a distributed electric circuit or transmission line. The voltage V on such a line obeys the differential equation

$ddz(1ZLdVdz)=YLV,$
where ZL is the line impedance per unit length, and YL is the line admittance per unit length. Equation (77) corresponds to the paraxial ray equation (30) if one puts
$jZL=k0(z)$
and
$YL=jk2(z).$
The transmission line is analogous to the lenslike medium if one interprets the voltage V as ray position x and relates the line current I to the ray slope by
$I(z)=x′ZL.$
For media with loss or gain, of course, the analogous transmission lines are lossy or show gain, and the line impedances and admittance are complex-valued.

The author would like to thank J. P. Gordon, J. W. Kluver, W. H. Louisell, and P. K. Tien for helpful comments and stimulating discussions.

## Figures Fig. 1 Wavefronts of spherical wave. Fig. 2 Ray path in lenslike medium and reference planes. Fig. 3 Gaussian beam propagating through lens system and corresponding circuit analogs.

## Footnotes

 * A treatment of the exact vector problem has to start from Maxwell’s equations, where one obtains this form of the wave equation for the electric field E in a medium of constant permeability u and space-dependent dielectric constant ∊ by neglecting the term grad div $E=−grad[1∊(grad∊.E)].$

1. S. A. Collins, Appl. Opt. 3,1263 (1964); [CrossRef]  T. Li, Appl. Opt. 3, 1315 (1964); [CrossRef]  J. P. Gordon, Bell System Tech. J. 43, 1826 (1964).

2. G. A. Deschamps and P. E. Mast, Quasi-Optics (Polytechnic Press, Brooklyn, 1964), pp. 379–395.

3. H. Kogelnik, Bell System Tech. J. 44, 455 (1965).

4. D. A. Kleinman, A. Ashkin, and G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).

5. M. A. Leontovich and V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).

6. P. K. Tien, J. P. Gordon, and J. R. Whinnery, Proc. IEEE 53, 129 (1965) [CrossRef]  .

7. D. W. Berreman, J. Opt. Soc. Am. 55, 239 (1965); [CrossRef]  D. Marcuse and S. E. Miller, Bell System Tech. J. 43, 1759 (1964);E. A. J. Marcatili, Bell System Tech. J. 43, 2887 (1964).

8. W. Brower, Matrix Methods in Optical Instrument Design (Benjamin, New York, 1964).

9. G. D. Boyd and J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau and F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv and J. P. Gordon, Proc. IEEE 51, 4 (1963); [CrossRef]  H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).

10. J. R. Pierce, Proc. Natl. Acad. Sci. (U.S.) 47, 1808 (1961) [CrossRef]  .

11. J. W. Kluver, submitted to J. Appl. Phys.

### References

• View by:

1. S. A. Collins, Appl. Opt. 3,1263 (1964);T. Li, Appl. Opt. 3, 1315 (1964);J. P. Gordon, Bell System Tech. J. 43, 1826 (1964).
[Crossref]
2. G. A. Deschamps, P. E. Mast, Quasi-Optics (Polytechnic Press, Brooklyn, 1964), pp. 379–395.
3. H. Kogelnik, Bell System Tech. J. 44, 455 (1965).
4. D. A. Kleinman, A. Ashkin, G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).
5. M. A. Leontovich, V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).
6. P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]
7. D. W. Berreman, J. Opt. Soc. Am. 55, 239 (1965);D. Marcuse, S. E. Miller, Bell System Tech. J. 43, 1759 (1964);E. A. J. Marcatili, Bell System Tech. J. 43, 2887 (1964).
[Crossref]
8. W. Brower, Matrix Methods in Optical Instrument Design (Benjamin, New York, 1964).
9. G. D. Boyd, J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau, F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv, J. P. Gordon, Proc. IEEE 51, 4 (1963);H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).
[Crossref]
10. J. R. Pierce, Proc. Natl. Acad. Sci. (U.S.) 47, 1808 (1961).
[Crossref]
11. J. W. Kluver, submitted to J. Appl. Phys.

#### 1965 (3)

P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]

H. Kogelnik, Bell System Tech. J. 44, 455 (1965).

#### 1961 (2)

G. D. Boyd, J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau, F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv, J. P. Gordon, Proc. IEEE 51, 4 (1963);H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).
[Crossref]

J. R. Pierce, Proc. Natl. Acad. Sci. (U.S.) 47, 1808 (1961).
[Crossref]

#### 1946 (1)

M. A. Leontovich, V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).

#### Ashkin, A.

D. A. Kleinman, A. Ashkin, G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).

#### Boyd, G. D.

G. D. Boyd, J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau, F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv, J. P. Gordon, Proc. IEEE 51, 4 (1963);H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).
[Crossref]

D. A. Kleinman, A. Ashkin, G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).

#### Brower, W.

W. Brower, Matrix Methods in Optical Instrument Design (Benjamin, New York, 1964).

#### Deschamps, G. A.

G. A. Deschamps, P. E. Mast, Quasi-Optics (Polytechnic Press, Brooklyn, 1964), pp. 379–395.

#### Fok, V. A.

M. A. Leontovich, V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).

#### Gordon, J. P.

P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]

G. D. Boyd, J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau, F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv, J. P. Gordon, Proc. IEEE 51, 4 (1963);H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).
[Crossref]

#### Kleinman, D. A.

D. A. Kleinman, A. Ashkin, G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).

#### Kluver, J. W.

J. W. Kluver, submitted to J. Appl. Phys.

#### Kogelnik, H.

H. Kogelnik, Bell System Tech. J. 44, 455 (1965).

#### Leontovich, M. A.

M. A. Leontovich, V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).

#### Mast, P. E.

G. A. Deschamps, P. E. Mast, Quasi-Optics (Polytechnic Press, Brooklyn, 1964), pp. 379–395.

#### Pierce, J. R.

J. R. Pierce, Proc. Natl. Acad. Sci. (U.S.) 47, 1808 (1961).
[Crossref]

#### Tien, P. K.

P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]

#### Whinnery, J. R.

P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]

#### Bell System Tech. J. (2)

H. Kogelnik, Bell System Tech. J. 44, 455 (1965).

G. D. Boyd, J. P. Gordon, Bell System Tech. J. 40, 489 (1961);G. Goubau, F. Schwering, Inst. Radio Engrs. Trans. AP-9, 248 (1961);A. Yariv, J. P. Gordon, Proc. IEEE 51, 4 (1963);H. Kogelnik in Advances in Lasers, A. K. Levine, ed. (Dekker, New York, 1965).
[Crossref]

#### Proc. IEEE (1)

P. K. Tien, J. P. Gordon, J. R. Whinnery, Proc. IEEE 53, 129 (1965).
[Crossref]

#### Proc. Natl. Acad. Sci. (U.S.) (1)

J. R. Pierce, Proc. Natl. Acad. Sci. (U.S.) 47, 1808 (1961).
[Crossref]

#### Zh. Eksperim. i Teor. Fiz. (1)

M. A. Leontovich, V. A. Fok, Zh. Eksperim. i Teor. Fiz. 16, 557 (1946);L. A. Vainshtein, Zh. Techn. Fiz. 34, 193 (1964);Soviet Phys.—Tech. Phys. 9, 157 (1964).

#### Other (4)

W. Brower, Matrix Methods in Optical Instrument Design (Benjamin, New York, 1964).

D. A. Kleinman, A. Ashkin, G. D. Boyd, to be published;H. Kogelnik, Proc. Symp. Quasi-Optics (Polytechnic Inst.Brooklyn, 1964), pp. 333–347;A. G. van Nie, Philips Res. Rept. 19, 378 (1964).

G. A. Deschamps, P. E. Mast, Quasi-Optics (Polytechnic Press, Brooklyn, 1964), pp. 379–395.

J. W. Kluver, submitted to J. Appl. Phys.

### Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

### Figures (3)

Fig. 2 Ray path in lenslike medium and reference planes.
Fig. 3 Gaussian beam propagating through lens system and corresponding circuit analogs.

### Equations (80)

$R 2 = R 1 + z .$
$1 R 2 = 1 R 1 − 1 f .$
$n ( r , z ) = n 0 ( z ) − 1 2 n 2 ( z ) · r 2 .$
$d d z ( n 0 d x d z ) = − n 2 x ,$
$( x 2 x 2 ′ ) = ( A B C D ) ( x 1 x 1 ′ ) .$
$R = x x ′ ·$
$R 2 = A R 1 + B C R 1 + D ·$
$E ( r , z ) = w 0 w exp { − j ( k z + φ ) − r 2 ( 1 w 2 + j k 2 R ) } ,$
$w 2 = w 0 2 [ 1 + ( λ z π w 0 2 ) 2 ] ,$
$R = z [ 1 + ( π w 0 2 λ z ) 2 ] ,$
$φ = tan − 1 ( λ z π w 0 2 ) .$
$1 q = 1 R − j λ π w 2 .$
$q 2 = q 0 + z ,$
$q 0 = j π w 0 2 λ .$
$q 2 = q 1 + z .$
$1 q 2 = 1 q 1 − 1 f .$
$q 2 = A q 1 + B C q 1 + D .$
$Δ E + k 2 E = 0 ,$
$k 2 = ω 2 ∊ μ ( 1 − j σ ∊ ω ) .$
$k 2 ( r , z ) = k 0 2 − k 0 k 2 r 2 .$
$E = ψ ( x , y , z ) · exp ( − j Φ ) ,$
$Φ ′ ( z ) = k 0 ( z ) ,$
$Δ t ψ − 2 j k 0 ψ ′ − j k 0 ′ ψ − k 0 k 2 r 2 ψ = 0 ,$
$Δ t = Δ − ∂ 2 ∂ z 2 = ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 .$
$ψ = exp − j ( p + 1 2 Q r 2 ) .$
$Q 2 r 2 + 2 j Q + k 0 ( 2 P ′ + Q ′ r 2 ) + j k 0 ′ + k 0 k 2 r 2 = 0 .$
$P ′ = − j ( Q k 0 + k 0 ′ 2 k 0 ) ,$
$Q 2 + k 0 Q ′ + k 0 k 2 = 0 .$
$Q = k 0 x ′ x .$
$( k 0 x ′ ) ′ + k 2 x = 0 .$
$Q m = − j k 0 k 2 .$
$x = a e j γ z + b e − j γ z ,$
$γ = k 2 k 0 .$
$Q 2 = Q m Q 1 + Q m + ( Q 1 − Q m ) e 2 j γ z Q 1 + Q m − ( Q 1 − Q m ) e 2 j γ z ,$
$k = β + j α ,$
$k 2 = β 2 + 2 j α β .$
$α = − σ 2 μ ∊ .$
$P 2 = P 1 : e 2 α z .$
$α = dB gain 8.686 ( z )$
$β = β 0 − 1 2 β 2 r 2 .$
$n = λ 2 π β = n 0 − 1 2 n 2 r 2 ,$
$α = α 0 − 1 2 α 2 r 2 .$
$k = k 0 − 1 2 k 2 r 2 .$
$k 0 = β 0 + j α 0 ,$
$k 2 = β 2 + j α 2 .$
$k 0 k 2 = β 0 β 2 + j ( α 0 β 2 + α 2 β 0 ) .$
$1 q = λ 2 π Q .$
$d q d z = 1 .$
$q 2 = q 1 + z ,$
$p 2 − p 1 = − j ∫ 0 z d z q 1 + z = j ln q 1 q 1 + z .$
$ln ( u + j υ ) = In u 2 + υ 2 + j tan − 1 υ u .$
$1 q 2 + n 0 ( 1 q ) ′ + n 0 n 2 = 0 .$
$1 q m = − j n 0 n 2 = − j λ π w m 2 .$
$1 R = n 0 w ′ w ,$
$n 0 w ( n 0 w ′ ) ′ = − n 0 n 2 + λ 2 π 2 w 4 .$
$k 0 k 2 = j 2 π λ α 2 .$
$1 q m = λ 2 π Q m = 1 2 ( 1 − j ) α 2 λ π .$
$R m = π w m 2 λ = 2 π α 2 λ ⋅$
$w m 2 = 2 λ π α 2 .$
$α = α 0 ( 1 − r 2 r 0 2 ) .$
$α 2 = 2 α 0 / r 0 2 ,$
$R m = r 0 2 π λ α , w m 2 = r 0 2 λ π α 0 .$
$α 0 = 11.5 m − 1 .$
$R m = 79 cm , w m = 0.94 mm .$
$p ′ = − j q m = − ( 1 + j ) z R m .$
$p 2 = p 1 − ( 1 + j ) 1 R m .$
$exp ( α 0 − 1 R m ) z .$
$γ = k 2 k 0 = j q m = 1 + j R m .$
$1 q 2 = 1 q m ( q 1 + q m ) − ( q 1 − q m ) exp 2 ( j − 1 ) z R m ( q 1 + q m ) + ( q 1 − q m ) exp 2 ( j − 1 ) z R m .$
$q 1 − q m q 1 + q m = exp 2 ( δ − j φ ) .$
$σ = 2 ( z R m − δ ) , ρ = 2 ( z R m − φ ) ,$
$1 q 2 = 1 ❘ q m e σ − e j ρ e σ + e j ρ = 1 q m sh σ − j sin ρ ch σ + cos ρ .$
$R 2 R m = c h σ + cos ρ s h σ − sin ρ .$
$w 2 2 w 2 2 = c h σ + cos ρ s h σ + sin ρ .$
$Y = j q = λ π w 2 + j R .$
$z = j q ,$
$d d z ( 1 Z L d V d z ) = Y L V ,$
$j Z L = k 0 ( z )$
$Y L = j k 2 ( z ) .$
$I ( z ) = x ′ Z L .$