Citation:

Ultrafast z-scanning for high-efficiency laser micro-machining


  • Light: Science & Applications  7, Article number: (2018)
More Information
  • Corresponding author:
    CB Arnold, Email: cbarnold@princeton.edu
  • Received: 28 June 2017
    Revised: 29 December 2017
    Accepted: 29 December 2017
    Published online: 20 April 2018

doi: https://doi.org/10.1038/lsa.2017.181

  • High-throughput laser micro-machining demands precise control of the laser beam position to achieve optimal efficiency, but existing methods can be both time-consuming and cost-prohibitive. In this paper, we demonstrate a new high-throughput micro-machining technique based on rapidly scanning the laser focal point along the optical axis using an acoustically driven variable focal length lens. Our results show that this scanning method enables higher machining rates over a range of defocus distances and that the effect becomes more significant as the laser energy is increased. In a specific example of silicon, we achieve a nearly threefold increase in the machining rate, while maintaining sharp side walls and a small spot size. This method has great potential for improving the micro-machining efficiency of conventional systems and also opens the door to applying laser machining to workpieces with uneven topography that have been traditionally difficult to process.
  • 加载中
  • [1] Malinauskas M, Žukauskas A, Hasegawa S, Hayasaki Y, Mizeikis V et al. Ultrafast laser processing of materials: from science to industry. Light Sci Appl 2016; 5: e16133. doi: 10.1038/lsa.2016.133
    [2] Shiu PP, Knopf GK, Ostojic M, Nikumb S. Rapid fabrication of tooling for microfluidic devices via laser micromachining and hot embossing. J Micromech Microeng 2008; 18: 025012. doi: 10.1088/0960-1317/18/2/025012
    [3] Saga T. Advances in crystalline silicon solar cell technology for industrial mass production. NPG Asia Mater 2010; 2: 96–102. doi: 10.1038/asiamat.2010.82
    [4] Kim TN, Campbell K, Groisman A, Kleinfeld D, Schaffer CB. Femtosecond laser-drilled capillary integrated into a microfluidic device. Appl Phys Lett 2005; 86: 201106. doi: 10.1063/1.1926423
    [5] Baker CA, Bulloch R, Roper MG. Comparison of separation performance of laser-ablated and wet-etched microfluidic devices. Anal Bioanal Chem 2011; 399: 1473–1479. doi: 10.1007/s00216-010-4144-3
    [6] Gu E, Jeon CW, Choi HW, Rice G, Dawson MD et al. Micromachining and dicing of sapphire, gallium nitride and micro LED devices with UV copper vapour laser. Thin Solid Films 2004; 453-454: 462–466. doi: 10.1016/j.tsf.2003.11.133
    [7] Wang DL, Cui HJ, Hou GJ, Zhu ZG, Yan QB et al. Highly efficient light management for perovskite solar cells. Sci Rep 2016; 6: 18922. doi: 10.1038/srep18922
    [8] Carey JE, Crouch CH, Shen M, Mazur E. Visible and near-infrared responsivity of femtosecond-laser microstructured silicon photodiodes. Opt Lett 2005; 30: 1773–1775. doi: 10.1364/OL.30.001773
    [9] Davis SP, Prausnitz MR, Allen MG Fabrication and Characterization of Laser Micromachined Hollow Microneedles. Proceedings of 12th International Conference on, 2003 TRANSDUCERS, Solid-State Sensors, Actuators and Microsystems; 8–12 June 2003; IEEE: Boston, MA, USA, 2003..
    [10] Brazzle J, Papautsky I, Frazier A. Micromachined needle arrays for drug delivery or fluid extraction. IEEE Eng Med Biol Mag 1999; 18: 53–58. doi: 10.1109/51.805145
    [11] Caliano G, Carotenuto R, Cianci E, Foglietti V, Caronti A et al. Design, fabrication and characterization of a capacitive micromachined ultrasonic probe for medical imaging. IEEE Trans Ultrason Ferroelectr Freq Control 2005; 52: 2259–2269. doi: 10.1109/TUFFC.2005.1563268
    [12] Oralkan O, Ergun AS, Johnson JA, Karaman M, Demirci U et al. Capacitive micromachined ultrasonic transducers: next-generation arrays for acoustic imaging? IEEE Trans Ultrason Ferroelectr Freq Control 2002; 49: 1596–1610. doi: 10.1109/TUFFC.2002.1049742
    [13] Nijs JF, Szlufcik J, Poortmans J, Sivoththaman S, Mertens RP. Advanced cost-effective crystalline silicon solar cell technologies. Solar Energy Mater Solar Cells 2001; 65: 249–259. doi: 10.1016/S0927-0248(00)00100-8
    [14] Gower MC. Industrial applications of laser micromachining. Opt Express 2000; 7: 56–67. doi: 10.1364/OE.7.000056
    [15] Neuhaus DH, Münzer A. Industrial silicon wafer solar cells. Adv OptoElectron 2007; 2007: 24521. doi: 10.1155/2007/24521
    [16] Bechtold P, Zimmermann M, Roth S, Alexeev I, Schmidt M. Ultrashort laser pulse technology. Opt Sci 2015; 3: 245.
    [17] Duocastella M, Arnold CB. Bessel and annular beams for materials processing. Laser Photon Rev 2012; 6: 607–621. doi: 10.1002/lpor.201100031
    [18] Kammel R, Ackermann R, Thomas J, Götte J, Skupin S et al. Enhancing precision in fs-laser material processing by simultaneous spatial and temporal focusing. Light: Sci Appl 2014; 3: e169. doi: 10.1038/lsa.2014.50
    [19] Sanner N, Huot N, Audouard N, Larat C, Huignard JP et al. Programmable focal spot shaping of amplified femtosecond laser pulses. Opt Lett 2005; 30: 1479–1481. doi: 10.1364/OL.30.001479
    [20] Lin JY, Huang RP, Tsai PS, Lee CH. Wide-field super-resolution optical sectioning microscopy using a single spatial light modulator. J Opt A 2008; 11: 015301. doi: 10.1088/1464-4258/11/1/015301
    [21] Bush K, German D, Klemme B, Marrs A, Schoen M Electrostatic membrane deformable mirror wavefront control systems: design and analysis. Proceedings of the SPIE 5553, Advanced Wavefront Control: Methods, Devices, and Applications Ⅱ; 12 October 2004, SPIE: Denver, CO, USA, 2004. 5553: 28–38..
    [22] Bifano TG, Perreault JA, Bierden PA, Dimas CE Micromachined Deformable Mirrors for Adaptive Optics. Proceedings of the SPIE 4825, High-Resolution Wavefront Control: Methods, Devices, and Applications Ⅳ; 6 November 2002, SPIE: Seattle, WA, USA, 2002. 4825: 10–13.
    [23] Ren H, Wu ST. Variable-focus liquid lens by changing aperture. Appl Phys Lett 2005; 86: 211107. doi: 10.1063/1.1935749
    [24] Ren HW, Fox D, Anderson PA, Wu B, Wu ST. Tunable-focus liquid lens controlled using a servo motor. Opt Express 2006; 14: 8031–8036. doi: 10.1364/OE.14.008031
    [25] Naumov AF, Loktev MY, Guralnik IR, Vdovin G. Liquid-crystal adaptive lenses with modal control. Opt Lett 1998; 23: 992–994. doi: 10.1364/OL.23.000992
    [26] Mermillod-Blondin A, McLeod E, Arnold CB. High-speed varifocal imaging with a tunable acoustic gradient index of refraction lens. Opt Lett 2008; 33: 2146–2148. doi: 10.1364/OL.33.002146
    [27] Chen TH, Ault JT, Stone HA, Arnold CB. High-speed axial-scanning wide-field microscopy for volumetric particle tracking velocimetry. Exp Fluids 2017; 58: 41. doi: 10.1007/s00348-017-2316-z
    [28] Duocastella M, Arnold CB. Enhanced depth of field laser processing using an ultra-high-speed axial scanner. Appl Phys Lett 2013; 102: 061113. doi: 10.1063/1.4791593
    [29] Lu QM, Mao SS, Mao XL, Russo RE. Delayed phase explosion during high-power nanosecond laser ablation of silicon. Appl Phys Lett 2002; 80: 3072–3074. doi: 10.1063/1.1473862
    [30] Yoo JH, Jeong SH, Greif R, Russo RE. Explosive change in crater properties during high power nanosecond laser ablation of silicon. J Appl Phys 2000; 88: 1638–1649. doi: 10.1063/1.373865
    [31] Schwarz-Selinger T, Cahill DG, Chen SC, Moon SJ, Grigoropoulos CP. Micron-scale modifications of Si surface morphology by pulsed-laser texturing. Phys Rev B Condens Matter Mater Phys 2001; 64: 155323. doi: 10.1103/PhysRevB.64.155323
    [32] Liu JM. Simple technique for measurements of pulsed Gaussian-beam spot sizes. Opt Lett 1982; 7: 196–198. doi: 10.1364/OL.7.000196
    [33] Hoel PG, Port SC, Stone CJ Introduction to Probability Theory. Boston, MA: Houghton Mifflin; 1971.
    [34] Srinivasan R, Braren B. Ablative photodecomposition of polymer films by pulsed far-ultraviolet (193 nm) laser radiation: dependence of etch depth on experimental conditions. J Polymer Sci 1984; 22: 2601–2609.
    [35] Srinivasan R, Braren B. Ultraviolet laser ablation of organic polymers. Chem Rev 1989; 89: 1303–1316. doi: 10.1021/cr00096a003
    [36] Wang WJ, Mei XS, Jiang GD, Lei ST, Yang CE. Effect of two typical focus positions on microstructure shape and morphology in femtosecond laser multi-pulse ablation of metals. Appl Surf Sci 2008; 255: 2303–2311. doi: 10.1016/j.apsusc.2008.07.100
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(8) / Tables(1)

