Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Methodology of selective metallic thin film ablation from susceptible polymer substrate using pulsed femtosecond laser

Open Access Open Access

Abstract

Electronic devices on flexible polymeric substrates allow new fields of applications. A maskless and flexible structuring process for such systems is offered by ablation using ultra-short pulse laser irradiation. Hereby, certain areas of a functional thin film coating (e.g. nickel-chromium) are locally removed from a substrate (e.g. polyimide) to yield the needed device structures. Micro laser patterning quality is influenced by the beam properties (beam profile, fluence) as well as by the pulse overlap, the substrate material and many other factors. A clear distinction must be made between the material ablation at the surface of a bulk material and the substrate selective removal of a thin metallic film. For the latter, general rules for the prediction of ablation results especially in the case of areal ablation, which were not known from the literature so far, are derived here in the form of mathematical criteria. A methodology for the parameter finding in different cases of ablation (dot, line, areal) is presented and exemplified using a practical example, but is also applicable to other flexible thin film based systems.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Nowadays, metallic thin films and coatings are an essential part of many advanced systems but also everyday devices [1]. Electronic devices, in former times solely based on solid printed circuit boards, are altered such that their functional structures can be fabricated on flexible substrates [2]. This allows for new fields of applications e.g. live science or medical products; as flexible systems can be conformably applied and integrated into all kinds of soft materials. Flexible substrates are often made from polymers. Most polymer foils are available in transparent or opaque form and with a wide range of mechanical properties. In addition, they show good dielectric properties and are relatively cheap to produce. For the fabrication of electronic devices mainly thermal and chemical robustness of the chosen polymer substrate have to be considered, because many microelectronic production processes are based on lithography and chemical etching [3]. The use of etchants, solvents or heat to cure photosensitive resins can reduce the selection of suitable substrate materials. As an alternative, the ablation by ultra-short pulse laser irradiation provides a maskless and very flexible process for structuring thin coatings on flexible polymeric substrates [4]. Due to its properties, this technology is used in industrial applications such as the production of touch screens, OLED systems, RFID tags, sensors and for EM shielding tasks [5]. Since very large areas can be processed, the structuring of conductive transistor electrodes for solar cell fabrication is another important area of application [6]. In general, highly energetic ultrashort pulses can be used to ablate material with marginal thermal input to the substrate, whereby the heat-affected volume in the vicinity of the ablation area is minimized [1,6].

The quality of micro laser patterning can be controlled via pulse overlap, fluence and scan repetitions [1]. For flexible electronic systems, usually certain areas of a conductive thin film coating are removed to yield the needed functional (sensor-) structures. Thus, complete galvanic insulation between separated structures has to be accomplished [7]. Precise movement of the laser beam along vertically stacked trajectory patterns allows areal ablation with pulsed laser systems. Despite the fact that single pulse ablation thresholds can be precisely determined with known techniques, the influence of scan parameters for line- and areal ablation is less investigated. Instead, trial-and-error parameter grids are mostly used in the field. Of course, these can also lead to a useful set of parameters whereby the ablation results are often evaluated according to rather subjective quality standards. However, the avoidance of substrate damage is crucial and is therefore regarded as the main optimization parameter in the context of this paper. For polyimide as a substrate, it can be shown that even minor substrate damage can be detected by purely optical measuring methods. The known methodology of ablation threshold determination evolved from rather bulky materials with the goal of cutting down into a substrate to a certain depth [8]. However, differences between thin film ablation and bulk material removal through laser induced evaporation of (bulk-) material cannot be ruled out [9]. Thus, methods and criteria allowing predictions for selective ablation results are still unclear or even missing.

Ideal selective ablation is achieved, when the thin film coating is removed without residues while the substrate stays untouched. This is especially challenging for metal coatings on susceptible polymer substrates. Because of this difficulty, authors in the past have considered small substrate damage (e.g. $1{\; }\mathrm{\mu}\textrm{m}$ abrasion from a $50{\; }\mathrm{\mu}\textrm{m}$ polymer foil after ablating a $150\;\textrm{nm}$ metallic thin film) inevitable to achieve proper electrical insulation [10]. As this is justified for some applications, it might limit the use of laser ablation processes when areal removal on transparent substrates should be accomplished (e.g. for solar cells). Even a slightly damaged surface may lead to a foggy, less transparent substrate here. Also, mechanical problems can arise from weakened substrates. To evade these limitations we will present a systematic methodology for laser parameter identification. As a result a selective ablation can be achieved, that can no longer be distinguished from the ideal selective case by optical methods. Our parameter selection approach is based on finding the sample/material specific ablation fluence threshold ${\phi _{\textrm{th}}}$ via multiple single spot D2-experiments. Two criteria will be introduced to predict ablation results such that the parameter search can be narrowed down beforehand.

2. Materials & methods

We chose the experimental setup, the inspection equipment, the sample fabrication method as well as the materials used such that results are application oriented and are easily reproducible at any time. To avoid influences caused by a wavy substrate foil, a spin coated polymer layer on a rigid carrier wafer was used.

2.1 Laser and inspection equipment

All laser based experiments were carried out using a pulsed Yb:KGW (Ytterbium-doped Potassium-Gadolinium Tungstate) solid-state laser system (microSTRUCT C, 3DMicromac), with a pulsed laser source ($\nu = 212\; \textrm{fs}$) emitting at a primary wavelength of $\mathrm{\lambda } = 1030\; \textrm{nm}$ in linear, horizontal polarization. Pulse frequency was 600 kHz for most ablation experiments. Beam quality is specified with $\textrm{M}^2\; \textrm{} < \; \textrm{}1,1$. The laser beam has a diameter of of $\; \textrm{}{\emptyset _{\textrm{Beam}}} = 4,7\; \textrm{mm}$ and can further be expanded using a Galilei lens setup. High speed laser beam positioning on the specimen is provided by a galvanometer-scanner (Scanlab intelliSCAN 14), while focusing is accomplished by a f-theta telecentric lens (Linos F-Theta Ronar) with an effective focal length of $100,1\; \textrm{mm}$ and a apodisation factor of $1,83$. The laser source was operated at a reduced power of 3W providing pulse energies between $0,28$ and $4,85\; \textrm{}\mathrm{\mu}\textrm{J}$. All testing was carried out under normal ambient conditions at room temperature.

Analysis of stationary spot experiments as well as inspection and evaluation of ablation quality was conducted with a confocal laser scanning microscope (CLSM, Keyence, VK-X260K). Its microscope objective provides a magnification of $150\textrm{x}$ at a numerical aperture of $0,95$ which enables morphological examination of the surface. In addition to the 3D height profile, a DIC (differential interference contrast) imaging was used to increase the contrast of smallest contours and edges on the surface. This technique creates a relief-like appearance to recognize otherwise invisible surface features.

2.2 Sample preparation

As mentioned before, we used polyimide substrates on rigid class carrier wafers for better handling. Polyimide foils were prepared by spin coating a liquid precursor (PI-2610/11, HD microsystems) onto the carrier as described by [11]. An adjacent thermal treatment was used for curing. Prior to sputter deposition, the foil surface was treated by inverse sputter etching to improve adhesion of the coating. A layer of $100\; \textrm{nm}$ of nickel-chromium (target composition: 80 wt% nickel and 20 wt% chromium) as piezoresistive metallic thin film alloy commonly used for strain gauge applications [12], was deposited on the substrate.

3. Calculation of fluence distributions

An areal ablation process with a pulsed laser source is accomplished by moving the beam in a scan pattern over the surface of the specimen. A combination of scanner speed v, pulse frequency f and spot radius ${\omega _0}$ determines the resulting pulse overlap parameter $\varphi $ (complete overlap at $\varphi = 0$ and no overlap for $\varphi > 1$):

