Advanced laser scanning for highly-efficient ablation and ultrafast surface structuring: experiment and model

Ultra-short laser pulses are frequently used for material removal (ablation) in science, technology and medicine. However, the laser energy is often used inefficiently, thus, leading to low ablation rates. For the efficient ablation of a rectangular shaped cavity, the numerous process parameters such as scanning speed, distance between scanned lines, and spot size on the sample, have to be optimized. Therefore, finding the optimal set of process parameters is always a time-demanding and challenging task. Clear theoretical understanding of the influence of the process parameters on the material removal rate can improve the efficiency of laser energy utilization and enhance the ablation rate. In this work, a new model of rectangular cavity ablation is introduced. The model takes into account the decrease in ablation threshold, as well as saturation of the ablation depth with increasing number of pulses per spot. Scanning electron microscopy and the stylus profilometry were employed to characterize the ablated depth and evaluate the material removal rate. The numerical modelling showed a good agreement with the experimental results. High speed mimicking of bio-inspired functional surfaces by laser irradiation has been demonstrated.

theoretical/experimental analysis of rectangular shaped cavity ablation and formation of multi-layer cavities with two-and-a-half-dimensional (2.5D) shape.
Here, we present the theoretical and experimental studies of the rectangular shaped cavity ablation taking into account the incubation phenomenon and the saturated ablation depth for multi-pulse treatment. The new model has been created for a scanned laser beam on the plane surface and ablation of the rectangular shape cavity. We applied ultra-short laser pulses to ablate the target material at various inter-pulse distances (beams scanning speeds) and laser fluences (spot sizes on sample). Proof-of-principle experiments on copper by various processing parameters sets demonstrate that an optimal point for highest ablation rate can be predicted by our theoretical model.

