English / Archive / FIRST ISSUE / akademik A. MARINČIĆ: Laser Beam Transformation through Space and Lenses

**Abstract**. *In transmission of laser beam through lossless space or through lenses or non-uniform media, the beam parameters change. If paraxial approximation is applicable, many problems in beam propagation through various media conditions can be reduced to closed form solutions. In this paper we consider transformation of circular and elliptical laser beam of Gaussian shape in passing through uniform medium or through lens or system of lenses.*

*Key words**: Laser circular and elliptical beams, lenses, beam parameters. *

Calculation of laser radiation field can be made using the aperture field method or using the spectrum of plane waves [1], or paraxial approximation of Helmholtz equation [2]. All these approaches lead to the same expression for the radiation field in the Fraunhofer region. In transmission through a homogeneous space, circular laser beam retains Gaussian shape regarding radial amplitude variation, but its phase front becomes curved and spherical at sufficiently far distance. Elliptical laser beam has to be considered independently along the main axes of the ellipse, and its amplitude field variation changes from elliptical to circular and then again to elliptical with the reversed axes. While the axial field of circular laser beam falls inversely with distance, it falls faster in the elliptical beam case.

Basic geometry used in calculation of circular laser beam is shown In Fig.1. In the aperture we assume the scalar field of the form

_{} (1)

where _{} is the axial field intensity at _{}, and _{} is the beam waist.

*Fig.1. Geometry used in defining laser beam*

Analytical expression for the beam at the plane_{}, can be expressed in the form

_{} (2)

where the beam radius and the radius of the phase front at _{} are

_{}, _{},

with _{}, _{}, _{} is the wavelength. At distances for which_{},

_{}, and_{}.

The far field is again Gaussian but with the complex parameters

_{} (3)

For a given laser _{}and_{} are constants. If we take into account that _{} we see that the far field beam is again circular and there is no_{}angle dependence. The axial field amplitude decreases inversely with the distance. At a fixed distance_{}, the contours of constant amplitude and phase are circles.

The transformation of Gaussian beam through a converging or diverging lens can be obtained by simple calculations based on the *q* parameter binomial transformation [3]. However, if the lens radius is comparable with the beam radius significant errors are encountered and we need to solve a complex electromagnetic problem. Here we will show steps to be taken if the lens radius is arbitrary.

_{} (4)

For the sake of simplicity we can now assume to have a lens of infinite size so that the aperture field behind the lens can be expressed as the Kirchhoff's-Huygens integral. If the size of lens is finite then the integral limits have to be reduced and no closed form solution is possible.

_{} (5)

where _{} and _{}. In the expression for _{}new axial coordinate is measured from the aperture to the observation point behind the lens. In this approximation we can put _{}in the denominator.

The integral (5) can be converted in the form suitable for integration. After some manipulations we reduce problem to solving the integral (shown only for the variable *x*’)

_{} (6)

where

_{}

After substitution _{}, the integral in Eq.(6) can be solved and the final result is

_{} (7)

In a similar way the integral with respect to *y*' leads to

_{} (8)

The radiated field is proportional to the product of (7) and (8) and it is also Gaussian. After further calculation the results for the field after an infinite lens is found to be identical with that obtained by the *q* parameters calculation. The effect of finite lens size can be found in ref. [1].

The electromagnetic field radiated from an aperture illuminated by the field:

_{} (9)

where _{} are the waists of the Gaussian elliptical beam and _{} is the electric field on the beam axis. The radiated field at some distance *z=z*_{1 }can be written in the form [5]:

_{}

_{} (10)

_{} (11)

For the sake of simplicity we can now assume to have a lens of infinite size so that the aperture field behind the lens can be expressed as:

_{}_{} (12)

where again we used _{} as the aperture coordinates, but they are different from those used in Eq. (9).

To find the radiation field we apply the Kirchhoff's-Huygens integral

_{} (13)

where _{}. Notice that the coordinate *z* is measured from the lens aperture so that we do not take into account the distance between the waists plane and the lens plane *z*_{1}.

Calculation of the double integral in (13) can be simplified if we assume that _{}in the denominator can be replaced by *z*, and in the exponential term by

_{} (14)

Now the radiation field can be calculated from two line integrals as the integrand in this case can be written in the form of product of functions _{}, where

_{} (15)

_{} (16)