$$\varphi = \frac{{{l_P}}}{{2{\omega _0}}}\quad {l_P} = \frac{v}{f}.$$
In literature, pulse overlaps of at least 90% ($\varphi = 0,1$) are commonly used to reduce edge line roughness. In this paper, we also work with much smaller overlaps as these will show to be advantageous when low average laser energy per ablation area is required. The real laser spot size can either be measured (by knife-edge method), or calculated from the specified beam width ${\emptyset _{\textrm{Beam}}}$ leaving the laser source. Our laser system provides a beam diameter of ${\emptyset _{\textrm{Beam}}} = 4,7\; \textrm{mm}$ measured with a beam profiler (Ophir Spiricon SP928) using the D4σ-standard. The collimated beam goes into the scanner unit and exits through a telecentric f-theta objective with an effective focal length ($EFL$) of $100,1\; \textrm{mm}$ and an apodizationfactor of 1,83. This leads to a spot size ${\emptyset _0}$ of:
$${\emptyset _{\textrm{0,calc}}} = 2{\omega _0} = \frac{{1,83 \times \lambda \times EFL}}{{{\emptyset _{\textrm{Beam}}}}} = 40,1\,\mathrm{\mu}\textrm{m}.$$
By measuring the laserpower $\textrm{P}$ with a thermal detector, the resulting pulse energy ${E_p}$ and maximum spot fluence in the beam center ${\phi _0}$ can be calculated as:
$${\phi _0} = \frac{{2 \times P}}{{\mathrm{\pi } \times \omega _0^2 \times f}} = \frac{{2 \times {E_P}}}{{\mathrm{\pi } \times \omega _0^2}}.$$
In our setup, the fluence profile along the radius r follows a spatial Gaussian distribution:
$$\phi (r )= {\phi _0} \times {e^{ - 2\;{{\left( {\frac{r}{{{\omega_0}}}} \right)}^2}}}.$$
Definitions such as FWHM, 4Dσ or 1/e2 coexist to specify the laser spotsize ${\emptyset _0} = 2{\omega _0}$. Here we use the latter two definitions which state that the border of the laserspot along the radius r in lateral direction is reached when the spot fluence has dropped to ${\phi _0}/\textrm{e}^2$. In the following ${\phi _0}$ will simply be called „fluence“. In Fig. 1, a gaussian fluence distribution is sketched.

 figure: Fig. 1.

Fig. 1. Gaussian fluence profile

Download Full Size | PDF

After scanning overlapping pulses along a straight line in $\textrm{X}$-direction, the resulting ablated cut width D can be thinner than the spotsize ${\emptyset _0}$ and depends on the specific material properties of the substrate.

In overlapping regions of two or more subsequent spots the fluence accumulates. For single laser lines the accumulated fluence profile can be gained by summing up the fluence distributions of individual pulses along its trajectory. The calculation of accumulated fluence profiles for areal ablation is more challenging, as the chosen ablation pattern has to be taken into account. We choose to work with a line pitch ${l_{P,y}}$ ($\textrm{Y}$-direction) of the stacked horizontal laser trajectories equal to the pulse spacing ${l_{P,x}}$ in the line scan direction $\textrm{X}$ (see Fig. 1: Gaussian fluence profile Fig. 2). This ablation pattern (${l_P} = {l_{P,x}} = {l_{P,y}}$) will subsequently be called “#-pattern”. It leads to an isotropic energy distribution.

 figure: Fig. 2.

Fig. 2. Pulse placement for an isotropic energy distribution. The red circles indicate the spot size of single pulses with their ${\omega _0}$ boundaries. Grey lines indicate the order in which ablation pattern is progressing (starting at the abscissa moving up, left to right).

Download Full Size | PDF

For illustration purposes, “uncovered areas” that are not irradiated (when ${\omega _0}$ is interpreted as setting a sharp irradiation boundary) are marked in Fig. 2: Gaussian fluence profile

Figure 2. In all our experiments the pulse overlap is chosen, such that full coverage of the ablation zone is achieved $0 < \varphi < 1\; \textrm{}$in order to avoid these uncovered regions.

3.1 Accumulated fluence profile (AFP) as an irradiation model

Eichstädt et al. have introduced an accumulated fluence profile (AFP) which is valid for laser engraved lines as well as spatially extended areas [13]. Thus, the fluence term includes x,y-coordinates shifted by multiples (${n_x},{n_y}$) of the distance between pulses (${l_P} = {l_{P,x}} = {l_{P,y}}$):

$$\begin{array}{l} \phi ({x,y,{n_x},{n_y}} )= {\phi _0} \times {\textrm{e}^{ - 2\left( {\frac{{{{({x - {l_P} \times {n_x}} )}^2} + {{({y - {l_P} \times {n_y}} )}^2}}}{{\omega_0^2}}} \right)}}\\ {n_x} \in [{0,\ldots ,{N_x}} ]\quad {n_y} \in [{0,\ldots ,{N_y}} ]\end{array}.$$
By summation of all successive pulse fluences the accumulated fluence profile $\mathrm{\Gamma }({x,y} )$ can be calculated as:
$$\Gamma ({x,y} )= \sum\limits_{{n_x}} {\sum\limits_{{n_y}} {\phi ({x,y,{n_x},{n_y}} )} }.$$
A line profile is obtained when one coordinate is not varied. The highest value in the set of data is considered as accumulated peak fluence ${\mathrm{\Gamma }_{\textrm{AP}}}$ while the lowest value is called the accumulated valley fluence ${\mathrm{\Gamma }_{\textrm{AV}}}$. We believe that ${\mathrm{\Gamma }_{\textrm{AV}}}$ is of particular interest for the optimization of ablation processes. In addition, the homogeneity of the distributed laser energy must be considered. As shown in Fig. 3 the accumulated fluence profile of a laser line varies along its $\textrm{X}$-trajectory.

 figure: Fig. 3.

Fig. 3. Accumulated fluence profile (red) obtained from a line of laser spot intensities.

Download Full Size | PDF

To avoid substrate damage, local high fluence maxima should be avoided, in particular if they exceed the substrate threshold fluence. To account for the fluence variance the overshoot variable ${\theta _P}$ is established, which relates ${\mathrm{\Gamma }_{\textrm{AV}}}$ to ${\mathrm{\Gamma }_{\textrm{AV}}}$:

$${\theta _P} = \frac{{{\mathrm{\Gamma }_{\textrm{AP}}} - {\mathrm{\Gamma }_{\textrm{AV}}}}}{{{\mathrm{\Gamma }_{\textrm{AV}}}}}\;[\%].$$
Ahmmed et al. considered a fluence profile constant or flat if the overshoot value is smaller than 1%, which can easily be achieved with large pulse overlaps [14]. This is illustrated in Fig. 4, which shows different line accumulated fluence profiles with increasing $\varphi $. Figure 4(a) represents the case $\varphi \ge 1$ where there is no overlap at all (spot size is specified by ${\phi _0}$/e2 border). In this case, ${\mathrm{\Gamma }_{\textrm{AP}}}$ equals ${\phi _0}$ but ${\theta _P}$ is large. Overlap can be increased to a level at which ${\theta _P}$ decreases but is still rather large while $\mathrm{\Gamma }(x )$ practically not yet exceeds ${\phi _0}\; $(see Fig. 4(b)). With even further decrease of $\varphi $ (as in Fig. 4(c)) a point can be reached where ${\mathrm{\Gamma }_{\textrm{AP}}}\; $exceeds ${\phi _0}$ by a non-negligible amount. This point can be marked by the condition ${\mathrm{\Gamma }_{\textrm{AV}}}$ = ${\phi _0}$. Beyond that, the accumulated fluence level can raise far beyond ${\phi _0}$ (in Fig. 4(d)) and eventually shows a very homogeneous fluence distribution, characterized by a very low ${\theta _P}$.

 figure: Fig. 4.

Fig. 4. Accumulated 1D fluence profiles (red) for line ablation with: a) No pulse overlap and b) Pulse overlap such that peak accumulated fluence is practically still equal to maximum spot fluence c) Pulse overlap such that valley accumulated fluence equals maximum spot fluence d) Large pulse overlap with accumulated fluence profile can be considered “flat” and overshoot indicator is small e.g. below a certain threshold.

Download Full Size | PDF