Article Metrics

Article views(373) PDF downloads(93) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

Ultrafast z-scanning for high-efficiency laser micro-machining

  • Department of Mechanical and Aerospace Engineering, and Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ 08544, USA
  • Corresponding author:

    Craig B Arnold, cbarnold@princeton.edu

doi: https://doi.org/10.1038/lsa.2017.181

Abstract: High-throughput laser micro-machining demands precise control of the laser beam position to achieve optimal efficiency, but existing methods can be both time-consuming and cost-prohibitive. In this paper, we demonstrate a new high-throughput micro-machining technique based on rapidly scanning the laser focal point along the optical axis using an acoustically driven variable focal length lens. Our results show that this scanning method enables higher machining rates over a range of defocus distances and that the effect becomes more significant as the laser energy is increased. In a specific example of silicon, we achieve a nearly threefold increase in the machining rate, while maintaining sharp side walls and a small spot size. This method has great potential for improving the micro-machining efficiency of conventional systems and also opens the door to applying laser machining to workpieces with uneven topography that have been traditionally difficult to process.

INTRODUCTION
  • The ability to laser machine materials with high resolution and high throughput is critical in advanced manufacturing for a vast array of applications, from photovoltaic cells to bio-compatible micro-components1-15. The precision of these manufacturing techniques relies on focusing a laser beam to a micron-sized spot onto the surface of the workpiece. This tight lateral focusing comes at a price since it narrows the depth of field (DOF) along the axial direction, causing a concomitant drop in machining efficiency outside the reduced range. Therefore, the surface of the workpiece has to be carefully maintained within the focal position of the laser beam to ensure a high efficiency; as material is removed, it is important to continue to adjust the location of the laser focus to follow the change in topography of the workpiece16.

    One common method used to mitigate a narrow vertical machining range consists of extending the DOF of the system by using, for example, a low focusing power lens or structured light. However, such attempts lead to a significant loss in lateral resolution17. Alternatively, the focused beam can be adjusted with respect to the surface location throughout the material removal process, but this requires real-time knowledge of the surface location as well as a fast focusing method. While some studies have attempted to acquire real-time surface location information to increase the machining efficiency, the complexity and inflexibility of the resulting machining system drastically reduces its potential for practical use18. In addition to the difficulty of real-time surface monitoring, most controlled focusing methods also suffer from low response rates compared to the repetition rate of the laser16, 18-25.

    Here, we demonstrate a new method for increasing the micro-machining efficiency by combining a high repetition rate laser with an ultrafast z-scanner. This creates an axial distribution for the focused beam that allows the pulses to be focused onto different surface locations of the material without real-time monitoring. More importantly, the extended range for focal positions also relaxes the constraints on surface flatness and positioning. The method shown in this paper represents a promising advance in high-efficiency micro-machining and provides great potential for laser machining in a broad range of materials.

