Strojniški vestnik - Journal of Mechanical Engineering 57(2011)1, 3-10 D01:10.5545/sv-jme.2010.120 Paper received: 27.05.2010 Paper accepted: 18.10.2010 Rotating-Mirror Q-Switched Er:YAG Laser for Optodynamic Studies Aleš Bizjak1* - Karolj Nemeš2 - Janez Možina3 1 I.H.S. d.o.o., Slovenia 2 Fotona d.d., Slovenia 3 University of Ljubljana, Faculty of Mechanical Engineering, Slovenia In this article a development of a simple, reliable and repeatable Q-switched Er:YAG laser with a pulse duration of 270 ns and energies up to 16 mJ in a TEM00 spatial beam mode is presented. The Q-switching is performed mechanically with a flat, highly reflective silicon rotating mirror. The laser pumping chamber is cooled with thermally stabilized de-ionized water. This laser is intended for optodynamic studies of the short Er:YAG laser pulses with various media, primarily those containing water, where the thermal diffusion effects are negligible. ©2011 Journal of Mechanical Engineering. All rights reserved. Keywords: Er:YAG laser, Q-switch, inverse population, optodynamics, Xe flash lamp, rotating mirror 0 INTRODUCTION The Er:YAG laser is a well-accepted tool in many medical applications, like the ablation of hard and soft dental tissue as well as in various areas of dermatology [1] to [3]. For all these applications the Er:YAG laser operates in a free-running mode with laser-pulse durations from several tens of ^s up to several ms. The interaction is based on the coincidence of the Er:YAG laser wavelength (2940 nm) with the absorption peak of water (absorption coefficient [ia = 13,300 cm-1 [4] and [5]. This high absorption leads to a very shallow (several ^m) penetration depth in the water as well as in all water-rich biological tissue. During the interaction of the laser pulses with the biological tissue, part of the laser-pulse energy is used up for ablation (depending on the instantaneous laser pulse power and the fluence), while the remaining part of the pulse energy during the laser pulse leaves the interaction region with diffusion that has a characteristic diffusion length of (4DT)1/2 (D is the diffusivity of the media and T is the laser-pulse duration). Recently, interactions in which the diffusion of heat from the interaction region during the laser pulse can be neglected became of some interest. Any study of this type of interaction obviously needs a very-short-pulse Er:YAG laser. Approaches for the reliable, Q-switched operation of an Er:YAG laser have already been described in the literature [6] to [8]. Most of these approaches need special materials and expensive components like undoped Nd:YAG crystals cut and polished at precise angles or specially designed electro-optic Q-Switches based on LiNbO3 crystals and high-voltage Q-switch drivers. The repeatability and long lifetime, as well as the easy set-up of the experimental laboratory, are important issues. Since the rotating-mirror and Q-switched lasers (Nd:YAG, Er:glass, and others) are already known for their reliability and repeatability, we decided to use this approach for the Q-switched Er:YAG laser as well. The Er:YAG laser is a complex, rare-earth, solid-state crystal laser in which 50% of the yttrium ions are replaced with Er3+ ions. The laser is exited with Xe flash lamps [9]. Since the laser operates at fairly high pumping energies, efficient cooling of both, the laser rod and the flash lamp is accomplished with de-ionized water. The Er:YAG laser is one of the rare lasers that operate at a wavelength near 3 ^m at room temperature. The basic Q-switching technique for the rotating-mirror laser is well known and described elsewhere [10], and only the main principles are shown in Fig. 6. One of the resonator mirrors rotates at a constant frequency m and passes the resonant alignment (with low losses) of the resonator once in a single rotation. The laser rod is excited with a Xe flash map until a steady inversion population n is built up during the non-resonant condition of the laser. When the rotating mirror reaches the resonant position at time t0, the resonator losses abruptly fall to a level where laser oscillation starts and a fast build-up of photons from the inverse population creates a Q-switched laser pulse (Fig. 1 photon flux 0). Fig. 1. Dynamics of the inverse population n and the photon flux according to the resonator losses following the time t0 In an analysis of the Q-switched laser it is very important to understand the laser crystal population dynamics, Q-switch operation itself as well as the laser resonator structure. The Er:YAG crystal's laser energy levels are shown in Fig. 2. The laser operates at room temperature, and the laser emission that we are interested in develops between the states I11/2 and I13/2 at a wavelength of 2940 nm. I11/2 is the laser upper level, while I13/2 is the laser lower level. The excitation with the Xe flash lamp directly pumps the levels 4S3/2 with a pumping rate r1 (cm-3s-1) as well as the upper laser level I11/2 with the rate r2 (cm-3s-1). The highly lying levels 4S3/2 are very quickly depleted with transitions to the upper laser level and therefore contribute to the creation of an inverse population. The up-conversion w2 from the upper laser level pumps the 4S3/2 level. The up-conversion w1 depletes the lower laser level in such a manner that it contributes to the inverse population twice. Both, the upper and lower laser levels are Stark split. The population distribution in the Stark split levels follows the Boltzmann distribution - i.e., the populations are determined by Boltzmann factors [11]. It is important to emphasize that both pumping rates r1 and r2 are tied with r = r1 + r2 . In this article we assumed the relationship between the rates r1 and r2 to be k = r2 / (r 1 + r2) [12]. Fig. 2. Simplified energy diagram with the most important transitions in an Er:YAG (50% at Er3+) laser crystal Fig. 3. a) 50% dopedEr:YAG crystal absorption spectrum; b) Xe flash lamp spectrum The Er:YAG laser is efficiently pumped by flash lamps filled with Xe since the pumping bands for the Er:YAG crystal are around 500, 650 and 750 to 800 nm and they overlap very well with the Xe flash lamp's spectrum (Fig. 3). 1 POPULATION DYNAMICS OF AN Er:YAG LASER The time evolution of pumping, the population of various levels as well as stimulated emission is described with rate equations. Due to the many processes that are present during the Er:YAG laser pumping these equations are fairly complex. For studying the population dynamics of the upper and lower laser levels of a flash-lamp-pumped Er:YAG laser as well as for the generation of laser emission, these equations can be simplified [12], and can be described using the equations below: dn3 , , n3 2 ~t = i(0 - +w2 n , dt t3 dn2 . . n3 2 2 —2 = r2(t) + — + wlnl -2■ w2n2 - dt t3 n —- - s2l (b2n2 - a7nj)0 dn n 2 n2 n It - 'wn +7"+ +s2l(b2n2 -a7nj)0 d® n ^ ° -= s21 (b2n2 - a7n )®--, dt t„ (1) (2) (3) (4) where is: n1 lower laser level population [cm-3] n2 upper laser level population [cm-3] n3 population of 4S3/2 state [cm-3] 0 photon flux in the resonator [cm-2] r t(t) pumping rate of 4S3/2 [cm-3s-1] r2(t) direct pumping rate of upper laser level [cm-3s-1] s21 cross-section for stimulated emission [cm2] b2 Boltzmann factor of Stark spited I11/2 upper laser level a7 Boltzmann factor of Stark spited I13/2 lower laser level t1 lower laser level lifetime [s] *2 t3 w2 tc upper laser level lifetime [s] 4S3/2 lifetime [s] coefficient of up-conversion process of upper laser level [cm3s-1] coefficient of up-conversion of lower laser level [cm3s-1] photon lifetime in laser resonator [s]. The inversion population that is important for the laser emission is: n = (b2n2 - a7ni\ (5) here, the degeneracy due to the Stark splitting of the upper and lower laser levels is taken into account. The pumping rates' time behaviors r1(t) and r2(t) are used in these equations as r1(t) = r1P(t) and r2(t) = r2P(t), where P(t) is the shape of the flash-lamp pumping pulse with an amplitude of 1. The relation between the pumping rates is a constant throughout, taken as k = r2/(r1+r2) = 0.6. In order to check the above rate equations, they were numerically solved for a free-running Er:YAG laser (no Q-switching) for a rectangular flash-lamp current pulse of 300 ^s duration, where for reasons of simplicity we used P(t) = 1-h(t-tp), where h(t) is a Heaviside step function and tp is the duration of the rectangular current pulse. The constants used in Eqs. (1) to (4) were taken from reference [12]. The photon lifetimes were calculated for a resonator with plane/plane mirrors, an out coupling mirror transmitivity of 85% at 2940 nm, a $ 4 x 90 mm laser rod and an overall resonator length of 178 mm. Typical results of the solutions are shown in Fig. 4. Since the laser in free-running mode operates in mixed-mode operation, without limiting apertures inside the resonator, except for the ones at the laser-rod ends, we can obtain the laser efficiency curve by determining the time integral of the photon flux (proportional to the laser pulse energy) as well as the pumping rate (proportional to the flash-lamp pulse energy) for various pumping rates. The curve in Fig. 5 predicts a linear relationship between the laser output energy on the flash-lamp pulse energy that is a typical result, obtained in all experiments for an Er:YAG laser operating in the free-running mode. Therefore, the 2 0 10 20 30 44 50 60 timers] Fig 4. a) lower laser level population - n1; b) upper laser level population - n2; c) inverse population n = (b2n2 - a7n1); d) photon flux in resonator - p; e) photon flux in resonator with higher temporal resolution (note the relaxation oscillation micro pulses with a time duration of the order of 1 /is, which are also observed experimentally) rate-equation model describes the behavior of a free-running Er:YAG laser well. The curve in Fig. 5 predicts a linear relationship between the laser output energy on the flash-lamp pulse energy that is a typical result, obtained in all experiments for an Er:YAG laser operating in the free-running mode. Therefore, the rate-equation model describes the behavior of a free-running Er:YAG laser well. 2 Q-SWITCHED OPERATION The simulation of the Q-switched operation was performed in two steps: Fig 5. Er:YAG laser output energy as a function of the flash-lamp input energy a) From Eqs. (1) and (2), as well as from Eqs. (3) and (4), the inverse population as a function of time is determined, where we assume no laser action (0 = 0) is assumed. Here we assume a rectangular time profile of the flash-lamp current and hence the pumping terms r1(t) and r2(t). The inverse population is calculated until time equals t0. b) From time t0 on, the rate Eqs. (1) to (3) and (6) are used, which were modified from Eq. (4) and presented below: d O d. O( L + s(t )) -= s21 • n —L O----— dt I t„ (6) Then, these equations, with the starting conditions from Eqs. (1) to (4) at time t0, are solved numerically. In Eq. 6 for a photon flux the term for the stimulated emission is multiplied by the ratio of the length of the active medium d2 to the overall resonator length lres. The photon lifetime tc is replaced with the resonator roundtrip time tR presented in Eq. 7: tR = 2 ' (d1 + nrodd2 + (7) where d1 is the separation between the out coupling mirror and the laser rod, nrod is the rod index of refraction, d3 is the separation from the rod to the rear mirror center. The losses in the resonator are expressed as: L = Ls + Ln (8) where Ls are the diffraction losses and Loc are the losses from the out coupling mirror with scattering and absorption. Ls was estimated to be Ls = 0.05, while Loc was taken as: L =-ln(R ), oc v oc / (9) where Roc is the reflectivity of the out coupling mirror. The term s(t) is the time-dependent losses of the rotating-mirror Q-switch. s(t), as a function of the rotating mirror's angular frequency m, was modeled as [10]: 3 ROTATING-MIRROR Q-SWITCHED LASER DESIGN The Q-switch in the form of a flat Si mirror (silver coated by II-VI Belgium) with a reflectivity of 99.7% at 2940 nm was fixed to a motor shaft. The position of the motor shaft at a defined angle a to the resonant position of the mirror was detected for triggering the flash-lamp discharge by optocouplers. The frequency of the rotating mirror was adjusted to 200 Hz. The Er:YAG laser-pumping chamber consisted of a BaS04 diffuse reflector with a Borofloat 33 flow tube. A 4 x 90 mm Er:YAG flat/flat rod (Sumitomo) doped with 50% at Er+3, AR coated at 2940 nm on both end faces, with fine ground sides, was used. The laser rod and the flash lamp were cooled with temperature-stabilized, de-ionized water (Twater = 18 ± 0.5 °C). The laser was pumped with a Xe flash lamp that had an inner diameter of 5 mm and an arc length of 89 mm with a cerium-doped envelope, made by Noblelight. The out coupling mirror was a flat/ flat mirror with a reflectivity of 85% at 2940 nm and an AR coating for 2940 nm manufactured by Layertec GmbH, Germany. The laser resonator is shown in Fig. 6. s(t) = a + b cos(©t) . (10) Fig. 6. Q-switchedEr:YAG laser resonator, where d1 = 52 mm, d2 =90 mm and d3 =36 mm The flash lamp was operated with a simmer current of approximately 150 mA, while the main discharge (pumping pulse) was generated from a 10 mF capacitor bank with an IGBT module. The IGBT module was activated with a programmable pulse generator, triggered by the rotating-mirror, Q-switch, position-sensor electronics. The repetition rate of the Q-switched laser pulses was set to 2.87 Hz. The flash-lamp triggering circuit as well as the typical flash-lamp current and voltage pulses are shown in Fig. 7. During the experiments for each laser pulse, the flash-lamp current and c the voltage pulse were measured and used for the calculation of energy delivered in the flash lamp. 3.1 Simulations of Q-Switched Operation For reasons of simplicity, an ideal rectangular flash-lamp pumping pulse with a duration of tp = 325 ^s was used in the simulations. The real pumping pulse is not rectangular due to the time needed for the flash-lamp Xe plasma to fill the flash lamp. Since we were only interested in the creation of a certain inverse population, the shape of the pumping pulse did not influence the results of the simulation. The values of the populations of all three levels n1, n2 and n3 were solved at the time t0 with a photon flux 0 = 0. For times t > t0 the photon flux was solved from Eqs. (1), (2), (3) and (6), where the time-dependent term for Q-switch losses e(t) was used. The simulations were made considering a Q-switch rotation frequency of 200 Hz. See Fig. 8a. Time [11s Fig. 7. a) flash-lamp triggering and main discharge control circuit and; b) typical flash-lamp voltage and current pulse shapes Fig. 8. Q-switched Er:YAG TEM00 laser pulse temporal shape; a) calculated from rate Eqs. (1) to (3) and (9); b) experimentally realized and measured 3.2 Measured Properties of Q-Switched Er:YAG Laser The Q-switched laser pulse shape, energy, efficiency, output-energy repeatability and laser-beam spatial profile were measured. The measured Q-switched pulse temporal shape is shown in Fig. 8b, together with the same laser pulse calculated from the rate equations (Fig. 8a). The temporal pulse shape was measured with a fast InAs detector (Teledyne Judson Technologies). Table 1. Calculated pulse-to-pulse repeatability Fig. 9. Q-switched Er:YAG laser pulse output energy as a function offlash-lamp input energy; the threshold energy for Q-switched operation was found to be 40.2 J The Q-switched laser-pulse temporal shape remained the same for input energies in the flash lamp of Ein = 52 J, with a pulse width (FWHM) of 270 ns. At higher pumping energies, the evolution of the second pulse was observed (typical in rotating-mirror Q-switching [10]) as well as the evolution of higher transverse modes, since no mode selection in the form of additional limiting apertures was used in the laser design. The laser output energy for single-mode, Q-switched operation as a function of the flash-lamp input energy is shown in Fig. 9. The repeatability of the laser pulses was measured at different laser-pulse energies (see Table 1) for a pulse-repetition frequency of 2.87 Hz. The repeatability is the worst at low energies (approximately 7%). For laser-pulse energies greater than 10 mJ the repeatability saturates to a remarkable value of 2%. Eflash [J] ^laar [mJ] St. dev. [mJ] Pulse-to-pulse repeatability [%] 45.9 7.1 0.5 6.9 47.3 9.6 0.4 4.6 48.6 11.4 0.2 2.1 50.8 13.8 0.3 2.1 52.1 15.9 0.3 1.9 The laser-beam spatial profile was measured at a position 341 mm from the laseroutput coupler. Due to a lack of equipment for better measurements, a technique of measuring the transmission through the apertures of various diameters was used and then the fluence as a function of aperture diameter and finally fitting the fluence to the Gaussian beam profile was calculated. The result is shown in Fig. 10. Fig. 10. Spatial profile of Er:YAG Q-switched laser beam at a distance of 341 mm from the laser out coupling mirror; the data points are from measured transmissions through apertures with various diameters; the laser-pulse energy was 15 mJ at repetition rate of 2.87 Hz; the measured data points were fitted to a Gaussian beam profile (full line); the 1/e-2points are at 2w =1.28 mm The measured spatial profile corresponds to a TEM00 mode with the 1/e-2 point at a radius of 0.64 mm. A higher certainty could be reached with measurements made at multiple distances from the laser out coupling mirror. The propagation of the fundamental laser mode created in the resonator was also calculated. Here, the thermal lensing of the laser rod under the conditions of the experiment (laser rod pumping power = 52.1 J x 3 Hz = 156.3 W) was taken into account. The obtained result for the 1/e-2 point of the spatial profile was 0.85 mm. For the reasons for the discrepancy between the measured and calculated value of the 1/ e-2 point in the spatial TEM00 beam profile may lie in the inaccuracy of the method of the beam-profile measurement used. The other possible reason could be that in the calculation of the TEM00 mode and its propagation we took an average value of thermal lensing of the Er:YAG rod. The actual thermal lensing of the laser rod during the creation of the Q-switched laser pulse could be different than the average value we used. 4 CONCLUSIONS We developed a simple, repeatable and reliable Q-switched Er:YAG laser for laboratory use. Q-switched, laser-pulse energies from 7 to 16 mJ with a pulse width of 270 ns and a TEM00 spatial beam profile were realized at a repetition rate of 2.87 Hz. The pulse-to-pulse repeatability reached values of approximately 2%. The laser is intended for various optodynamic laboratory studies of 2940 nm, Q-switched, laser pulses' interaction with various media that contain water. In the future, we intend to run experiments to detect and define the thickness of the water condensing on various surfaces, like aluminum, copper and glass. 5 REFERENCES [1] Grad, L., Možina, J., Šušterčič, D., Funduk, N., Skalerič, U., Lukač, M., Cenčič, S., Nemeš, K. (1994). Optoacoustic studies of Er:YAG laser ablation in hard dental tissue. SPIEProc, vol. 2128, p. 456-465. [2] Jelinkova, H., Dostalova, T., Hamal, K., Krejsa, O., Kubelka, J., Prochazka, S. (1998). Er:YAG Laser in dentistry. Laser Physics, vol. 8, no. 1, p. 176-181. [3] Paghdiwala, A.F., Vaidyanathan, T.K., Paghdiwala, M.F. (1993). Evaluation of erbium:YAG laser radiation of hard dental tissues: analysis of temperature changes, depth of cuts and structural effects. Scanning Microsc., vol. 7, no. 3, p.989-997. [4] Shori, R.K., Walston, A.A., Stafsudd, O.M., Fried, D., Walsh, J.T., Jr. (2001). Quantification and modeling of the dynamic changes in the absorption coefficient of water at X = 2.94 ^m. IEEE Journal on Selected Topics in Quantum Electronics, vol. 7, no. 6, p. 959-970. [5] Samokhin, A.A., Vovchenko, V.I., Ilichev, N.N. (2008). Photoacoustic and evaporation pressure signals in water irradiated with erbium laser pulses. Phys. Wave Phenomena, vol 16, no. 4, p. 275-282. [6] Livingston, F.E., George, S.M., Shori, R.K. (2002). Optimization of rotary Q-switched Er:YAG laser. Rev. Sci. Instrum, vol. 73, no. 7, p. 2526-2532. [7] Vodopyanov, K.L., Shori R.K., Stafsudd, O.M. (1998). Generation of Q-switched Er:YAG pulses using evanescent wave absorption in ethanol. Appl. Phys. Lett., vol. 72, no. 18, p. 2211-2213. [8] Korada, P., Nemec, M., Jelinkova, H., Sulc, J., Cech, M., Shi, Yi-Wei, Matsuura, Y., Miyagi, M. (2005). Electro-optically Q-switched Er:YAG laser and its application. Proceedings of the SPIE, XV International Symposium of Gas Flow, Chemical Lasers, and High Power Lasers, vol. 5777, p. 384389. [9] Lukač, M., Nemeš, K., Cenčič, S. (1993). Influence of direct and cross-relaxation pumping processes on the output energy and thermal load of an Er:YAG laser. Appl. Opt., vol. 32, p. 7399-7401. [10] Koechner, W. (1999). Solid state laser engineering, 4th ed. Springer-Verlag, Berlin. [11] Pollack, S.A., Chang, D.B., Birnbaum, M., Kokta, M. (1991). Upconversion-pumped 2.8-2.9 ^m lasing of Er3+ ion in garnets. J. Appl. Phys., vol. 70, p. 7227-7239. [12] Majaron, B., Rupnik, T., Lukač, M. (1996). Temperature and gain dynamics in flash lamp-pumped Er:YAG. IEEE J. Quantum Electron, vol. 32, no. 9, p 1636-1644.