The fluence distribution in the case of 2D areal ablation patterns is illustrated in Fig. 5. Note that values for the accumulated fluence profile differ between line (1D)- and areal (2D)-ablation for the same value of $\varphi $. Assuming a standard areal ablation pattern in which laser lines (written in $\textrm{X}$-direction) are stacked above each other (in $\textrm{Y}$-direction), the accumulated fluence profile sums up by the contribution of successive pulses along the laser line as well as by pulses from vertically neighboring lines. Thus, $\mathrm{\Gamma }$ will always be higher for 2D accumulated profiles than for 1D accumulated line profiles assuming the same value of $\varphi $. A homogeneously spread accumulated fluence profile (small overshoot variable ${\theta _P}$) is required to determine well-defined process parameters allowing a controlled areal ablation.

 figure: Fig. 5.

Fig. 5. 2D areal accumulated fluence distributions (${\phi _0}$ = 1 J/cm2, ${\omega _0}$ = 16 µm) with: a) Large gaps between individual spots such that fluence level locally drops to almost zero b) Overlapping pulses resulting in an accumulation of pulse fluences c) Pulse spacing similar to the pulse radius resulting in a practically homogeneous (flat) accumulated fluence profile far above ${\phi _0}$.

Download Full Size | PDF

3.2 Ablation pattern

Process parameters are considered well suited, if they allow complete homogenous ablation of a metallic thin film. At the same time, the underlying polyimide substrate must stay transparent indicating that surface damage is avoided. This is possible only if the accumulated fluence is homogeneously spread. To obtain a texture free homogeneous areal ablation, an isotropic and practically homogeneous 2D fluence distribution is required. The isotropic fluence distribution requirement excludes some of the ablation patterns that are in practical use. Here we use the #-pattern as shown in Fig. 2: Gaussian fluence profile Figure 2 where horizontal lines are vertically stacked at a pitch of ${l_{P,y}} = {l_P,x}$.

As seen in Fig. 6(a), the maximal pulse distance ${l_{P,\textrm{max}}}$ to completely cover an area with the "#-pattern” is given by ${\approx} 0,71\cdot D$. With increasing pulse overlap ($\varphi \to 0$) more area on the substrate is irradiated multiple times (e.g. Figure 6(b)). This makes the comparison of single pulse thresholds with ablation thresholds for line or areal ablation a challenging task as the incubation effect must be considered (see example of areal fluence distribution in Fig. 7).

 figure: Fig. 6.

Fig. 6. Illustration of #-pattern: a) Explaining the maximal pulse distance ${l_{P,\textrm{max}}}$ to prevent uncovered areas of the specimen b) Providing an example pattern where the borders of overlapping pulse areas and spot centers are marked. The white area is receiving fluence contributions from four single pulses.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. 2D fluence distribution obtained with a “#-pattern” at ${l_P} = 20\; \textrm{}\mathrm{\mu}\textrm{m} = 0,63\; \textrm{}\cdot D\; \textrm{}$. As this is smaller than ${l_{P,\textrm{max}}}$, all areas between pulses are irradiated multiple times.

Download Full Size | PDF

3.3 Further relevant ablation parameters

Liu et al. stated that any laser ablation mechanism starts with the absorption of laser energy in the target [15]. Because metals show high absorptivity, light is fully absorbed well within the first micrometer below the surface [7]. Typical coating thicknesses are between 50 and a couple of hundred nanometers. For a pure nickel coating, penetration depth was found to be approximately 50 nm using Lambert-Beer’s law [7,16]. The coating on our samples was fabricated with a sputter target made with 80 wt% nickel and 20 wt% chromium, such that the penetration depth of nickel gives a good approximation for the NiCr alloy.

Pulse radius as well as threshold fluence are important factors for laser induced ablation processes. Liu et al. have found a linear relation between the square of the optically measurable diameter of the laser modified zone D and the fluence plotted on a logarithmic scale [17]:

$${D^2} = 2\omega _0^2 \times \ln \left( {\frac{{{\phi_0}}}{{{\phi_{\textrm{th}}}}}} \right).$$
This relation is commonly used to calculate (by squared error regression with $y = a\ln (x )+ b$ [7]) effective pulse radius and threshold fluence for bulk material ablation or laser-induced surface modifications but has some limitations when it comes to thin film ablation.

Jee et al. have shown that material specific threshold fluence diminishes with increasing number of incident pulses [18]. This holds true for coating and substrate, which is why large pulse overlap ablations often result in substrate damage. The law of incubation can be stated as a formula, using the so-called incubation factor $\xi$:

$${\phi _{\textrm{th}}}({{N_P}} )= {\phi _{\textrm{th}}}(1 )\times N_P^{\xi - 1}.$$
A value of $\xi = 1$ means that no incubation influence can be observed. Typical values of $\xi $ for metals are between 0,8 and 0,96 [7,1].

4. Spot ablation

To describe the selective ablation criteria in case of thin film ablation the minimum level of the accumulated fluence profile ${\mathrm{\Gamma }_{\textrm{AV}}}$ will be related to the previously determined accumulated single spot threshold fluence. As the threshold fluence is dependent on the number of incident pulses, the so called incubation factor must also be determined. Thus, the material specific threshold fluence ${\phi _{\textrm{th}}}$ and the incubation factor $\xi $ of the coating were identified at first. This was accomplished with the aid of multiple D2-experiments where the coating was irradiated with different numbers of single laser pulses ${N_P}$ at increasing energy. In order to gain the material specific incubation factor, the threshold values for different numbers of pulses was plotted such that a line fit could be conducted. For these experiments the laser energy was configured, such that a single shot ${N_P} = 1$ can fully trepan the full coating thickness at the highest fluence level used. We used confocal laser scanning microscopy to measure diameter and ablation depth at the irradiated spot. The energy per pulse was varied in the range of $0,29\; \textrm{}\mathrm{\mu}\textrm{J}$ to $4,83\; \textrm{}\mathrm{\mu}\textrm{J}$, giving fluence values between $70\; \textrm{mJ}/\textrm{cm}^2$ and $1,11\; \textrm{J}/\textrm{cm}^2$.

4.1 Ablation zones

An ablation spot on a thin film coated polymeric substrate is not clearly confined and multiple characteristic zones can be observed. Figure 8 shows three distinguishable ablation zones. Outer areas with low fluences only show surface modifications and cracks in the coating. We can define the size of this influence zone by the expansion of the cracks. Higher fluences yield a ring-shaped area of selective ablation in which the substrate is unveiled but not damaged. This will be called selective ablation zone. For the highest fluence around the spot center a destruction zone evolves, which is recognized optically by a milky and/or rough appearance of the polymeric substrate. Fluence thresholds for each of the three zones can be determined.

 figure: Fig. 8.

Fig. 8. LSM micrograph (x150) showing three distinguishable ablation zones for a NiCr-coated sample (100 nm) with two different imaging technics: a) 3D CLSM image b) Differential interference contrast-image (DIC)

Download Full Size | PDF

4.2 Single pulse experiment

The development of the ablation zones at increasing laser fluence is illustrated in the subsequent micrographs of Fig. 9. In Fig. 9(a) the fluence is low so that only some surface modification of the coating was observed. With increasing fluence a clear influence zone with cracks started to appear (Fig. 9(b)). Close to the selective ablation threshold fluence the coating lifted up and formed a bladder (Fig. 9(c)). Also, the coating was sometimes completely removed in this area due to the strong particle suction system inside the laser machine. Further increase in fluence allowed to generate a clear selective ablation zone as shown in Fig. 9(d). The step in height between the coating (green area) and the unveiled substrate (light blue) was in the range of 100 nm which is the thin film thickness.

 figure: Fig. 9.

Fig. 9. Results of single pulse experiments with increasing fluence (left: 3D CLSM image / right: direct optical view) a) First material reaction b) Cracks form inside the influence zone c) The coating lifts up d) Selective ablation showing no substrate damage

Download Full Size | PDF

4.3 Multiple pulse experiments

Incubation experiments with multiple pulse repetitions (at $\varphi = 0$) were carried out with gradually increasing pulse energy where the increasing size of the influence, the ablation and the destruction zone was observed. For small fluences a cracked influence zone with a small selective ablation area in the spot center became visible on the surface (Fig. 10(a)). Note that laser fluence is equal to the one used for the single pulse experiment in Fig. 9(a). With increased fluence the diameter of the influence zone and selective thin film ablation zone in the spot center became larger (Fig. 10(b)). Solely some dust particles sometimes remained on the substrate. These could be cleaned using an ultrasonic bath. With even higher fluence, substrate damage started to appear (Fig. 10(c)). In Fig. 10(d) the destruction zone had expanded drastically such that the zone of selective areal ablation is narrowed.

 figure: Fig. 10.

