## 1. Introduction

Since the first laser was invented in 1960, many different types of lasers have been developed due to rapidly increasing areas of applications. In the selection of a laser for a particular application, the following parameters are usually considered: pulse duration, wavelength, fluence, pulse repetition rate, beam uniformity and stability, bandwidth, weight and size of laser system, lifetime and reliability, hardware design, and cost. Among them, the first five are more crucial for material processing. Sufficient laser pulse energy and proper repetition rate (if a pulse laser is selected) are needed for processing the material, wavelength that determines laser energy absorption, and pulse duration that controls the heat affected zone. The interaction among the above four effects, together with the beam quality, controls the material process quality.

Ultrashort pulsed lasers (USPLs) are a particular class of lasers whose pulse duration (*t _{p}*) is less than a few picoseconds (ps). This classification is based on the fact that for most solids the time to establish thermal equilibrium between electrons and phonons is of the order of picoseconds. For a 100 femtoseconds (fs) and 0.1 J laser pulse, for example, the peak power is one terawatt (1.0 × 10

^{12}W). If the laser pulse is focused onto an area of 0.1 mm

^{2}, its peak intensity can reach 10

^{15}W/cm

^{2}. In view of the extremely short pulse duration and high intensity, USPLs have been explored and demonstrated great potential for a large variety of applications in physics, chemistry, life science, materials, and engineering science in the past 20 years.

Ultrashort time and high peak power are the two unique features of USPLs such that a thin layer of material can be ablated before the absorbed laser energy diffuses into the surrounding bulk. Because of that, a material can be processed very precisely with minimal or even without collateral damage. Such concepts and applications to micro/nano precision processing have been demonstrated in laboratories since the mid-1990s. Furthermore, with the possibility of reaching an extremely high intensity by focusing ultrashort pulses in a tiny volume or area, multi-photon absorption is another niche that makes USPLs well suited for 3D micro and nano processing and fabrication. The advantages of USPL material processing include: (i) high precision [1, 2], (ii) minimal collateral damage, (iii) capability of processing practically any material, (iv) 3D fabrication [3–5], and (v) a single step for creating different surface structures, such as laser-induced periodic surface structures (LIPSS) [6]. Larger material removal rates than those with conventional pulsed lasers by ultrafast bursts of pulses have also been reported recently [7].

In this chapter, three areas of laser material processing are considered: (a) mechanisms of laser-material interaction and ablation, (b) modeling of laser material ablation, and (c) applications. For the mechanisms, the characteristic times of materials and the difference in material interactions with long, short and ultrashort laser pulses are discussed. The modeling is limited to laser-material interaction for metal and semiconductor materials, including ablation for metals and damage for semiconductors. Examples of the USPL application focus on modification of material properties below the surface of transparent materials and LIPSS for sensors, precision molds, and medical devices.

## 2. Mechanisms of laser-material interaction and ablation

The mechanisms of laser-material interaction and ablation mainly depend on the characteristic times and thermophysical properties of materials as well as laser parameters such as wavelength, pulse duration and fluence. This section will focus on laser interactions with industrial materials.

### 2.1. Laser interaction with metals

Laser irradiation to a metal solid is a two-step heating process [8]. The incident photons collide with free electrons that are confined within the laser light propagation path. The absorbed laser energy is then converted to kinetic energy of the excited electrons, allowing them to travel relatively long distance (tens or hundreds of nanometer (nm), depending on the materials) before collisions with other electrons. This is referred to as ballistic electron transport in metals. After the electron-electron collisions, their kinetic energy is spread among the electrons into the Fermi-Dirac distribution. This phase is referred to as thermalization, and is completed in a few tens of femtoseconds. The electron temperature (*T _{e}*) is measureable at this time. Afterwards, the electron thermal energy diffuses, through electrons, into a deeper part of the material. In the meantime, a part of the electron thermal energy is transferred to lattice (phonons) through electron-phonon collision. It takes few picoseconds (ps) for the excited phonons to become thermalized so that the lattice temperature (

*T*) can be measured. The energy exchange between the electrons and the lattice lasts for tens or hundreds of picoseconds, leading to thermal equilibrium (

_{l}*T*=

_{e}*T*) between the two subsystems (electrons and lattice), then common thermal diffusion drives the heat dissipation in the bulk.