Materials and Methods
Experimental setup and procedures. The industrial-grade diode-pumped solid-state laser (Atlantic, Ekspla) with a pulse duration of 10 ps emitting at a wavelength of 1064 nm was used in the experiments. The laser provided light pulses with pulse energy up to 130 μJ at a repetition rate of 100 kHz with an average laser power of 13.0 W. The Gaussian beam quality factors M 2 declared by the laser manufacturer was of 1.062 and 1.043 in transverse directions x and y, respectively. The scheme of experimental setup is given in Fig. 1a.
The beam position on the sample was controlled by using a galvanometer scanner (Scangine 14, Scanlab). Translation of the laser spot on the target material at a controllable speed up to 300 mm/s provided the adjustable distance between the transverse irradiation spots Δx. The telecentric f-theta objective lens with the focal length of 80 mm was used to focus the beam on the surface of the target material. The array of rectangular cavities with the transverse spatial dimensions of 2.5 mm × 1.0 mm and the depth-dependant on the processing parameters was ablated (Fig. 1b). The scanning speed of the beam was changed in a horizontal axis. The laser spot size on the sample was changed in the vertical axis of the array by controlling the elevation position z of the sample. The variable spot sizes provided different laser fluences. The path of the scanned beam on the copper sample consisted of parallel lines of overlapped laser pulses as shown in Fig. 1c. The distance between laser pulses was chosen to be identical with the distance between scanned lines as Δx = Δy in order to reduce the number of processing parameters. The ablation rate measurement experiments were conducted by using maximal available laser power of 13.0 W in order to have the maximal material removal rate.
Laser beam characterization. The power meter (Nova ΙΙ, Ophir) equipped with a thermal power sensor (30A-BB-18, Ophir) was used to measure the average laser power. The profiles of the laser beam distribution measured by using the knife-edge method 48,49 are given in Fig. 2.
The experimental data for the measured transmitted laser power P(x) was fitted using the equation 50 , which assumes a Gaussian intensity profile of the beam: x 0 c Δx is the interpulse distance between laser spots equal to the ratio of the beam scanning speed and pulse repetition rate; Δy is the distance between two adjacent scanned lines; w 0 is the laser spot radius on the sample. where P 0 = 13.0 W is the laser power, x c the centre coordinate of the Gaussian intensity profile, and erf is a standard error function. The result of experimental data points fit by Equation (1) are shown in Fig. 2a. The beam radiuses of w x = 2.26 ± 0.02 mm and w y = 2.30 ± 0.03 mm at 1/e 2 level were received for x and y directions, respectively. The first derivatives of the measured transmitted laser power provide the Gaussian beam intensity distributions in the x and y directions as shown in Fig. 2b. It is clearly seen from fits of experimental data points by first derivative of Equation (1) in Fig. 2b that the beam has well defined Gaussian distribution in transverse plane. By performing goodness of fit testing (comparing how well experimental intensity profiles correlates with Gaussian function fits) the reduced chi-squared values of 0.012 and 0.011 were achieved for x and y directions, respectively. Those two values show that experimental data points deviate from Gaussian distribution fits only by 1.2% and 1.1%. Such Gaussian profiles of the beam confirms beam quality factor declared by laser manufacturer. A two-dimensional Gaussian beam at the entrance of the galvanometer scanner exhibits a Gaussian fluence distribution F(x,y) on the target material surface according to: where F 0 is the maximum fluence in the centre, x 0 and y 0 are the coordinates of the central point, w 0 is the Gaussian beam radius. At radius w 0 , by definition F(x, y) decreases to F 0 /e 2 . The dependence of the peak fluence on the single-pulse energy can be written as where E p is the pulse energy. The radius of the transverse focal spot on the target sample surface was measured by the technique described in 51 . Measuring technique provides information about the actual laser beam spot size despite the beam quality factor M 2 . It is valid for a Gaussian beam with a good beam qualify factor of M 2 = 1.1 52 , moderate beam quality factor of M 2 = 1.3 [53][54][55] , and even for multimode beam with a poor beam quality factor of M 2 = 95 56 . Also, this measuring technique is valid for a highly elliptical Gaussian beam 57 . Assuming that a laser beam has the Gaussian spatial beam profile, the relation between the crater diameter D and the peak laser fluence F 0 in the centre of the Gaussian beam can be written as 51 : where w 0 is the radius of the Gaussian beam on the sample, F th (N) is the ablation threshold, N is the number of laser pulses applied to a single spot. The optical microscopy was employed for the characterization of ablated craters (dimples). Diameters and depths of craters ablated by the single and multi-pulse regimes in copper for determination of ablation threshold, incubation parameter and penetration depth have been performed by using an optical microscope (Eclipse LV100, Nikon) equipped with high-definition 5-megapixel CCD digital camera (DS-Fi1, Nikon) with a resolution of 2560 × 1920. The camera was equipped with a controller (Digital Sight DS-U2, Nikon) and image processing software (NIS-Elements D, Nikon). The objective (LU Plan Fluor 20x, Nikon) with a numerical aperture (NA) = 0.5 and a magnification factor of 20X was used in bright field. The specimen was illuminated by a halogen lamp (LV-HL50PC, Nikon). The pixel size of the digital microscope image was ≈0.2 µm, which was much smaller than the variation of crater size from different craters ablated under the same experimental conditions. Therefore, the standard deviation of five ablated craters was taken as a measurement error. By having Equation (4) with the pulse energy values, an ablation threshold pulse energy E th is calculated from a semi-logarithmic plot of the diameter squared of the ablated area D 2 versus pulse energy E p . The slope  (3) and (4). The spot radius dependence on the sample vertical position is given in Fig. 3b. The beam propagation equation along the z-direction for spot radius w(z) of non-ideal Gaussian beam can be written as 58 : where w 0 is the spot size in focus, λ = 1064 nm is the wavelength of irradiation, z 0 is the focal position, M 2 is the beam quality parameter. From the experimental data fit by Equation (5) in Fig. 3b the Gaussian beam quality factor of M 2 = 1.06 ± 0.01 was retrieved. The received beam quality parameter is in good agreement with the data provided by the laser manufacturer.
Ablation rate characterization. The solid copper (CW004A, Ekstremalė) target plate with the dimensions of 50 mm × 50 mm × 5 mm, purity of 99.9%, and surface roughness of R a < 0.1 µm (mirror finish) was used in the ablation tests. The surface topography was studied by using a scanning electron microscope (SEM) (JSM-6490LV, JEOL). The SEM micrograph of an array of laser ablated rectangular shaped cavities is shown in Fig. 4a. The profiles of laser-ablated cavities were measured by a stylus profiler (Dektak 150, Veeco) with the measurement resolution of 0.1 μm, 1 μm and 1 nm in x, y and z directions, respectively. The typical profile is shown  in Fig. 4b. The ablation rate was evaluated from the volume of the ablated cavities using data from the profiles. The measuring errors of ablation depth were related to the surface roughness of the cavity bottom. Therefore, the average relative error of ablation rate was less than 2.5%.
Ablation model of rectangular shaped cavity. The laser ablation model of rectangular shaped cavity including the incubation phenomenon and the saturated ablation depth for multi-pulse processing has been proposed in this work. The graphical representation of a theoretical model in the case of rectangular shaped cavity ablation is given in Fig. 5. The laser beam is scanned at a certain speed in parallel lines over the copper target, and the rectangular shape cavity is ablated (Fig. 5a). It can be assumed that a single laser pulse ablates the primitive cell marked in a black dashed line ( Fig. 5a) with the volume lateral dimensions of Δx and Δy and height of h: The primitive cell is an imaginative concept of the laser ablation volume per single pulse. The volume of the primitive cell is equal to the volume of the rectangular cavity divided by the total number of pulses applied. However, the depth of the primitive cell is influenced by many surrounding pulses within the cut-off radius. The beam scanning path consisting of the parallel scanned lines of the laser beam is shown (Fig. 5b). The three main areas are marked by Roman numerals, and they represent the irradiation points of a different origin: (I) is a selected central point where ablation volume per single pulse is calculated; (II) is any point in the active area pulses have influence in the ablation depth at the (I) position with a non-zero ablation depth; (III) the irradiation points outside the active area laser pulses have zero influence in the ablation depth at the selected point of interest (I).
The three-dimensional (3D) representation of a primitive cell of ablated volume per pulse is shown in Fig. 5c. The symbols Δh I , Δh II and Δh III represent the ablated depths of three points of interests in Fig. 5b (I), (II) and (III), respectively. The overall ablated depth consists of the ablated depths from pulses from the active area. The total ablation depth per pulse is evaluated by summing all the ablation depths of surrounding irradiated spots in x and y directions which are in the range of a critical radius r c : x y The ablation rate can be assumed as p where f p = 100 kHz is the laser pulse repetition rate. The main purpose of the proposed model is to calculate the ablated volume per pulse of the scanned laser beam by taking into account the incubation phenomenon and saturation of ablation depth for increasing number of pulses per spot. The results of ablation using picosecond laser treatment at different numbers of pulses per spot (N = 1, 10, 100, and 1000) are presented in Fig. 6a. The slope of the straight lines provides information about the Gaussian beam radius w 0 = 41.8 ± 2.1 µm on the target material from Equation (4). With the known beam radius, the peak laser fluence values were calculated by using Equation (3). The parallel lines represent the fixed slopes because of the same spot size on the sample. However, different modification thresholds are observed due to the incubation phenomena. The extrapolation of lines at D 2 = 0 μm gives the ablation threshold fluences: F th (1) = 2.0 ± 0.35 J/cm 2 ; F th (10) = 0.87 ± 0.21 J/cm 2 ; F th (100) = 0.42 ± 0.08 J/cm 2 ; F th (1000) = 0.19 ± 0.01 J/cm 2 . The copper ablation thresholds correspond well to the multi-pulse ablation thresholds F th (N) = 1.73, 0.74, 0.50, and 0.33 J/cm 2 for N = 1, 10, 100, and 1000 obtained in our previous work (1064 nm, 10 ps, 50 kHz) 27 . The incubation model describes the relation between the single-pulse ablation threshold F th (1) and the multi-pulse ablation threshold F th (N) in form 39 : S th th 1 with a material-dependent incubation parameter S and the laser fluence sufficient to ablate the material for a fixed number of laser pulses per spot N is named the ablation threshold fluence F th (N) for a multi-pulse treatment. The accumulated laser fluence N × F th (N) versus a number of pulses N for copper are given in Fig. 6b. A calculation by Equation (9) 60 . The accumulation coefficient measured in our work coincides well with the results from other groups. For photochemical as well as for photo-thermal ablation and many materials, a relation of the ablation depth per pulse can be described as 2 : where δ is the energy penetration depth. If many pulses N are applied per single spot, ablation depth h(N) increases up to a material and laser dependent pulse number N 0 and then reaches saturation 39,[43][44][45]47 . However, an analytical expression is not found in the scientific literature for a prediction of ablation depth growth and saturation with increasing number of pulses per spot. In this work, the analytical function of ablated dimple depth for multi-pulse irradiation has been empirically constructed:   (11) has an advantage of having two parameters with defined physical meanings required to describe the saturation phenomenon. Also, it has an analytical form which is advantageous when used in numerical simulations and it is differentiable. The first derivative of softplus is the sigmoid function. The function with initial linear increase has been chosen because straight-line increase in ablation depth versus number of pulses has been experimentally observed in numerous scientific works [62][63][64][65][66][67] . The saturation after certain number of laser pulses has been also observed experimentally 43,45,46 . Therefore, the empirically selected function (11) fulfilled our needs of a fully analytical expression to be later used for approximation of experimental data, derivation of saturation parameters, and calculation of ablation rates in our numerical modelling.
The dimples on a copper target were ablated by using a Gaussian beam by controlling the number of incident pulses per spot. The dimple depth has been experimentally measured by using an optical microscope (Eclipse LV100, Nikon) and the focus knob technique described in 68 . The method is based by the position difference measurement of the top and bottom of the object, using the markings on the fine-focus knob of the microscope. The optical microscope was equipped with xyz stage (LV-S64, Nikon) and a 100X microscope objective (LU Plan Fluor 100X, Nikon) with NA = 0.9 and depth of field of 0.19 µm have been used. The z accuracy of the fine-focus knob was 1 µm. The depth of field of the objective was much smaller than the accuracy of the fine-focus knob. Therefore, measurement errors of dimple depth were determined by the accuracy of the fine-focus knob. The depth of ablated dimple dependence on the laser pulses per spot is given in Fig. 7a.
The ablated dimple depth increases with the number of pulses applied have been tested using three values of peak laser fluences and similar saturation behaviour has been observed for all fluences (Fig. 7a). The fit of experimental data points by Equation (11) provides: saturation value of N 0 ≈ (1.0 ± 0.1) × 10 3 , the saturation softness of ΔN ≈ (2.0 ± 0.2) × 10 2 pulses, and penetration depth for single pulse irradiation of δ(1) ≈ 50 ± 5 nm (Fig. 7a). The fit has a good agreement with the experimental data (Fig. 7a). For small irradiation doses <10 4 J/ cm 2 , the depth of the dimple increases linearly with an increasing number of laser pulses. However, after reaching the saturation value of >10 3 pulses, the depth of the dimples stops growing and saturates. By differentiating the Equation (11) for a given laser fluence by the number of pulses, the ablated depth per pulse dependence on pulse number is achieved: 1 exp ln , The function (12) has constant value for the initial laser pulses N < N 0 − 3ΔN, followed by the decrease to zero in the range N 0 − 3ΔN ≤ N ≤ N 0 + 3ΔN, and a zero value for the pulse number exceeding N > N 0 + 3ΔN. This function is known as a sigmoid function 69 or logistic function and becomes a step or Heaviside function with ΔN = 0. The sigmoid functions primitive is the softplus function 61 . The function coincides with scientific works, where constant ablation depth per pulse has been observed experimentally for an initial pulse number 64,70 .
By calculating the first derivative of data points of Fig. 7a, the depth per pulse dependence on the pulse number can be achieved (Fig. 7b). For a large number of laser pulses >10 3 , the ablated depth per pulse has a value close to zero. The data points are scattered from ~50 to ~125 µm. However, its average value stays constant at around Δh(300 < N) ≈ 60 µm. The reason for the high scattering of the data points is the large relative errors of the focus-knob method for small depths of <20 µm 71 (Fig. 7a, N < 300) and consequently amplified scattering of its first derivative (Fig. 7b, N < 300). However, despite a large scattering of the experimental data points and large measurement errors, the fit of experimental data points by Equation (12) has a good agreement (Fig. 7b). The penetration depth reduction and saturation behaviour with an increase of number of laser pulses is mostly related to two physical aspects in the keyhole evolution: recast layer formation and multiple reflections. The redeposition of the ablated material and recast layer formation in laser drilling has already been explained in great detail in ref. 72 . Geometrical aspects of deep dimple ablation and multiple reflections influence in keyhole evolution has already been investigated in refs 73,74 and refs [75][76][77][78] , respectively. Therefore, we did not want to go into the physics of already well-known aspects of multi-pulse ablation. Our main goal of this research was to show the responsible combination of physical effects related to the laser ablation and their influence on the material removal rate of the rectangular shape cavity formation. Our approach was to directly compare experimental data with simulation results of the proposed model. However, usage of saturation equation (11) for rectangular cavity ablation with a scanned laser beam, which was derived for dimple ablated with fixed beam position, is not straight forward. Therefore, we provide a graphical scheme and discussion with emphasised similarities between the fixed and scanned beam ablation.
The validity conditions of saturation function (11) usage for scanned laser beam which was derived for fixed position ablation is discussed below. There are several known physical aspects related to the saturation behaviour of depth using fixed laser beam position in multi-pulse irradiation, but the most important one is the multiple-reflections of the laser beam from the internal walls of the ablated dimple 46,75,76,[79][80][81] . The Fresnel reflection is drastically increased because of an increase in the incidence angle as shown in Fig. 8a.
The increased reflectivity results in a decrease of absorbed energy. Also, an internal surface area of the dimple is larger than an irradiated disk on a flat surface in the initial stage of ablation. The increased beam-matter interaction surface area causes the decreased energy density on the internal surface walls of a dimple (Fig. 8a). As the result, the dimple depth stops increasing for the following laser pulses and the saturation starts. With the fixed laser beam position, the beam interacts with interior walls of the parabolic dimple (Fig. 8a). If compared to the scanned laser beam and grove ablation with the shape of the parabolic cylinder depicted in Fig. 8b, the direct beam interaction area is about a half of the internal paraboloid area shown in Fig. 8a. In the rectangular shaped cavity formation by bidirectional scanning and hatching the beam interaction area is even smaller as illustrated in Fig. 8c. However, the same effect of increase in incidence angle, its related increase of Fresnel reflection, and decrease of absorption are observed in all three situations. The effect of decrease of energy density on the ablated surface is common for the fixed and scanned beam ablation. Only the magnitude of the effect is different, which is related to the geometrical aspects of ablated cavities, for all three cases. The similarities presented in Fig. 8 justifies that equation (11) derived for a fixed position is valid for a scanned laser beam, with some tolerance. Therefore, having in mind that fixed position and scanned beam ablation are different procedures, only having similarities in saturation behaviour, the formula (11) will be used for numerical simulation in our model.
The model incorporated a numerical calculation of an equation system with the input of the copper ablation parameters. All laser ablation parameters required for the proposed ablation model of the rectangular shape cavity are given in Table 1. For the selected point of interest (I) (Fig. 5b) the active circle area πr c 2 with the cut-off radius r c . All the laser pulses within this radius have non-zero ablation depths in position (I) (Fig. 5b). Ablation will not be possible by the laser pulses, which centre positions x and y are outside of cut-off radius as (x − x 0 ) 2 + (y − y 0 ) 2 > r c 2 , because laser fluence from the Gaussian wings will be below the threshold for multi-pulse ablation. The number of all those pulses N within the cut-off radius can be evaluated as a ratio between the active circle area and the primitive cell base area: where Δx and Δy are the distance of the irradiation spots in the transverse directions x and y, respectively. On the other hand, the cut-off r c of the ablated area evaluated by taking into account that in the processed area N pulses are absorbed, and the ablation threshold is decreased. From Equations (4) and (9): S c 0 0 th 1 By numerically solving Equations (13) and (14) the cut-off radius and the number of laser pulses per spot can be evaluated. The inter-pulse distance Δx and the distance between scanned lines Δy have been chosen with equal values in order to reduce the number of laser processing parameters. Moreover, they can be replaced by an experimentally controllable processing parameter -a beam scanning speed v on the sample: p where f p is the pulse repetition rate. Therefore, the main well known equations, defining the laser ablation by a Gaussian laser beam were included into this model: diameter squared of the ablated crater dependence on the peak laser fluence by Equation (4); incubation behaviour for multi-pulse treatment by Equation (9); ablation depth dependence on the laser fluence by Equation (10); saturation of ablation depth by Equation (11).