Fig. 10. Results of multiple (${N_P} = 6$) pulse experiments with increasing fluence (left: 3D CLSM image / right: direct optical view) a) Large influence zone with cracks b) Selective ablation zone, dust particles remain but no substrate damage can be detected c) Substrate damage starts to appear d) Selective ablation zone ring narrows

Download Full Size | PDF

The measured diameters of different experiments were used to draw the D2-plot in Fig. 11. The plot shows only the data of the selective ablation zone. Selective threshold fluence for ${N_P} = 80$ is found to be ${\phi _{\textrm{th},80}} = 0,07\; \textrm{J}/\textrm{cm}^2$.

 figure: Fig. 11.

Fig. 11. D2-plot for a thin film selective ablation of 100 nm NiCr on polyimide. The diagram shows lines according to Liu's formula. They were obtained for each pulse repetition series by linear regression. Differences in the slopes result from measurement uncertainties and small instabilities of the ablation setup as the spot radius ${\omega _0}$ stayed constant during all experiments.

Download Full Size | PDF

The accumulated fluence ${\mathrm{\Gamma }_{\textrm{AP}}}$ that achieves selective ablation was calculated for stationary spot experiments by multiplying the selective threshold fluence with the number of irradiated pulses:

$${\Gamma _{\textrm{AP}}} = {\phi _{\textrm{th}}}({{N_P}} )\times {N_P}.$$
If the number of irradiated pulses and the chosen fluence were well picked, full selective ablation with no substrate damage or dust generation could be achieved (see Fig. 12).

 figure: Fig. 12.

Fig. 12. 3D CLSM image proving that full selective ablation is achievable

Download Full Size | PDF

In section 3 we already calculated the gaussian spot size to be ${\emptyset _{0,\textrm{calc}}} = 40,1\; \textrm{}\mathrm{\mu}\textrm{m}$. In comparison, the evaluation of the D2-experiment in the selective ablation regime lead to a spot radius of only ${\omega _{0,\textrm{D}^2}} \approx 16\; \textrm{}\mathrm{\mu}\textrm{m}\; \Rightarrow {\emptyset _{0,\textrm{D}^2}} = 32\; \textrm{}\mathrm{\mu}\textrm{m}$. In Fig. 13, it can be seen that the D2-evaluation tends to yield higher values if the number of incident pulses increases. If plotted data is averaged for the three different zones, the radius of the influence zone ${\omega _{0,\textrm{D}^2}} \approx 20\; \textrm{}\mathrm{\mu}\textrm{m} \Rightarrow {\emptyset _{0,\textrm{D}^2}} = 40\; \textrm{}\mathrm{\mu}\textrm{m}$ correlates best with the calculated size.

 figure: Fig. 13.

Fig. 13. Plot of different zone radii evaluated from D2-experiments

Download Full Size | PDF

This means that influence zone data is most valuable to define the spot size, whereas selective ablation zone data must be used to define the threshold fluence.

4.4 Incubation effect

From our D2-experiments of NiCr we found the incubation factor for selective ablation to be ${\xi _{\textrm{NiCr},\; \textrm{sel}.\; \; \textrm{abl}}} = 0,53$. This is smaller than typical values for bulk material ablation of metals (see section 3.3). Therefore, we conclude that selective ablation of thin films is not only dependent on material evaporation but also majorly influenced by stress assisted interface destruction between coating and substrate. Figure 14 summarizes the threshold fluences for the three ablation zones defined in section 4.1 at increasing pulse numbers. It can be seen, that the threshold fluence of ablation and destruction zone are identical for number of pulses ≥ 40. Even though selective ablation can also be achieved in a narrow window at higher pulse numbers, we interpret this as a suitable upper boundary condition.

 figure: Fig. 14.

Fig. 14. Plot of threshold fluences for different ablation zones with increasing number of pulses at 600 kHz. Incubation factor $\xi $ can be derived from these curves.

Download Full Size | PDF

Thus, laser frequency, mark speed and pulse divider must be chosen in such a way, that the number of pulses hitting the very same position on the specimen stays below the upper boundary condition (${N_P} \le 40$ in this case). Otherwise, substrate damage is likely to occur.

4.5 Influence of laser frequency

The laser pulse frequency is not considered to influence the ablation threshold. Our experiments however revealed, that ablation behavior changes for higher numbers of pulses depending on laser frequency. This demands for an adjustment of the current well known methodology for ablation threshold determination developed by Liu et al. [2] which is not considering frequency influence. In Fig. 15, it can be seen, that the experimental selective threshold fluence for smaller number of pulses (1-10) was indeed frequency independent. For higher pulse numbers (≥20) however, threshold values obtained using equation [Eq. (8)] for 16 Hz and 600 kHz differed.

 figure: Fig. 15.

Fig. 15. Direct comparison of D2-experimental data confirms that ablation thresholds are influenced by laser frequency only at higher number of pulses ≥ 20.

Download Full Size | PDF

An incubation factor ${\xi _{\textrm{NiCr},\; \textrm{sel}.\; \textrm{abl},16\textrm{Hz}}} = 0,67$ for the selective ablation zone was obtained at 16 Hz pulse frequency. Figure 16 clarifies this phenomenon. As the threshold calculation is based on the measurable ablated zone diameter, the diameter values gained from the two upper rows of the figure were very similar. In contrast, in case of higher pulse numbers (the two lower rows) the affected zone is always larger for the higher pulse frequency. Also, the selective ablation and destruction zone have joined and cannot be perceived separately. Thus, the purely selective ablation is absent at higher fluences. This means that the transition from purely selective ablation to substrate damage becomes more sudden at higher pulse numbers.

 figure: Fig. 16.

Fig. 16. Comparison of multiple pulse experiments at increasing single pulse fluence. In case of ${N_P} = 8$ (two upper rows) the ablated diameter is similar and independent of laser frequency. At the bottom rows ${N_P} = 20$ the selective zone ablated at 600 kHz is always larger than at 16 Hz. Black areas indicate severe substrate damage where the crater is deeper than 200 nm.

Download Full Size | PDF

In addition, optically verifiable ablation happened at a fluence below the selective ablation threshold when pulse frequency is high (second row), even though the calculated threshold values matched for both frequencies. We assume this to be the results of different delamination effects. In case of low frequency ablation, the coating was “slowly” pulverized leading to an increased amount of residual particles and a poorly definable ablation zone border. For high frequency ablation though, the coating experienced a sudden flat lift and was removed at once. This is why the coating tends to bulge upwards at the boundary of the ablated zone at 600 kHz.

5. Selective ablation in scan patterns

From theory, there are two ways to achieve homogeneous line and areal ablation. On the one hand, the substrate can be irradiated with low fluence and a high pulse overlap resulting in a smooth accumulated fluence profile. In this case, the accumulated fluence level is much higher than the maximum spot fluence value ${\phi _0}$. On the other hand, lower pulse overlap values allow for higher single pulse fluence levels (which showed to be beneficial for coating removal) while the overall fluence variance (${\theta _P} < \; \textrm{e}.\textrm{g}.\; \textrm{}1\; \textrm{\%}$) is still acceptably low.

Suitable material ablation process conditions are typically found for a fixed scan strategy. The ablation threshold is then experimentally obtained using a parameter grid. In the following, we show that selective ablation for line- and/or areal ablation can be achieved based on the consideration of stationary spot threshold values only. This will be shown for scan patterns with practically isotropic and homogeneous fluence distribution. In the following, we will derive the according rules and show how they apply for the exemplary case of thin film NiCr ablation. We do not discuss multiple repetitions of the same ablation pattern as they do not improve the surface quality obtained after ablation. They are only necessary in cases of thicker layers.

5.1 Line ablation