_{l}During pulsed laser irradiation, since the heat capacity of electrons is about two orders of magnitude lower than that of the metal lattice, the electron temperature can shoot up to a very high temperature (e.g., 10^{3}–10^{4} K) while the lattice mainly remains in low temperature state. Hence, the electron and lattice temperatures can be quite different before thermal equilibrium is established. This two-step heating is particularly true for laser pulses that are shorter than 100 ps [8, 9]. **Figure 1** shows the time histories of the calculated electron and lattice temperatures at the irradiated surface of a copper foil by a single femtosecond laser pulse (duration 120 fs, wavelength 800 nm) of two fluences (0.5 J/cm^{2} and 3.0 J/cm^{2}). The smooth evolutions of both *T _{e}* and

*T*with the time shown in the case of a laser fluence of 0.5 J/cm

_{l}^{2}are the same as most simulation results of USPL heating reported previously. For heating by nanosecond (ns) or longer laser pulses, since the characteristic times of energy transfer between the excited electrons and lattice are much shorter than the laser pulse duration, thermal equilibrium between the two subsystems is nearly established while lasing. This is particularly true for laser pulses longer than 10 ns [10]. Therefore, one-step heating (

*T*=

_{e}*T*) is usually considered for the case of conventional, long-pulse laser heating.

_{l}### 2.2. Laser ablation of metals

The laser energy is primarily deposited within the optical penetration depth (a few tens of nm for metals) with an exponential decay. Because the penetration depth is so shallow, the gradient of the induced electron temperature during and shortly after USPL irradiation can be immense. This non-uniform electron temperature can generate a considerably large hot-electron blast force that could cause severe deformation in the cold lattice [9–12]. The non-uniform lattice temperature that develops later is another cause for the ultrafast deformation. Due to the small size of the affected thermal zone, the gradient of lattice temperature could be very sharp. As a result, severe thermal stresses are induced in the lattice. In the case of a relatively low laser fluence, the hot-electron blast force and thermal stress are the vital power that destroys the bulk material underneath the irradiated surface before the thermal energy is significantly conducted into the bulk [11, 13]. For high laser fluences, the lattice, on the other hand, can be superheated and subsequently undergo a metastable phase transformation from solid to liquid and then, further to vapor if the fluence is sufficiently high [14–17].

Using a simple criteria for material removal, Momma et al. [2] and Nolte et al. [18] deduced two well-known logarithmic functions for the USPL material ablation depth: for low fluences and for high fluences, where η is the optical penetration depth; ζ is the effective electron heat diffusion length; is the laser fluence; and and are the ablation threshold for low and high laser fluences, respectively. The failure of these two simple functions to describe the ablation rate for higher fluences [19] suggests that more mechanisms may be involved in the ablation process.

For USPL heating at laser fluences slightly exceeding the ablation threshold, only a thin layer of material is ablated by a single pulse due to the effects of the hot-electron blast and/or thermal expansion. Material ablation can take place during the lasing or after the laser pulse is off but before the next pulse arrives. When the material is removed, its associated thermal energy is also eradicated. Because very little or no molten liquid remains in the bulk material, hydrodynamic motion is negligible or never even occurs. As a result of a very thin layer ablated with little or no thermal damage to the bulk material by each laser pulse, the hole or crater resulting from this non-thermal ablation process can be very clean and precise [2].

Although non-thermal ablation can precisely process materials with minimal damage, it suffers from its very slow processing due to insufficient laser energy. This, therefore, creates a need for developing a more efficient means of processing. One potential approach to efficiently process a material target, for example, is to use a laser beam that irradiates high-fluence pulses (for high ablation rate) and then, is followed by relatively low-fluence pulses (for high precision). One can consider the low-fluence finishing as an “integrated” processing step, which results in high quality without the need for a separate post-processing. **Figure 2** shows the excellent quality of a hole produced by a femtosecond laser beam of 120 fs and 800 nm with low-fluence finishing.