Results
Experimental and modelling results of rectangular cavity ablation. The ablation rate has been tested experimentally and theoretically at various beam scanning speeds (inter-pulse delays) and spot radiuses on the sample (peak laser fluences) (for details see Materials and Methods). The ultrafast laser with pulse duration of 10 ps, irradiation wavelength of 1064 nm, pulse energy of 130 μJ, and pulse repetition rate of 100 kHz was used to ablate an array of rectangular areas in the copper plate (for details see Materials and Methods), and ablation rate was measured (Fig. 9a). Equations (7-9 and 12-15) have been solved numerically by using the measured ablation characteristics of copper from Table 1 and symbolic/numeric computing software (Maple, Maplesoft). The calculated ablation rate versus laser processing parameters is given in (Fig. 9b). The experimental data of measured ablation rates was quantitatively compared to the numerical results of the model equations in Fig. 9c. The series of cross sections from the experiment and model with variable peak laser fluence (from 1.20 J/cm 2 to 18.4 J/cm 2 ) and laser spot radius (from 21.2 μ to 83.0 μm) are shown as a function of inter-pulse distance. Our simulation results at each value of the laser fluence were in good agreement with the ablation rate measurements. The experimental data points can be well described by the modelling curves that take into account the incubation effect and saturation of ablation depth for multi-pulse irradiation. At a low scan speed v < 100 mm/s, a high number of pulses per spot is achieved because of the high pulse overlap. Therefore, the ablation threshold decreases drastically, and its related total ablated depth increases. However, the saturation of ablation depth is achieved, and the ablated depth is not increased any more. The ablation rate reaches moderate values and decreases even more by reducing the scanning speed. Meanwhile, with high scanning speed v > 280 mm/s, the number of laser pulses per spot becomes small because of the low pulse overlap. Therefore, the ablation threshold increases due to the accumulation effect and it causes a small ablation depth. Contrary, the large ablated transverse area for a single pulse is caused by a large inter-pulse distance. The product of a small total ablation depth and large transverse area gives only moderate ablation rate.
For small laser spot sizes w 0 < 40 μm, the laser fluence high above the ablation threshold is achieved, which should result in large ablation depths. However, low numbers of pulses per spot are achieved even for small inter-pulse distances. Therefore, low ablation rates are achieved because of high ablation thresholds. For high laser spot sizes w 0 > 75 μm, it is difficult to reach fluence above the ablation threshold, so low ablation depths are achieved, as well. However, the number of laser pulses is sufficiently high even for large scanning speed because of the high pulse overlap. Thus, low and moderate ablation rates are achieved. The optimal laser parameter set exists for maximum available ablation rate: scanning speed 180 mm/s < v < 220 mm/s and beam spot radius 60 μm < w 0 < 70 μm (Fig. 9). In the optimal set up the maximal ablated volume per pulse is reached because of a combination of all the factors influencing material removal rate. The moderate scanning speeds and its related distance between adjacent pulses give sufficiently large transverse ablation area for a single pulse. Also, the moderated spot size gives the sufficiently high number of laser pulses and its related appreciable decrease of ablation threshold due to the accumulation effects. Volume pulses per unit point are less, so the result is a high ablation threshold, which decreases increasing the number of pulses. The saturation depth is just reached in this processing region. Therefore, the ablation rate has a maximum at scanning speed of v max ≈ 200 mm/s and a laser spot radius of w 0 max ≈ 62.4 μm (Fig. 9). The experiment data results have a good agreement with the numerical calculation results of the new model.
By having laser pulse energy E p = 130 μJ and spot radius w 0 = 62.4 μm the calculated peak laser fluence Equation (3) with laser pulse energy is F 0 = 2.13 J/cm 2 . By solving Equations (13) and (14) using optimal inter-pulse distances between scanned spots and scanned lines of Δx = Δy = 2.0 μm, one can calculate the cut-off radius r c ≈ 62.7 μm and a number of pulses per spot N ≈ 2.5 × 10 3 pulses. The threshold for that number of pulses The ratio between the optimal laser fluence and the threshold is 7.6. This number is close to e 2 ≈ 7.4 which is theoretically predicted 26,27 . The experimental measurement of ablation rate was based on the rectangular shape cavity profiles that were achieved by a stylus profiler and average depth evaluation. The measurement errors were related to the bottom surface roughness of an ablated cavity. The average relative error of a cavity depth and its related ablation rate measurement was of 2.5%. The peak laser fluence can be calculated by Equation (3) and is proportional to the laser power and inversely proportional to the spot radius squared. The laser power was measured by a power measurement unit (for details see Laser beam characterization, Materials and Methods) with a relative error of 0.7%. The beam radiuses on the sample at each height of the specimen were measured by a technique described in 51 (Fig. 3). The average relative error of the beam radius measurement was 2.7%. The scanning speed (inter-pulse distance) and its related pulse overlap were controlled by a galvanometer scanner with an accuracy of 1.0% 82,83 . The ablation rate model with a separate variation of laser fluence, beam radius and scanning speed (inter-pulse distance) by 0.7%, 2.7% and 1.0% gave the highest relative errors of 0.4%, 1.5% and 0.8%, respectively. The ablation model errors were always comparable with the experimental errors. The largest relative errors of the model calculations of ablation rate of 1.5% were achieved when spot radius was varied. The spot size variation influenced the calculations mostly because it has the direct influence on several physical parameters: peak laser fluence by Equation (3) and cut-off radius by Equation (14). Moreover, peak laser fluence is directly involved in the ablation depth per pulse calculations by Equation (12). The largest relative variation of the calculated ablation rate of 1.5% is smaller than the experimental error of 2.5%. Efficient deep laser engraving and structuring. One of the aims of the study was to show that the laser ablation is capable of 2.5D structure formation and replication of bio-inspired structure at high speed. The copper was structured by using pulsed laser irradiation with picosecond laser duration (for details see Experimental setup and procedures, Materials and Methods). The examples of deep laser engraving and 2.5D structuring are given in Fig. 10.
The laser ablated cavity with the removal rate of 330 μm 3 /pulse has been performed by using 130 μJ laser pulse energy Fig. 10a. The current world record of ablation cooled process by bursts of laser pulses is ~600 μm 3 /pulse for 130 μJ total incident energy of 200 pulses in the burst at 864 MHz 3 . The achieved removal rate of efficient ablation in our experiments for copper is ~2 times slower than the record obtained with a specially designed GHz burst laser 3 . The efficient ablation laser parameter set has been used to create a rectangular cavity (Fig. 10a), the logo of Center for Physical Sciences and Technology (Fig. 10b), a bio-inspired shark skin structure (Fig. 10c) and a drag reducing blade-riblet structure (Fig. 10d) in the copper sample. The ablation speed of 2.0 mm 3 /min has been achieved at an average laser power of 13.0 W.