Single laser lines structured into a metallic thin film coating can be used to electrically insulate certain regions of a coating. This can yield sensor structures or help to separate conductive areas. Because very little material must be removed, line ablation, which we call contour cut [19] is a time and cost effective technique. For line ablation laser pulses are placed along a line – each pulse overlapping its predecessor as defined by the parameter $\varphi $.

In line ablation, the average number of pulses ${\tilde{N}_P}\; $that hit a single position on the specimen can be defined as the product of pulse frequency f and the spot diameter ${\emptyset _0}$ divided by the mark speed v and pulse divider value $PD$. The pulse divider periodically guides a certain amount of pulses into a beam dump. The pulse frequency which is effective on the substrate can thereby be changed without changing the pulse generation frequency (600 kHz for our machine) and the single pulse fluence. $PD$=1 means that all pulses pass through the device such that they leave the source at their generation frequency:

$${\tilde{N}_P} = \left\lceil {\frac{{f \times {\emptyset_0}}}{{v \times PD}}} \right\rceil.$$
Figure 17(a) shows a micrograph with transmitted light illumination of a large line pattern test grid. Single pulse fluence was increased from left to right while mark speed was increased from the bottom to the top of the grid. Each “box” (the boxframe is solely a position marker) contains two scan lines having the same pulse spacing. The upper line was written at a pulse frequency of 600 kHz. The lower line consists of pulses placed on the specimen at a lower frequency of only 16 Hz. The pulse fluence and average number of pulses ${\tilde{N}_P}$ is equal for both lines. In Fig. 17(b) the experimental results of the optically measured cut width are compared with the calculable theoretical values.

 figure: Fig. 17.

Fig. 17. a) Line ablation grid. Inside each box frame, there are two ablated lines: Upper line is marked at 600 kHz, lower line is marked by placing individual pulses in a row yielding a slow frequency of 16 Hz (pulse fluence ${\phi _0}$ and the number of pulses on a single spot ${\tilde{N}_P}$ are the same for both lines). White numbers indicate cut width of lines [µm] ablated at 600 kHz. Side pictures: i) A fluence of 0,22 J/cm2 and a mark speed of 2500 mm/s allowed for a selective line ablation using the standard 600 kHz provided by the laser source. In contrast, the line ablated at 16 Hz is intermittent ii) The upper line shows severe substrate damage while the lower line is still in the selective region. Edgeline quality of lower line is low iii) At these parameters the upper line is selectively ablated whereas no ablation is observed in the case of slow pulse frequency. b) Comparison of cut width D calculated for mark speed v at increasing fluence ${\phi _0}$ using equations [Eq. (8),(9)] and [Eq. (11)] with cut width measured optically (see Visualization 1 for high-resolution image of the entire figure).

Download Full Size | PDF

Three different areas are separated (indicated by no, yellow and red shading). The area on the left (no color) corresponds to conditions not leading to any ablation. The middle region (yellow) is the selective ablation area. White numbers inside the fields represent the measured cut-width of the upper line in µm. The red zone on the right marks those fields, where substrate damage is detectable by optical means. Only the lines written at 600 kHz determine to which area (no color, yellow or red) the corresponding parameter set of the box is assigned to. From the grid the single pulse fluence values ${\phi _0}$ at which selective ablation starts to occur with respect to the mark speed (e.g. ${\phi _0} = 0,16\; \textrm{J}/\textrm{c}{\textrm{m}^2}$ and mark speed of $1500\; \textrm{mm}/\textrm{s}$) can be read off. The border between either no or selective ablation is given by the following criterion (red line in the plot).

Selective ablation criterion for line-ablation:

$$\boxed{{{\textrm{1}^{\textrm{st}}}\textrm{criterion:}\quad {\phi _0} \ge {\phi _{\textrm{th}}}({{{\tilde{N}}_P}} )}}.$$
The single pulse ablation thresholds (see Fig. 15) show up to ${N_P} = 10$ overlapping pulses no dependency on pulse frequency. We saw a pulse frequency dependence in our line patterns (the two lines in most box frames appear differently) which cannot be explained without assuming a difference in delamination effects, which likely have to do with progressing interfacial damage. Fast successive pulses beyond the threshold for selective ablation cause a coating delamination and result in flaking. However, this effect did not occur for low pulse frequencies. This is important to note as this interferes with any attempt to forecast ablation behavior at very low frequencies based on the first criterion. In our data we noticed this frequency dependency only, if pulse overlap was quiet large (as for the grid in Fig. 17). For reduced pulse overlap ablation results for high and low pulse frequencies converge.

The graph in Fig. 18 contains the threshold fluence values for selective ablation and substrate destruction determined in the previous stationary spot D2-experiments in dependence of the number of incident pulses. They define the window for selective ablation with ${\tilde{N}_P}$ as the decisive parameter in line pattern ablation.

 figure: Fig. 18.

Fig. 18. Logarithmic plot of fluence threshold values as obtained by stationary spot experiments (multiple pulses) which define the selective ablation window (yellow) for selective line ablation.

Download Full Size | PDF

The blue diamonds mark the fluence level ${\phi _0}$ in the plot where visible ablation started to occur in dependence of averaged number of overlapping pulses ${\tilde{N}_P}$. All these data points lie inside the proposed selective ablation zone. The knowledge of the pulse fluence ${\phi _0}$, the spot size and the threshold fluence corrected by the incubation effect allow the calculation of the resulting cut width D (see Fig. 1). The plot in Fig. 18 proves that experimental data and theory align.

5.2 Areal ablation using #-pattern