For laser fluences well exceeding the ablation threshold, non-thermal ablation only occurs in the early time of the USPL heating [11]. Because the lattice temperature rises so drastically, the irradiated material rapidly undergoes phase transformation. Due to the fact that time is too short to allow necessary heterogeneous nuclei to form, the melted material is unable to boil [14, 20]. Instead, it is superheated past the normal boiling point to more or less the thermodynamic equilibrium critical temperature (*T _{tc}*). At that state, the tensile strength of the superheated liquid falls to zero and fluctuation in volume becomes dramatic, leading to a tremendous number of homogeneous nuclei being formed. As a result of the bubbles formed at such an extremely high rate, the subsurface layer of metastable liquid relaxes explosively into a mixture of vapor and equilibrium liquid droplets, which are immediately ejected from the bulk material. This thermal ablation mechanism is referred to as phase explosion. The phase explosion and material removal of metals caused by nanosecond [15] and USPL [16, 17, 21–23] heating have widely been investigated experimentally and theoretically. However, most studies investigated the thermal ablation within a low laser fluence regime (e.g., <10 J/cm

^{2}), and a high laser fluence regime was rarely investigated.

As discussed above, high stress (non-thermal) and phase explosion (thermal) are the two main mechanisms for metal ablation. Non-thermal ablation provides a precision controllability which is crucial to micro- and nano-scale material processing. Another advantage is that little or no post processing is needed. Therefore, USPLs with a fluence slightly exceeding the ablation threshold are a good choice for microfabrication for which precision and small sizes are required. On the other hand, USPLs with fluences well beyond the ablation threshold may not offer an advantage for precise material processing due to the strong thermal (superheating and bubble forming) and mechanical (fluid ejection) effects that accompany them. However, there could be a trade-off between precision and efficiency of laser material processing [18].

The LIPSS created on a metal surface by multiple ultrashort laser pulses with a fluence near the ablation threshold are another interesting result of laser matter interaction [6]. Formation of micro- and nano-sized periodic ripples and column arrays are dependent on the number of laser pulses, pulse duration, fluence, wavelength, radiation polarization, and the ambient conditions. The formation of periodic-like surface structures on metal can be attributed to the interference between an incident femtosecond laser beam and the surface scattered wave [24–27].

### 2.3. Laser interaction with semiconductor and dielectric materials

The difference between semiconductor and dielectric materials is the number of free electrons in the conduction band. When a semiconductor or dielectric material is irradiated by a laser, the electrons in the valence band absorb the photon energy and then transit to the conduction band via single- or multi-photon absorption, depending on the photon energy (*hv*, where *h* is the Planck constant and *v* is the laser-light frequency) and the band-gap energy (*E _{g}*) of the material. The interband transition of the electrons creates holes in the valence band. The excess energy of the created electron-hole pairs,

*khv*−

*E*(

_{g}*k*= 1 for single-photon absorption and

*k*> 1 for

*k*-photon absorption), is the kinetic energy of the excited carriers that determines the carrier temperature. As the electrons and holes undergo temporal and spatial evolution, some of them recombine, for example, through the three-body Auger process. Meanwhile, additional electron-hole pairs could be generated via impact ionization from those free carriers with kinetic energy that is equal to or greater than the band gap. The excited electron-hole pairs thermalize to the Fermi-Dirac distribution via carrier-carrier collisions on a <100 fs timescale after photon absorption. In the meantime, thermalization between the carriers and phonons proceeds until the thermal equilibrium state is inevitably established.

### 2.4. Laser ablation of semiconductor and dielectric materials

Three main factors that affect damage and ablation of semiconductor and dielectric materials by lasers are: (a) free-electron density, (b) free-electron kinetic energy, and (c) carrier and phonon temperature. Increasing of the carrier and phonon temperature leads to the hot-electron blast force, thermal expansion, superheating, and phase explosion that are basically the same as those in metal materials described previously.