Discussion and Conclusions
Scientific background of this research was two previous works 26, 27 that emerged a decade ago attempting to predict efficient laser ablation. However, the accumulation effect and saturation of ablation depth for multi-pulse ablation was not taken into account. In this research, the complete model of laser ablation incorporating the accumulation effect and the saturation of ablation depth, if many pulses are applied to a single spot, is presented. Laser irradiation with picosecond pulse duration has been used for the ablation of the target material of copper. The processing parameters: scanning speed (inter-pulse distance, distance between scanned lines) and spot size on the sample (peak laser fluence) has been varied in order to reach maximum removal rate. The new proposed model of rectangular shape cavity ablation has been created which takes into account the ablation threshold and the penetration depth decrease with increasing number of laser pulses per spot. The model also assumes that laser beam has a Gaussian intensity profile, the ablation depth per pulse is proportional to the logarithm of fluence applied with a certain ablation threshold, the beam is scanned on a flat target material surface by parallel overlapping lines, the inter-pulse distance is equal to the distance between scanned lines, and the ablation volume per pulse is calculated by summing ablated depths from surrounding pulses within cut-off radius. The experiments have been conducted for the ablation rate measurements and confirmation of modelling results. The experimental evidence of the model is clearly seen in Fig. 9 as the numerical simulation of the model equations is in good agreement with experimental results. Thus, the modelling results are confirmed by the experimental data. The maximum removal rate of 2.0 mm 3 /min (330 μm 3 /pulse) has been achieved for the copper target by using the laser processing parameters: wavelength of irradiation 1064 nm; mean laser power 13.0 W; pulse repetition rate 100 kHz; pulse duration 10 ps; inter-pulse distance 2.0 µm (scanning speed 200 mm/s); peak laser fluence 2.13 J/ cm 2 (spot radius w 0 = 62.4 μm). From a scientific point of view, it is worth mentioning that the incorporation of the accumulation effect and saturation of ablation depth for multi-pulse ablation to laser ablation model has been demonstrated for the first time. Comparing our model to previous attempts of predicting the efficient ablation point, where only logarithm dependence on the fluence of ablation depth and diameter squared for ablation by Gaussian beam has been included, our model is much more complete. From an engineering point of view, it is essential to have a theoretical tool and laser scanning procedure to find the most efficient ablation spot. In this work, theoretical model and laser scanning procedure are presented enabling high texturing speed. The possibility to mimic the bio-inspired functional surfaces at high processing rate has been demonstrated. The results of laser application contribute to technological conclusions, that controlled laser irradiation is a reliable tool to mimic functional surfaces found in nature.