Multi-focus non-periodic scanning method for femtosecond lasers based on DMD and galvanometer scanners [Invited]

Huaming Li; Yu Wang; Qinglei Hu; Zhuoyu Zhang; Xiaohua Lü; Shaoqun Zeng

doi:10.3788/COL202422.051701

1. Introduction

Femtosecond lasers are widely used in microfabrication^[1,2], biological imaging^[3-7], and optical storage^[8,9] due to their nonlinear optical effects. Generally, the femtosecond laser beam is focused to achieve a high peak power density for nonlinear optical effects, which means that it is only suitable for point scanning. However, the low scanning throughput of point scanning limits the application of the femtosecond laser. One of the methods to improve scanning throughput is multi-focus parallel scanning. Conventional optical devices that produce multipoints include the micro lens array (MLA), cascaded beam splitters (CBSs), diffraction optical elements (DOEs), acousto-optic deflectors (AODs), and spatial light modulators (SLMs). The MLA and CBS can simply generate multiple foci^[10,11], but the uniformity of the focus intensity is poor. The multi-focus generated by the DOE can possess superior intensity uniformity^[12], but because it is a diffractive element, it will have an angular dispersion effect when applied to the femtosecond laser. In addition, the MLA, CBS, and DOE are passive components, and the number and distribution of foci generated by them cannot be changed, so they lack flexibility in application. The AOD can change the number and distribution of foci by adjusting its loaded actuation signal^[13,14], which is more flexible than passive components but needs to solve the problems of time and spatial dispersion. At present, a more common method for generating multi-focus is to display a computer-generated hologram (CGH) on SLM^[15,16], which has high diffraction efficiency and optical modulation ability. However, due to the refresh rate of SLM being around 100 Hz, it is usually necessary to use it in conjunction with fast scanning components such as a galvanometer scanner, and the scanning trajectory can only be of a periodic structure, resulting in a loss of flexibility. Recently, research on multi-focus non-periodic scanning for femtosecond lasers has also developed. Geng et al. used the digital micromirror device (DMD), which can be seen as a special SLM with a refresh rate of up to 32 kHz, as both a scanner and a multi-focus generator. By switching the holograms displayed on the DMD, multi-focus parallel aperiodic scanning can be achieved^[17], but the scanning throughput is lower compared to the method of SLM combined with the galvo scanners. To solve the scanning flux problem, Jiao et al. proposed a method that combines scanning and projection methods, using the AOD as the scanner, the DOE as the multi-focus generator, and the DMD as the digital mask^[18]. Yang et al. advanced a method of combining the DMD with the MLA, where the DMD acts as a digital mask to change the number and distribution of multiple points generated by the MLA^[19].

In this Letter, we propose a new multi-focus non-periodic trajectory scanning method for femtosecond lasers based on DMD and galvo scanners. By loading different binary holograms, the DMD can generate and control one to dozens of foci, simultaneously. During the scanning process of the galvo scanners, the number of foci is constantly changing during the scanning process by rapid switching of the displayed hologram on the DMD, enabling the scanning of complex aperiodic structures. However, the regular method of changing the number of foci based on diffraction results in a change in intensity of focus, which leads to uneven processing quality. Here, we present a method to change the foci number by controlling aberration. It will significantly improve the intensity uniformity of foci arrays with different numbers.

2. Methods and Experimental Results

The principle of multi-focus non-periodic trajectory scanning for femtosecond lasers based on DMD and galvo scanners is shown in Fig. 1. Take scanning a triangle by four foci as an example. First, four foci distributed in a linear array are generated through the DMD, and the galvo scanners cause these foci to be moved parallel in the same direction. Second, when the red focus is scanned outside the triangular area, quickly switch the hologram loaded on the DMD that will change the number of foci to 3 (the red focus is “off”). By analogy, as shown in Fig. 1(a), once the foci are out of the target scanning area, the hologram will be switched to make these foci “off.” Similarly, when the foci enter the target scanning area, switch the hologram to make these foci “on,” as shown in Fig. 1(b). Finally, the scanning trajectory of the triangle area as shown in Fig. 1(c) is realized.