At the very beginning of USPL excitation, the density rise of the electron-hole pairs generated grows exponentially after a very short transient [28]. Depending on the peak laser intensity, for example, 12 TW/cm^{2}, this could result in a high density (10^{21}–10^{22} cm^{−3}) of the electron-hole pairs. When the kinetic energy of an excited electron exceeds the surface barrier and the momentum component normal to the surface of the material is positive, the excited electrons can escape from the irradiated material into the surrounding air or vacuum. Consequently, the irradiated surface gains highly positive charges, leading to a repulsive force. If the repulsion force between ions is greater than the lattice binding strength, the atomic bonds are broken and a subsurface layer of the material is disintegrated. Damage and removal of material by this repulsive force is referred to as Coulomb explosion [29]. This ablation mechanism is non-thermal because Coulomb explosion occurs prior to significant heating of the phonon subsystem. Like the non-thermal ablation in metals, this is the mechanism that can lead to high precision in processing semiconductor or dielectric materials with USPLs.

It is quite natural to use focused USPLs for 3D micromachining of semiconductor and dielectric materials because of their extremely high peak intensity. With the advantages of superior localization and great optical penetration depth (longer laser wavelength), multi-photon absorption, in general, is a better mechanism for 3D micromachining, especially for processing inside bulk material. One-photon absorption may only be suitable for surface processing due to linear absorption and smaller penetration depth (shorter laser wavelength).

When a USPL beam is tightly focused inside a transparent material, multi-photon ionization, tunneling ionization, and avalanche ionization could occur in the focal volume. Like metals, the excited electrons transfer their energy to the ions via collision. With some time delay, the electrons and ions eventually reach thermal equilibrium. When the laser intensity is below a certain threshold, electrons recombine with holes in a non-radiative way. For laser intensities exceeding the certain level (referred to as plasma threshold), a plasma spark, resulting from the recombination of high-density electron plasma with holes, can be observed. The spark is associated with an optical breakdown of the transparent materials, having a characteristic duration of the order of 10 ns. After recombination, the thermal energy diffuses away from the focal volume on a microsecond timescale. At even higher intensities, laser-induced permanent modification of material properties is observed [4]. If the laser intensity is above the modification threshold, the hot electron-hole plasma and ions explosively expand from the focal volume into the surrounding material.

With tight focusing of ultrashort laser pulses, different kinds of permanent structure changes in a transparent material can be made by adjusting the incident pulse energy, in principle. Fabrication/machining/processing can be performed on the surface or within the bulk of the material by moving the laser focus along the desired paths. More detailed information can be found in the papers cited in the review paper [30].

## 3. Modeling of laser material ablation

In this section, we focus on two-step heating models for ultrafast thermal transport, ablation models for metals, and damage models for semiconductor materials.

### 3.1. Two-temperature model for metals

The two-temperature (2T) model was pioneered by Anisimov et al. [31] in 1974 to describe a two-step heating process of metals subjected to short-pulse laser irradiation. It was not until the early 1990s that USPL material interactions received considerable attention. Since then, numerous modified versions of 2T have been proposed [9, 32]. Based on the Boltzmann transport approximation, a semi-classical 2T model is derived for the dynamics of electron concentration (*n*), mean (drift) velocity , and average energy for metal materials [9]:

In the above equations, *m _{e}* is the electron mass, is the heat flux vector in the electron subsystem,

*k*is the Boltzmann constant,

_{B}*e*is the electron charge, μ

*represents the mobility of the electrons, β*

_{o}*is the temperature-dependent parameter for free electrons,*

_{e}*C*is the heat capacity of the electron,

_{e}*G*represents the electron-phonon coupling factor,

*S*is the volumetric laser heat source, and is the divergence operator. Since the electron relaxation time of metals is around ten femtoseconds, the constitutive relation for and is given by [32]:

where τ* _{e}* is the electron relaxation time (the mean time for electrons to change their states), and

*K*is the electron thermal conductivity.

_{e}For the lattice subsystem, the thermal transport equation includes an energy exchange with the electrons and a thermal relaxation effect in a general case [33]:

In Eqs. (5) and (6), is the heat flux vector in the lattice, τ* _{l}* is the relaxation time in phonon collisions, and

*C*and

_{l}*K*are the heat capacity and thermal conductivity of lattice, respectively.

_{l}By neglecting the electron drift velocity, energy equation (3) is simplified to:

Thus, the semi-classical 2T model, seen in Eqs. (1)–(6), is reduced to the dual-hyperbolic 2T model (HH2T) [33].

For pure metals, *K _{l}* is much smaller than

*K*. The hyperbolic 2T model (H2T) [32] can be retrieved from the HH2T model by neglecting the heat conduction in the lattice, i.e., in Eq. (5). Further neglecting the electron relaxation effect in the H2T model leads to the parabolic 2T model (P2T) [8, 34].

_{e}After the thermal equilibrium, and . Combining Eqs. (5) and (7) yields the macroscopic energy balance equation:

The above energy equation (8) together with Fourier's heat conduction law is a well-known, parabolic one-temperature heat conduction model (P1T). When the Cattaneo equation [35] is considered, the result becomes a hyperbolic (thermal wave) one-temperature model (H1T). The above 1T models represent one-step heating and have been widely used for long-pulse laser heating since the electrons and lattice are assumed to be in thermal equilibrium instantaneously when a medium is heated.

The thermophysical properties, *C*, *K* and *G*, control thermal response in laser-irradiated material. The expressions for these three thermophysical properties can be found in Refs. [10, 36]. Recently, Lin et al. [37] presented data *C _{e}* and

*G*for different metals over a wide range of electron temperatures from room temperature to 5 × 10

^{4}K. A formula for

*K*up to Fermi temperature was given by Anisimov and Rethfeld [38].

_{e}The mechanical stresses caused by the hot-electron blast force and the thermal stresses induced by non-uniform temperatures in the lattice can be solved by the equation of motion of lattice:

where *u _{i}* (

*i*=

*x*,

*y*, and

*z*) is the displacement vector of a material (lattice) point, σ

*is the stress tensor at the point, and the subscript “,” denotes the spatial derivative. The last term on the right-hand side of Eq. (9) is the hot-electron blast force. The thermal strains resulting from a non-uniform lattice temperature are included in the stress-strain relations. The deformation can be linearly elastic or non-linearly plastic.*

_{ij}The exchange of thermal and mechanical energy should not be neglected due to the extremely high strain rate (~10^{9} s^{−1}). Thus, the energy equation of lattice, Eq. (5), is re-written as:

where is the Lamé constant, is the shear modulus, is the thermal expansion coefficient, and is the volume dilatation. In addition, the strain rate effects on plastic deformation and failure could be substantial, and the existing stress-strain relations and fracture strength may be insufficient for ultrafast deformation. To accurately predict non-thermal material ablation by stresses, the strain rate effects on material behavior should be investigated further and quantified.

### 3.2. Laser heat source

The most popular laser beams are Gaussian in both space and time. For simplicity, let us consider an incident laser beam that is normal to the material target. The volumetric laser heat source *S* in Eqs. (3) and (7) is given as [39]:

where *r* is the radial distance from the center of the beam, *z* is the coordinate in the direction of laser beam propagation, with denoting the pulse energy, *r _{o}* being the beam radius defined at the

*e*

^{−2}distance,

*R*(

*r*, 0;

*t*) is the surface reflectivity of the material, α is the absorptivity coefficient,

*t*is the full width half maximum (FWHM) of the Gaussian temporal pulse,

_{p}*m*is an integer, and

*β*= 4ln(2). The lasing, assumed to start at

*t*= 0, reaches its peak power at

*t*=

*mt*, and ends at

_{p}*t*= 2

*mt*. The value of

_{p}*m*can be set at 2 or 3 since the laser energy outside this period of time is inconsequential. If the effects of the hot electrons ballistic motion [40] are taken into account, the laser heat source is modified to:

where is the temperature-dependent optical penetration depth and δ* _{b}* is the ballistic electron penetration depth. The temperature-dependent optical properties,

*R*and α, in Eqs. (11) and (12) can be determined from optical properties and the Fresnel function.

Recently, an extended Drude model [41] was proposed to characterize the temperature-dependent *R* and α, and a critical point model with three Lorentzian terms for interband transition [42] for copper. The numerical results show that when a metal is irradiated by a USPL at high fluence, the dynamic changes of *R* and α during the laser irradiation could significantly alter laser energy absorption and the distribution of laser heat density. Although it is only true for low electron temperatures, the constant *R* and α at room temperature have been widely employed in 2T modeling. In this case, the integral term becomes a constant, *z/δ* in Eq. (11) and *z/(δ + δ _{b})* in Eq. (12).