MATERIALS AND METHODS
  • The experimental setup is shown in Figure 1a. A Nd:YVO4 laser (Coherent, Inc., Santa Clara, CA, USA, 355 nm, 15 ns pulse duration) is guided through a Tunable Acoustic Gradient-Index (TAG) lens (TAG Optics, Inc.) and focused onto the substrate, which is in our case a 500-μm-thick silicon wafer, using a microscope objective (5 ×, N.A. 0.13). The TAG lens is an ultrafast vari-focal device that has a lens power ($\frac{1}{{{f_{\rm{T}}}}}$, fT is the focal length of TAG lens) that varies sinusoidally, with the oscillation frequency set to 140 kHz26. This results in an adjustable focus position after the objective, which is indicated as the scanning range, RT, in Figure 1a. The scanning range is simulated in Zemax (Zemax LLC) using the TAG lens model detailed in Ref. 27. Table 1 shows a few relevant scanning ranges.

    Fig. 1   

    (a) Experimental setup of the system. The laser beam is guided through the TAG lens and objective before being focused onto the substrate. The TAG lens is an ultrafast vari-focal device that generates sine wave oscillations in lens power (a quantity that is linearly dependent on the driving voltage). These oscillations result in variations in the focal position below the objective, which we call the laser scanning range, or RT. (b) In this paper, we adopt the convention that the laser beam travels in the negative z direction. The beam waist is shown as a curved black line; the position at which the beam waist is narrowest defines the focal point. Ablation occurs when the fluence is larger than the threshold fluence. The blue circle indicates the region in which ablation occurs: the ablation range. The bottom diagram shows the line scan path during the milling experiment, where the blue circles represent the focused beam location and the red lines indicate the beam path along the substrate.
    Lens power (m−1) RT (5 ×) (mm) RT(15 ×) (mm)
    0.45 0.73 0.09
    0.61 0.99 0.14
    0.76 1.24 0.16
    0.91 1.51 0.19
    1.06 1.80 0.23
    1.21 2.09 0.27

    Table 1.  The scanning range corresponding to the TAG lens power is simulated in Zemax with 5 × and 15 × objectives

    For the silicon experiment, we machine a nominal 200 μm by 200 μm square hole with and without z-scanning. The silicon thickness is 500 μm. The spacing between pulses is 1 μm, and the laser repetition rate is 1 kHz. Each pass is accomplished by bi-directional line scanning over the entire square, as in Figure 1b, with 175 to 200 pulses per line, depending on the spot size.

    A similar setup is utilized in experiments carried out for Kapton with a thickness of 135 μm. A Nd:YAG laser (EKSPLA, 355 nm, 30 ps) and a 15 × microscope objective lens (N.A. 0.32) is used in this case. The different materials machined under various scanning ranges are examined and characterized under a scanning electron microscope and a laser scanning confocal microscope (Olympus LEXT).

