CC BY
38
LASER BEAM CHARACTERIZATION BY MEANS OF DIFFRACTIVE OPTICAL CORRELATION FILTERS
V.S. Pavelyev, V.A. Soifer, M. Duparre*, B. Luedge* Image Processing Systems Institute, Russian Academy of Sciences, Samara, Russia *Institute for Applied Optics, Univ. of Jena, D-07743 Jena, Germany
Analyzing of amplitude-phase characteristics of laser beam is topical in experimental physics and in a great number of laser applications, such as, for example, laser material treatment. The task of analyzing the amplitude-phase beam structure may be treated as that of analyzing the modal composition, if this is thought of as both analyzing individual modal powers and intermode phase shifts. In this paper the problem is tackled using a special diffractive optical element (DOE), called MODAN, matched to a group of laser radiation modes and their special combinations. The experimental results reported indicate that such an approach shows promise.
Key words: laser beam, Gaussian modes, intermode power distribution, intermode phase shifts.
1. Introduction
The analysis of amplitude-phase characteristics of the laser beam is a topic in experimental physics and in a great number of laser applications, such as laser material treatment [1,2,3]. In this paper the problem is tackled using a special diffractive optical element (DOE) -MODAN - matched to a group of transversal laser modes. Laser radiation modes are beams whose cross-section complex amplitude distributions are described by eigenfunctions of the operator of light propagation
[3]. Most familiar are the Bessel or free space modes
[4], and the Gaussian (Gauss-Hermite and Gauss-Laguerre) modes [3,5]. A fundamental property of the Gaussian modes is their ability to retain structure and mutual orthogonality while propagating in free space or lens-like medium [3,5]. The Gaussian modes describe the modes of an open passive resonator with spherical mirrors, and the modes of a gradient-index fiber with parabolic profile [3,5]. Coherent light modes may be interpreted as "generalized spectral components" since they propagate in their medium with modal eigenvalue and individual transverse amplitude-phase structure, acquiring only a phase shift [3] proportional to the optical path passed. Having in mind the above mode concept, we can represent any given amplitude-phase light distribution W (u) as a superposition of a number of modes
V (u) with definite complex-valued coefficients C n [3]:
arg C n, = arg C n,_! +A^n
(2)
W (u)=ÎC n^n (u),
(1)
where n = (n, m ), u = (u ,v ). Thus, analyzing the amplitude-phase structure of the beam reduces to analyzing the modal composition understood as measurement of
both individual modal powers, pn |2, and intermodal phase shifts, Apn = argC n _ arg C n-1. Then, based on
the measurements made for N
mod
modes we can re-
construct the initial amplitude-phase distribution taking any arbitrary constant value as phase p0 = arg C 0:
W (u) = P o|exp(i-Po \o (u) +
N mod -1
+ £ P n(u)exp[i(Apn, + afg C nM)]'
The accuracy of the estimate in Eq. (2) of the series in Eq. (1) have been discussed in [3]. We just note that in most practical tasks it is the maximal value of the mode Ncut with non-zero coefficient CN , ^ 0 in
cut
Eq. (1) which is known, so in practical applications
Nmod = f (Ncut). Using the measured individual
mode powers and intermodal phase shifts for modes ranging from fundamental one to that cut-off allows the amplitude-phase distribution in the laser beam cross-section to be reconstructed. Numerical methods discussed in [6,7] make it possible to reconstruct the power distribution between the modes, while the methods in Refs. [8,9] provide the direct reconstruction of the amplitude-phase distribution in the beam cross-section by numerical processing of the results of measuring the beam intensity in various planes. Note, however, that such an approach requires a high-precision measurement of intensity, analog-to-digital transformation and further numerical processing of intensity distributions or interferograms. Besides, the iterative phase reconstruction [8,9] from the results of intensity measurements is fundamentally unable to guarantee the convergence of the iterative procedure. In [3,10] it has been shown that the use of fundamental properties of laser light modes makes possible the synthesis of diffractive optical elements - MODANS - for analysis of modal composition. Rapid progress in technologies for DOE fabrication, including those for high-power lasers [11], open wide prospects to the use of such elements in the studies of the laser beam cross-section structure.
2. Basic formalism: DOE design for analysis of the beam modal spectrum
In [3] they report on the designing of phase mo-dans intended for the analysis of the coherent beam transversal modal spectrum. To change from the complex-valued transmission function (in the general case) to the purely-phase function both analytical coding methods [3] (with a carrier introduced into the element's phase) and iterative coding methods [9] have been developed.
n=0
T(u)
ra
W(u)
f
f
Fig. 1. Formulation of the problem: analysis of the beam modal composition
To this end, we need to seek for the DOE complex transmission function as a superposition of the complex transmission functions of "virtual" subapertures with corresponding weight coefficients.
"Virtual" apertures correspond to optical filters matched to the modes (or their combinations) and intended for deriving information about the modal powers, |Cn' |2 , and/or intermodal phase shifts Apn'.
Let us consider the calculation of the complex transmission function for such an element in more detail. In [3] they have shown that if a beam of cross-section complex-valued distribution W (u) interacts with a
filter of complex transmission function /*n (u) (where /n (u) is the mode numbered n, and * denotes complex
conjugation) at the Fourier-stage input, the resulting complex distribution at the center of the Fourier-stage output plane will be given by
w (0) = jjW (u)/n(u)d2u = Cn , (3)
D
where C n is the modal coefficient of the n-th mode in
the decomposition in Eq. (1) and D is the element aperture.
Measuring with a point-wise pick-off the intensity C n | at the output plane center allows one to derive information about the energy contribution of the mode number n into the beam W (u). To measure the N mod modal powers
simultaneously, in [3] they consider an element with the complex transmission function in the form
N mod */\ i \
T (u)= Z ZnW n'(u)exp(n,u).
(4)
n' =1
The exponential additions exp(ivn.u) are introduced to spatially separate channels corresponding to different modes in the element's plane, and the weight coefficients Zn ' > 0 are introduced to produce the desired energy redistribution between the subapertures (if needed). The interaction of the beam with a filter of complex transmission function in Eq. (4) at the input of a Fourier-stage of focal length f will result in the emergence of complex distributions w(/Vn') = Zn ' Cn ' in the output plane (Fig. 1). In [3], to measure the intermode phase shifts Apn ' = argCn ' - argCn ^ by modal optical correlation filters they introduce into the modan phase of extra virtual subapertures with the transmission functions given by
(5)
,kW n '-il
-i (u)+w*n ' (u)]exp (ivk u)
Zs W*n '-i (u)exp (iA< )+ y*n ' (u)-]exp(ivsu), (6) where A< is a certain phase shift introduced. It can be demonstrated [3] that the intensity values w k|2 = W (fk ) and w s|2 = W (fys ) are defined as follows (lets assume Z = 1,Z = 1 ):
Wk s w„
(7)
+ 2 • |wn'-l| • Iwn'|cos(A<n')
Iws|2 = Iwn'-l|2 + Iwn'|2 + 2 • Iwn'-l| •
• |wn'|[cos(A<n)cos A< - sin(A<n')sin A<]
(8)
Then, as it has been shown in [3], the values of the intermode phase shifts A<n< can be derived from Eqs. (7) and (8). Knowing the maximal admissible mode number in the beam under analysis we can write the relation for the complex transmission function of a filter intended for the analysis of the beam amplitude-phase composition:
N mod *
T (u)= Z Zn'W n (u)exp(ivu) +
n'=1
N
Z Zk [[**k (u) +W*k+1 (u)]exp(ivku) +
k=1
N „
i (u)]>
(9)
+ Z ~s W*s (u)exp(iA< )+W*s+1
s=1
x exp (ivs u)
Taking into account that certain mode(s) | (u) can have zero-valued coefficient C k = 0 , Eq.(9) should be changed to
N mod *
T (u)= Z Zn'W n'(u)exp(ivn'u) +
n'=1
N mod -1
Z Zk [[ (u)+W*k+1 (u)]exp(ivku)
k=0,k^r-1
Z Zs W*r (u)exp(iA< )+W*s+1 (u^]>
+
s=0,s^r-1
x exp (ivs u)
(10)
where r is the number of the "reference" mode, which has always nonzero valued coefficient C ^ 0 (for the
most of real laser resonators this can be the fundamental mode n=(0,0) ). Now, based on the measurements made
taken for N mod modes we can estimate the initial amplitude-phase distribution except a constant phase shift of arg C r = cpr :
W (u)sCr|exp(i•<|r(u) +
'|Wn'(u)exp [i(A<n' +<r )]'
n'=1,n'^r
A< = argC n ' - <r.
(11)
L
2
2
2
+ wn' +
+
+
3. Design and manufacturing of a diffractive optical element for mode power distribution measurement
The derived complex transmission function in Eq. (4) (which is certainly not phase-only in general case) can be reduced to the purely phase function using well-known encoding techniques [3] based on the introduction of a carrier into the phase function, for example, by generalization of well-known Kirk-Jones method [12] described in [3]. Thus, the modan for measurement of modal power distribution for 10 Gaussian-Laguerre modes was calculated by Eq. (4) and manufactured and investigated in the frame of present paper. The analyzing MODAN used for the investigations presented in this paper was designed to detect 10 Gauss-Laguerre(GH)-modes which are below the chosen cutoff condition 2n + \m\ < 3: (n,m)=(0,0), (1,0), (0,1),
(1,1), (0,2), (0,3), (0,-1), (1,-1), (0,-2), (0,-3). In accordance with our experimental conditions, the design value of the beam size parameter for these ten modes was selected as a= 0.8 mm. The element was calculated for a wavelength of 2 = 632.8 nm, and should work together with a Fourier lens of f = 452 mm. The corresponding „correlation points" for each mode, respectively, should develop in the focal plane of this lens in a distance of 1.5 mm from each other (see Fig. 3). The element was coded as a grating with 55.5 lines/mm, with a rectangular shaped carrier function „slowly" modulated across the aperture („synthetic coefficients' method"). The grating was calculated with a resolution of 2048 x 2048 pixels, with a pixel size of 3.0 ^m x 3.0 ^m. Figure 2 depicts the calculated phase distribution at the aperture of modan. Fig. 3 depicts the amplitude distribution in the used (zeroth) diffraction order of the element under the suggestion that it was illuminated by a pure (0,0) gaussian mode with a= 0.8 mm.
The calculated phase would have to be transferred into a corresponding surface depth, taking into consideration the intended transmission mode of the element, and the refractive index of the resist. In a next step, this two-dimensional „continuous" surface pixel map should be copied out into a corresponding electron dose map, which in its turn could have been written subsequently by the e-beam lithograph. However, thermal load resulting from such a one-shot-per-pixel procedure negatively influences final resist- and surface quality. As a consequence, the „continuous" surface pixel map was approximated by a staircase structure with 16 permitted levels in depth (=15 steps of equal height in the staircase). Each of these 15 steps had to be generated by an individual „etching" process. Therefore, 15 single binary masks had to be generated (as data fields) by software, starting from the „continuous" surface pixel map. Manufacturing started from a substrate from quartz glass, covered with a thin film of PMMA (thickness nearly 2.5^m), into which the spatially varying dose of electrons was implanted. The 15 dose levels each corresponding to one of the future new surface levels were realized by 15 times application of a binary electron beam writing process, using a commercial ZBA 23H
system. After finishing the development procedure of the PMMA film, the profile depth over all 15 steps should be 1200nm.
I > • 4 '
It ■»-*■' 1
Ill^y
■
Fig.2. Calculated phase distribution at the aperture of
modan
GL(0.0) GL(1,0) GL(0,1)
GL(0,-1) GL(1,-1) GL(0,-2| GL|0,-3) GL(l.l) GL(0,2) GL(0,3)
Fig. 3. Result of simulation of detection of the (0,0)-mode by the designed DOE
4. Optical set-up
The manufactured analyzing MODAN had to demonstrate its performance in a series of optical experiments. For this purpose, the set-up schematically shown in Fig. 4 has been used. The 10-channel modan (a=0.8 mm) has been illuminated by gaussian beam with a=0.525 mm (see Fig. 5,6). The decomposition of gaussian distribution with a=0.525 mm by using basis of Gaussian-Laguerre modes with a=0.8 mm giving us non-zero values of (0,0) mode coefficient and (1,0) mode coefficient (nearly 25% of (0,0) coefficient). Another modes with nonzero coefficients in this decomposition are of order higher then cutoff number was used.
He-Ne-Laser Modan Lens Camera
Fig.4. Optical set-up for modan investigation
Fig. 5. Result of simulation of illumination ofmodan by (0,0)-mode with o=0.525 mm
Fig. 6. Result of measurement of intensity in the output plane of modan illuminated by (0,0)-mode with <o=0.525 mm
Note good agreement between results of computational and natural experiments. Note also parasitical intensity peak in the central diffractive order (at the place of (1,-1) mode correlation peak) which is connected with technological problem - calculated surface profile was not perfectly achieved.
3. Conclusion
Designing a diffractive optical element intended for the analysis of the amplitude-phase composition of a laser beam has been reported in this paper. The results of computer simulation reported confirm the potentialities of the developed approach.
4. Acknowledgements
The authors would like to thank B. Kley from IAP for assistance in manufacturing the 3d resist profile.
References
1. Ward B.A., "In-process sampling and analysis of the far-field of a 21kW CO2 laser beam" in Laser Beam Characterization, P. M. Mejias, H. Weber, R. Martinez-Herrero and A. Gonzales-Urena, editors, (SEDO, Madrid, 1993), 53-64.
2. Teppo E. A., "Diagnostic tools for laser beam characterization" in Laser Beam Characterization, P. M. Mejias, H. Weber, R. Martinez-Herrero and A. Gonzales-Urena, editors, (SEDO, Madrid, 1993), 23-30.
3. Soifer V.A., Golub M.A., Laser Beam Mode Selection by Computer Generated Holograms .- 1994.-CRC Press. 215 p.
4. Durnin J., Miceli J.J., Eberly Jr. And J.H., "Diffraction-Free Beams", Physical Review Letters, 58 (15), 1499-1501.
5. Yariv A., Optical electronics, Holt, Rinehart, and Winston, New York, 1985.
6. Cutolo A., Pierri T., Zeni L., "Measurement of the mode content of a laser beam" in Laser Beam Characterization, P. M. Mejias, H. Weber, R. Martinez-Herrero and A. Gonzales-Urena, editors, (SEDO, Madrid, 1993), 263-273.
7. Du K.-M., Herziger G., Loosen P., Ruhl F., "Measurement of the mode coherence coefficients", Optical and Quantum Electronics, 24, S1119- S1127, (1992).
8. Stark H., Image Recovery: Theory and Application, Academic Press Inc., 1987.
9. Soifer V.A., Kotlyar V.V., Doscolovich L.L., Iterative Methods for Diffractive Optical Elements Computation, Taylor&Francis Ltd., 1997
10. Golub M.A., Prokhorov A.M., Sisakian I.N., Soifer V.A., "Synthesis of spatial filters for investigation of transversal modal content of coherent radiation", Quantum Electronics, 9, N 9, 1866-1868, (1982).
11. Kononenko V.V., Konov V. I., Pimenov S. M., Prokhorov A. M., Pavelyev V. S., Soifer V. A., „CVD diamond transmissive diffractive optics for CO2 lasers", New Diamond and Frontier Carbon Technology, 10, N 2, (2000) (accepted for publication).
12. Kirk J.P., Jones A.L., "Phase-only complex valued spatial filter", JOSA, 61, N 8 , 1023-1028, (1971).
13. Duparre M., Pavelyev V.S., Ludge B., Soifer V.A., Kowarschik R., "Laser beam characterization by means of optical correlation filters", Proceedings of EOS Topical Meeting in Diffractive Optics DO99, Jena , Germany 23-25 August 1999, 284-285.