### 3.3. Comparison of different thermal models

**Figure 3** compares the predictions of four thermal models, HH2T, H2T, H1T, and P1T, for lattice temperatures at the front surface of a 200 nm gold film heated by a laser pulse of *F*_{o} = 0.5 J/cm^{2} and *t _{p}* = 100 ps and 1 ns [10]. Evidently, both of the one-temperature models significantly overestimated the bulk temperature for the 100 ps pulse; hence, they are inadequate for simulating USPLs heating. The results, shown in

**Figure 3(b)**, confirm that the Fourier heat conduction model is sufficient for long-pulse laser heating. In fact, the difference among the four models becomes indistinct for a 10 ns laser pulse [10].

### 3.4. Ablation models for metals

The 2T models discussed above only characterize thermal transport in a solid metal. For material ablation, other physical processes such as phase transitions and material removal need to be considered. Different approaches have been proposed to simulate laser material ablation, including ultrafast thermoelasticity [11], dynamics of thermal ablation [17], hydrodynamic modeling [43], molecular dynamics [22], etc. Although those numerical analyses [17] show a certain degree of success in comparison with experimental measurements, none of the aforementioned models have been validated thoroughly with experimental data for USPL ablation over a wide pulse duration and fluence.

A comprehensive USPL ablation model is needed to be able to accurately describe the entire ablation process and predict the ablation depth and rate. It should be developed based on universal regulations and methods that allow for the description of non-equilibrium responses and cover the complex physical phenomena for the entire laser process. Such phenomena include photon-electron interaction, laser-pulse propagation and ionization, phase transitions, superheating through a metastable liquid phase, rapid nuclei formation, explosion and dynamics of homogeneous nuclei (bubbles), generation of recoil pressure, hydrodynamic motion, formation and dynamics of the plasma plume, laser-plasma interaction, radiation from the resulting plasmas, condensation and re-solidification of the liquid/vapor, etc.

To accurately simulate thermal ablation by USPLs, a semi-classical 2T model integrated with models of phase transformations for ultrafast melting, evaporation and re-solidification and a phase explosion model for ejection of metastable liquid and vapor was attempted recently [44]. When superheated liquid temperature reaches 0.9*T _{tc}*, phase explosion is assumed to take place [20] and that material, including both electrons and lattice, is then removed under the assumption of phase explosion. Once phase explosion no longer occurs, vaporization could continue until the lattice temperature drops significantly.

**Figure 4**shows the time history of the ablation depth for copper foil irradiated by a single femtosecond laser pulse (duration 120 fs, wavelength 800 nm) of different fluences. The steep occurrences of material ablation result mainly from phase explosion, while the sloping parts result from vaporization. As shown in

**Figure 4(b)**for the case of laser fluence 6.1 J/cm

^{2}, the ablation depth by phase explosion and vaporization are 258.4 nm (up to 310 ps) and 5.9 nm (from 310 to 600 ps), respectively. The results shown here indicate that for high fluences, phase explosion is the dominating material ablation mechanism in USPL material ablation. It is noted that this ablation model is not yet fully comprehensive.

### 3.5. Ultrafast transport dynamics models for semiconductors

To model the transport process of a large number of hot electrons, holes, and phonons in a semiconductor, the formalism must be based on the principle of statistical mechanics. A self-consistent model for transport dynamics in semiconductors caused by USPL heating can be found in [46]. Based on the relaxation-time approximation of the Boltzmann equation, the rate equations are derived for the dynamics of electron-hole carrier number density and current, ambipolar energy current, carrier energy, and lattice energy.

#### 3.5.1. Rate equation for carrier pairs

The balance equation for the electron-hole pairs number density generated by a laser pulse for two-photon absorption is:

where β_{1} is the avalanche coefficient, β_{2} is the 2-photon absorption coefficient, *γ* is the Auger recombination coefficient, *θ* is the impact ionization coefficient, and *J* is the carrier current. The magnitude of β_{1} is much greater than that of β* _{k}* (

*k*> 1). For silicon at room temperature, for instance, and . The last three terms on the right-hand side of Eq. (13) represent Auger recombination, impact ionization, and the loss due to carrier current, respectively. The Auger recombination is described as the time when an electron-hole pair recombines, giving up its energy to a third electron in the conduction band; this reduces the carrier number density. The reverse effect is impact ionization. Due to the nonlinearity nature (

*I*), multi-photon absorption could be more efficient and occur faster than single-photon absorption when laser intensities are sufficiently high.

^{k}#### 3.5.2. Rate equation for carrier energy

The balance equation for the electron-hole pairs energy produced by a laser pulse for two-photon absorption is:

where is the heat capacity of electron-hole pairs, is the free carrier absorption cross-section, and *W* is the ambipolar energy current which is the sum of the carrier energy currents in electrons and holes.

#### 3.5.3. Rate equation for lattice energy

For semiconductor materials the thermal conductivity of the lattice is comparable with the bulk value. It is necessary to consider the thermal transfer in the lattice. The lattice energy rate is expressed as:

The three rate equations (13)–(15), together with the constitutive equations for the carrier pair current and the ambipolar energy current [46], compose a complete self-consistent model for ultrafast transport dynamics in semiconductors subjected to USPL irradiation.

**Figure 5** shows the time histories of the density and temperature of carriers and the lattice temperature at the incident surface of a silicon sample irradiated by a 500 fs laser pulse with fluence 0.005 J/cm^{2} [46]. The peak laser power is at *t* = 1.5 ps (*m* = 3). As shown in **Figure 5**, the carrier temperature reaches its maximum at about 0.68 ps and the number density reaches its maximum at about 2.21 ps. It should be noted that when the peak carrier temperature occurs, the laser power is only at about 0.06% of its maximum. The contradictory intuition of this high temperature is attributed to a very small amount of electron-hole pairs that are created during this early time. Consequently, the carrier heat capacity is very small, thereby leading to a rapid, noticeable rise in the carrier temperature although the net carrier thermal energy is quite low. As time is prolonged, the carrier density increases drastically since much more laser energy has been absorbed. As a result, the carrier heat capacity becomes greater and greater. Meanwhile, the energy loss, due to the change in carrier energy density (the second last term on the right hand side of Eq. (14)), becomes pronounced, and so does the loss due to thermal transfer from the carriers to the lattice. Those changes make the rate of change of the net energy unable to remain positive at some point and thereafter, even though considerable laser energy is absorbed. This explains why the carrier temperature quickly reaches its peak very early in the laser irradiation and then falls. The numerical result in Ref. [46] also shows that one-photon absorption and Auger recombination are two crucial factors that alter the carrier's number density.

### 3.6. Damage models for semiconductors

For pulse durations longer than the carrier–lattice energy relaxation time (a few picoseconds), it has generally been accepted that the solid–liquid phase transition by high laser fluences is a thermal (melting) process. On the other hand, experiments with femtosecond laser pulses of high fluence have demonstrated an ultrafast phase transformation on a subpicosecond time scale. The mechanism of this ultrafast phase transformation is nonthermal melting, differing from the above thermal melting.

Among the continuum models, there are two approaches that are frequently used in evaluation of the damage threshold for semiconductor materials. One approach employs a single rate equation to evaluate electron density in the conduction band [47]. The equation includes ionization rates for single and multi-photon absorption, avalanche ionization, and other relaxation terms. Damage is assumed when the calculated electron density exceeds a critical value, which is often determined semi-empirically by matching with the experimental data. In most cases, the critical value used differs from that of the critical density of electron-hole plasma that makes the plasma opaque (optical breakdown). The other approach uses a self-consistent model in which the rate equation for the carrier number density is coupled to the energy balance equations for both the carriers and phonons [46]. The onset of damage is determined by which condition, carrier number density, or lattice temperature is first met. **Figure 6** shows the comparison between the theoretical damage fluence thresholds with measured values for Si [9, 21]. It appears that the self-consistent model agrees fairly well with the experimental data for laser pulse durations up to several picoseconds; on the other hand, the carrier density model is only in good agreement for laser pulse durations of subpicoseconds [46]. In view of both approaches failing to predict the damage fluence thresholds for longer laser pulses, a more robust model or approach is suggested.