Figure 1.Schematic for multi-focus non-periodic trajectory scanning for a femtosecond laser based on DMD and galvo scanners. (a) The blue area is the target scanning area. The galvo scanners cause foci to achieve raster scanning, and the DMD controls the “on or off” of each of the four foci by quickly transforming the hologram. (b) Scanning control in the other direction of raster scanning. When foci enter the target area, the foci are turned “on.” (c) Final scanning trajectory.

Download full size

View all figures

High-intensity uniformity foci array generation is a key issue of the multi-focus scanning system, which determines the quality of high-precision fabrication. Using the same CGH method to generate multiple foci, the total diffraction efficiency is fixed, which means that the change in the foci amount will lead to a change in energy allocated to each focus. For example, each focal energy of two foci generated by the hologram is twice that of four foci. To solve this problem, a common way is to add energy-adjusting devices to the optical system, such as an acoustic-optic modulator (AOM) and electro-optic modulator (EOM), to quickly adjust the energy of the laser beam. However, due to the addition of these devices, the system complexity and system control difficulty increase. In this paper, we propose a method to independently control the “on and off” of each focus through aberration. Aberrations of the foci that need to be “on” are corrected to increase the focus energy density, and aberrations are introduced into the foci that need to be “off” to reduce the focus energy density. Because of the threshold effect and nonlinear effect of the femtosecond laser, the energy density of “off” focus cannot reach the processing threshold energy density, so it has no processing effect. In this method, the number of foci generated by each hologram does not change, so there are unnecessary additional energy modulators for adjusting the diffraction energy of different holograms.

Generally, the holograms that generate different numbers and distributions of focal points can be directly calculated by the weighted Gerchberg–Saxton (GSW) algorithm^[20]. As shown in Eq. (1), the phase value corresponding to the $j$ th pixel of the hologram generating $K$ foci after the $n$ th iteration is $Ø_{j}^{(n)} = \arg {\sum_{k = 1}^{K} w_{k}^{(n)} e^{i [φ_{k} (x_{j}, y_{j}) + θ_{k}^{(n)}]}},$ (1)where $φ_{k} (x_{j}, y_{j})$ is the phase of the $k$ th focus. The iteration proceeds of weight $w_{k}^{(n)}$ and phase $θ_{k}^{(n)}$ are as follows: $w_{k}^{(n)} = w_{k}^{(n - 1)} \frac{\frac{1}{K} \sum_{k = 1}^{K} | E_{k}^{(n - 1)} |}{| E_{k}^{(n - 1)} |},$ (2) $θ_{k}^{(n)} = \arg (E_{k}^{(n - 1)}) .$ (3)