In order to predict areal ablation results the incubation effect plays an important role because the 2D-irradiation pattern causes exponentially more pulse overlaps than for line ablation. Depending on the ablation pattern (here we focus on the isotropic #-pattern as described in section 3.2), the number of laser pulse hits on a single position must be determined. This value can be calculated knowing the pulse spacings used in the experiment. In order to do so, for each laser pulse of the ablation pattern, a cylinder on the sample surface with a height of one is modelled. By adding up the heights of all distributed cylinders on the surface, the number of laser impacts at any position $Z({x,y} )$ can be determined as:

$$\begin{array}{l} Z({x,y} )= \sum\limits_{{n_x}} {\sum\limits_{{n_y}} {\Delta Z({x,y,{n_x},{n_y}} )} } \\ {n_x} \in [{0,\ldots ,{N_x}} ]\quad {n_y} \in [{0,\ldots ,{N_y}} ]\end{array}$$
$$\textrm{with}\quad \Delta Z = \left\{ {\begin{array}{{c}} {0\quad \textrm{if}\,\sqrt {{{({x - {l_P} \times {n_x}} )}^2} + {{({y - {l_P} \times {n_y}} )}^2}} \ge {\omega_0}}\\ {1\quad \textrm{if}\,\sqrt {{{({x - {l_P} \times {n_x}} )}^2} + {{({y - {l_P} \times {n_y}} )}^2}} < {\omega_0}} \end{array}}. \right.$$
A Z distribution calculation is exemplified in Fig. 19. Meanvalue of incident pulses is two for the shown parameter set.

 figure: Fig. 19.

Fig. 19. Map and histogramm of the specimen. In this case, a #-pattern with ${\omega _0}$=16 µm and a pulse spacing of ${l_P}$=20 µm (results from a mark speed of 3000 mm/s, a laser frequency of 600 kHz and a pulse divider of 4) was evaluated. Those regions that were hit by most pulses are marked yellow, whereas those areas struck only once remain blue. The histogram on the right shows how large the surface fraction is (in arbitrary units) that is hit by a certain number of laser pulses

Download Full Size | PDF

From a histogram as shown in Fig. 19 the average number of pulse hits ${\tilde{N}_P}$ on the surface can be calculated (weighing them with their corresponding area). Depending on the pulse overlap, the mean number of pulses per spot always differs assuming a constant spot diameter.

We created an areal ablation parameter grid to test ablation thresholds determined by stationary spot experiments with practical data gained from areal ablation. Figure 20 shows this parameter grid using the #-pattern each field having a size of 1000 × 1000 µm. Fluence and mark speed were adjusted in the same manner as for line ablation (look back to Fig. 17). It can be seen, that areal ablation with a single pulse fluence of ${\phi _0} = 0,11\; \textrm{J}/\textrm{cm}^2$ caused interface destruction between coating and underlying substrate (field 3/3). This ablation mechanism was denoted “flaking”. The polymer substrate appeared very clear and bright after flaking, which indicates a damage free surface. Even though selective areal ablation in this regime seemed a promising method for areal ablation, it has to be taken into account that for geometric shapes more complicated than squares we observed flakes to adhere much stronger. In our experiments, flakes generated by areal ablation contours with straight edges broke away easily whereas for any pointy features (in larger contours), such as triangles, tips or thin lines, residual flakes usually remained on the surface. The narrower the gap between the untouched coating sidewalls, the more the ablated flakes in the middle adhered to the surface.

 figure: Fig. 20.

Fig. 20. Areal ablation grid with pulse divider 1 at 600 kHz pulse frequency (see Visualization 2 for high-resolution image). Fluence is increased left to right. Values on the left indicate the mark speed dependent mean number of pulses ${\tilde{N}_P}$ for each row. Blue area marks all those fields where the coating has lifted in flakes from the substrate. All fields inside the selective ablation zone show some unwanted greyish substrate color. This is because of residues and dust. Nevertheless, they are electrically isolating. The red area defines those fields where severe substrate damage has occurred. The five micrographs on the bottom enlarge some of the fields inside the grid.

Download Full Size | PDF

Right of the flaking region the ablation zone appears. Laser parameters for field 9/5 have yielded a complete coating removal, even though unveiled substrate is not as bright as for field 3/3 and 4/3 (flaking ablation). Laser parameters of field 11/5 caused a distinct substrate damage noticeable by the color change. Generally speaking, substrate damage came along with surface roughness scattering the through light. Thus, the darker the substrate appeared the more severe was the damage. Ablation by flaking could provide ideally preserved substrate surfaces but appeared to be inappropriate for many applications as flake residues could cause electrical short circuits. In the ablation area to the right of the flaking regime, however, the substrate appeared comparatively dark. The explanation can be found in Fig. 14, where differences in the threshold fluences for selective ablation and damage almost vanish when the mean number of incident pulses is equal or greater than 40.

A simple way to reduce the number of incident pulses per area is the use of a pulse divider. With this, the pulse overlap can be decreased without affecting the fluence. In addition, laser position accuracy is not influenced, as this would be the case using very fast mark speeds alternatively (also technical limitations for fastest galvo mirror speed must be considered). Figure 21 shows the corresponding parameter grid with a pulse divider of four. As the level of the accumulated fluence profile was much lower as compared to the case of pulse divider one (in Fig. 20) the single spot fluence could be much larger without causing substrate defects (therefor the grid was extended to the right). In this case five regions could be identified. Right of the flaking zone appeared a dotting zone (green area). “Dotting” means, that selective ablation was achieved only in the center of the pulses, but pulse spacing was such, that unveiled circles did not intersect and coating residues remained between them. There was a rather smooth transition from the dotting zone to the selective ablation zone. Laser parameters in the selective ablation area allowed flake free ablation of coatings without damaging the substrate. If one compares the substrate coloration in the flaking region (4/1) with the ablation results in the selective region (e.g. 11/2, 4/3 and 7/4), there is hardly any difference. No selective ablation results can be detected in the bottom line of the grid because the mean number of pulses (in this case ${N_P} = 73$) was too large, such that selective ablation threshold of the coating and destruction threshold of the substrate overlap in the incubation plot (see Fig. 14). As a result, flaking and destruction region are direct neighbors in the bottom row.

 figure: Fig. 21.

Fig. 21. Areal ablation grid with pulse divider 4 at 600 kHz pulse frequency (see Visualization 3 for high-resolution image). Fluence is increased left to right. Mean number of pulses ${\tilde{N}_P}$, pulse spacing, pulse overlap and mark speed are indicated on the left. Five ablation zones can be differentiated: no ablation, flaking, dotting, selective ablation and destruction regime. In the dotting regime the coating shows a uniform pattern of holes but retains areal conductivity. Two criterion boarder lines are marked inside the grid. The ten micrographs at the bottom enlarge fields inside the grid which are of particular interest.

Download Full Size | PDF

The first criterion, validated for the line experiments, directly compares the fluence with the pulse number dependent selective ablation threshold obtained in stationary beam experiments. The boundary given by this criterion is indicated in the areal parameter grid (Fig. 21) by a pink line. However, it may be seen that this criterion does not allow an exact delimitation of the selective ablation range. One reason amongst others for this deviation is a simplification in the calculation of ${\tilde{N}_P}$. The Gaussian profile of the laser was disregarded so far. Therefore, the accumulated valley fluence ${\mathrm{\Gamma }_{\textrm{AV}}}$ as defined in section 3.1 is considered to account for the Gaussian beam profile and will be used for a more precise criterion in the following.

Figure 22 shows a plot where the product of ablation threshold and the number of pulses is related to the accumulated valley fluence. The ablation threshold ${\phi _{\textrm{th}}}({{{\tilde{N}}_P}} )$ defines the minimal fluence which is needed to achieve selective ablation with a number of pulses ${\tilde{N}_P}$. Thus the product of the experimentally obtained selective ablation threshold and the mean number of pulses is equal to the minimal accumulated fluence level needed to achieve this ablation. This value can therefore be compared to the calculated accumulated valley fluence ${\mathrm{\Gamma }_{\textrm{AV}}}$. This calculation yields the second criterion in a rearranged form:

$$\boxed{{{\textrm{2}^{\textrm{nd}}}\textrm{criterion:}\quad \frac{{{\Gamma _{\textrm{AV}}}}}{{{{\tilde{N}}_P}}} \ge {\phi _{\textrm{th}}}({{{\tilde{N}}_P}} )}}.$$
Only if the accumulated valley fluence devided by the average pulse number is larger than the threshold, defect-free areal ablation will occur.

 figure: Fig. 22.

Fig. 22. Accumulated fluence plot. The accumulated valley fluence is depicted for different mark speeds in relation to spot fluence (blue dotted lines). Yellow shaded area marks the selective ablation zone. If accumulated valley fluence is larger than the product of mark speed dependent pulse number and the corresponding ablation threshold then selective ablation is achieved. (Refer to Fig. 19 for correlation between pulse number and mark speed e.g. ${N_P} = 18$ for 1000 mm/s. The corresponding threshold fluence can be drawn from the plot in Fig. 18.)

Download Full Size | PDF

The second criterion is added to the areal ablation grid indicated by a red line. As expected, it precisely meanders along the transition region between dotting and selective ablation. Thus, it can be used as a boundary condition to find suitable selective ablation parameters without the need for a large parameter grid. Furthermore, a set of parameters close to this calculated boundary can be considered to keep a safe distance from substrate damage. This provides a numerical value (distance to the boundary condition), which can be used to evaluate the ablation results in parameter optimizations when optically hardly any differences are perceptible.

In conclusion, we have experimental evidence that proves functionality of two simple criteria to predict ablation behavior in selective metallic film removal from polymeric substrates. The well-known first criterion, which compares the laser fluence with the pulse number dependent ablation threshold is valuable to forecast the ablation results of single spot and line ablation. In case of areal ablation this criterion is not piercing enough though. To predict at which laser fluence level areal selective ablation can be achieved we formulated a second, sharper criterion.

The described method for areal selective ablation promotes a combination of high-energy pulses and rather large pulse distancing to reduce the risk of substrate damage. However, a small pulse overlap can lead to noticeably rippled edges (as can also be seen in Fig. 23). Therefore, keeping the overall accumulated fluence low with this method is an appropriate way to achieve selective areal ablation on susceptible substrates but a subsequent smoothing of the contour with low fluence but high pulse overlap might be necessary.

 figure: Fig. 23.

Fig. 23. Micrograph of a pattern generated by successful areal ablation. Left side shows the sample directly after laser irradiation with parameters chosen according to the second criterion. Right hand side picture was taken after additional ultrasonic cleaning (to remove coating dust) at through light illumination only (thus areas with remaining opaque coating appear black). Any substrate damage would appear dark on the light yellowish polyimide substrate, but such destruction has clearly not occurred here. Because of the relatively small pulse overlap, minor positioning errors of the laser spot left some tiny coating flakes behind. Edgeline quality can be enhanced if areal ablation is followed by a contour line ablation with a larger pulse overlap.

Download Full Size | PDF

6. Conclusion

A method has been established to provide laser parameters for full selective ablation of a metallic thin film (NiCr) on a susceptible polymeric substrate (polyimide) in an efficient and reproducible manner. Thus, it can improve process optimization in research as well as for industrial applications whenever laser induced ablation of thin films is used for system fabrication. In particular, pulse overlap and pulse fluence have to be carefully controlled. The methodology is based on experimentally derived ablation thresholds and a calculated accumulated fluence profile to minimize overall energy input into the substrate. This may greatly reduce the risk of substrate damage. In practice, it may be sufficient to calculate ${\tilde{N}_P}$ for simple line ablations and work with the incubation dependent ablation threshold. In case of areal ablation, the second criteria must be considered. This allows to predict suitable parameters for an areal selective ablation (for methodology overview see Table 1).

Tables Icon

Table 1. Methodolody of selective metallic thin film ablation

Funding

Deutsche Forschungsgemeinschaft (BR 2178/36-1, DI 1934/8-1).

Disclosures

The authors declare no conflicts of interest.

References

1. J. F. Düsing, T. Temme, and R. Kling, “Micro patterning of metal thin films on polyimide foils using high-repetition picosecond laser,” in Laser Assisted Net Shape Engineering 5, Proceedings of the 5th LANE 2007, M. Geiger, A. Otto, and M. Schmidt, eds. (Meisenbach-Verlag, 2007), pp. 1157–1166.

2. S. ONLINE, Hamburg, and Germany, “Galaxy Fold 5G: Samsung kündigt erneut Verkaufsstart für sein Falthandy an - SPIEGEL ONLINE - Netzwelt,” https://www.spiegel.de/netzwelt/gadgets/galaxy-fold-5g-samsung-kuendigt-verkaufsstart-seines-falt-handys-an-a-1285328.html.

3. J. van den Brand, J. de Baets, T. van Mol, and A. Dietzel, “Systems-in-foil – Devices, fabrication processes and reliability issues,” Microelectron. Reliab. 48(8-9), 1123–1128 (2008). [CrossRef]  

4. R. Mandamparambil, H. Fledderus, G. van Steenberge, and A. Dietzel, “Patterning of Flexible Organic Light Emitting Diode (FOLED) stack using an ultrafast laser,” Opt. Express 18(8), 7575–7583 (2010). [CrossRef]  

5. M. Kasischke, E. Subaşı, C. Bock, D.-V. Pham, E. L. Gurevich, U. Kunze, and A. Ostendorf, “Femtosecond laser patterning of graphene electrodes for thin-film transistors,” Appl. Surf. Sci. 478, 299–303 (2019). [CrossRef]  

6. N. Farid, H. Chan, D. Milne, A. Brunton, and G. M. O’Connor, “Stress assisted selective ablation of ITO thin film by picosecond laser,” Appl. Surf. Sci. 427, 499–504 (2018). [CrossRef]  

7. O. Suttmann, M. Gosselin, U. Klug, and R. Kling, “Picosecond laser patterning of NiCr thin film strain gages,” A. Heisterkamp, J. Neev, S. Nolte, and R. P. Trebino, eds. (SPIE, 2010), p. 758914.

8. N. Bärsch, K. Körber, A. Ostendorf, and K. H. Tönshoff, “Ablation and cutting of planar silicon devices using femtosecond laser pulses,” Appl. Phys. A 77(2), 237–242 (2003). [CrossRef]  

9. J. Mur, J. Petelin, J. Schille, U. Loeschner, and R. Petkovšek, “Ultra-fast laser-based surface engineering of conductive thin films,” Appl. Surf. Sci. 509, 144911 (2020). [CrossRef]  

10. D. Vollberg, A.-C. Probst, M. Langosch, A. Landes, D. Göttel, M. Cerino, A. Lellig, O. Freitag-Weber, and G. Schultes, “Hochempfindliche Folien-Dehnungsmessstreifen auf dem Weg zur technologischen Reife,” tm - Technisches Messen 82 (2015).

11. M. Schwerter, L. Hecht, E. V. Koch, M. Leester-Schädel, S. Büttgenbach, and A. Dietzel, “Liquid polyimide as a substrate for aeronautical sensor systems,” (International Society for Optics and Photonics, 2015), 94352Y.

12. E. Koch and A. Dietzel, “Skin attachable flexible sensor array for respiratory monitoring,” Sens. Actuators, A 250, 138–144 (2016). [CrossRef]  

13. J. Eichstädt, G. R. B. E. Römer, and A. J. Huis in’t Veld, “Determination of irradiation parameters for laser-induced periodic surface structures,” Appl. Surf. Sci. 264, 79–87 (2013). [CrossRef]  

14. K. T. Ahmmed, E. J. Y. Ling, P. Servio, and A.-M. Kietzig, “Introducing a new optimization tool for femtosecond laser-induced surface texturing on titanium, stainless steel, aluminum and copper,” Opt. Lasers Eng. 66, 258–268 (2015). [CrossRef]  

15. X. Liu, D. Du, and G. Mourou, “Laser ablation and micromachining with ultrashort laser pulses,” IEEE J. Quantum Electron. 33(10), 1706–1716 (1997). [CrossRef]  

16. S. Amoruso, R. Bruzzese, X. Wang, N. N. Nedialkov, and P. A. Atanasov, “Femtosecond laser ablation of nickel in vacuum,” J. Phys. D: Appl. Phys. 40(2), 331–340 (2007). [CrossRef]  

17. J. M. Liu, “Simple technique for measurements of pulsed Gaussian-beam spot sizes,” Opt. Lett. 7(5), 196–198 (1982). [CrossRef]  

18. Y. Jee, M. F. Becker, and R. M. Walser, “Laser-induced damage on single-crystal metal surfaces,” J. Opt. Soc. Am. B 5(3), 648 (1988). [CrossRef]  

19. C. von der Heide, M. Grein, and A. Dietzel, “Femtosecond laser-contoured micro-strain gages,” Microelectron. Eng. 214, 81–86 (2019). [CrossRef]  

Supplementary Material (3)

NameDescription
Visualization 1       High-resolution image of line ablation grid without colour shading for different ablation zones.
Visualization 2       High-resolution image of areal ablation grid with pulse divider 1 without colour shading for different ablation zones.
Visualization 3       High-resolution image of areal ablation grid with pulse divider 4 without colour shading for different ablation zones.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (23)

Fig. 1.
Fig. 1. Gaussian fluence profile
Fig. 2.
Fig. 2. Pulse placement for an isotropic energy distribution. The red circles indicate the spot size of single pulses with their ${\omega _0}$ boundaries. Grey lines indicate the order in which ablation pattern is progressing (starting at the abscissa moving up, left to right).
Fig. 3.
Fig. 3. Accumulated fluence profile (red) obtained from a line of laser spot intensities.
Fig. 4.
Fig. 4. Accumulated 1D fluence profiles (red) for line ablation with: a) No pulse overlap and b) Pulse overlap such that peak accumulated fluence is practically still equal to maximum spot fluence c) Pulse overlap such that valley accumulated fluence equals maximum spot fluence d) Large pulse overlap with accumulated fluence profile can be considered “flat” and overshoot indicator is small e.g. below a certain threshold.
Fig. 5.
Fig. 5. 2D areal accumulated fluence distributions (${\phi _0}$ = 1 J/cm2, ${\omega _0}$ = 16 µm) with: a) Large gaps between individual spots such that fluence level locally drops to almost zero b) Overlapping pulses resulting in an accumulation of pulse fluences c) Pulse spacing similar to the pulse radius resulting in a practically homogeneous (flat) accumulated fluence profile far above ${\phi _0}$.
Fig. 6.
Fig. 6. Illustration of #-pattern: a) Explaining the maximal pulse distance ${l_{P,\textrm{max}}}$ to prevent uncovered areas of the specimen b) Providing an example pattern where the borders of overlapping pulse areas and spot centers are marked. The white area is receiving fluence contributions from four single pulses.
Fig. 7.
Fig. 7. 2D fluence distribution obtained with a “#-pattern” at ${l_P} = 20\; \textrm{}\mathrm{\mu}\textrm{m} = 0,63\; \textrm{}\cdot D\; \textrm{}$. As this is smaller than ${l_{P,\textrm{max}}}$, all areas between pulses are irradiated multiple times.
Fig. 8.
Fig. 8. LSM micrograph (x150) showing three distinguishable ablation zones for a NiCr-coated sample (100 nm) with two different imaging technics: a) 3D CLSM image b) Differential interference contrast-image (DIC)
Fig. 9.
Fig. 9. Results of single pulse experiments with increasing fluence (left: 3D CLSM image / right: direct optical view) a) First material reaction b) Cracks form inside the influence zone c) The coating lifts up d) Selective ablation showing no substrate damage
Fig. 10.
Fig. 10. Results of multiple (${N_P} = 6$) pulse experiments with increasing fluence (left: 3D CLSM image / right: direct optical view) a) Large influence zone with cracks b) Selective ablation zone, dust particles remain but no substrate damage can be detected c) Substrate damage starts to appear d) Selective ablation zone ring narrows
Fig. 11.
Fig. 11. D2-plot for a thin film selective ablation of 100 nm NiCr on polyimide. The diagram shows lines according to Liu's formula. They were obtained for each pulse repetition series by linear regression. Differences in the slopes result from measurement uncertainties and small instabilities of the ablation setup as the spot radius ${\omega _0}$ stayed constant during all experiments.
Fig. 12.
Fig. 12. 3D CLSM image proving that full selective ablation is achievable
Fig. 13.
Fig. 13. Plot of different zone radii evaluated from D2-experiments
Fig. 14.
Fig. 14. Plot of threshold fluences for different ablation zones with increasing number of pulses at 600 kHz. Incubation factor $\xi $ can be derived from these curves.
Fig. 15.
Fig. 15. Direct comparison of D2-experimental data confirms that ablation thresholds are influenced by laser frequency only at higher number of pulses ≥ 20.
Fig. 16.
Fig. 16. Comparison of multiple pulse experiments at increasing single pulse fluence. In case of ${N_P} = 8$ (two upper rows) the ablated diameter is similar and independent of laser frequency. At the bottom rows ${N_P} = 20$ the selective zone ablated at 600 kHz is always larger than at 16 Hz. Black areas indicate severe substrate damage where the crater is deeper than 200 nm.
Fig. 17.
Fig. 17. a) Line ablation grid. Inside each box frame, there are two ablated lines: Upper line is marked at 600 kHz, lower line is marked by placing individual pulses in a row yielding a slow frequency of 16 Hz (pulse fluence ${\phi _0}$ and the number of pulses on a single spot ${\tilde{N}_P}$ are the same for both lines). White numbers indicate cut width of lines [µm] ablated at 600 kHz. Side pictures: i) A fluence of 0,22 J/cm2 and a mark speed of 2500 mm/s allowed for a selective line ablation using the standard 600 kHz provided by the laser source. In contrast, the line ablated at 16 Hz is intermittent ii) The upper line shows severe substrate damage while the lower line is still in the selective region. Edgeline quality of lower line is low iii) At these parameters the upper line is selectively ablated whereas no ablation is observed in the case of slow pulse frequency. b) Comparison of cut width D calculated for mark speed v at increasing fluence ${\phi _0}$ using equations [Eq. (8),(9)] and [Eq. (11)] with cut width measured optically (see Visualization 1 for high-resolution image of the entire figure).
Fig. 18.
Fig. 18. Logarithmic plot of fluence threshold values as obtained by stationary spot experiments (multiple pulses) which define the selective ablation window (yellow) for selective line ablation.
Fig. 19.
Fig. 19. Map and histogramm of the specimen. In this case, a #-pattern with ${\omega _0}$=16 µm and a pulse spacing of ${l_P}$=20 µm (results from a mark speed of 3000 mm/s, a laser frequency of 600 kHz and a pulse divider of 4) was evaluated. Those regions that were hit by most pulses are marked yellow, whereas those areas struck only once remain blue. The histogram on the right shows how large the surface fraction is (in arbitrary units) that is hit by a certain number of laser pulses
Fig. 20.
Fig. 20. Areal ablation grid with pulse divider 1 at 600 kHz pulse frequency (see Visualization 2 for high-resolution image). Fluence is increased left to right. Values on the left indicate the mark speed dependent mean number of pulses ${\tilde{N}_P}$ for each row. Blue area marks all those fields where the coating has lifted in flakes from the substrate. All fields inside the selective ablation zone show some unwanted greyish substrate color. This is because of residues and dust. Nevertheless, they are electrically isolating. The red area defines those fields where severe substrate damage has occurred. The five micrographs on the bottom enlarge some of the fields inside the grid.
Fig. 21.
Fig. 21. Areal ablation grid with pulse divider 4 at 600 kHz pulse frequency (see Visualization 3 for high-resolution image). Fluence is increased left to right. Mean number of pulses ${\tilde{N}_P}$, pulse spacing, pulse overlap and mark speed are indicated on the left. Five ablation zones can be differentiated: no ablation, flaking, dotting, selective ablation and destruction regime. In the dotting regime the coating shows a uniform pattern of holes but retains areal conductivity. Two criterion boarder lines are marked inside the grid. The ten micrographs at the bottom enlarge fields inside the grid which are of particular interest.
Fig. 22.
Fig. 22. Accumulated fluence plot. The accumulated valley fluence is depicted for different mark speeds in relation to spot fluence (blue dotted lines). Yellow shaded area marks the selective ablation zone. If accumulated valley fluence is larger than the product of mark speed dependent pulse number and the corresponding ablation threshold then selective ablation is achieved. (Refer to Fig. 19 for correlation between pulse number and mark speed e.g. ${N_P} = 18$ for 1000 mm/s. The corresponding threshold fluence can be drawn from the plot in Fig. 18.)
Fig. 23.
Fig. 23. Micrograph of a pattern generated by successful areal ablation. Left side shows the sample directly after laser irradiation with parameters chosen according to the second criterion. Right hand side picture was taken after additional ultrasonic cleaning (to remove coating dust) at through light illumination only (thus areas with remaining opaque coating appear black). Any substrate damage would appear dark on the light yellowish polyimide substrate, but such destruction has clearly not occurred here. Because of the relatively small pulse overlap, minor positioning errors of the laser spot left some tiny coating flakes behind. Edgeline quality can be enhanced if areal ablation is followed by a contour line ablation with a larger pulse overlap.