We can convert the above expressions in the form suitable for integration. After some manipulations that lead to the integral (shown only for the *x*' )

_{} (17)

where

_{}

After substitution _{}, the integral in Eq.(17) can be solved and the final results is

_{} (18)

In a similar way the integral with respect to *y*' leads to

_{}_{} (19)

where

_{}

Finally

_{} (20)

The radiated field is also Gaussian and the 1/e half width at some distance z can be found after separating the real part of the term

_{} (21)

The required square of the 1/e half-widths are then

_{} (22)

_{} (23)

For *z*=0 we have _{} and _{} as it should be. The position of the beam waist after the lens is obtained from the derivative of (22) or (23). For the *x*-component it is

_{} (24)

After substitution of *w*_{z1x} and *R*_{z1x} from Ref.[5] into Eq.(24), identical expression as the one obtained by the *q*-parameter calculation is found:

For the *y*-component similar calculation leads to

_{} (26)

Upon substitution of the beam waist position from Eq.(24) into Eq.(22) it is found that the waist is

_{} (27)

for the amplitude at the axis one obtains which is identical to the one derived on the basis of *q*-parameter for the *x*- component. The y- component is obtained in a similar way and the only difference is the index in (27) “*x*” which has to be replaced by “*y*”. After further elaborate calculations we find

_{} (28)

where _{} and _{} are obtained from Eq.(22) and (23), respectively.

Lens systems consisting of more than one lens are used to modify the waist of a Gaussian beam. Systems consisting of two converging lenses or diverging and converging lenses are called *the beam expander*. The scheme of two converging lenses is shown in the Fig.2.

*Fig.2 Two lens system*

** **

In investigation of two or more lens systems, the matrix method of forming transfer matrix is advisable. Matrix of the system is formed by the matrix product of the individual sections as explained earlier. For example, the matrix of the two-lens system is obtained from the product of three matrices:

_{} (29)

where _{} and _{} are the focal lengths of the first and the second lens to which the beam is incident. Distance between the lenses is _{}

If the beam waist is at the distance _{}from the first lens then the *q _{ }*factors are:

_{}, at the beam waist cross-section plane;

_{}, at the cross-section plane in front of the first lens;

_{}, at the cross-section after the first lens;

_{}, at the cross-section of the beam waist after the first lens;

_{}, at the cross-section before the second lens, and _{}, at the cross-section of the output plane of the second lens.

The parameter _{} is obtained from

_{} (30)

_{Fig. 3 Beam transformation through a lens
}

The beam waist and its position after the first lens are obtained from _{}, where _{}. Hence,

_{}, _{} (31)

After substitution of defined parameters, we have also:

_{}, _{} (32)

To find beam waist and its position after the two lenses we have to use the matrix of the whole system as given by (1) so that

_{} (33)

To find the beam parameter _{} after the second lens we have simply:_{}. The new beam waist is obtained for a specific _{}, where _{}, i.e. the new beam waist after the system of two lenses is located at the plane where _{} is imaginary. Since _{}, the new beam waist is obtained from:

_{} (34)

_{} (35)

General solution of (34) and (35) is complex and will be given here in closed form only for effective magnification ratio:

_{} (36)

where_{}, is the maximum magnification, and _{}.

In a special case when _{} and _{} the solution for (34) and (35) is simple:

_{} , and _{} (37)

This case is referred to as the astronomical beam expander if_{}.

In transmission through various conditions in propagation, laser beam changes the beam waist and its position along the propagation axis. In paraxial cases or in Fraunhofer region closed form expressions for the electromagnetic field can be found. If infinite space beam distribution is not applicable, or if finite sizes of lenses are present, more precise approach is needed and the solutions can be found by numerical means.

- A. Marinčić,”
*Gausova funkcija u nekim problemima foton-elektronskih telekomunikacija*“, Glas CCCXCV Srpske akademije nauka i umetnosti, Odeljenje tehničkih nauka, knj. 34-2003. - M. Born, E. Wolf, “
*Principles of Optics*”, Pergamon Press, Oxford, Third Ed., 1965. - H. A. Haus, ”
*Waves and Fields in Optoelectronics*”, Prentice Hall, Oxford, New Jersey, 1984. - A. Marinčić, B. Milovanović, “
*The three-dimensional Gaussian beam field in the focal region*”, PES, 2005, Niš. - A. Marinčić,: ‘Huygens-Kirchhoff’s Theory in Calculation of Elliptical Gaussian Beam Propagation Through a Lens’, XXXVII Intern.Scientific Conference ICEST 2002, Oct.2002, Niš.