The initial value of iteration $w_{k}^{(0)} = 1 / K$ , and $θ_{k}^{(0)}$ is a random phase. $E_{k}^{(n - 1)}$ is the complex amplitude of the $k$ th focus after the $(n - 1)$ th iteration. In pure phase modulation, the calculation expression for $E_{k}^{(n)}$ is $E_{k}^{(n)} = \sum_{j = 1}^{N} \frac{1}{N} e^{i [Ø_{j}^{(n)} - φ_{k} (x_{j}, y_{j})]} .$ (4)The DMD is a binary amplitude modulator; thus, the phase modulation $Ø_{j}$ needs to be transformed into amplitude modulation via Lee holography^[21]. The DMD-generated amplitude modulation $H_{j} (x_{j}, y_{j})$ is expressed as $H_{j} (x_{j}, y_{j}) = {\begin{matrix} 1, & | \frac{Ø_{j} (x_{j}, y_{j})}{2 π} + b | < \frac{1}{4} \\ 0, & otherwise \end{matrix},$ (5)where $b$ is an integer. Therefore, the calculation expression for $E_{k}^{(n)}$ is modified to $E_{k}^{(n)} = \sum_{j = 1}^{N} \frac{1}{N} H_{j} (x_{j}, y_{j}) e^{i [- φ_{k} (x_{j}, y_{j})]} .$ (6)

In this paper, in order to achieve focus “on or off” through aberration control, the focus phase is the sum of the scanning phase and the aberration correction phase, namely, $φ_{k} (x_{j}, y_{j}) = φ_{scan, k} (x_{j}, y_{j}) + φ_{AO, k} (x_{j}, y_{j}),$ (7)where $φ_{scan, k} (x_{j}, y_{j}) = \frac{2 π (x_{j} x_{k} + y_{j} y_{k})}{λ f} + \frac{z_{k} π (x_{j}^{2} + y_{j}^{2})}{λ f^{2}}$ (8)is the scanning phase of the $k$ th focus, and $φ_{AO, k} (x_{j}, y_{j}) = \sum_{a = 4}^{15} c_{a} Z_{a} (\frac{ρ}{ρ_{0}}, θ)$ (9)is the aberration correction phase of the corresponding focus, where ${\begin{cases} ρ = \sqrt{x_{j}^{2} + y_{j}^{2}}, \\ θ = \arctan (y_{j} / x_{j}), & x_{j} \geq 0 \\ θ = \arctan (y_{j} / x_{j}) + π, & x_{j} < 0 \end{cases} .$ (10) $(x_{k}, y_{k}, z_{k})$ represents the $k$ th focus coordinates. $c_{a}$ is the Zernike coefficient of the $a$ th Zernike mode, and $ρ_{0}$ is the radius of the objective lens pupil. Previously, we have demonstrated the obtained Zernike coefficient $c_{a, best}$ for aberration correction through iterative algorithms^[22]. In this way, when we make $c_{a} = - c_{a, best}$ in Eq. (9), the focus will introduce a huge aberration, achieving focus “off.”

The GSW algorithm, which is improved and applied to generate binary holograms in this article, is called the GSW for binary (GSW-B) algorithm. The intensity uniformity of multi-focus generated by binary holograms that are directly calculated using the GSW-B is low. To improve the uniformity of intensity, the hologram is adjusted through true focal intensity captured by the camera. As shown in Fig. 2, the actual light intensity is fed back to calculate the weight scaling factor $g_{k}^{(m)}$ . The expression for calculating the phase of a hologram is $Ø_{j}^{(n)} = \arg {\sum_{k = 1}^{K} g_{k}^{(m)} w_{k}^{(n)} e^{i [φ_{k} (x_{j}, y_{j}) + θ_{k}^{(n)}]}},$ (11)where $m$ is the number of times that the camera has feedback, and the initial value $g_{k}^{(0)}$ of $g_{k}^{(m)}$ is 1. During the iteration process, $g_{k}^{(m)} = g_{k}^{(m - 1)} \frac{{\bar{I}}^{(m - 1)}}{I_{k}^{(m - 1)}},$ (12)where $I_{k}^{(m - 1)}$ is the intensity of the $k$ th focus captured by the camera, and ${\bar{I}}^{(m - 1)}$ is the average intensity of all foci.

Figure 2.Schematic for aberration control and intensity homogenization process of the foci array. Scanning phase φ_scan,k determines the position of the foci; correction phase φ_AO,k determines the aberration (i.e., “on or off”) of the foci; the weight w_k can change the intensity of the foci; θ_k starts with a random phase. g_k is the weight scaling factor calculated by the actual focus intensity captured by the camera. Iteratively correct g_k, w_k, and θ_k to improve the intensity uniformity of the foci.

Download full size

View all figures

Figure 3(a) shows the generated focal point arrays of 1 to 36 recorded by the camera. It can be seen that all foci arrays actually have 36 foci, but the intensity of the “off” focus drops sharply by introducing aberration. As shown in Fig. 3(b), the intensity uniformity of the “on” focus of all foci arrays can reach 98%, and the intensity of the “off” focus is less than 20% of the intensity of the “on” foci. The intensity uniformity of the foci arrays is calculated as $U = 1 - \frac{\max (I) - \min (I)}{\max (I) + \min (I)},$ (13)where $U$ is the intensity uniformity of the foci and $I$ is the intensity of the foci.

Figure 3.Results of intensity uniformity for foci arrays with different numbers. (a) shows the foci arrays with a number of foci from 1 to 36, with each row being a focus array. The focus within the green triangle area is the “off” focus, in which the gray value is multiplied by 8. (b) Normalized intensity of all foci. The blue data points represent the “on” foci, and the red data points represent the “off” foci. The intensity uniformity of the “on” foci reaches 98%, and the intensity of the “off” foci is less than 20% of the intensity of the “on” foci.

Download full size

View all figures

The optical configuration of the scanning system is schematized in Fig. 4. The light source is a fiber femtosecond laser (HR-Femto-IR-50-30, Huaray) with an average power of 30 W at 1035 nm, which is polychromatic. The laser beam is expanded and collimated by a beam-expanding system composed of L1 and L2 ( $f_{L 1} = 50 mm$ , $f_{L 2} = 75 mm$ ). The DMD (DLP650L NIR, Texas Instruments) acts as a blazed grating in this system, with spatial dispersion that needs to be compensated^[23].Therefore, the blazed grating ( $1158_30 \times 13_3_N$ , Gitterwerk) with a period of 1 µm and a blazed angle of 31° is used to compensate for the spatial dispersion in the center direction of the scanning range. The lenses L3 and L4 ( $f_{L 3} = 75 mm$ , $f_{L 4} = 200 mm$ ) constitute a 4-f system to match the spatial dispersion produced by the blazed grating and DMD. After the DMD, a Keplerian dispersion compensation module (KDCM), which is a 4-f system with different magnifications at different wavelengths, is configured to compensate for residual dispersion for which the blazed grating cannot compensate^[24]. Then, a spatial filter is placed at the Fourier plane of the KDCM to block the unwanted diffraction order. Lastly, the 4-f system composed of lenses L5 and L6 (self-designed, with a magnification of $1 \times$ ) relays the galvo scanners (MINISCAN-II-10, RAYLASE) to the back-pupil of the objective (CFI Apo NIR $40 \times$ W, $NA = 0.8$ , Nikon). To expand the scanning range, the sample is mounted on a motorized precision XYZ stage (ANT130-160-XY-AVL125, Aerotech). The objective, tube lens ( $f_{TL} = 180 mm$ ), and CCD form the microscope system to monitor the cutting process in situ. An LED light source (XHP70, XLamp) is coupled into the system through the condensing lens and beam splitter (50:50) to provide illumination for the microscope.

Figure 4.Optical configuration of the scanning system. FS laser, fiber femtosecond laser; M1, reflectivity mirror; L1-L6, lens; DMD, digital micromirror device; KDCM, Keplerian dispersion compensation module; SF, spatial filter; Galvo, galvanometer scanners; DM, dichroic mirror; OBJ, objective lens; BS, beam splitter; CL, condensing lens; TL, tube lens.

Download full size

View all figures

Using our system, we make corneal flaps, which is a step of femtosecond laser-assisted in situ keratomileusis (FS-LASIK) in vitro porcine eyes. In FS-LASIK as shown in Fig. 5, a circular surface in the cornea with a diameter of about 8 mm needs to be cut by femtosecond laser. It is necessary to make a cylindrical side-cut with the same radius as the circular surface, and its height is about 120 µm. At the same time, the corneal flap pedicle with a width of about 4 mm was reserved for cornea and corneal flap connection. The scanning mode of mosaic splicing is used for circular surface scanning. As shown in Fig. 5(a), the pale blue area is the target scanning surface, each red square (i.e., mosaic) corresponds to an independent scanning field of view, and a circular scanning surface is composed of the mosaics (the mosaic gap is used to distinguish different scanning fields of view, which do not exist in the actual scanning process). There are two types of mosaics, namely mosaic 1 and mosaic 2, as shown in Fig. 5(a). Mosaic 1 is a square scanning surface, which can be scanned with multi-focus parallel raster scanning by the galvo scanners, as shown in Fig. 5(b). Mosaic 2 is located in the area of the edge of the circular surface, which is a typical non-periodic structure that cannot be scanned by the conventional multi-focus parallel raster scanning method. Thus, the scanning strategy described in this paper is adopted in mosaic 2. Figures 5(d)–5(g) show the results of making corneal flaps on the porcine cornea in vitro by the multi-focus scanning method described in this paper and the conventional multi-focus scanning method, respectively. The corneal flap made by our system does not yield serrated edges, indicating that our system can achieve multi-focus non-periodic structure scanning. It is also proved that the “on or off” of each focus can be controlled independently by aberration adjustment.

Figure 5.(a) The scanning mode of mosaic splicing is used to scan a circular surface. The light blue area is the target scanning surface. (b) The mosaic 1, which is completely inside the circle, is a square scanning surface that can be fully scanned with multi-focus parallel raster scanning. (c) The target scanning area of the mosaic 2, which is partly outside the circle, is non-periodic, and thus needs to use the scanning method in Fig. 1. (d) The result of making the corneal flap porcine cornea in vitro by the proposed multi-focus scanning method. (e) The result of making the corneal flap by the conventional multi-focus scanning method. (f), (g) The enlarged view of the corresponding part in (d), (e). (h) The result of lifting the corneal flap for (d).

Download full size

View all figures

Table 1 compares the time consumption for non-periodic structure scanning (a circular surface with a diameter of 8 mm) by using three scanning methods with the following parameters: the width of mosaic $l_{mosaic} = 320 μm$ , the height of mosaic $h_{mosaic} = 400 μm$ , the number of mosaics $N_{mosaic} = 416$ , scanning points pitch $d = 1.6 μm$ , the number of foci $N_{focus} = 8$ , the scan frequency of galvo scanners $f_{g} = 180 Hz$ , and the refresh rate of DMD $f_{DMD} = 10.75 kHz$ . The mosaic scanning time for the three scanning methods can be calculated as follows: $T_{single} = \frac{l_{mosaic}}{2 d f_{g}} \cdot N_{mosaic}, T_{multi, DMD} = \frac{l_{mosaic} h_{mosaic}}{d^{2} N_{focus} f_{DMD}} \cdot N_{mosaic}, T_{ourwork} = \frac{l_{mosaic}}{2 d N_{focus} f_{g}} \cdot N_{mosaic} .$ (14)

Table 1. Fabrication Time for a Circular Surface with a Diameter of 8 mm under Different Scanning Methods

View table
View all Tables
Table 1. Fabrication Time for a Circular Surface with a Diameter of 8 mm under Different Scanning Methods

Scanning Method Single-Focus Based on Galvanometers Multi-Focus Based on DMD Our Work
Mosaic scanning time 231 s 242 s 29 s
Splicing movement time 15 s 15 s 15 s
Total time 246 s 257 s 44 s

As a scanner, the scanning speed of the DMD is lower than that of the galvo scanner, because the DMD only scans one point when switching one hologram. The line scanning time of the DMD is given by $T_{line, DMD} = h_{mosaic} / (d f_{DMD})$ . However, the line scanning time of the galvo scanner is $T_{line, galvo} = 1 / (2 f_{g})$ . Under the scanning parameters used in this work, the scanning time of the multi-focus scanning method based on the DMD is better than the single-focus scanning method based on galvo only when $N_{focus} > 8$ . Our work combines the fast scanning speed of galvo scanners with the flexibility of DMD multi-focus scanning, while achieving non-periodic scanning and improving scanning speed. Compared to the other two methods, the total time-cost of our solution has been reduced by 82% and 83%, respectively.

3. Conclusion

In summary, we present a multi-focus parallel scanning method based on DMDs and galvanometers to realize high-throughput non-periodic trajectory scanning. Meanwhile, the GSW algorithm is improved to realize the independent control of the “on or off” of the foci by introducing and correcting aberration. The method described in this paper can be applied to the field of laser processing and laser medical treatment.

Category: Biophotonics

Received: Oct. 9, 2023

Accepted: Feb. 5, 2024

Published Online: May. 20, 2024

The Author Email: Shaoqun Zeng (sqzeng@mail.hust.edu.cn)

DOI:10.3788/COL202422.051701

CSTR:32184.14.COL202422.051701

Table 1. Fabrication Time for a Circular Surface with a Diameter of 8 mm under Different Scanning Methods

Table 1. Fabrication Time for a Circular Surface with a Diameter of 8 mm under Different Scanning Methods