CC BY
26
M.A.Golub, L.L.Doskolovich, V. V.Kotlyar, I. V.Nikolskv. V.A.Soifer
ITERATIVE-PHASE METHOD FOR DIFFRACTIVELY LEVELLING THE GAUSS BEAM
INTENSITY
1. Introduction
The demand for optical elements capable of levelling the Gaussian intensity distribution exists in such problems of optical data processing as laser-based superficial strengthening and laser projection printing. The diffraction redestribution of the light beam intensity can be performed by means of so called focusators of laser radiation (1,2|. Methods of geometrical optics make it possible to derive 12,31 analytical relationships for a focusator phase which is the smooth function of two variables. However, such an approach ignores diffraction effects that can result in low accurace in formation of the required focusator-produced intensity distribution. On the other hand, in computing a phase of kinoforms, an iterative algorithm [4-61 is employed, which takes into account diffraction effects taking place under the light propagation. Making use of the adaptive correction procedure 17,81 allows this algorithm to be successfully applied to compute the phase of focusators.
This paper deals with numerical comparison of performances of focusators from the Gaussian beam into the uniform intensity square computed by geometrical-optical with those computed by iterative methods.
2. Focusing from the Gaussian beam into a rectangle
As it has been shown |9], for the focusator that form the rectangular intensity-uniform distribution as
1,
0, |x|>dv \y\>d2
from the plane beam with the Gaussian profile of the intensity distribution
l0(u,v)-l0ex^-u2 + v*
I a2
we can (given a lens located immediately adjacent to the focusator) find its phase function in the form of the sum
<t>(iy.^)=<|»1(i/)-<j>2(v) ,
with its terms satisfying the set of two simultaneous equations
dx_W du /(x)
x,u+flr*d№ du
where/is the focal length of the lens in whose focal plane the light rectangle is formed, k is the wavenumber of light, (w,v) and (.v,y) are the coordinates in the planes of the focusator and of the Fourier spectrum, respectively, 2d„ and 2d2are the measures of the rectangle and a is the parameter of the Gaussian beam.
The first equation in the system (1) sets the equality of the density of light energy of a 1D focusator to the density in the corresponding areas of the straight-light segment. The second equation in the system (1) describes the stationary points in the paraxial approximation.
The discrete variant of the solution of the system (1) can be written as follows
Mx{6n/V-1erf(6nA/"V
(2)
<JW>n(1 0)
+7T Jexp(-36n2/V *)UM;6mA/ \
xerf(6/77A/->7r 2exp(-36mzA/ z)
where erf(x)=2n 1 |exp(-f2)df ,
w=0,± 1 ,±2,...,±°'V/2, «=0,±1,±2.....±N/2,
M^Ntffiinoy, My=TtN(6v/InT0)-\ Mx and Mv are the number of the minimum diffraction spots that fall into the rectangle with the measures NxxNr In this case, the square aperture of the focusator measured NxN equals to the measure of the square, 6a. and the Gaussian collimated beam that illuminates the focusator produces the light distribution as
/0w,=exp[-36/V>2+m2)] (3)
In section 4 we employ the phase (2) to numerically simulate the operation of a geometrical-optical focusator.
3. Adaptive-iterative computation of focusators
Following 17,81, let us briefly describe how the iterative algorithms can apply to computing the
focusators as kinoforms. Let us proceed from the preset complex amplitude A{u) of the illuminating beam and the required intensity distribution I(x) in a focal plane of the lens. The complex amplitude F(x) in a focal plane is related to the complex amplitude
/'(a)=a(i;)exp[/<t)(iy)]
immediately behind the focusator through the Fourier transform
b
Psx)=j f(u)exp(ikxu/f)dx, -b
where 2b is the size of the focusators aperture. The focusator's equation takes the form
№)|2=/(x) <4>
To find iterative solution of Eq.(4) with respect to the phase <p(u), one should perform some preliminary estimate of the phase 4>0(u) followed by the computation of the complex light amplitude in a focal plane. In this case the complex amplitude Fn(x) calculated in the n-th step of iterations is replaced by the function F°(x) according to the rule
F„(x) |Fn(x) |, |*| id
(5)
**Fn{x)
\x\>d
where In(x)=( l + a)I(x)-a J/n(.*) |2, /(a) is the required intensity distribution within the interval \-d,d\ of a focal plane, a is the parameter that controls the rate of convergence of the calculated intensity to the required one. For a=0, the replacement (5) changes to the standard replacement in the Gerchberg-Saxton algorithm [5|.
The amplitude of light in the plane of focusator./„(//) is calculated with the help of the inverse Fourier transform and is replaced by the function f„°(x) according to the rule
f°(x)-l A(u)fn(umu)r , \u\sb (6) n( ) { 0 , \u\>b .
In contrast to the algorithms of the conditional gradient |10|, the a parameter is here introduced directly into the intensity function.
The rate of convergence of the intensity |Fn(x)to the required one /(x) is checked by the root-mean-square deviation
' fm-\F„(x)l']2à*V * ■d
0 =
fl2(x)dx
Besides, a new parameter e characterizing energy efficiency of focusing is introduced:
(7)
the
f\Fn(x)\2dx
e
(8)
Fn(x)\2dx
In the next section we apply the phase derived from Eqs. (4)-(6) to numerically simulating the operation of diffractive focusators.
4. Numerical results
We examine the focusator into a square that comprises 32x32 pixels and measures 10 minimum diffraction spots. For a focusator that focuses from the Gaussian collimated beam into a square, the phase deduced from Eq.(2) on the net of pixels 256x256 and taken to the modulus 2n represents the set of rings (lines of equal phase) changing to the lines of the square perimeter (F-ig.l). Figure 2 illustrates the light intensity distribution in the lens focal plane calculated as the Fourier transform of the amplitude
Fig. I Geometrical -opticalphase eg the focusator into a square
fig.2 The intensity distribution obtained Iron the geometrical-where /0mn is taken «P'^l focusator in a leas focal
from Eq. (3) and <t>nin
from Eq. (2). The root-mean-square deviation of the obtained distribution from the uniform one amounted to 5% and the efficiency was 91.6%.
To calculate the focusator that focuses the Gaussian beam into the uniform intensity square, as an initial phase estimate we have chosen the geometrical-optical phase function described above (Fig.l). The subsequent iterative calculation has been conducted using three techniques. The first approach to the calculation of the phase based on the standard variant of the Gerchberg-Saxton algorithm with the replacement (5). for ar^O. yields, at first, the increase of the error <5 in the first iteration steps and then it gives the slow decrease of the error (Fig.3a, curve 1). In this case during 10 iterations the error does not diminish below 13% though the efficiency e grows up to 98.9% (Fig.3b, curve 1).
Fig. J TIk root-mean-squarc deviation (a) at the energy c 111 -ciency (b) agains the number of iterations for different variants of an iterative algorithm Gerchberg-Saxton (1), combined (2). and adaptive (3).
The second approach is combined: the first three iterations are conducted with the replacement (5) for a=0, the remaining seven iterations for <r=l (Fig.3a, curve 2). This method after 10 iterations results in formation of the square characterized by the uniform intensity, with the error of 2.4% and the efficiency of 96.4% (Fig.3b. curve 2). This way is seen to be more effective as compared with the first approach, since it yields the decrease of the error from 13% to 2% (by 6-fold) without essential decrease in the efficiency.
The third method is purely adaptive, this means that the replacement (5), with df— 1, is performed in each step of iteration (Fig 3a. curve 3). As one can see in this case, the error decreases monotonously and during 10 iterations it becomes equal to 01%, but the efficiency falls to 92.2% (Fig.3b, curve 3). Figure 4 illustrates the phase of the focusator to the modulus 2* that has been calculated during 10 iterations on the basis of the third approach from the geometrical-optical phase (2). First, one can see that the use of the iterative algorithm does not result in the essential change of the initial phase The pronounced changes of the phase
Fig 4 The phase of the focusator into a square obtained after 10 iterations from the initial geo-metro-optical phase
Fig.5 The radius-random phase chosen as the initial estimate for iterations.
Fig.6 The phase of the focusator into a square obtained alter 10 iterations from the initial random phase
(Fig. I ) take place only on the edges of the focusator (Fig.4). Second, we can draw the conclusion that the adaptive iterative algorithm based on the replacement (5) for a=\ makes it possible to improve the geome-trical-optical phase in such a manner that the error in the formation of a light square reduces more than by an order, with the efficiency almost unchanged.
For the aim of f^ comparison, there was conducted another numerical simulation where we made use a radius-random phase function (Fig.5) as the initial approximation. The phase shown in Fig.6 was calculated during 10 iterations. The intensity distribution formed by the focusator characterized by such a phase is presented in Fig.7. The error and the efficiency were equal (5=6.4% and 6=80.8%, respectively. Further iterations did not produce the essential change in these values. Data given above show that the iterative algorithm that begins with the random phase estimate, leads to the non-regular structure of the focusator's zones and results in the accuracy and efficiency which are somewhat less than those for a geometrical-optical focusator.
The table summarizes all the focusators discussed in this paper. One can see the advantages in terms of the uniformity, the efficiency, and the regular phase structure that we achieve when using physically warranted initial approximation in the form of the geometrical-optical phase (compare rows 1-4 with row 5). The comparison between rows 1 and 2 in the table shows that the Gerchberg--Saxton method yields the increase of 7% in the energy efficiency of the geometrical-optical approximation (row I) but produces considerable non-
Fig.7 The intensity distribution in a lens focus obtained from the focusator with a phase shown in Fig 6
Table
Type of a focusator and a phase The energy efl-ici-cncy f. % ' The relative Г 111 s. error 5. %
(ieometiical-optical phase 91.6 5.0
10 iterations of the Gerchberg-Saxion algorithm on the geometrical-optical phase 98.9 13.0
10 iterations of the adaptive-iterative algorithm on the gcomctrical-optical phase 92.2 0.1
10 combined iterations on the geometrical-optical phase 96.4 2 4
10 iterations on the radius-random phase 80.8 6 4
uniformity of the intensity over the square (13% instead of 5%).
On the contrary, the adaptive method ensures the high degree of the intensity uniformity (0.1%) but it practically does not improve the efficiency as compared with the geometrical-optical approximation (row 3). The combined method enables us to obtain rather high efficiency in combination with the quite satisfactory intensity uniformity (row 4 in the table). Row 5 shows that not knowing the geometrical-optical solution but based upon the random initial phase in the estimate of the focusator's phase, and using the iterative method we can achieve suitable results in terms of the accuracy and the efficiency in the formation of the required intensity distribution.
5. Conclusions
In the present work the authors have numerically shown that an optimal approach to a problem of computation of phase optical elements focusing the laser light into the small areas of the spatial spectrum plane consists in the solving of the inverse geometrical-optical task and finding of the phase function which is then chosen as the initial approximation for the adaptive-iterative procedure of obtaining of tn final phase In this case the iterative procedure of the correction for the initial phase with the regular zone structure does not result in its essential change.
References
|1| Golub, M.A.. Karpeyev, S.V., Prokhorov, IN.. Sisakyan. I.N.. and Soifer, V.A., 1981, Lett. J. Techn. Phys.. 7, 618. |2| Golub, M.A., Sisakyan, I N., and Soifer V.A.. 1991, Opt. & Lasers in Eng.. 15, 297. 131 Han, C.--Y., Ishii. Y.. and Murata. K..
1983, Appl. Opt., 22, 3644. |4| Hirsch, P.M., Jordan. J A., and Lesem.
L.B., 1971, Patent US 3619022. [5] Gerhberg, R.W., and Saxton. W.D., 1972,
Optik, 35, 237. |6| Fienup. J.R., 1980. Opt. Eng., 19, 297. |7| Kotlyar, V.V., Nikolsky, I.V., and Soifer,
V.A., 1991, Optik. 88. 17. |SI Kotlyar, V.V., and Nikolsky, I.V., 1991, Opt
«Sc Lasers in Eng., 15, 323. |9| Golub. M.A., Doskolovich, L.L., Kazansky. N.L., Sisakyan, I N., Soifer, V.A., and Kharitonov, S.I., 1990, Computer Optics (Moscow), 7, 42. [I0( Vorontzov, M.A., and Shmalgauzen, V.l., 1985, Principles of Adaptive Optics (Moscow: Science).
Вниманию читателей
Подготавливается к изданию сборник
Компьютерная оптика
выпуск 14
телефон для справок: 198-72-10