Tables (1)

Tables Icon

Table 1. Methodolody of selective metallic thin film ablation

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

φ = l P 2 ω 0 l P = v f .
0,calc = 2 ω 0 = 1 , 83 × λ × E F L Beam = 40 , 1 μ m .
ϕ 0 = 2 × P π × ω 0 2 × f = 2 × E P π × ω 0 2 .
ϕ ( r ) = ϕ 0 × e 2 ( r ω 0 ) 2 .
ϕ ( x , y , n x , n y ) = ϕ 0 × e 2 ( ( x l P × n x ) 2 + ( y l P × n y ) 2 ω 0 2 ) n x [ 0 , , N x ] n y [ 0 , , N y ] .
Γ ( x , y ) = n x n y ϕ ( x , y , n x , n y ) .
θ P = Γ AP Γ AV Γ AV [ % ] .
D 2 = 2 ω 0 2 × ln ( ϕ 0 ϕ th ) .
ϕ th ( N P ) = ϕ th ( 1 ) × N P ξ 1 .
Γ AP = ϕ th ( N P ) × N P .
N ~ P = f × 0 v × P D .
1 st criterion: ϕ 0 ϕ th ( N ~ P ) .
Z ( x , y ) = n x n y Δ Z ( x , y , n x , n y ) n x [ 0 , , N x ] n y [ 0 , , N y ]
with Δ Z = { 0 if ( x l P × n x ) 2 + ( y l P × n y ) 2 ω 0 1 if ( x l P × n x ) 2 + ( y l P × n y ) 2 < ω 0 .
2 nd criterion: Γ AV N ~ P ϕ th ( N ~ P ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.