We present an experimental demonstration of adaptive control of modal properties of optical beams. The control is achieved via heat-induced photothermal actuation of transmissive optical elements. We apply the heat using four electrical heaters in thermal contact with the element. The system is capable of controlling both symmetrical and astigmatic aberrations providing a powerful means for in situ correction and control of thermal aberrations in high power laser systems. We demonstrate a tunable lens with a focusing power varying from minus infinity to −10 m along two axes using SF57 optical glass. Applications of the proposed system include laser material processing, thermal compensation of high laser power radiation, and optical beam steering.
© 2010 OSA
The control of laser beam spatial profiles has become more important with the increased use of lasers in industry, research, and everyday life. Absorption of high-power laser radiation in optical materials deposits heat; the resulting changes in shape and refractive index distort the beam shape as the radiation is transmitted or reflected by the optical element. These distortions act as a lens, producing aberrations that depend on the absorbed laser power and incident spatial profile [1–4]. Consequently, the modal properties of the beam change as a function of incident power. An adaptive optics system that can compensate the distortion or otherwise control the modal properties of the laser beams is highly desirable. Possible applications of adaptive optics includes high power laser interferometers [5–7], laser material processing , image processing [9,10], and optical displays . For those applications where mechanical repositioning of lenses is not convenient or even possible, it is desirable to control adaptively the laser beam shape by electrical or optical means. There are a number of solutions available including thermal compensation using negative thermo-optic materials , CO2 induced heating and electrical heating of the elements [6,7], deformable mirrors , tunable liquid crystals , and MEMS mirrors . For high optical power applications, electrical heating from the perimeter of the optic, using “ring heaters,” can easily be used for adaptive compensation by controlling the electrical power delivered to one or more heating elements . These systems have been proposed and are planned to be used in high power laser systems such as Laser Interferometer and Gravitational-wave Observatory .
We can extend the concept of symmetric ring heaters to include angular segmentation, allowing the application of heat from localized points on the barrel of the optic. This allows for correction of non-cylindrically symmetric aberrations. The use of quadrant ring heaters to compensate for the astigmatic effects resulting from non-normal incidence mirrors, non-uniform absorption, and birefringence of the materials was introduced by us in Refs. 16 and 17. In this paper, we describe in detail the theory and experimental proof of the concept. Issues related to thermal beam steering, response time, and expected focal range are also discussed.
2. Figure of Merit
We begin by developing a figure of merit that allows us to compare different optical materials. The system consists of an optical element heated from the barrel as shown in Fig. 1 . We consider a cylinder of radius a and thickness h in a cylindrical co-ordinate system such that the radial direction of the optical element coincides with the radial coordinate r, with and the axial direction of the optical element is along the coordinate z, with The heater (or a system of radially segmented heaters) is placed on the barrel of the optic in good thermal contact with the optical element.
The thermal diffusion equation describing the temperature distribution of the optical element is:Equation (1) can be solved for both steady-state and transient solutions. In the steady state, Eq. (1) becomes:
In addition to the generation (or absorption) of heat within the optical element, there will be a heat flux incident on (or radiated from) the surfaces of the optical element. We break the heat flux q (W m−2) into two parts:Eq. (2) becomes a non-linear partial differential equation. The problem-specific boundary conditions are:
The simultaneous solution of Eqs. (2) and (5–7) will provide the steady state temperature distribution in the optical element heated by external heaters in thermal contact with the barrel of the optical element. Similarly, a transient solution can be obtained from Eq. (1).
However, we can simplify the solution considerably by assuming that the temperature distribution depends only on the radial coordinate. Then, the temperature distribution of a cylinder under constant heat load at the barrel is :Eq. (8) is approximate as it ignores the temperature dependence in the axial direction. Nevertheless, Eq. (8) shows how i) the temperature depends inversely upon the thermal conductivity of the material, and ii) a thermal gradient exists in the radial direction. This thermal gradient will induce through the temperature dependence of the refractive index a lens-like profile in the optical element. The power of this lens depends upon the optical properties of the material. The induced optical path length change is:Eq. (9) signifies that the barrel of the optical element is constrained and cannot expand radially. Equation (9) also represents a gradient index lens equation that can be described as a lens of focal length f:
The quantities used in Eqs. (9) and 10 are listed in Table 1 for a number of materials. Equation (10) provides a figure of merit (FOM) representing the strength of the lens produced for a given geometry and incident heat flux. This FOM, which gives an easy way to compare optical materials, is:
It should be noted that the actual temperature dependence requires a solution of Eq. (2) satisfying the boundary conditions. Equations (8-10) are stated here for the purpose of comparing material properties and not for the absolute determination of focal power. Based upon FOM values, we have chosen to use SF57 glass in our device. We use COMSOLTM, a finite element code, to model our specific geometry, which we compare with experiment below.
3. Experimental demonstration
The experimental arrangement is shown in Fig. 2 and 3 . The sample is a 1.0 cm long by 2.5 cm diameter SF57 glass disk (a flint glass produced by SchottTM USA) with four symmetrically positioned independent heating elements in thermal contact with its barrel (Fig. 2). To minimize conductive heating losses and maximize efficiency, the holder is made of Teflon®, a very low conductor of heat. The SF57 is placed in a vacuum chamber to eliminate convection cooling effects. The dominant heat transport mechanism is radiation; in addition, a small amount of heat is lost via conduction from the holder on which the SF57 is mounted. The heater pads are obtained from Minco Inc. and consist of nichrome wire sealed in Kapton. These heaters can be heated to 200° C, and are attached to the SF57 barrel by copper tape. The heaters are connected via copper wires brought outside the vacuum using electrical feedthroughs. There we use four independent power supplies to control the heating. The heaters have a resistance of 25 Ω and thus can be connected to fairly ordinary power supplies.
To characterize the heat-induced thermal lensing in SF57, we use a probe beam that passes through the SF57 via windows in the vacuum chamber. This probe beam is a 10 mW, 1.064 μm Nd:YAG laser, expanded to provide a nearly collimated (large Rayleigh range) 2 mm radius beam waist (1/e2 intensity) at the SF57 sample. The probe beam, after passing through the SF57 and exiting from the vacuum chamber, is divided into two beams. One part passes through a 1 m focal length spherical lens. This lens reduces the beam radius at its focal plane to about 100 μm. A CCD camera (WinCAM) installed on a rail is used to scan the beam in the axial (z) direction. The beam passes through a waist whose location and size depends upon the combination of the heat-induced SF57 thermal lens and the f = 1 m static lens. As the focal length of the thermal lens in SF57 changes, the beam waist size and location change. This change is measured by performing successive beam scans using the CCD camera. The analysis software associated with the camera provides Gaussian fits to the beam size at each axial (z) location. The thermal lens is computed by experimentally performing an axial beam scan of this beam and using standard ABCD matrices to calculate the thermal lens as a function of heating power.
The beam splitter sends the other portion of the beam to a 4-quadrant photodetector (QPD). This beam does not pass through any lens. The QPD is capable of measuring small beam angular deviations by measuring the differential power in the right and left halves for horizontal displacements and in the top and bottom halves for vertical displacements. Measurements of the distance between the QPD and SF57 and the displacement in the QPD plane provide the angular displacement induced by the heat-induced lens in the SF57.
3.1 Thermal lensing dynamic range test
The first series of tests were performed to check the linearity of the system. We varied the power of the heaters from 0 to 2.4 W. This power produces a 200°C temperature at the barrel of the optic. This temperature is the upper rating of the heaters, so the heating input was not increased further. This limitation does not apply to the SF57, which could be used at higher temperatures with different heaters. The temperature of the system was monitored by a thermocouple in situ as well as by an external IR thermal viewer for temperature profiling of the SF57 surface outside the vacuum. The scan data are shown in Fig. 4 for the X (horizontal) axis at each axial (Z-axis) location. Because the same power level was used for the vertical axis (top and bottom heaters), similar behavior was observed for y. The fitted M2 values do not change with power .
The extracted thermal lens from these measurements is shown in Fig. 5 . Here we have also plotted the COMSOL generated cumulative OPL change on the left axis in blue. These results were obtained by using the 3-D temperature distribution and integrating along the axial direction to obtain the OPL. The temperature distribution was obtained by solving Eqs. (1–7) numerically using COMSOL. The OPL is then converted to focal length and plotted on the right y-axis along with the measured values. The measured data is in good agreement with the simulation. The only free parameter used in the calculation is the temperature dependence of the thermal conductivity of SF57, which is not known to us. A value close to that for fused silica has been used. At the full heating power of 2.4 W per heater, a lens of −10 m focal length is obtained. Thus, the system focal length can be tuned from infinity to −10 m for 9.6 W of total heating power.
To verify further the numerical solution and the data, we obtained a thermal image of SF57 under full heating power with the thermal imaging camera (window limitations required this measurement to be done in air). The heat input was 2.4 W applied to each heater. The temperature data and the COMSOL simulation agree well, as shown in Fig. 6 . Also, the isotherms do indeed show a circularly symmetric temperature profile near the center, where the laser beam passes.
3.2 Astigmatic lens creation
One of the most powerful capabilities of this adaptive lens system is its ability to correct astigmatism (or to provide different focal lengths in the horizontal and vertical axes). This effect is achieved by differentially heating the SF57 in the x- or y-axes, resulting in the creation of two cylindrical lenses with their axes perpendicular to each other. The focal length of the two lenses can be controlled over a range of infinity to −8 m and −9 m in x and y axes respectively. Images taken at a fixed distance from the SF57 sample location are shown in Fig. 7 . The inferred thermal lens along the horizontal and vertical axes are given in Table 2 .
An important feature that can be observed from Table 2 is the cross coupling between the horizontal and vertical axes; i.e., a positive thermal lens is observed in the axis orthogonal to the heat axis. As an example, when a y-oriented thermal lens of −9 m is generated via heaters in the vertical direction, a 27 m focal length positive lens is measured in the x orientation. This behavior is due to the finite extent of the heaters. The heater geometry shown in Fig. 2 has been modeled in COMSOL and the result of applying heat in vertical direction only is shown in Fig. 8 . The data clearly show that, due to the finite heater size, a positive gradient (lower temperature at the edges) appears in the horizontal. However, this lens is weaker than the thermal lens generated on the axis to which heat is applied. This effect also manifests itself in the difference between thermal lensing in a particular axis when heat is applied only in one axis versus the case when heat is applied in both axes. For example, when the maximum heat is applied in both axes, the thermal lens has a focal length of −13 m whereas when the heat is applied in one direction only the thermal lens in that direction is −9 m.
This cross coupling can be reduced by ensuring that the heaters cover the maximum possible area on the barrel of the optics. The maximum possible area for a single heater is limited by the fact that independent control of each heater is desirable. We simulated the effect of the gap between two adjacent heaters on the cross coupling between x and y directions. We applied heat in the y-direction and measured the integrated OPL in both directions. If there were no cross coupling, then the OPL along x-direction would be zero. Figure 9 shows three data sets for 0, 10, and 20 degree gaps between the adjacent heaters. Heat input to the y-direction heaters was adjusted to provide same OPL in the y-direction. The x-direction OPL is also shown. It is evident that the cross coupling cannot be made zero but it can be reduced somewhat by decreasing the gap between the heaters. Going from a 20 degree gap to 0 degree gap, reduces the coupling from–0.44 to −0.33 where the coupling coefficient is defined as. We selected 4 mm radial distance as more than 99% of the energy in the probe beam is contained in this area.
3.3 Response time
We characterize the response time of the system by measuring the beam size as a function of time. These data are shown in Fig. 10 , where we plot the beam size at a fixed location as a function of time after the application of heat. The data cover the whole set of heating values shown in Fig. 4, i.e., for the cold case and for four heating powers. The thermal lens generated in the SF57 depends upon the measured beam difference defined as where and are the final and initial beam sizes. In the steady state, becomes . We define the response time as the time beam takes to change size from 10% to 90% of , excluding the first occurrence of the 90% beam size change caused by the overshoot. The response time of the adaptive device is found to be approximately 500 s.
The response time is also reproduced in the COMSOL simulation. We evaluated the optical path length change evolution at two points; at the center and at a point 2 mm away from the center. Because the beam size at the SF57 sample is 2 mm in radius, this position corresponds to OPL change across the beam spot. We can determine the response time of the system from the transient response of the OPL change. This response is shown in Fig. 11 where the green and the blue lines represent the simulated OPL change as a function of time for positions 2 mm off center and at the center of the SF57 sample respectively. Note that the thermal lensing only depends upon the radial OPL change and not on the absolute value of the OPL. The red curve shows the OPL difference (multiplied by 100X) between the green and the blue curves. The response time is approximately 600 s, in adequate agreement with measurement. The reason for the large overshoot is due to the low thermal conductivity of the SF57. When the heat is applied to the edge of the optic, the temperature of the edge starts rising instantly because it is in thermal contact with the heater. However, the heat takes some time before it reaches the center. Therefore a larger gradient between the edge and the center is developed initially. However, as the heat diffuses toward the center of the optic, the gradient decreases.
3.4 Thermal beam steering
The application of heat may produce some power-dependent beam steering from either off-axis beam propagation or asymmetries due to differences among the heating segments. If the centering of the beam is not perfect, the steering can be corrected through differential heating of the various parts of the SF57. The QPD shown in Fig. 3 is used to measure the thermally-induced beam steering. The beam size at the QPD was approximately 2.5 mm in diameter. In Fig. 12 , we plot the measured beam displacement (normalized to the beam size) as a function of time as heat is applied to various combinations of the segmented heaters. Here, both horizontal and vertical thermal beam drift has been measured. Each block in Fig. 12 represents a specific state of SF57 heating as indicated in the plot. We start by turning all heaters to the ‘ON’ state (full heating power). The bottom heater is turned ON at a slightly later stage. This delay does show up as beam drift in the vertical direction. Once all heaters are turned ON, the beam drift is less than 0.01% of the beam size. This level of beam drift corresponds to 0.25 μradian angular deviation. This small amount observed could actually be due to a slight misalignment in the vertical direction. For the horizontal beam drift, there is almost no change from the cold to the hot case.
Figure 12 illustrates another use for the system: high resolution beam steering. This effect demonstrated by selectively turning the heat ON as shown in separate grid blocks. For example, state 2 represents the case when the beam is steered diagonally by turning the left and bottom heaters to the OFF state. Various combinations of heating control show that the beam can be effectively steered in all directions. The range of the available steering is low; moreover there is also deviation of about 0.15% on the beam size if 2.4 W of heating is applied to one of the heaters. Thus, one should be careful in using this mechanism in steering applications. Coupled with the steering will be second-order phase front changes.
3.5 Beam quality
The quality of the thermal lens is also important to investigate. We have shown beam scan data in Fig. 4 that shows virtually no change in the M2 value as compared with the original beam. The M2 value remains at 1.20 ± 0.01 for all cases of heating as well as for no heating. In addition, we took beam profile measurements and monitored the TEM00 Gaussian content of the beam. There is no noticeable change in the TEM00 Gaussian content of the beam. This process was done for various beam sizes starting from 1 to 5 mm, with no noticeable degradation observed after the beam passed through SF57 thermal lens.
One possible concern with the experimental realization of the adaptive lens comes from the lack of rotational heating symmetry and possible non-uniformity. Figure 2 shows the configuration, with a gap between two adjacent heaters, making the temperature at the perimeter of the SF57 vary as a function of angle. This behavior is evident from measured and simulated data in Fig. 6, where the temperature at the barrel is lower in the gaps between the heaters, i.e., at 45° to the heater central axis. To show this more clearly, we have plotted in Fig. 13 the temperature profile along the heating axis and along the 45° line. (See Fig. 2) This plot represents the worst case of unintended astigmatism generated via SF57 heating. There is no difference in the temperature profile within the central 4.5 mm radius. Thus, this is not a fundamental limitation of the experiment, although it could be rectified by using heaters that have more angular extent.
The lens created by heating the barrel in a particular radial direction is very close to an ideal lens. This fact can be shown by taking a slice through a particular axis and fitting a quadratic function to the profile. Figure 14 shows a cross-section of the measured temperature distribution along the y-axis using a thermal IR camera. Here we have also plotted the residuals to the quadratic fitting equation. It can be seen that the residues are well within the ± 1°C resolution of the thermal camera.
There could be concerns about non-linearity of optical properties of materials at elevated temperature. However, it should be noted that the temperature differential across the optic is much smaller compared to the increase in over-all temperature of the optic. We simulated the temperature rise as well as the resultant integrated OPL in COMSOL. The results are shown in Fig. 15 . Here it is clear that the temperature profile and the OPL both have same regression coefficient. Therefore, it can be inferred that if the temperature profile is quadratic, the resultant OPL will also be quadratic as well.
3.6 Thermal depolarization
Optical birefringence can arise from heating-induced photoelastic stresses in glasses, resulting in thermal depolarization. This anisotropy can be an issue for polarization-sensitive experiments where depolarization is undesirable. To investigate this issue, we placed the SF57 between two crossed polarizers and measured the transmission of a specific polarization through the SF57 with no heating and with full heating. The results are summarized in Table 3 , which shows that there is negligible thermal depolarization. The theoretical justification of this experimental result would be obtained by doing a photoelastic analysis of SF57. However, because the temperature-dependent stress-optical coefficient data for SF57 is not known, we are unable to complete such an analysis. In view of the experimental results above, we can say that, to the precision of the experiment, SF57 does not show any thermal depolarization.
3.7 Coating reflectivity change
Another potential issue is the possible change in the reflectivity of SF57 when heated to high temperature. This change would be directly related to the quality of the anti-reflection (AR) coatings available. In these experiments, we used commercial off-the-shelf AR coating on both SF57 surfaces. These coatings were determined to have a reflectance of approximately 500 ppm per surface at room temperature. The reflectivity increased to 1000 ppm per surface at 200°C. Potentially better coatings could be obtained if optical losses are an issue in a particular application or the coating could be designed for a specific working temperature.
In conclusion, we have demonstrated an adaptive beam shaping element that can change the beam shape effectively independently along two orthogonal axes. A focal length dynamic range of minus infinity to −10 m has been achieved by heating from room temperature to 200°C. We have investigated beam steering, beam shape, beam quality, thermal depolarization, and reflectivity change as a function of temperature. The beam quality does not change significantly over this range of temperature. The proposed system was also tested for astigmatic lens creation. A differential lens of −10 m can be obtained by applying heating to one axis only. Another possible application of the system is its ability to be used for high-resolution, slow-speed beam steering.
This segmented heating system provides a powerful way to shape optical beams without the need for moving elements. It can be used in vacuum, where accessing and moving lenses proves difficult. The proposed system is particularly useful with high laser power systems where other alternatives like deformable mirrors, digital micro-mirror devices, and liquid crystals may not work. Another application of the proposed system is adaptive correction of high power thermal effects and an ability to correct astigmatism.
The authors acknowledge the encouragement of the LIGO Science Collaboration. Discussions with Volker Questchke, Rich Abbott, Peter Fritschel, and Dennis Coyne have been very helpful. This work is supported by the National Science Foundation under grant PHY-0855313.
References and links
1. J. D. Foster and L. M. Osterink, “Thermal effects in Nd:YAG Laser,” Appl. Opt. 41, 3656–3663 (1970).
2. E. Greninger, “Thermally induced wave-front distortions in laser windows,” Appl. Opts., 41, 549–552, (1986).A. E. Siegman, Lasers, University Science books, Sausalito, CA (1984).
3. P. Hello and J. Vinet, “Analytical models of thermal aberrations in massive mirrors heated by high power laser beams,” J. Phys. (France) 51, 1267–1282 (1990). [CrossRef]
4. P. Hello and J. Vinet, “Analytical models of transient thermoelastic deformations of mirrors heated by high power CW laser beams,” J. Phys. (France) 51, 2243–2261 (1990). [CrossRef]
5. G. Mueller, R. S. Amin, D. Guagliardo, D. McFeron, R. Lundock, D. H. Reitze, and D. B. Tanner, “Method for compensation of thermally induced modal distortions in the input optical components of gravitational wave interferometers,” Class. Quantum Gravity 19(7), 1793–1801 (2002). [CrossRef]
7. M. A. Arain, V. Quetschke, J. Gleason, L. F. Williams, M. Rakhmanov, J. Lee, R. J. Cruz, G. Mueller, D. B. Tanner, and D. H. Reitze, “Adaptive beam shaping by controlled thermal lensing in optical elements,” Appl. Opt. 46(12), 2153–2165 (2007). [CrossRef]
8. R. Schmiedl, “Adaptive optics for CO2 laser material processing, ” in 2nd International Workshop on Adaptive Optics for Industry and Medicine, G. D. Love, ed. (World Scientific Publishing Co Pte Ltd, 2000), pp. 32–36.
9. S. Sato, “Liquid-crystal lens-cells with variable focal length,” Jpn. J. Appl. Phys. 18(9), 1679–1684 (1979). [CrossRef]
10. W. L. IJzerman, S. T. de Zwart, and T. Dekker, “Design of 2D/3D switchable displays,” SID Symposium Digest, 36, 98–101 (2005).
11. T. L. Kelly, A. F. Naumov, M. Yu. Loktev, and M. A. Rakhmatulin, “Focusing of astigmatic laser diode beam by combination of adaptive liquid crystal lenses,” Opt. Commun. 181(4-6), 295–301 (2000). [CrossRef]
12. I. Kanno, T. Kunisawa, T. Suzuki, and H. Kotera, “Development of deformable mirror composed of piezoelectric thin films for adaptive optics,” IEEE J. Sel. Top. Quantum Electron. 13(2), 155–161 (2007). [CrossRef]
13. H. Ren and S.-T. Wu, “Adaptive liquid crystal lens with large focal length tunability,” Opt. Express 14(23), 11292–11298 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-23-11292. [CrossRef]
14. T.-Y. Chen, C.-H. Li, J.-L. Wang, C. E. Chiu, and G. J. Su, “A MEMS-based Organic Deformable Mirror with Tunable Focal Length,”2007 IEEE/LEOS International Conference on Optical MEMS and Nanophotonics, July 16 (2007), pp. 103–104.
15. M. Smith, and P. Willems, Auxiliary Optics Support System Conceptual Design Document, Volume 1 Thermal Compensation System, LIGO-T060083–00-D, http://docuserv.ligo.caltech.edu/docs/public/T/T060083-00/T060083-00.pdf
16. M. A. Arain, V. Quetschke, L. F. Williams, G. Mueller, D. B. Tanner, and D. H. Reitze, “Elements for Future Gravitational Wave Interferometers,’ Frontiers in Optics, OSA meeting, San Jose, CA, September 2007.
17. M. A. Arain, V. Quetschke, L. F. Williams, R. Martin, G. Mueller, D. B. Tanner, and D. H. Reitze, “Adaptive Optical Elements for Laser Beam Shaping,” US provisional patent, Serial No. 61/086,661, US patent pending.
18. F. P. Incropera, and D. P. DeWitt, Fundamentals of Heat and Mass Transfer, 5th ed., John Wiley & Sons, USA, (2002).
19. A. E. Siegman, “Defining, measuring, and optimizing laser beam quality,” Proc. SPIE 1868, 2 (1993). [CrossRef]