RESULTS AND DISCUSSION
  • We machine a 200 μm by 200 μm square hole with and without z-scanning while maintaining a constant laser peak fluence of 500 J cm−2. Without scanning, four passes are required to completely machine the square hole, as shown in Figure 2a. This results in an average machining rate of 9.04 × 108 mm3 per pulse. After enabling the z-scanner at a lens power of 0.91 m−1, the machining is ~88 % complete after two passes and 100 % complete after three passes, as shown in Figure 2b and 2c. Increasing the lens power to 1.06 m−1 or 1.21 m−1, we find that 90 % of the machining is complete after two passes and is fully complete after three, as shown in Figure 2d and 2e. The measured ablated volume versus pulse number is shown in Figure 2f. This corresponds to an increased rate of 2.50 × 107 mm3 per pulse at a lens power of 1.06 m−1; therefore, the machining rate is more than doubled with the assistance of the ultrafast z-scanner.

    Fig. 2   

    (a) SEM images of a square hole ablated without a z-scanner (lens power of 0 m−1). The hole is a 200 μm by 200 μm square with a depth of 500 μm. Without z-scanning, we require four passes over the substrate to ablate the desired volume. However, we require only three passes with lens powers of (b) 0.76 m−1 and (c) 0.91 m−1. When we further increase the lens power to (d) 1.06 m−1 and (e) 1.21 m−1, 90% of the desired ablation is completed in only two passes. The measured ablated volume versus pulse number is shown in (f). The scale bar shown corresponds to 100 μm.

    This increase in machining rate is a highly promising result, but such gains are often accompanied by degradation of the machining quality, i.e. a poor lateral resolution and wall angle. Variations in the lateral resolution are directly related to the smallest spot size that the laser can produce on the substrate. With z-scanning enabled, the spot size is predicted to be larger: in this specific case, three times larger than that achieved without scanning due to the fluence being 500 times larger than the threshold fluence28. Such a minor degradation in resolution is acceptable in most of the aforementioned applications.

    To investigate the wall angle produced by our system, we machine a 200 μm by 200 μm square on silicon at different lens powers. Figure 3a and 3b shows the machining results obtained for three passes at lens powers of 0 and 0.91 m−1. The horizontal and vertical cross-sectional profiles of the square holes are shown in the second and third rows of Figure 3a and 3b. The difference between the horizontal and vertical cross-sectional profiles of the square hole can be caused by the re-deposition of materials29-31. The wall angles corresponding to each side of the square are measured separately and are shown in Figure 3c as angular deviations from the desired right angle. The data show that this discrepancy increases slightly at higher lens powers; however, the discrepancies are more uniform compared to the case for zero lens power and never exceed five degrees. Therefore, we doubled the machining efficiency of the system with a minimal loss in the quality of the lateral resolution and wall angle.

    Fig. 3   

    (a) To compare machining quality, we ablate a 200 μm by 200 μm square on a 500-μm-thick silicon wafer with three passes at a lens power of (a) 0 m−1 (no z-scanning) and (b) 0.91 m−1. The second row of (a and b) shows the horizontal cross-sectional profiles, as indicated by red lines in the first row. The third row similarly shows the vertical cross-sectional profiles, as indicated by green lines in the first row. In (c), we plot the wall angles as a function of the lens power, where the wall angle is a measure of the deviation from the desired right angle. Clearly, for a fixed number of passes, the addition of z-scanning produces a much more uniform wall. In addition, deviations grow slowly with increasing lens power.

    The benefits of z-scanning extend beyond an increased machining efficiency. As our system scans through different focal positions, the laser pulse reaches different axial positions at the surface, which enables uniform machining over a range of axial substrate positions. This suggests that similar machining efficiencies can be achieved even for a non-flat surface. To test our hypothesis, we place the sample at a fixed defocus distance, ${z'_0} = {z_{{\rm{s}}0}} - {z_{{\rm{f}}0}}$, which is the distance between the focal position, zf0, and the surface, zs0, without z-scanning, and machine a square, as before, at three different lens powers. The process is illustrated in Figure 4a. We then measure the normalized exit-hole area, noting that a value of one indicates ideal performance. A plot of this area over a range of ${z'_0}$ values and at three different lens power levels is shown in Figure 4b. Our hypothesis implies that at a higher lens power, we should obtain similar machining results over a range of ${z'_0}$ values. This behavior is confirmed in Figure 4b. This means that z-scanning effectively relaxes tight focusing constraints, which is particularly beneficial for samples with rough surfaces, such as those found in most real-world materials.

    Fig. 4   

    (a) The surface of the substrate at zs0 is placed at a fixed distance from the focal point at zf0 throughout the machining process. We refer to the difference in these quantities as the defocus distance and measure the exit-hole area at three different lens powers and over a range of ${z'_0}$ values (three passes were used in all cases). The results shown in (b) are normalized by the designed hole area and are shown with error bars determined from ten different experiments.

    So far, we have compared the machining results from fast scanning to those from a fixed focus. However, one can hypothesize that a stepwise focal adjustment after each machining pass could further improve the material removal rates. We consider this method in the Supplementary Information (Supplementary Fig. S1) and find that it is not as efficient as the fast scanning method. The reason for this inefficiency is that stepwise machining assumes perfect knowledge of the surface location, which is not necessarily uniform because of material re-deposition and other phenomena. Use of our new method results in significantly increased machining rates at all defocus distances, with a negligible loss of product quality, and further enables efficient machining on non-flat surfaces that would have otherwise required expensive and slow focus control systems.

  • Non-scanning. To explain how focus scanning improves the machining efficiency and enables an extended machining range in the z direction, we calculate the machining efficiency as a function of the defocus position, z'. The efficiency of a micro-machining system is directly related to the amount of material removed by a single laser pulse, which we call the single-shot ablated volume, or Vss. This value can be expressed as the integral of the ablation depth, D, over the ablated area, A:

    $$ {V_{{\rm{ss}}}} = \int_A D \text{d}A $$ (1)

    The ablation depth can be derived based on knowledge of specific material and laser properties. In our case, we assume that the Beer–Lambert model holds, so

    $$ D = L{\rm{ln}}\left( {\frac{F}{{{F_{{\rm{th}}}}}}} \right) $$ (2)

    where F is the laser fluence, L is the effective laser penetration depth, and Fth is the threshold fluence of the material. Although L and Fth are measured experimentally, as detailed in the Supplementary Information, we must still determine an expression for F. The laser is a Gaussian beam, which for a fixed focal position, zf, causes the laser fluence to take on the following functional form:

    $$ F(E,r,{z^\prime }) = \frac{{2E}}{{\pi {\omega ^2}({z^\prime })}}{e^{ - \frac{{2{r^2}}}{{{\omega ^2}({z^\prime })}}}} $$ (3)

    where $z' = z - {z_{\rm{f}}}$ is the defocus distance measured from the focal position, $\omega \left({z'} \right) = {\omega _0}\sqrt {1 + {{\left({\frac{{z'}}{{{z_{\rm{R}}}}}} \right)}^2}} $ is the beam diameter at a particular z position, ${z_{\rm{R}}} = \frac{{\pi \omega _0^2}}{\lambda }$ is the Rayleigh length, E is the pulse energy, and r is the radial coordinate in the xy plane. The beam waist, ω0, is measured by experiments detailed in the Supplemental Information32. By substituting Equation (3) into Equation (2), we find

    $$ D(E,r,{z^\prime }) = L{\rm{ln}}\left( {\frac{{2E}}{{\pi {\omega ^2}({z^\prime }){F_{{\rm{th}}}}}}} \right) - \frac{{2L{r^2}}}{{\omega {{({z^\prime })}^2}}} $$ (4)

    The ablation depth is a function of r2, indicating that the ablation hole will form a paraboloid. As the volume of a paraboloid has a simple analytic expression, we can rewrite Equation (1) as

    $$ {V_{{\rm{ss}}}} = \frac{1}{2}{D_{{\rm{ab}}}}{A_{{\rm{ab}}}} $$ (5)

    where Dab is the depth of the hole and Aab is the base area. We find Dab by setting r=0, in Equation (4),

    $$ {D_{{\rm{ab}}}}(\varepsilon ,{z^\prime }) = L{\rm{ln}} \left( {\frac{{2E}}{{\pi {\omega ^2}({z^\prime }){F_{{\rm{th}}}}}}} \right) = L{\rm{ln}} \left( {\frac{\varepsilon }{{1 + {{\left( {\frac{{{z^\prime }}}{{{z_R}}}} \right)}^2}}}} \right) $$ (6)

    where we have introduced the non-dimensionalized pulse energy $\varepsilon \equiv \frac{E}{{{E_{{\rm{th}}}}}}$, where ${E_{{\rm{th}}}} = \frac{{\pi \omega _0^2{F_{{\rm{th}}}}}}{2}$ is the threshold energy required for ablation at the focal point. Turning to Aab in Equation (5), we have

    $$ {A_{{\rm{ab}}}} = \pi r_{{\rm{ab}}}^2 $$ (7)

    where the radius rab can be found by solving $D\left({E, {r_{{\rm{ab}}}}, z'} \right) = 0$. The latter gives

    $$ {r_{{\rm{ab}}}}(E,{z^\prime }) = \frac{{\omega ({z^\prime })}}{{\sqrt 2 }}{\left( {{\rm{ln}}\left( {\frac{{2E}}{{\pi {\omega ^2}({z^\prime }){F_{{\rm{th}}}}}}} \right)} \right)^{1/2}}. $$ (8)

    Thus, the base area is

    $$ \begin{array}{*{20}{l}} {{A_{{\rm{ab}}}}(\varepsilon ,{z^\prime }) = \frac{{{\omega ^2}({z^\prime })}}{2}{\rm{ln}}\left( {\frac{{2E}}{{\pi {\omega ^2}({z^\prime }){F_{{\rm{th}}}}}}} \right)}\\ \ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} = \frac{{\omega _0^2}}{2}{{\left( {1 + \left( {\frac{{{z^\prime }}}{{{z_R}}}} \right)} \right)}^2}{\rm{ln}}\left( {\frac{\varepsilon }{{1 + {{\left( {\frac{{{z^\prime }}}{{{z_R}}}} \right)}^2}}}} \right)} \end{array} $$ (9)

    Finally, combining Equations (6) and (9) in Equation (5) we find,

    $$ \begin{array}{*{20}{l}} {{V_{{\rm{ss}}}}(\varepsilon ,{z^\prime }) = \frac{1}{2}{A_{{\rm{ab}}}}{D_{{\rm{ab}}}}}\\ \ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} = \frac{{\pi \omega _0^2L}}{4}\left( {1 + {{\left( {\frac{{{z^\prime }}}{{{z_R}}}} \right)}^2}} \right){{\left( {{\rm{ln}}\left( {\frac{\varepsilon }{{1 + {{(\frac{{{z^\prime }}}{{{z_R}}})}^2}}}} \right)} \right)}^2}} \end{array} $$ (10)

    In Figure 5a–5c, we plot, from left to right Dab, Aab, and Vss as functions of $\mathcal E$ and ${z'}$. We can see from Figure 5c that the maximum Vss is located at the origin of ${z'}$ for smaller values of ɛ, while it splits at higher energy levels. This implies that the maximum ablation rate can be achieved by defocusing at some higher energy level. We find the defocus distance at which the ablated volume Vss has a maximum by taking the derivative with respect to the defocus distance ${z'}$. We provide detailed calculation in the Supplementary Information. The maximum Vss, Vm, occurs at ${{z'}_{\rm{m}}}$, such that

    $$ z_{\rm{m}}^\prime = \left\{ {\begin{array}{*{20}{l}} {0,}&{{\rm{ if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \varepsilon \le {e^2}}\\ { \pm {z_\text{R}}\sqrt {\frac{\varepsilon }{{{e^2}}} - 1} ,}&{{\rm{ if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \varepsilon > {e^2}} \end{array}} \right. $$ (11)

    Fig. 5   

    (a) The central depth of the single-shot ablated hole is plotted as a function of the energy level, ɛ, and defocus distance, ${z'}$. Similarly, the base area of the ablated hole, Aab, is plotted in (b). The volume of the ablated hole, Vss, can then be calculated with knowledge of the central depth and the base area of the ablated hole, as shown in (c). To illustrate the ablated volume profile change, we select four ɛ to represent different cases and plot the corresponding Dab, Aab and Vss in (df), respectively.

    leading to the expressions,

    $$ {V_{\rm{m}}} = \left\{ {\begin{array}{*{20}{l}} {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \frac{{\pi \omega _0^2L}}{4}{{({\rm{ln}}(\varepsilon ))}^2},}&{{\rm{ if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \varepsilon \le {e^2}}\\ {\frac{{\pi \omega _0^2L\varepsilon }}{{{e^2}}},}&{{\rm{ if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \varepsilon > {e^2}} \end{array}} \right. $$ (12)

    We pick four different energy levels for Dab, Aab and Vss and plot them as a function of ${z'}$ in Figure 5d–5f, respectively. We can see that Dab monotonically increases as the energy level increases, while the profile of Aab changes from a unimodal profile to a bimodal profile at energy level ɛ=e, as discussed in our previous work28. Interestingly, when these curves are combined to form Vss, we find that the transition point shifts to ɛ=e2. Above this value, the maximum material removal rate is not located at ${z'}=0$, implying again that a greater ablated volume can be achieved by purposefully defocusing the laser.

    Scanning. So far, we have discussed the case for a fixed focus system (zf=zf0). We now incorporate into our analysis of the effect of the variable focal position introduced by the TAG lens. The focal position varies sinusoidally, as described in Materials and Methods:

    $$ {z_{\rm{f}}} = {z_{{\rm{f0}}}} + {R_{\rm{T}}}{\kern 1pt} {\kern 1pt} {\rm{cos}}({\omega _{\rm{T}}}{t_{\rm{n}}} + \phi ) = {z_{{\rm{f0}}}} + {R_{\rm{T}}}{\kern 1pt} {\kern 1pt} {\rm{cos}}(\theta ) $$ (13)

    tn is the reciprocal of the laser repetition rate, which is constant during the experiment. ωT is the driving frequency and varies by ~1% throughout the experiment, and ϕ is the phase, which takes a new value for each new line. The equation is simplified by grouping terms into θ=ωTtn+ϕ, which accounts for the variations in focal positions.

    Since the ablation rate is a function of ${z'}$, i.e. the distance between the focal position and surface, as discussed in Equation (10), the ablation rate can be calculated by considering the probability distribution (PDF) of the focal position, p(zf), by assuming a fixed surface position. For a given set of parameters, RT and zf0, the PDF of the focal position can be determined by knowing the probability distribution of θ, p(θ), and basic formulae for the probability transformations33. This gives

    $$ p({z_{\rm{f}}};{R_{\rm{T}}},{z_{{\rm{f0}}}}) = \sum\limits_{\theta {\kern 1pt} \in {\kern 1pt} {\kern 1pt} {\rm{co}}{{\rm{s}}^{ - 1}}\frac{{{z_{\rm{f}}} - {z_{{\rm{f0}}}}}}{{{R_{\rm{T}}}}}} {p(\theta )\left| {\frac{{{\rm{d}}\theta }}{{{\rm{d}}{z_{\rm{f}}}}}} \right|} . $$ (14)

    In the case of a statistically representative number of pulses, such as in our experiments, θ will assume all possible values between 0 to 2π with equal probability. Thus, we express the probability of θ as

    $$ p(\theta ) = \left\{ {\begin{array}{*{20}{c}} {\frac{1}{{2\pi }},{\rm{ if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \theta {\kern 1pt} {\kern 1pt} {\kern 1pt} \in {\kern 1pt} {\kern 1pt} {\kern 1pt} [0,2\pi ]}\\ {0,{\rm{ otherwise }}} \end{array}} \right. $$ (15)

    In the range of [−RT, RT], there are two solutions to the equation zf=zf0+RTcos(θ) for θ∈[0, π]: one in [0, π] and the other in [π, 2π]. Thus

    $$ \begin{array}{l} p({z_{\rm{f}}};{R_{\rm{T}}},{z_{{\rm{f0}}}}) = \sum\limits_{\theta {\kern 1pt} \in {\kern 1pt} {\kern 1pt} {\rm{co}}{{\rm{s}}^{ - 1}}\frac{{{z_{\rm{f}}} - {z_{{\rm{f0}}}}}}{{{R_{\rm{T}}}}}} p (\theta )\frac{1}{{{R_{\rm{T}}}\sqrt {1 - {{\left( {\frac{{{z_{\rm{f}}} - {z_{{\rm{f0}}}}}}{{{R_{\rm{T}}}}}} \right)}^2}} }}\\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} = \left\{ {\begin{array}{*{20}{c}} {\frac{2}{{2\pi {\kern 1pt} }}\frac{{{\kern 1pt} 1}}{{{R_{\rm{T}}}\sqrt {1 - {{\left( {\frac{{{z_{\rm{f}}} - {z_{{\rm{f0}}}}}}{{{R_{\rm{T}}}}}} \right)}^2}} }}{\rm{,if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {z_{\rm{f}}} \in {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ { - {R_{\rm{T}}},{R_{\rm{T}}}} \right]}\\ {0,{\rm{otherwise}}} \end{array}} \right. \end{array} $$ (16)

    The PDF of the focal position is plotted in Figure 6a. From Equation (10), we see that the material removal rate is dependent on the defocus distance, $z' = z - {z_{\rm{f}}}$, which is the distance between the focal position and surface. In our model, the focal position is a probability function of the scanning range, as described in Equation (16). Thus, we can calculate the expected material removal rate for a given scanning range through

    $$ {\langle {V_{{\rm{ss}}}}(\varepsilon ,{z^\prime })\rangle _{{R_{\rm{T}}}}} = \int_{ - {R_{\rm{T}}}}^{{R_{\rm{T}}}} p ({z_{\rm{f}}};{R_{\rm{T}}}){V_{{\rm{ss}}}}(\varepsilon ,z - {z_{\rm{f}}}){\rm{d}}{z_{\rm{f}}} $$ (17)

    Fig. 6   

    The PDF of the pulse positions is normalized for each RT and plotted in (a). We then calculate the expected ablated volume, 〈Vss〉, by using Vss and the PDF of the focal position for each scanning range, as plotted in (b). The scanning range that maximizes the ablated volume, optimal RT, is then calculated and plotted in (c). To compare with our silicon experiment, we plot the expected ablated volume as a function of scanning range, which shows the optimal scanning range is at 1.74 mm, as marked by the blue dashed line. This result is close to that used in our experimental conditions.

    If we assume that the sample is placed at the focal point of the fixed focus system (z=zf0, and, therefore zzf=−RTcos(θ)), the resulting plot of 〈Vss〉 at different scanning ranges is given in Figure 6b. The maximum 〈Vss〉 is located at RT=0 in the small energy regime. However, the maximum 〈Vss〉 is shifted to higher scanning range in the large energy regime. This implies that in the small energy regime, z-scanning is unfavorable; however, in the large energy regime, z-scanning should serve to increase the maximum material removal rate.

    We then calculate the scanning range that maximizes the ablated volume at each energy level and plot the results in Figure 6c. Here, we can clearly see that a scanning range of zero is favored when ɛe2, where Vss has a unimodal profile. When ɛ > e2, Vss takes on a bimodal profile, and one can improve the material removal rate through the addition of z-scanning at a well-defined scanning range, as shown in Figure 6c. Interestingly, the optimal scanning range takes the same value as the defocus distance that maximizes the ablated volume, ${{z'}_{\rm{m}}}$, as described in Equation (11). This finding explains why extending the z-scanning range to ${{z'}_{\rm{m}}}$ enables one to achieve an optimal ablation rate at the surface. Since the maximum ablation rate increases linearly with laser energy, the enhancement in the machining rate provided by z-scanning is expected to increase with increasing pulse energy, as described in Equation (12). However, this gain in the machining rate would be unavoidably accompanied by an increase in spot size. Duocastella et al. has reported a loss in lateral resolution when the energy level is larger than e, which is always the case for this method, with this loss increasing with the energy level28. Therefore, one should carefully balance these factors according to the demands of each particular application.

    We now compare the predicted optimal RT with the silicon machining experimental conditions provided in Table 1. The expected ablated volume at different scanning ranges is plotted in Figure 6d, where we can see that the maximum removal rate is achieved at a scanning range of 1.74 mm, as marked by a dashed line in the figure. This is close to the conditions used for which we observed the best experimental results, which shows that our model provides a useful framework for understanding the general concepts at play.

  • So far, we have qualitatively shown that scanning provides a benefit for improving the machining rate. However, quantitative prediction of the ablated volume in silicon is difficult due to the complex ablation mechanism involved29-31. To make a quantitative comparison between theory and experiments, we focus on a model system, Kapton, whose ablation mechanism follows a simple photochemical desorption model34, 35.

    We machine a pocket into Kapton at different defocus positions and scanning ranges, as explained in Materials and methods section and Figure 4a. The profile of the machined pockets at ${z'}=0$ is shown in Figure 7. The laser energy level used is ɛ=e2. We measure the ablated volume over a 100 μm by 100 μm area for each defocus setting. The fixed focus ablated volume as a function of ${z'}$, $V(z')$ is obtained by interpolating experimental data and is shown in Figure 8a. The asymmetric shape of the ablation depth over the defocus distance could be caused by the intersecting beam divergence, as discussed in the previous literature36. The PDF of different scanning ranges is calculated and shown in Figure 8b. Similarly, with knowledge of the ablated volume function V and PDF of the focal positions, we can derive the expected ablated volume value, as described in Equation (17). In each subfigure of Figure 8c, we fix the scanning range and vary the defocus distance, calculating the expected ablated volume for each new experiment (scattered points). We compare this to theoretical predictions derived from Equations (16) and (17) (smooth curve) and find excellent quantitative agreement. Figure 8c also confirms our hypothesis from Silicon machining results section that similar machining efficiencies can be attained over a wide range of ${z'}$ values, especially at a higher lens power. This highlights the other significant benefit of this system in that using a z-scanner not only increases the machining efficiency but also enables an extended machining range in the z direction. We achieve a uniform machining rate over a range of ${z'}$ values by using ultrafast z-scanning and therefore eliminate the need for real-time focus control. Figure 8d compares experimental results with theoretical predictions for multiple passes, which again shows strong quantitative agreement. Here, our predicted curves for multiple passes are found simply by scaling the single pass curve by the number of passes. Since the incremental increase in ablation depth is nearly linear for small numbers of pulses, we expect a simple scaling approximation to be valid for our settings.

    Fig. 7   

    The depth profile of the pocket machined in Kapton using different lens powers. All of these profiles are machined at a defocus distance ${z'}=0$. The unit for the color bar is μm and for the scale bar is 100 μm.

    Fig. 8   

    We machine a pocket into Kapton at different defocus positions. The ablated volume as a function of ${z'}$ is shown in (a). The PDF of the focal positions with lens powers of 0 (m−1), 0.45 (m−1), 0.61 (m−1) and 0.76 (m−1) is normalized and shown in the same order from top to bottom in (b). The calculated and experimental results for the ablated volume function with different lens powers are shown in (c). The calculated data (smooth curve) show good agreement with experiment (scattered points). We further apply linear extrapolation for multiple passes, and find the calculation still shows good agreement with the experimental data, as shown in (d).
CONCLUSION
  • To meet the growing demand for micro-machined products, such as photovoltaic cells, electronic devices, and medical micro-elements, we demonstrate a new high-efficiency laser machining method enabled by an ultrafast z-scanner. Machining efficiency can be derived as the material removal rate function, which is heavily dependent on the relative positions of the focal point and substrate surface. To date, efficiency gains have been limited by the difficulty of pairing real-time surface monitoring with requisite tight focus control. Here, we eliminate the need for this type of control by combining rapid focal position scanning with rapid laser firing. This creates a system in which the focal point is no longer at a fixed position, but is rather described by a probability distribution whose values are dependent on the scanning range. Our model predicts that focus scanning will result in an increased machining efficiency when the laser energy is higher than e2 times the threshold energy and that such improvements will grow as the laser energy further increases. We verify this using the model substrate of Kapton. However, increasing the laser energy for a Gaussian beam can also lead to an increase in the minimum resolution; therefore, one should carefully consider such a trade-off in light of the desired product features.

    We further demonstrate our method using the technologically important substrate of silicon. Based on our experimental conditions, our model predicts an optimal scanning range of 1.74 mm, which is close to the optimal experimental conditions. At these settings, we find a nearly threefold increase in machining efficiency compared to a system without z-scanning. The extended range of focal positions not only increases machining efficiency but also relaxes the requirement of an accurate focus control. Thus, the method presented in this paper significantly enhances the efficiency of existing systems, while opening the door to the machining of non-flat materials found in real-world applications.

CONFLICT OF INTEREST
  • The authors declare no conflict of interest.

ACKNOWLEDGEMENTS
  • We thank TAG Optics Inc. for technical support, and gratefully acknowledge financial support from the NSF (Grant No. CMMI-1235291) and Taiwan Authority of Education. We sincerely thank Alexander Holiday for valuable discussions and feedback.

Reference (36)

Catalog

    /

    DownLoad:  Full-Size Img PowerPoint
    Return
    Return