Modeling of the non-thermal ablation of semiconductor materials via Coulomb explosion can be found in the study by Stoian et al. [29]. On the other hand, a comprehensive modeling of thermal ablation via phase explosion, to the authors' knowledge, has not yet been reported in open literature.

## 4. Applications

As aforementioned, the great applications of USPLs are attributed to two unique features: ultrashort pulse duration and extremely high intensity. The benefits of USPL material processing include high precision with a minimal heat affected zone, the ability to locally modify below the surface of transparent materials, the single step for creating different surface structures, etc. A great number of applications of USPLs have been proposed, for example, to industrial technologies and bio/medicine [48, 49]. Recently, applications of LIPSS have been of particular interest.

The ability to locally modify material properties below the surface of transparent materials is an important area of USPL application. One interesting example is the repair of a packaged.

TFT-LCD panel. In this repair, a femtosecond laser beam passes through the polarizer, glass substrate, and then, focuses on the color filter (color photoresist) corresponding to the hot pixels such that the phenomenon of nonlinear multi-photon absorption can be induced within the localized area. Then, the property of the color photoresist layer is changed. The hot pixels become dead pixels, and any defective pixels inside the panel of the packaged LCD can be repaired directly. **Figure 7** shows a red colored filter after laser repair. The image is captured from different focal positions. As shown in **Figure 7(a)**, the red colored filter (marked by an arrow) has been blackened. In **Figure 7(b)**, however, the polarizer corresponding to the red colored filter (marked by an arrow) is not damaged after laser irradiation. Accordingly, the defective pixel is eliminated by the transformation of the bright one into a dark one.

LIPSS by USPLs allows for a new range of functionalized surface processing that permits novel applications for sensors, medical devices, precision molds, etc. For sensor applications, fabrication of Ag nanostructure-covered surface structures can be used for surface-enhanced Raman scattering (SERS) application. The SERS intensities of rhodamine 6G (R6G) at LIPSS-treated Ag substrates are 15 times greater than those at non-treated Ag substrates [50].

For medical devices applications, LIPSS by USPLs can be employed to fabricate electrosurgical blades [51, 52]. The average temperatures of tissue cut by conventional and LIPSS-treated electrosurgical blades are 145°C and 116°C, respectively. This corresponds to an improvement of about 29°C in temperature. As a result, there is less damage to tissue that is cut by the proposed electrosurgical blade.

LIPSS by USPLs is also employed to fabricate medical stainless steel [53]. As shown in **Figure 8**, the attachment assay with 72 h in a NIH3T3 culture demonstrates that significantly more cells are attached to the LIPSS-treated sample. This means that the cells spread over the treated surface more rapidly than over the non-treated surface. LIPSS-treated stainless steel is believed to possess better biocompatibility than stainless steel without surface modification.

For precision mold applications, for example, LIPSS by USPLs can be employed to fabricate structured molds and fast replication of large-area hierarchical micro/nano structures (lotus-leaf-like patterns) on a plastic part by injection molding [54]. Compared with an ordinary plastic surface, the contact angle of the structured plastic parts is increased by 38%, from 97 to 134°.

## 5. Conclusion

In reviewing laser micro- and nano-structuring of materials through pulse laser ablation, three main areas have been covered, including laser system development, experimentation and modeling, and applications. In this single chapter, we focused on laser material ablation mechanisms, modeling, and some applications of material properties modification and LIPSS that have been developed or demonstrated. Four areas in the modeling of laser matter interaction have been presented, including temperature response with temperature-dependent material properties and optical properties, ablation models for metals, ultrafast transport dynamics models for semiconductors, and damage models for semiconductors. Understanding the mechanisms and modeling of laser material interaction with materials can provide an insight for optimizing the processing parameters for precisely machining a variety of materials and can allow the development of high-value and innovative laser process techniques for a variety of applications.