**Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics**

Roman V. Dyukin, George A. Martsinovskiy, Olga N. Sergaeva, Galina D. Shandybina, Vera V. Svirina and Eugeny B. Yakovlev

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/46237

## **1. Introduction**

196 Laser Pulses – Theory, Technology, and Applications

*Elelctron.,* 24, 449 (1994)].

*Elelctron.,* 25, 799 (1995)].

(1994)].

(2011)].

*[Soy. J. Quantum Elelctron.,* 21, 643 (1991)].

[6] Zaykin A.E., Levin A.V., Petrov AL., Stranin S.A. *Kvantovaya Elektron.,* 18, 708 (1991)

[7] Zaykin A.E., Levin A.V., Petrov AL. *Kvantovaya Elektron.,* 21, 486 (1994) *[Quantum* 

[8] Vasil'ev S.V., Ivanov A.Yu., Lyalikov AM. *Kvantovaya Elektron.,* 22, 830 (1995) *[Quantum* 

[11] Rabinovich MI., Trubetskov DI. *Vvedenie v teoriyu kolebanii i voln* (Introduction to the

[12] Goncharov V.K., Kozadayev K.V. *Inzhenerno-fizicheskii zhurnal,* 83, 80 (2010) *[Journal of* 

[13] Chumakov AN., Bereza NA., Hu Dz.D., Bosak NA., Guo Z.H., Hie K.K. *Jnzhenernofizicheskiizhurnal,* 84, 524 (2011)*[Journalof Engineering Physics and Thermophysics,* 84, 567

[15] Goncharov V.K., Kozadayev K.V., Shchegrikovich DV. *Inzhenerno-flzicheskiizhurnal,* 84,

[16] Ashmarin I.I, Bykovskii Yu A, Gridin V A, et al. *KvantovayaElektron. (Moscow)* 6 1730

781 (2011) *[Journalof Engineering Physics and Thermophysics,* 84, (2011)].

Theory of Oscillations and Waves) (Moscow: Nauka, 1984).

*Engineering Physics and Thermophysics,* 83, 90 (2010)].

[14] Klimentov SM. et al. *Laser Phys.,* 8 (6), 1(2008).

(1979) *[Sov. J. Quantum Electron.* 9 1019(1990)].

[9] Dorofeev IA., Libenson MN. *Opt. Spektrosk.,* 76,73(1994) *[Opt. Spectrosc.,* 76, 66 (1994)]. [10] Chumakov AN. et al. *Kvantovaya Elektron.,* 21, 773 (1994) *[Quantum Electron.,* 24, 718

> Femtosecond lasers bring new opportunities in a variety of technological applications [1] in micro- and nanotechnologies, including electronics, mechanics, medicine and biology. Technologies, based on femtosecond effects, are used, for example, to make light absorbers for solar energy devices [2], for direct fabrication of integrated optical components [3], enhancing performance of photo-electronic devices [4], friction reduction and improvement of mechanical wear resistance [5], surface conditioning of medical implants [6], etc. Further development of the above technologies requires deeper understanding of the physical processes occurring under the ultrashort laser pulse action on different materials.

> Changes in the material optical properties under the action of intense radiation represent the key feature of the interaction of laser radiation with condensed media. Dynamics of optical properties of solids under the action of femtosecond laser pulse determines a number of physical effects which are of the great interest for both fundamental science and new applications. In particular, the feedbacks, which are being formed in this case, fundamentally change the properties of condensed matter [7].

> During the action of femtosecond pulse on solids the electronic subsystem undergoes intensive photoexcitation while the lattice stays cold. The processes of excitation of the electrons and release of the absorbed energy are spaced in time. High intensity of the laser radiation results in modification of the state of the electron subsystem thus significantly changing the optical properties of the medium [8]. Studies of the femtosecond pulses effects on semiconductors and insulators [9, 10] showed that the concentration of nonequilibrium carriers generated by laser radiation is so high that the surface layer acquires properties of the metals during the pulse.

© 2012 Svirina et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2012 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

It was observed experimentally that semiconductors can be disintegrated during the femtosecond laser pulse action [11]. In [11] it was proposed a mechanism of destruction based on the crystal lattice destruction by the electric field resulting from the violation of quasi-neutrality in the irradiated area due to the external electron emission - Coulomb explosion. Conditions for Coulomb explosion occurrence in metals were not found.

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 199

It is known that the light is absorbed by the conduction electrons in metals. Depending on the concentration of conduction electrons *n* and the wavelength of the incident light *λ* the electrons can be considered as free provided *n*>>1/*λ*3, and the free electron model cannot be applied if *n*<<1/*λ*3. In typical conditions in metals, where *n*≈1022 cm-3, and for the light wavelength of *λ*≈1μm the free electron model works well. It is also well known that the electron gas in metals is degenerated at practically all temperature range, and the distribution function of the electron gas just slightly differs from the distribution function at

The Fermi energy *εF* for metals is very high. For example, for copper *εF*=7.1 eV, for silver *εF*=5.5 eV. For this reason heat effects engage the electrons, whose energy lies in a narrow

The concentration of electrons, which absorb the incident radiation in metal, can be estimated by using the relation *n*=*hνN*/*εF*, where *ν* − the frequency of the incident light. For example, *n/N=*25% of the conduction electrons are affected in copper under *hν=*1.7 eV.

The electron transfer part of excessive energy, which they receive due to light absorption, though their collision to other electrons, ions and lattice defects (dislocations, grain boundaries, etc.). Heating of the metal is determined by efficiency of the collisions, which depends on the particles that exchange with energy. Typically the following relation takes place: *νee*>*νei*>*νep*, where *νee* − the frequency of electron - electron collisions, *νei* − the frequency of collisions between electrons and phonons, *νep* − the frequency of electrons collisions with impurities and defects in the metal. The electron-phonon mechanism is quite feasible, since the Fermi energy of electrons in metals is high and is essentially represented by the energy of translational motion of free electrons. The Fermi velocity is *υF*=(2*εF/me*)1/2=1,5·108 cm/s is

also high compared to the speed of sound in metals, that is typically: *υ<sup>F</sup>* ~ 105 cm/s.

the absorbed energy transfer in metal is

the conditions:

A single collision may not be enough for the electron, which absorbed a photon, to release the excessive energy, i.e. energy relaxation is a multi-stage process of a diffusion character. For this reason, the energy redistribution occurs not only at the skin layer (*δ<sup>s</sup>* ~ 10-6 cm), but in the deeper layer of *ld*=(*D*/*νee*)1/2=10-5 cm, where *D* is the electron diffusion coefficient.

Thus there is a heating of the metal. The heat, which is released in the layer *lD*, is further transferred depthward into the material through heat conduction. The characteristic time of

*τei* = 1/*νei* = 10-12-10-11 s. The electron gas and the lattice of the metal are two weakly interacting subsystems. Under

> , *ee t* >>τ

, *ii ei* ν <<ν

(*τee* − time to establish an equilibrium energy distribution in the electron gas) and

(1)

(2)

energy range 2 *<sup>B</sup>* ≈ *k T* (*kB -* Boltzmann constant) near the Fermi level.

the absolute zero.

The currently used experimental approaches such as femtosecond pump-probe technology [12, 13] and mass-spectroscopy (for example, [14]) provide measurement of integral characteristics, but have limited capability for retrieval of dynamics of the processes. The limitations of the experiment approach are being compensated by extensive use of mathematical modeling, where the fast non-linear processes are simulated in a wide range of the initial data.

Two-temperature model, which is traditionally used to describe the ultrashort laser pulse interaction with matter, has proved its validity in various conditions. The phenomenological two-temperature model of parabolic type was proposed in the 50's of last century by M.I. Kaganov, I.M. Lifshitz, L.V. Tanatarov [15]. It has been used by S.I. Anisimov to describe transient phenomena in a nonequilibrium electron gas and lattice under the submicrosecond laser action [16]. The model represents the primary approach to mathematical description of the nonequilibrium heating of the condensed medium by the action of short- and ultrashortpulse laser radiation. According to the model it is assumed that the energy absorbed by free electrons increases their temperature, then the interaction of the heated electrons with the lattice results in increasing lattice temperature. Heat transfer takes place through the heat conduction mechanisms.

When analyzing the effect of femtosecond laser pulses on matter one has to consider the following: applicability of two-temperature model to description of the electron temperature, which is determined by the electron equilibrium and applicability of the notion of temperature; and taking into account multi-quantum effects in description of electron emission.

**The applicability of a two-temperature model for description of action of femtosecond laser pulse on metals.** Femtosecond pulse action on the metal can be generally described as follows. Absorption of photons by free electrons in metals results in increase in the electron kinetic energy and the energy distribution becomes nonequilibrium. This well-known feature determines the behavior of metals in a wide spectral range. One can use the diffusion approximation for obtaining qualitative characterization of photoexcitation of solids during the action of femtosecond pulse. With this approach the distribution of free electrons is described by the integral concentration *n*(*z*, *t*), which varies in time and space (along the axis *z*, directed into the material) due to photo- and thermionic emission from the surface layers and electron diffusion. It is assumed that the thermalization of the electron gas occurs so fast that the notion of electron temperature can be immediately applied. "Hot" electrons contribute to the photo- and thermal emission: they withdraw a part of the energy stored in the electronic subsystem, thus reducing its temperature and, eventually, the temperature of solid as a whole. At the same time, change in the electrons concentration in the surface region results in a change in the optical characteristics of the material.

It is known that the light is absorbed by the conduction electrons in metals. Depending on the concentration of conduction electrons *n* and the wavelength of the incident light *λ* the electrons can be considered as free provided *n*>>1/*λ*3, and the free electron model cannot be applied if *n*<<1/*λ*3. In typical conditions in metals, where *n*≈1022 cm-3, and for the light wavelength of *λ*≈1μm the free electron model works well. It is also well known that the electron gas in metals is degenerated at practically all temperature range, and the distribution function of the electron gas just slightly differs from the distribution function at the absolute zero.

198 Laser Pulses – Theory, Technology, and Applications

of the initial data.

conduction mechanisms.

emission.

It was observed experimentally that semiconductors can be disintegrated during the femtosecond laser pulse action [11]. In [11] it was proposed a mechanism of destruction based on the crystal lattice destruction by the electric field resulting from the violation of quasi-neutrality in the irradiated area due to the external electron emission - Coulomb

The currently used experimental approaches such as femtosecond pump-probe technology [12, 13] and mass-spectroscopy (for example, [14]) provide measurement of integral characteristics, but have limited capability for retrieval of dynamics of the processes. The limitations of the experiment approach are being compensated by extensive use of mathematical modeling, where the fast non-linear processes are simulated in a wide range

Two-temperature model, which is traditionally used to describe the ultrashort laser pulse interaction with matter, has proved its validity in various conditions. The phenomenological two-temperature model of parabolic type was proposed in the 50's of last century by M.I. Kaganov, I.M. Lifshitz, L.V. Tanatarov [15]. It has been used by S.I. Anisimov to describe transient phenomena in a nonequilibrium electron gas and lattice under the submicrosecond laser action [16]. The model represents the primary approach to mathematical description of the nonequilibrium heating of the condensed medium by the action of short- and ultrashortpulse laser radiation. According to the model it is assumed that the energy absorbed by free electrons increases their temperature, then the interaction of the heated electrons with the lattice results in increasing lattice temperature. Heat transfer takes place through the heat

When analyzing the effect of femtosecond laser pulses on matter one has to consider the following: applicability of two-temperature model to description of the electron temperature, which is determined by the electron equilibrium and applicability of the notion of temperature; and taking into account multi-quantum effects in description of electron

**The applicability of a two-temperature model for description of action of femtosecond laser pulse on metals.** Femtosecond pulse action on the metal can be generally described as follows. Absorption of photons by free electrons in metals results in increase in the electron kinetic energy and the energy distribution becomes nonequilibrium. This well-known feature determines the behavior of metals in a wide spectral range. One can use the diffusion approximation for obtaining qualitative characterization of photoexcitation of solids during the action of femtosecond pulse. With this approach the distribution of free electrons is described by the integral concentration *n*(*z*, *t*), which varies in time and space (along the axis *z*, directed into the material) due to photo- and thermionic emission from the surface layers and electron diffusion. It is assumed that the thermalization of the electron gas occurs so fast that the notion of electron temperature can be immediately applied. "Hot" electrons contribute to the photo- and thermal emission: they withdraw a part of the energy stored in the electronic subsystem, thus reducing its temperature and, eventually, the temperature of solid as a whole. At the same time, change in the electrons concentration in

the surface region results in a change in the optical characteristics of the material.

explosion. Conditions for Coulomb explosion occurrence in metals were not found.

The Fermi energy *εF* for metals is very high. For example, for copper *εF*=7.1 eV, for silver *εF*=5.5 eV. For this reason heat effects engage the electrons, whose energy lies in a narrow energy range 2 *<sup>B</sup>* ≈ *k T* (*kB -* Boltzmann constant) near the Fermi level.

The concentration of electrons, which absorb the incident radiation in metal, can be estimated by using the relation *n*=*hνN*/*εF*, where *ν* − the frequency of the incident light. For example, *n/N=*25% of the conduction electrons are affected in copper under *hν=*1.7 eV.

The electron transfer part of excessive energy, which they receive due to light absorption, though their collision to other electrons, ions and lattice defects (dislocations, grain boundaries, etc.). Heating of the metal is determined by efficiency of the collisions, which depends on the particles that exchange with energy. Typically the following relation takes place: *νee*>*νei*>*νep*, where *νee* − the frequency of electron - electron collisions, *νei* − the frequency of collisions between electrons and phonons, *νep* − the frequency of electrons collisions with impurities and defects in the metal. The electron-phonon mechanism is quite feasible, since the Fermi energy of electrons in metals is high and is essentially represented by the energy of translational motion of free electrons. The Fermi velocity is *υF*=(2*εF/me*)1/2=1,5·108 cm/s is also high compared to the speed of sound in metals, that is typically: *υ<sup>F</sup>* ~ 105 cm/s.

A single collision may not be enough for the electron, which absorbed a photon, to release the excessive energy, i.e. energy relaxation is a multi-stage process of a diffusion character. For this reason, the energy redistribution occurs not only at the skin layer (*δ<sup>s</sup>* ~ 10-6 cm), but in the deeper layer of *ld*=(*D*/*νee*)1/2=10-5 cm, where *D* is the electron diffusion coefficient.

Thus there is a heating of the metal. The heat, which is released in the layer *lD*, is further transferred depthward into the material through heat conduction. The characteristic time of the absorbed energy transfer in metal is

$$\tau\_{ei} = 1/\nu\_{ei} = 10^{-12} \text{-} 10^{-11} \text{ s.} \quad $$

The electron gas and the lattice of the metal are two weakly interacting subsystems. Under the conditions:

$$t \gg \mathfrak{r}\_{\mathfrak{e}\prime} \tag{1}$$

(*τee* − time to establish an equilibrium energy distribution in the electron gas) and

$$\nu\_{\dot{\imath}i} \ll \nu\_{\dot{\imath}i\prime} \tag{2}$$

(*νii* − frequency of ion-ion collisions) the electron gas and the lattice can be described separately with electron *Te* and lattice *Ti* temperatures. Condition (1) infers rapid redistribution of the absorbed energy between the conduction electrons, and (2) means that the energy transferred to the lattice by electrons, rapidly redistributes between the ions.

Let us excessively consider the relaxation rates *νef*, *νee*, *νei*, *νii.* The relations between these relaxation rates substantially determine the processes in metals during the absorption of radiation.

The frequency of collisions between electrons and photons *νef* is proportional to the power density of laser radiation *q* absorbed by the metal. It can be estimated using the relation:

$$\nu\_{cf} = \frac{\alpha q}{h\nu \, m} \,\, \,\tag{3}$$

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 201

10-14–10-13 s, which is comparable with the pulse duration (tens of femtoseconds). However, further studies [16-19] have shown that the electron-electron relaxation time can be reduced up to 10-16 s at the electron gas temperatures of ~ 100 000 K, achieved by the action of ultrashort laser pulse. In other words, two-temperature model could be applied to analyze

**Multiphoton absorption.** Although the theory of multiphoton absorption is pretty well developed [17], there is a certain difficulty regarding definition of multi-photon absorption cross-sections of real media, when this theory is used for the analysis of multiphoton

In terms of quantum mechanics multiphoton process can be represented as a series of successive transitions of an electron to the virtual states [18]. Only the initial and final states are real in this case. The energy conservation law is valid with an accuracy of a natural width of the energy level only for the initial and final states. For virtual states the energy conservation law takes place with the accuracy, that is determined by the energy-time uncertainty relation *δE·δt*≥*ћ*. At each virtual state the quantum system lives for the time:

> δ

It means that the quantum system can potentially absorb a photon of any energy, however the lifetime of the quantum system, absorbing a photon, would differ. If the quantum system absorbs a photon with energy *hν*=*ΔEmn*, the system transits into a real state, where the lifetime δ*t* is determined by the probability of spontaneous decay of this state. If the quantum system absorbs photon with energy *hν*≠*ΔEmn*, the system transits into a virtual state, where the lifetime is determined by the energy-time uncertainty relation: *δt*≥*ћ*/*δE*,

With this approach one can use the following method for estimation of multiphoton absorption cross sections for various medium. One-photon absorption cross section

metals can be determined by the known absorption coefficient α and the concentration of

photon energy of about 1 eV τ0 ~ 10-16 s and assuming the absorption cross section of excited

0. Similar consideration provides

σ

**2. Effect of electron emission on metal heating and destruction by** 

σ*m+*1*=*σ1*m+1n*τ

2=10-28cm s,

σ

3=10-61cm3s2.

Let us analyze the influence of hot electron emission on heating and destruction of metals based on a two-temperature model and using as example the numerical simulation of

*/n*. Estimating the lifetime of an electron on a virtual level with the

σ3*=*σ13*n*τ

0*<sup>m</sup>*. For example, for

*t E* ≥ . (6)

1, one obtains two-photon absorption cross

α

02. Consequently for *m*+1-

=105cm-1 and *n*=1022cm-3

σ1 for

δ

Absorption of another photon makes the system transit to the next state.

the effects of femtosecond laser pulse.

*δE*=│*hν*-*ΔEmn*│.

free electrons *n*:

section to be

one can obtain

Coulomb explosion.

σ1*=*α

σ2*=*σ12*n*τ

σ

**femtosecond laser pulse** 

photon absorption cross section is

electrons in all virtual levels to be equal to

1=10-17cm2,

σ

absorption during femtosecond laser action.

where *α* is the absorption coefficient in metals: *α* ~ 105 cm-1.

The frequency of electron-electron collisions *νee* in the metal is mainly determined by the number of electrons in the Fermi smearing and is calculated using the relation:

$$\nu\_{\rm ce} = \nu\_{\rm f} \sigma\_{\rm ce} n\_e \left(\frac{k\_B T\_e}{\varepsilon\_{\rm f}}\right)^2,\tag{4}$$

where *υ<sup>F</sup>* − the electron velocity on the Fermi surface ~ 108 cm/s; *σee* − electron - electron interaction cross-section, *σee* ≈ 5⋅10<sup>−</sup>16 cm2; *kBTe* − Fermi smearing region. Then, at *Te* ~ 103 K, the value of *νee* ≈ 1014 s<sup>−</sup>1, and the time to establish an equilibrium distribution of the electron gas *τee* ~ 1/*νee* ~ 10<sup>−</sup>14 s. The rate of energy transfer from the electron gas to the lattice and the lattice temperature is determined by the heat source and heat transfer coefficient of the electrons with the lattice *βei*.

The energy, which the lattice obtains from the electron gas per unit volume per unit time is ~ *βei*(*Te*-*Ti*), *βei* ~ 1010 W/(cm3K).

The frequency of electron-ion relaxation can be expressed in terms of heat transfer coefficient

$$\nu\_{ei} = \frac{\beta\_{ei}}{\rho\_i c\_i} \,\prime \tag{5}$$

where *ρici* − the volumetric heat capacity of the lattice [10<sup>−</sup>7 W⋅s/(cm3K)]. Substituting numerical values one obtains *νei*~1011 s<sup>−</sup>1 and *τei*~10<sup>−</sup>11 s.

Comparing the expressions (3) and (4), one can show that (1) is always satisfied with the flux densities *q*<sup>0</sup> ≤ 109 W/cm2 and the electron gas in a metal is described by temperature *Te*. Similarly, the condition (2) is also satisfied, therefore temperature *Ti* can be introduced to describe the thermal state of the lattice.

Thus, the electron relaxation time is estimated as *τr*~1/*νee*, where *νee* − the frequency of electron-electron collisions. According to the estimates given above, its magnitude is about 10-14–10-13 s, which is comparable with the pulse duration (tens of femtoseconds). However, further studies [16-19] have shown that the electron-electron relaxation time can be reduced up to 10-16 s at the electron gas temperatures of ~ 100 000 K, achieved by the action of ultrashort laser pulse. In other words, two-temperature model could be applied to analyze the effects of femtosecond laser pulse.

200 Laser Pulses – Theory, Technology, and Applications

radiation.

electrons with the lattice *βei*.

*βei*(*Te*-*Ti*), *βei* ~ 1010 W/(cm3K).

(*νii* − frequency of ion-ion collisions) the electron gas and the lattice can be described separately with electron *Te* and lattice *Ti* temperatures. Condition (1) infers rapid redistribution of the absorbed energy between the conduction electrons, and (2) means that the energy transferred to the lattice by electrons, rapidly redistributes between the ions.

Let us excessively consider the relaxation rates *νef*, *νee*, *νei*, *νii.* The relations between these relaxation rates substantially determine the processes in metals during the absorption of

The frequency of collisions between electrons and photons *νef* is proportional to the power density of laser radiation *q* absorbed by the metal. It can be estimated using the relation:

> , *ef q h n* α

The frequency of electron-electron collisions *νee* in the metal is mainly determined by the

<sup>=</sup>

where *υ<sup>F</sup>* − the electron velocity on the Fermi surface ~ 108 cm/s; *σee* − electron - electron interaction cross-section, *σee* ≈ 5⋅10<sup>−</sup>16 cm2; *kBTe* − Fermi smearing region. Then, at *Te* ~ 103 K, the value of *νee* ≈ 1014 s<sup>−</sup>1, and the time to establish an equilibrium distribution of the electron gas *τee* ~ 1/*νee* ~ 10<sup>−</sup>14 s. The rate of energy transfer from the electron gas to the lattice and the lattice temperature is determined by the heat source and heat transfer coefficient of the

The energy, which the lattice obtains from the electron gas per unit volume per unit time is ~

The frequency of electron-ion relaxation can be expressed in terms of heat transfer coefficient

*ei*

ν

, *ei*

*i i c* β

ρ

where *ρici* − the volumetric heat capacity of the lattice [10<sup>−</sup>7 W⋅s/(cm3K)]. Substituting

Comparing the expressions (3) and (4), one can show that (1) is always satisfied with the flux densities *q*<sup>0</sup> ≤ 109 W/cm2 and the electron gas in a metal is described by temperature *Te*. Similarly, the condition (2) is also satisfied, therefore temperature *Ti* can be introduced to

Thus, the electron relaxation time is estimated as *τr*~1/*νee*, where *νee* − the frequency of electron-electron collisions. According to the estimates given above, its magnitude is about

2 , *B e*

*F k T*

ε

 

ν

<sup>=</sup> (3)

<sup>=</sup> (5)

(4)

ν

number of electrons in the Fermi smearing and is calculated using the relation:

ν υσ*n*

*ee F ee e*

where *α* is the absorption coefficient in metals: *α* ~ 105 cm-1.

numerical values one obtains *νei*~1011 s<sup>−</sup>1 and *τei*~10<sup>−</sup>11 s.

describe the thermal state of the lattice.

**Multiphoton absorption.** Although the theory of multiphoton absorption is pretty well developed [17], there is a certain difficulty regarding definition of multi-photon absorption cross-sections of real media, when this theory is used for the analysis of multiphoton absorption during femtosecond laser action.

In terms of quantum mechanics multiphoton process can be represented as a series of successive transitions of an electron to the virtual states [18]. Only the initial and final states are real in this case. The energy conservation law is valid with an accuracy of a natural width of the energy level only for the initial and final states. For virtual states the energy conservation law takes place with the accuracy, that is determined by the energy-time uncertainty relation *δE·δt*≥*ћ*. At each virtual state the quantum system lives for the time:

$$
\delta t \geq \hbar / \delta E.\tag{6}
$$

Absorption of another photon makes the system transit to the next state.

It means that the quantum system can potentially absorb a photon of any energy, however the lifetime of the quantum system, absorbing a photon, would differ. If the quantum system absorbs a photon with energy *hν*=*ΔEmn*, the system transits into a real state, where the lifetime δ*t* is determined by the probability of spontaneous decay of this state. If the quantum system absorbs photon with energy *hν*≠*ΔEmn*, the system transits into a virtual state, where the lifetime is determined by the energy-time uncertainty relation: *δt*≥*ћ*/*δE*, *δE*=│*hν*-*ΔEmn*│.

With this approach one can use the following method for estimation of multiphoton absorption cross sections for various medium. One-photon absorption cross section σ1 for metals can be determined by the known absorption coefficient α and the concentration of free electrons *n*: σ1*=*α*/n*. Estimating the lifetime of an electron on a virtual level with the photon energy of about 1 eV τ0 ~ 10-16 s and assuming the absorption cross section of excited electrons in all virtual levels to be equal to σ1, one obtains two-photon absorption cross section to be σ2*=*σ12*n*τ0. Similar consideration provides σ3*=*σ13*n*τ02. Consequently for *m*+1 photon absorption cross section is σ*m+*1*=*σ1*m+1n*τ0*<sup>m</sup>*. For example, for α=105cm-1 and *n*=1022cm-3 one can obtain σ1=10-17cm2, σ2=10-28cm s, σ3=10-61cm3s2.

## **2. Effect of electron emission on metal heating and destruction by femtosecond laser pulse**

Let us analyze the influence of hot electron emission on heating and destruction of metals based on a two-temperature model and using as example the numerical simulation of Coulomb explosion.

In theoretical studies dealing with the different variants of two-temperature model, the most important aspect is to define the nonlinear optical and thermal properties in a wide temperature range, as well as quantitative characteristics of the electron-electron and electron-phonon interaction, controlling the temperature of the electron gas and the energy exchange between electrons and lattice, respectively.

In terms of the numerical analysis classical version of a two-temperature model is a system of interrelated differential equations of heat conduction, the accuracy of whose solution depends strongly on adaptation of the computational grid to the required numerical solution.

Let us consider the two-temperature model of metal heating by a femtosecond laser pulses on metals [19-22]. The model consists of a system of heat-conduction equations for the electrons and the phonons (lattice) subsystems, where thermophysical properties depend on electron temperature *Te* and electron concentration *n*, and the equation describing the temporal evolution of the electron concentration.

$$\mathcal{L}\_{\varepsilon} \left( T\_{\varepsilon}, n \right) \frac{\partial T\_{\varepsilon}}{\partial t} - \frac{\partial}{\partial z} \left[ \lambda\_{\varepsilon} \left( T\_{\varepsilon}, n \right) \frac{\partial T\_{\varepsilon}}{\partial z} \right] = -\mathcal{B}\_{\varepsilon i} \left( T\_{\varepsilon}, n \right) \left( T\_{\varepsilon} - T\_i \right) + q\_{\text{v}} \tag{7}$$

$$\left[c\_i \frac{\partial T\_i}{\partial t} - \frac{\partial}{\partial z}\right] \lambda\_i \frac{\partial T\_i}{\partial z}\bigg| = \mathcal{J}\_{ei}\left(T\_{e'}n\right)\left(T\_e - T\_i\right) \tag{8}$$

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 203

with initial and boundary conditions:

photoemission for the *m*-photon absorption:

estimation of the critical electric field [24]:

emission *Fe*=*Ft*+*Fmpho*.

0

*z <sup>n</sup> D F z* <sup>=</sup>

The thermionic emission is determined by the law of Richardson:

2 0

=

0

0

*J* is the absorbed photons flow, *J=q/hν*, *σm* is the multi-photon absorption cross section.

The emission of electrons leads to accumulation of positive charge on the metal surface and, therefore, to generation of the electric field. The electric field resulted from breaking of quasi-neutrality of the irradiated area can be calculated from the following equation [24]:

> 0 . εε *e i <sup>E</sup> <sup>q</sup> n n*

This electric field, which is induced by the charge separation, can reach extremely high magnitude and exceed the energy of atomic bonds resulting in Coulomb explosion. To determine the conditions for initiation of the Coulomb explosion the electric field (10) is compared to the threshold magnitude required for removal of an atom from the target. The

0

<sup>2</sup> , *th z <sup>n</sup> <sup>E</sup>*

*z*

( )

0

εε <sup>=</sup> <sup>Λ</sup> <sup>=</sup>

<sup>∂</sup> = − ∂ (10)

∞ <sup>=</sup> =− −

*mpho z <sup>m</sup> <sup>e</sup> F* σ

, *e*

<sup>∂</sup> = − <sup>∂</sup> <sup>0</sup> <sup>0</sup> , *<sup>z</sup> n n* <sup>=</sup> <sup>=</sup> <sup>0</sup> , *<sup>z</sup> n n* =∞ <sup>=</sup>

where *n*0=1022 cm-3 is the initial concentration of electrons evenly distributed over the volume, *Fe* is the flow of the electrons resulted from thermionic emission and photo-

> *t ez e e bez F BT zl q k T* ϕ

External photoelectric effect implies that the energy of absorbed photons is used by the electrons to overcome the work function, i.e. the minimum energy required for the electron to escape from the surface. The work function in metals is several eV (for silver *φe*=4,28 eV). Therefore, for *hν*=1.55 eV one should expect the three-photon absorption. For calculating the photoemission let us assume that the free electrons, which are involved into the multiphoton process, reach the surface without energy loss and leave the metal. The emissive layer thickness is limited by the electron mean free path and also depends on the electron concentration and the emission coefficient. The flow of electrons (cm-2/s) caused by

= − <sup>−</sup> <sup>−</sup>

here *B* is the Richardson coefficient, *φ<sup>e</sup>* is the work function, *qe* is the electron charge.

0 exp exp( / ) / , *<sup>e</sup>*

exp( / ) , *<sup>m</sup>*

*J z l dz*

=

where *λe*, *λi*, *c*e, *ci* are the electron and lattice heat conductivity and heat capacity, β*ei* – electron-ion energy transfer coefficient, *qv* – the absorbed power density released in the electron subsystem.

Dependence of optical and thermo-physical properties of metal [23] on the electron temperature and concentration is taken into account for electron heat capacity: 2 2 ( ,) , <sup>2</sup> *b e e F knztT c* π ε <sup>=</sup> electron heat conductivity: *<sup>λ</sup>e=*ν*e* 2τ*eece*/3, electron-ion energy transfer coefficient: *βei=ce/tei*, electron velocity: 3 /, *e be* ν = *kT m* the electron mean free path: *l*e*=*1*/*(*nσ√*2). The absorption coefficient of metals depends only on the concentration of free electrons: α(*n*)*=*(α*/n*0)*n*, *n*0 is the initial concentration of free electrons.

The boundary conditions:

$$\left. \lambda\_{\epsilon} \frac{\partial T\_{\epsilon}}{\partial z} \right|\_{z=0} = -j\_{\epsilon \prime} \quad \lambda\_i \left. \frac{\partial T\_i}{\partial z} \right|\_{z=0} = 0, \quad T\_{\epsilon} \Big|\_{z=\infty} = T\_i \Big|\_{z=\infty} = T\_{n\prime}$$

where *je* – is the heat flow carried away by the emitted electrons, *Tn* – the initial temperature. Evolution of the electron density distribution is described by the following diffusion equation:

$$\frac{\partial n(z,t)}{\partial t} = D \frac{\partial^2 n(z,t)}{\partial z^2} \,, \tag{9}$$

with initial and boundary conditions:

202 Laser Pulses – Theory, Technology, and Applications

exchange between electrons and lattice, respectively.

temporal evolution of the electron concentration.

<sup>=</sup> electron heat conductivity: *<sup>λ</sup>e=*

coefficient: *βei=ce/tei*, electron velocity: 3 /, *e be*

0 λ , *<sup>e</sup> e e z <sup>T</sup> <sup>j</sup> z* <sup>=</sup>

<sup>∂</sup> = − <sup>∂</sup> <sup>0</sup>

*i*

electron subsystem.

2 2 ( ,) , <sup>2</sup> *b e*

*F knztT*

ε

α(*n*)*=*(α

The boundary conditions:

*e*

π

electrons:

*c*

In theoretical studies dealing with the different variants of two-temperature model, the most important aspect is to define the nonlinear optical and thermal properties in a wide temperature range, as well as quantitative characteristics of the electron-electron and electron-phonon interaction, controlling the temperature of the electron gas and the energy

In terms of the numerical analysis classical version of a two-temperature model is a system of interrelated differential equations of heat conduction, the accuracy of whose solution depends strongly on adaptation of the computational grid to the required numerical solution.

Let us consider the two-temperature model of metal heating by a femtosecond laser pulses on metals [19-22]. The model consists of a system of heat-conduction equations for the electrons and the phonons (lattice) subsystems, where thermophysical properties depend on electron temperature *Te* and electron concentration *n*, and the equation describing the

> ( ) ( ) ( )( ) <sup>ν</sup> , λ , ,, *e e e e e e ei e e i*

∂ ∂ <sup>∂</sup> <sup>−</sup> = − − +

λ ( )( ) , , *i i i i ei e e i*

β∂ ∂ <sup>∂</sup> −= − ∂∂ ∂

β

ν*e* 2τ

<sup>∂</sup> <sup>=</sup> <sup>∂</sup> , *ein z z TTT* =∞ =∞ = =

(7)

β*ei* –

(8)

*eece*/3, electron-ion energy transfer

= *kT m* the electron mean free path:

∂ ∂ <sup>=</sup> <sup>∂</sup> <sup>∂</sup> (9)

*c Tn T n TnT T q tz z*

*c TnT T*

where *λe*, *λi*, *c*e, *ci* are the electron and lattice heat conductivity and heat capacity,

ν

*/n*0)*n*, *n*0 is the initial concentration of free electrons.

λ 0, *<sup>i</sup>*

*T z* <sup>=</sup>

*z*

where *je* – is the heat flow carried away by the emitted electrons, *Tn* – the initial temperature. Evolution of the electron density distribution is described by the following diffusion equation:

> (,) (,) , *nzt nzt <sup>D</sup> t z*

2

2

electron-ion energy transfer coefficient, *qv* – the absorbed power density released in the

Dependence of optical and thermo-physical properties of metal [23] on the electron temperature and concentration is taken into account for electron heat capacity:

*l*e*=*1*/*(*nσ√*2). The absorption coefficient of metals depends only on the concentration of free

*T T*

∂∂ ∂

*T T*

*tz z*

$$D\frac{\partial \mathfrak{n}}{\partial z}\bigg|\_{z=0} = -F\_{\mathfrak{e}'} \cdot \mathfrak{n}\Big|\_{z=0} = \mathfrak{n}\_{0'} \cdot \mathfrak{n}\Big|\_{z=\infty} = \mathfrak{n}\_{0'} \cdot \mathfrak{n}$$

where *n*0=1022 cm-3 is the initial concentration of electrons evenly distributed over the volume, *Fe* is the flow of the electrons resulted from thermionic emission and photoemission *Fe*=*Ft*+*Fmpho*.

The thermionic emission is determined by the law of Richardson:

$$F\_t = -BT\_e^2 \Big|\_{z=0} \exp\left(-\frac{\varphi\_e}{k\_b T\_e \Big|\_{z=0}}\right) \exp(-z \Big/ l\_e) / q\_{e'} $$

here *B* is the Richardson coefficient, *φ<sup>e</sup>* is the work function, *qe* is the electron charge.

External photoelectric effect implies that the energy of absorbed photons is used by the electrons to overcome the work function, i.e. the minimum energy required for the electron to escape from the surface. The work function in metals is several eV (for silver *φe*=4,28 eV). Therefore, for *hν*=1.55 eV one should expect the three-photon absorption. For calculating the photoemission let us assume that the free electrons, which are involved into the multiphoton process, reach the surface without energy loss and leave the metal. The emissive layer thickness is limited by the electron mean free path and also depends on the electron concentration and the emission coefficient. The flow of electrons (cm-2/s) caused by photoemission for the *m*-photon absorption:

$$\left.F\_{mpho}\right|\_{z=0} = -\bigwedge\_{0}^{\infty} \sigma\_m f^m \exp(-z \left.l\_e\right|) dz, \star$$

*J* is the absorbed photons flow, *J=q/hν*, *σm* is the multi-photon absorption cross section.

The emission of electrons leads to accumulation of positive charge on the metal surface and, therefore, to generation of the electric field. The electric field resulted from breaking of quasi-neutrality of the irradiated area can be calculated from the following equation [24]:

$$\frac{\partial E}{\partial z} = \frac{q\_e}{\varepsilon \varepsilon\_0} (n\_i - n). \tag{10}$$

This electric field, which is induced by the charge separation, can reach extremely high magnitude and exceed the energy of atomic bonds resulting in Coulomb explosion. To determine the conditions for initiation of the Coulomb explosion the electric field (10) is compared to the threshold magnitude required for removal of an atom from the target. The estimation of the critical electric field [24]:

$$E\_{th}|\_{z=0} = \sqrt{\frac{2\Lambda n}{\varepsilon \varepsilon\_0}},$$

where *n*0- the concentration of atoms (cm-3), Λ = 2.951 J/atom - the heat of sublimation, *ε* = 4.9 - relative permeability of silver, *ε*0=8.854·10-14 F/cm - dielectric constant.

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 205

*ph* (3) and thermionic emission *ne*

*ph* (3) and thermionic emission *ne*

*th* (4). Laser pulse

*th* (4). Laser pulse

**Figure 2.** Laser pulse shape (1), the temporal evolution of the net surface electron density *N* (2) and the

**Figure 3.** Laser pulse shape (1), the temporal evolution of the net surface electron density *N* (2) and the

The gradient of the electron density results in an electric field, which grows and reaches the

Fig. 5 illustrates the nonlinear dependence of the moment when Coulomb explosion starts on the laser power density. It is seen that the dependence is nonlinear, and Coulomb

We can draw the following conclusions from the results of numerical simulation of the influence of electron emission on heating and destruction of metals irradiated by

explosion for a pulse duration ~ 100 fs can occur when *q >* 1015 W/cm2.

density of emitted electrons due to photoemission *ne*

density of emitted electrons due to photoemission *ne*

shape *q=qm*(*t/tm*)exp(*-t/tm*), *tm*=100 fs

femtosecond laser pulse.

Coulomb explosion threshold (Fig. 4).

shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs

For the numerical solution of the heat conduction equation and the equation describing the electron density, which are non-stationary partial differential equations, the finite difference method was applied. To calculate the values of the temperature and the electron concentration explicit difference scheme was used. Though explicit scheme provides a relatively high speed of calculation, it has a serious disadvantage related to the need for satisfying the stability conditions, which impose limits on the amount of steps partitions with respect to coordinate and time.

Fig. 1 - 4 show the results for silver for the laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs. The calculation was performed until the beginning of the Coulomb explosion, when the electric field resulted from the charge separation exceeds the threshold required for the removal of atoms. For comparison, the simulation was performed without taking the emission into account, but preserving the dependence of material properties on temperature.

Fig. 1 shows that temperature of the electrons and the lattice increase, but does not reach its maximum during the pulse. The difference between the temperature, calculated taking into account the emission and without it at the beginning of Coulomb explosion *t*=0.15 fs is Δ*T*≈290 K for temperature of the electrons (Fig. 1b) and Δ*T* ≈60 K for temperature of the lattice.

**Figure 1.** a – Transient behavior of electron and lattice temperature: (1) and (2) – with taking into account emission, (3) and (4) – disregarding emission. (5) – laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs, maximum power density *qm*=5·1014 W/cm2. b – an enlarged part of the temporal dependence of the electron temperature: (1) – with taking into account emission, (2) – disregarding emission

Fig. 2, 3 illustrate temporal behavior of the free electrons concentration, and the number of electrons emitted by photo-and thermionic emission for different pulse shapes. At the initial stage of pulse action photoemission dominates, but thermionic emission increases with the electron temperature rapidly, so thermionic emission begins to prevail over the photoemission. The maximum of the photoemission rate maximum corresponds to the maximum of laser power density. This indicates that the laser pulse shape significantly influences the dynamics of the processes.

with respect to coordinate and time.

influences the dynamics of the processes.

where *n*0- the concentration of atoms (cm-3), Λ = 2.951 J/atom - the heat of sublimation, *ε* = 4.9

For the numerical solution of the heat conduction equation and the equation describing the electron density, which are non-stationary partial differential equations, the finite difference method was applied. To calculate the values of the temperature and the electron concentration explicit difference scheme was used. Though explicit scheme provides a relatively high speed of calculation, it has a serious disadvantage related to the need for satisfying the stability conditions, which impose limits on the amount of steps partitions

Fig. 1 - 4 show the results for silver for the laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs. The calculation was performed until the beginning of the Coulomb explosion, when the electric field resulted from the charge separation exceeds the threshold required for the removal of atoms. For comparison, the simulation was performed without taking the emission

Fig. 1 shows that temperature of the electrons and the lattice increase, but does not reach its maximum during the pulse. The difference between the temperature, calculated taking into account the emission and without it at the beginning of Coulomb explosion *t*=0.15 fs is Δ*T*≈290 K for temperature of the electrons (Fig. 1b) and Δ*T* ≈60 K for temperature of the lattice.

a) b)

**Figure 1.** a – Transient behavior of electron and lattice temperature: (1) and (2) – with taking into account emission, (3) and (4) – disregarding emission. (5) – laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs, maximum power density *qm*=5·1014 W/cm2. b – an enlarged part of the temporal dependence of the electron temperature: (1) – with taking into account emission, (2) – disregarding emission

Fig. 2, 3 illustrate temporal behavior of the free electrons concentration, and the number of electrons emitted by photo-and thermionic emission for different pulse shapes. At the initial stage of pulse action photoemission dominates, but thermionic emission increases with the electron temperature rapidly, so thermionic emission begins to prevail over the photoemission. The maximum of the photoemission rate maximum corresponds to the maximum of laser power density. This indicates that the laser pulse shape significantly

into account, but preserving the dependence of material properties on temperature.


**Figure 2.** Laser pulse shape (1), the temporal evolution of the net surface electron density *N* (2) and the density of emitted electrons due to photoemission *ne ph* (3) and thermionic emission *ne th* (4). Laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs

**Figure 3.** Laser pulse shape (1), the temporal evolution of the net surface electron density *N* (2) and the density of emitted electrons due to photoemission *ne ph* (3) and thermionic emission *ne th* (4). Laser pulse shape *q=qm*(*t/tm*)exp(*-t/tm*), *tm*=100 fs

The gradient of the electron density results in an electric field, which grows and reaches the Coulomb explosion threshold (Fig. 4).

Fig. 5 illustrates the nonlinear dependence of the moment when Coulomb explosion starts on the laser power density. It is seen that the dependence is nonlinear, and Coulomb explosion for a pulse duration ~ 100 fs can occur when *q >* 1015 W/cm2.

We can draw the following conclusions from the results of numerical simulation of the influence of electron emission on heating and destruction of metals irradiated by femtosecond laser pulse.

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 207

Emission processes have a significant impact on the processes of heating and destruction of the semiconductors, because the initial concentration of conductivity electrons in semiconductors can be below concentration of free electrons generated by the action of femtosecond laser radiation in contrast to metals. Let us consider the emission impact on the

The formation of periodic surface structures (PSS) is a perfect evidence of an induced change in the surface optical properties under the femtosecond pulse action on semiconductor and wide-band dielectrics (Fig. 6). Being formed under different conditions PSS exhibit the same formation regularities: the structures orientation depends on the direction of polarization vector of the laser radiation; the structures period depends on the wavelength, incidence angle of the radiation and dielectric permeability of the medium. The observed regularities suggest that the mechanism of PPS formation is determined by electromagnetic field, which is the result of interference of the incident wave and the excited

**Figure 6.** SEM image of the monocrystalline silicon surface irradiated by the laser pulse. Arrow indicates the direction of laser radiation polarization. Left: after irradiation by 1200 pulses, the energy

Microstructures shown in Fig. 6, were obtained on silicon irradiated by 80 fs laser pulses at a wavelength of 1.25 μm [25, 26]. Surface structures turned during rotation of the laser radiation polarization vector, and depending on the density of laser flux were oriented either perpendicular (*Q*0 ~ 1 J/cm2) or parallel to the direction of polarization (*Q*<sup>0</sup> ~

It should be noted that low concentration of free electrons in semiconductor in the initial state does not provide the surface optical properties, which are necessary for excitation of

density is ~ 1 J/cm2. Right: after exposure by 300 pulses, the energy density is ~ 2 J/cm2

example of a femtosecond microstructuring of the silicon surfaces.

surface electromagnetic waves (SEW).

2 J/cm2).

**3. Effect of electron emission on changes in optical properties of** 

**semiconductors under the femtosecond laser pulse action** 

**Figure 4.** Temporal dependence of the electric field *E* (1), arising due to emission of electrons and the threshold field value *E*th (2) corresponding to the beginning of Coulomb explosion. Laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs (3) and the temporal evolution of the net surface electron density (4)

**Figure 5.** Time of Coulomb explosion onset as function of laser density for laser pulse shape *q=qm*exp(- (*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs

Various types of emission have different impacts on the concentration of emitted electrons (see Fig. 2-3): thermionic emission dominates over the photoemission and increases with the electron temperature increasing. The pulse shape significantly affects on the dynamics of the emission processes. However, according to the calculations the effect of emission processes on the electron gas temperature and the lattice temperature (Fig. 1a, b) is negligible. Also, numerical experiment showed that the occurrence of Coulomb explosion (caused by the emission processes during the pulse) in metals requires high-power incident radiation, which is impossible in the real exposure modes.

Emission processes have a significant impact on the processes of heating and destruction of the semiconductors, because the initial concentration of conductivity electrons in semiconductors can be below concentration of free electrons generated by the action of femtosecond laser radiation in contrast to metals. Let us consider the emission impact on the example of a femtosecond microstructuring of the silicon surfaces.

## **3. Effect of electron emission on changes in optical properties of semiconductors under the femtosecond laser pulse action**

206 Laser Pulses – Theory, Technology, and Applications

density (4)

(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs

which is impossible in the real exposure modes.

**Figure 4.** Temporal dependence of the electric field *E* (1), arising due to emission of electrons and the threshold field value *E*th (2) corresponding to the beginning of Coulomb explosion. Laser pulse shape *q=qm*exp(-(*t-tm*)2/*tm*12), *tm*=100 fs, *tm*1=50 fs (3) and the temporal evolution of the net surface electron

**Figure 5.** Time of Coulomb explosion onset as function of laser density for laser pulse shape *q=qm*exp(-

Various types of emission have different impacts on the concentration of emitted electrons (see Fig. 2-3): thermionic emission dominates over the photoemission and increases with the electron temperature increasing. The pulse shape significantly affects on the dynamics of the emission processes. However, according to the calculations the effect of emission processes on the electron gas temperature and the lattice temperature (Fig. 1a, b) is negligible. Also, numerical experiment showed that the occurrence of Coulomb explosion (caused by the emission processes during the pulse) in metals requires high-power incident radiation, The formation of periodic surface structures (PSS) is a perfect evidence of an induced change in the surface optical properties under the femtosecond pulse action on semiconductor and wide-band dielectrics (Fig. 6). Being formed under different conditions PSS exhibit the same formation regularities: the structures orientation depends on the direction of polarization vector of the laser radiation; the structures period depends on the wavelength, incidence angle of the radiation and dielectric permeability of the medium. The observed regularities suggest that the mechanism of PPS formation is determined by electromagnetic field, which is the result of interference of the incident wave and the excited surface electromagnetic waves (SEW).

**Figure 6.** SEM image of the monocrystalline silicon surface irradiated by the laser pulse. Arrow indicates the direction of laser radiation polarization. Left: after irradiation by 1200 pulses, the energy density is ~ 1 J/cm2. Right: after exposure by 300 pulses, the energy density is ~ 2 J/cm2

Microstructures shown in Fig. 6, were obtained on silicon irradiated by 80 fs laser pulses at a wavelength of 1.25 μm [25, 26]. Surface structures turned during rotation of the laser radiation polarization vector, and depending on the density of laser flux were oriented either perpendicular (*Q*0 ~ 1 J/cm2) or parallel to the direction of polarization (*Q*<sup>0</sup> ~ 2 J/cm2).

It should be noted that low concentration of free electrons in semiconductor in the initial state does not provide the surface optical properties, which are necessary for excitation of the SEW. In case of longer laser pulses (longer than tens of picoseconds) experimentally observed excitation of SEW on semiconductors is related to properties of the melt formed on the surface due to laser heating. This explanation can't be used for femtosecond action because photoexcitation and thermal processes are separated in time and the surface does not melt during the laser pulse. The conditions for excitation of SEW during ultrashort laser pulse result from high concentration of nonequilibrium carriers, which are generated in the semiconductor by the light.

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 209

*e*

<sup>∂</sup> <sup>∂</sup> (11)

∂ (12)

<sup>=</sup> =− − (13)

(14)

(15)

τ

2

2 2 ( ,) ( ,) ( ,) , *phi*

*nzt nzt nzt J D <sup>t</sup> <sup>z</sup>*

=+ −

where *D* is the diffusion coefficient of electrons in a solid, *τe* is the time of electron

Similarly to the above considerations for metals, the losses of electrons at the surface caused

*nzt D FF*

− =+

∞

20 2 0 exp( / ) . *pho z <sup>e</sup> F* σ

02 0 0 (,) , *z pho z t z*

*J z l dz*

0 exp exp( / ) / *<sup>e</sup>*

/ ,

/ / ,

<sup>=</sup> (16)

α ν

=

*t ez e e bez F BT zl q k T* ϕ

Temperature of electrons gas (*Te*) is determined by heat conduction equations (15), where

( )

*T T a T T c Jh c <sup>t</sup> <sup>z</sup>*

*e ei e i e e e*

2 2 ( ,) , <sup>2</sup> *b e*

*F knztT*

ε

The Bouguer–Lambert differential law of Eq. (17) determines the intensity distribution *J*(*z, t*) inside the solid (the *z* axis is directed depth ward). In this model absorptance of the material

( )

*i ei e i i*

= −− +

*a TT c*

π

β

β

= − <sup>−</sup> <sup>−</sup>

*z* = ==

2

Expression (13) describes two-photon emission. Single-photon absorption cross-section *σ*1 is estimated assuming the maximum absorption coefficient, that corresponds to absorption by free carriers in metal *σ*<sup>1</sup> ≈ *αmax*/*nmax* ~ 10–17 cm2 (*nmax* is electron concentration in the metal). Photo-ionization of atoms is considered in terms of electron transition through virtual states. The electron life time (*τ*0) at these virtual states is determined basing on the uncertainty relation between energy and time: *τ*0~10-16 s. In this case two-photon cross-section is given by *σ*2≈*σ*12*n*(*z*, *t*)*τ*0 ~ 10–28 cm s. Expression (14) describes external thermo-emission

α

by external emission are taken into account in the model as follows:

∂

2 0

=

heat capacity of electron gas is a function of temperature and concentration.

∂ ∂ = −− <sup>∂</sup> <sup>∂</sup>

*e*

where *F2pho* is determined from (13) for two-photon photoeffect, *Ft* - from (14).

*c*

*e e*

∂ ∂

∂ ∂

*i i*

*T T*

*t z*

collisions.

(Richardson low).

(*A*) is assumed to be constant.

∂ ∂

Dynamics of the optical properties change at the surface of semiconductors under femtosecond laser action is related to change of the non-equilibrium carriers plasma frequency. In order to analyze the behavior of the optical properties let us first consider the basic mechanisms of light absorption and recombination of the absorbed energy.

Total absorption coefficient in semiconductors can be considered as a sum of absorption coefficients associated with different mechanisms 1 *i i* α = α = . First of all these are fundamental

band-to-band absorption of light (*hν*>*E*g, *E*g is the band-gap) and intraband absorption, i.e. absorption by free carriers - electrons and holes (from now on we consider only the electrons for simplicity). Rate of relaxation of the crystal electron system from the excited state to the equilibrium state is determined by the recombination mechanisms with characteristic times ~ 10-12-10-10s.

During femtosecond laser pulse the light intensity can achieve very high level without destruction of the matter, therefore initiating multiphoton processes. In the wavelength range from IR to near UV the energy of a single photon is not sufficient for electron transition from the valence band to the conduction band (*hν*<*E*g). Such transition takes place as a result of simultaneous absorption of several photons. In this case the rate (or probability) of multiphoton ionization is highly dependent on the laser power. The multiphoton ionization rate is proportional to *σI*m, where *I* is the laser radiation intensity, *σ<sup>m</sup>* - *m*-photon absorption cross-section. The required number of photons is determined by the lowest value of *m*, satisfying the relation *mhν*>*E*g.

During femtosecond laser pulse only photoexcitation and fast electronic processes are observed, while recombination and lattice heating can be neglected, because the characteristic times of these processes are much higher than the pulse duration. This is one of the major differences between action of ultrashort and longer laser pulses.

**Model.** According to the model the dynamics of distribution of nonequilibrium electron concentration *n*(*z*, *t*) during the femtosecond pulse is determined by generation of the nonequilibrium electrons due to two-photon absorption (*hν*<*E*g) followed by their diffusion and participation in collision processes as described in the diffusion equation (11).

In this consideration nonequilibrium electrons are the electrons that transit from the valence band to the conduction band under photo-excitation, and then contribute to increasing electron gas temperature, emission, and finally recombination after the end of the pulse. Below by "electrons" we will understand "non-equilibrium electrons" and by "electron gas" we will understand "nonequilibrium electrons gas".

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 209

$$\frac{\partial n(z,t)}{\partial t} = \alpha\_{2\text{-}phi}I + D\frac{\partial^2 n(z,t)}{\partial z^2} - \frac{n(z,t)}{\tau\_\ell},\tag{11}$$

where *D* is the diffusion coefficient of electrons in a solid, *τe* is the time of electron collisions.

208 Laser Pulses – Theory, Technology, and Applications

coefficients associated with different mechanisms

lowest value of *m*, satisfying the relation *mhν*>*E*g.

we will understand "nonequilibrium electrons gas".

semiconductor by the light.

~ 10-12-10-10s.

the SEW. In case of longer laser pulses (longer than tens of picoseconds) experimentally observed excitation of SEW on semiconductors is related to properties of the melt formed on the surface due to laser heating. This explanation can't be used for femtosecond action because photoexcitation and thermal processes are separated in time and the surface does not melt during the laser pulse. The conditions for excitation of SEW during ultrashort laser pulse result from high concentration of nonequilibrium carriers, which are generated in the

Dynamics of the optical properties change at the surface of semiconductors under femtosecond laser action is related to change of the non-equilibrium carriers plasma frequency. In order to analyze the behavior of the optical properties let us first consider the

Total absorption coefficient in semiconductors can be considered as a sum of absorption

band-to-band absorption of light (*hν*>*E*g, *E*g is the band-gap) and intraband absorption, i.e. absorption by free carriers - electrons and holes (from now on we consider only the electrons for simplicity). Rate of relaxation of the crystal electron system from the excited state to the equilibrium state is determined by the recombination mechanisms with characteristic times

During femtosecond laser pulse the light intensity can achieve very high level without destruction of the matter, therefore initiating multiphoton processes. In the wavelength range from IR to near UV the energy of a single photon is not sufficient for electron transition from the valence band to the conduction band (*hν*<*E*g). Such transition takes place as a result of simultaneous absorption of several photons. In this case the rate (or probability) of multiphoton ionization is highly dependent on the laser power. The multiphoton ionization rate is proportional to *σI*m, where *I* is the laser radiation intensity, *σ<sup>m</sup>* - *m*-photon absorption cross-section. The required number of photons is determined by the

During femtosecond laser pulse only photoexcitation and fast electronic processes are observed, while recombination and lattice heating can be neglected, because the characteristic times of these processes are much higher than the pulse duration. This is one

**Model.** According to the model the dynamics of distribution of nonequilibrium electron concentration *n*(*z*, *t*) during the femtosecond pulse is determined by generation of the nonequilibrium electrons due to two-photon absorption (*hν*<*E*g) followed by their diffusion

In this consideration nonequilibrium electrons are the electrons that transit from the valence band to the conduction band under photo-excitation, and then contribute to increasing electron gas temperature, emission, and finally recombination after the end of the pulse. Below by "electrons" we will understand "non-equilibrium electrons" and by "electron gas"

of the major differences between action of ultrashort and longer laser pulses.

and participation in collision processes as described in the diffusion equation (11).

1 *i i* α=

=

α

. First of all these are fundamental

basic mechanisms of light absorption and recombination of the absorbed energy.

Similarly to the above considerations for metals, the losses of electrons at the surface caused by external emission are taken into account in the model as follows:

$$\left. \frac{\partial n(z,t)}{\partial z} \right|\_{z=0} = F\_{2\,pho} \Big|\_{z=0} + F\_t \Big|\_{z=0} \tag{12}$$

$$\left.F\_{2\,ph\nu}\right|\_{z=0} = -\int\_0^\alpha \sigma\_2 f^2 \exp(-z \left.\left.l\_e\right|\_e) dz.\tag{13}$$

Expression (13) describes two-photon emission. Single-photon absorption cross-section *σ*1 is estimated assuming the maximum absorption coefficient, that corresponds to absorption by free carriers in metal *σ*<sup>1</sup> ≈ *αmax*/*nmax* ~ 10–17 cm2 (*nmax* is electron concentration in the metal). Photo-ionization of atoms is considered in terms of electron transition through virtual states. The electron life time (*τ*0) at these virtual states is determined basing on the uncertainty relation between energy and time: *τ*0~10-16 s. In this case two-photon cross-section is given by *σ*2≈*σ*12*n*(*z*, *t*)*τ*0 ~ 10–28 cm s. Expression (14) describes external thermo-emission (Richardson low).

$$F\_t = -BT\_e^2 \Big|\_{z=0} \exp\left(-\frac{\varphi\_e}{k\_b T\_e \Big|\_{z=0}}\right) \exp(-z \Big/ l\_e) / q\_e \tag{14}$$

Temperature of electrons gas (*Te*) is determined by heat conduction equations (15), where heat capacity of electron gas is a function of temperature and concentration.

$$\begin{aligned} \frac{\partial T\_e}{\partial t} &= a\_e \frac{\partial^2 T\_e}{\partial z^2} - \beta\_{ei} \left(T\_e - T\_i\right) / c\_e + \alpha\_e l l \nu / c\_{e\prime} \\ \frac{\partial T\_i}{\partial t} &= a\_i \frac{\partial^2 T\_i}{\partial z^2} - \beta\_{ei} \left(T\_e - T\_i\right) / c\_{i\prime} \end{aligned} \tag{15}$$

$$c\_e = \frac{\pi^2 k\_b^2 n(z, t) T\_e}{2\varepsilon\_F} \tag{16}$$

where *F2pho* is determined from (13) for two-photon photoeffect, *Ft* - from (14).

The Bouguer–Lambert differential law of Eq. (17) determines the intensity distribution *J*(*z, t*) inside the solid (the *z* axis is directed depth ward). In this model absorptance of the material (*A*) is assumed to be constant.

$$\begin{aligned} \frac{\partial f(z,t)}{\partial z} &= - (\alpha\_{2phi} + \alpha\_c + \alpha\_{2pho}) f(z,t), \\ f(0,t) &= A f\_0(t); f\_0(t) = q\_0(t) / \langle h\nu \rangle, \\ q\_0(t) &= q\_m(t \,/\ t\_m) \exp(-t \,/\ t\_m) \end{aligned} \tag{17}$$

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 211

= (21)

(20)

( )

( ) <sup>2</sup> 4 , , *<sup>e</sup>*

*e n zt q m*

+

 γ

*p*

ω γ

ω

*p*

γ

ω

,

Re ,

= − <sup>+</sup>

*n*

 ε*i*

= +′ ′′

εε

ε ε

ε

*p*

ω

imaginary parts of the permeability.

**Figure 7.** The dynamics of silicon photoexcitation

above. The calculation results are shown on Fig. 8-10.

=

Im ,

ωω

π

where *εn* is the initial value of the semiconductor permittivity, and *ε*΄ and *ε*΄΄ is the real and

Let us consider dynamics of the dielectric permeability in the surface layer of semiconductor taking into account the change in the plasma frequency of non-equilibrium carriers by using a mathematical model (11-21). This will help us to identify the role of different emission

Numerical simulation was performed for silicon (*εn* = 12, *γ* = 1014 s-1) with initial data given

When there is no emission surface the real part of dielectric permeability at the surface quickly decreases and becomes negative stabilizing by the end of the first quarter of the pulse (Fig. 10, curve 1). If external photo-emission is taken into account, the character of the dynamics of Reε does not change, but the magnitude at which the permeability stabilizes is increased (Fig. 8, curve 2). Thermionic emission strongly affects the dynamics of the

processes in evolution of the optical properties during a femtosecond pulse.

where *AJ*0(*t*) is the density of the absorbed photon flux, with a bell-shaped temporal distribution of intensity of the laser radiation, *α*2*phi* is two-photon absorption coefficient of the inner photoeffect, *α*2*pho* is two-photon absorption coefficient of the extrinsic photoeffect, and *αe* is coefficient for absorption by free electrons, which in turn are defined as

$$\begin{aligned} \alpha\_{2\text{phi}} &= \sigma\_1^2 n\_{\text{max}} \mathfrak{r}\_\nu \mathfrak{l}\_\nu\\ \alpha\_{2\text{pho}} &= \sigma\_1^2 n(z, t) \mathfrak{r}\_\nu \mathfrak{l}\_\nu\\ \alpha\_\varepsilon &= \sigma\_1 n(z, t) .\end{aligned} \tag{18}$$

Expressions (11-18) along with initial and boundary conditions at *t* = 0 and *z* = ∞

$$\begin{aligned} \left. m \right|\_{z=\omega} &= n \Big|\_{t=\omega} = 0, \\ \left. T\_e \right|\_{z=\omega} &= T\_i \Big|\_{z=\omega} = T\_e \Big|\_{t=0} = T\_i \Big|\_{t=0} = 0 \end{aligned} \tag{19}$$

allow one to obtain spatial and temporal distribution of the electron concentration in semiconductor *n*(*z*, *t*).

The first qualitative estimates were made for a simplified model, in which the emission mechanism is generalized and the emission flux in the expression (12) is taken into account using the emission factor *μ* without separation of photo- and thermal emission

$$-D\frac{\partial n(z,t)}{\partial z}\Big|\_{z=0} = -\mu m.$$

Fig. 7 shows obtained variation of the electron concentration with depth during the laser pulse action on silicon. The following initial data were used for the calculation: *D =* 80 cm2/s*, τе =*10-14s*, μ =* 2·108 cm/*s, Q*0 *=* 2 J/cm2, pulse duration is 80 fs, wavelength is 1.25 μm. The calculation results indicate that the semiconductor surface acquires properties of a metal during the laser pulse action. The maximum of electron density is located at some distance from the surface. It shifts from the surface into the bulk of material and its value increases during the pulse action. By the middle of the pulse the distribution of electron concentration stabilizes following the shape of the laser pulse.

The above estimates showed that dynamics of the optical properties of semiconductor under action of ultrashort pulse can be described within the same approach used for metals. According to Drude dispersion theory the dielectric permeability of photoexcited semiconductor can be determined by the plasma frequency of the electron gas (*ωp*), incident radiation frequency (*ω*) and the frequency of electron collisions (*γ*), according to following expressions:

$$\begin{aligned} \varepsilon &= \varepsilon' + i\varepsilon'', \\\\ \text{Re}\,\varepsilon &= \varepsilon\_n - \frac{o\rho\_p^2}{o^2 + \gamma^2}, \\\\ \text{Im}\,\varepsilon &= \frac{o\rho\_p^2 \gamma}{o\left(o^2 + \gamma^2\right)}, \\\\ o\rho\_p &= \sqrt{\frac{4\pi n(z,t)q\_e^2}{m\_e}} \end{aligned} \tag{21}$$

where *εn* is the initial value of the semiconductor permittivity, and *ε*΄ and *ε*΄΄ is the real and imaginary parts of the permeability.

**Figure 7.** The dynamics of silicon photoexcitation

210 Laser Pulses – Theory, Technology, and Applications

semiconductor *n*(*z*, *t*).

expressions:

stabilizes following the shape of the laser pulse.

0

∂

*z*

∂

2 2

*phi e pho*

ν

(17)

(18)

(,) ( ) ( , ),

 αα

*Jzt <sup>J</sup> z t*

00 0

*mm m*

where *AJ*0(*t*) is the density of the absorbed photon flux, with a bell-shaped temporal distribution of intensity of the laser radiation, *α*2*phi* is two-photon absorption coefficient of the inner photoeffect, *α*2*pho* is two-photon absorption coefficient of the extrinsic photoeffect,

> 2 2 1 max 2 2 1 1

αστ

= = =

( , ).

*nzt*

0,

=∞ =∞ = =

allow one to obtain spatial and temporal distribution of the electron concentration in

The first qualitative estimates were made for a simplified model, in which the emission mechanism is generalized and the emission flux in the expression (12) is taken into account

<sup>∂</sup> <sup>−</sup> = − <sup>∂</sup>

Fig. 7 shows obtained variation of the electron concentration with depth during the laser pulse action on silicon. The following initial data were used for the calculation: *D =* 80 cm2/s*, τе =*10-14s*, μ =* 2·108 cm/*s, Q*0 *=* 2 J/cm2, pulse duration is 80 fs, wavelength is 1.25 μm. The calculation results indicate that the semiconductor surface acquires properties of a metal during the laser pulse action. The maximum of electron density is located at some distance from the surface. It shifts from the surface into the bulk of material and its value increases during the pulse action. By the middle of the pulse the distribution of electron concentration

The above estimates showed that dynamics of the optical properties of semiconductor under action of ultrashort pulse can be described within the same approach used for metals. According to Drude dispersion theory the dielectric permeability of photoexcited semiconductor can be determined by the plasma frequency of the electron gas (*ωp*), incident radiation frequency (*ω*) and the frequency of electron collisions (*γ*), according to following

0 (,) . *<sup>z</sup> nzt D n z*

μ=

*e i ei z z tt*

*T T TT*

,

ν

0 0

0

= === (19)

( ,) ,

 τ

ν

*n J nzt J*

(0, ) ( ); ( ) ( ) / ( ), ( ) ( / )exp( / ),

*J t AJ t J t q t h qt q tt tt*

= = = −

=− + +

α

and *αe* is coefficient for absorption by free electrons, which in turn are defined as

*phi pho e*

Expressions (11-18) along with initial and boundary conditions at *t* = 0 and *z* = ∞

= =

using the emission factor *μ* without separation of photo- and thermal emission

α σ

*z t*

=∞ =∞

*n n*

ασ

Let us consider dynamics of the dielectric permeability in the surface layer of semiconductor taking into account the change in the plasma frequency of non-equilibrium carriers by using a mathematical model (11-21). This will help us to identify the role of different emission processes in evolution of the optical properties during a femtosecond pulse.

Numerical simulation was performed for silicon (*εn* = 12, *γ* = 1014 s-1) with initial data given above. The calculation results are shown on Fig. 8-10.

When there is no emission surface the real part of dielectric permeability at the surface quickly decreases and becomes negative stabilizing by the end of the first quarter of the pulse (Fig. 10, curve 1). If external photo-emission is taken into account, the character of the dynamics of Reε does not change, but the magnitude at which the permeability stabilizes is increased (Fig. 8, curve 2). Thermionic emission strongly affects the dynamics of the

permeability. Influence of the thermionic emission is small in the beginning of the pulse and value of Re*ε* abruptly decreases. Few femtosecond later contribution of the thermionic emission grows and value of Reε, returned to its initial level after several damping oscillations (Fig. 8, curve 3). The observed inertia is typical for the thermionic emission mechanism, and the return of the dielectric permeability to the initial value means that all the "hot" electrons leave the surface as a result of thermal emission. If photo-emission is also taken into account escape of the electrons speeds up (Fig. 8, curve 4).

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 213

**Figure 9.** Dynamics of the real part of the permeability on the surface of silicon for various values of the light-flux energy density. (a) Photoemission only, (b) Combined action of photoemission and thermal

Figure 10 shows the calculated distribution of permeability Re*ε*(*z*), at the end of the pulse (*Q*<sup>0</sup> *=* 2 J/cm2). If there is no emission, a very thin metal-like layer appears on the surface (curve 1). If external photo-emission is taken into account, a metal-like layer (several tens of nanometers thick, curve 2) is formed. The loss of the electrons due to thermionic emission qualitatively changes the permeability distribution depth wards (curve 3). A layer with smaller permeability than its initial value is formed at a distance of tens of nanometers from the surface. Combined action of photo-emission and thermionic emission shifts this layer

Let us consider again the conditions necessary for the SEW excitation in order to compare the results of numerical simulation with the above experimental data on femtosecond silicon

It is known that the formation of PSS oriented perpendicular to the polarization of laser radiation usually results from excitation of surface plasmon-polaritons under the laser pulse action on the metal. Surface plasmon-polaritons are partially longitudinal electromagnetic waves of the TM-type propagating along the interface between two media with the wave electro-magnetic field being localized near the interface. Excitation of plasmon-polaritons is

emission. 1 – *Q*0= 1 J/cm2, 2 – *Q*0 = 2 J/cm2, 3 – *Q*0 = 3 J/cm2

microstructuring associated with the SEW excitation.

deeper (curve 4).

**Figure 8.** Dynamics of the real part of the permeability at the surface of silicon (*Q*0 *=* 2 J/cm2). 1 – no external emission, 2 – only photo-emission is taken into account, 3 – only thermionic emission is taken into account, 4 – combined action of photo-emission and thermal emission

Figure 9 shows dynamics of the real part of the permeability at the silicon surface for different values of the light flux density. If only two-photon photo-emission is taken into account, the surface acquires metal-like properties at the first femtoseconds of the pulse, and the permeability remains negative during the entire pulse. Increasing the radiation-flux density reduces the time to reach the steady-state level without changing the character of the dependence (Fig. 9a). When both photo-emission and thermionic emission are taken into account the picture changes (Fig. 9b). If *Q*<sup>0</sup> ≤ 1 J/cm2, Re*ε* smoothly returns to its initial value (curve 1). When the energy density increases the permeability oscillates during transition to its initial value (curves 2 and 3). These oscillations result from dependence of heat capacity of the electron gas on electron concentration. The higher light flux induces higher temperature of the electron gas and contribution of thermionic emission increases. As a result the permeability returns to its initial value faster. At the same time increasing flux density raises the electron concentration and, accordingly, the heat capacity of the electron gas increases thus reducing its temperature and the contribution of the thermal emission. This results in oscillating of dielectric permeability.

permeability. Influence of the thermionic emission is small in the beginning of the pulse and value of Re*ε* abruptly decreases. Few femtosecond later contribution of the thermionic emission grows and value of Reε, returned to its initial level after several damping oscillations (Fig. 8, curve 3). The observed inertia is typical for the thermionic emission mechanism, and the return of the dielectric permeability to the initial value means that all the "hot" electrons leave the surface as a result of thermal emission. If photo-emission is also

**Figure 8.** Dynamics of the real part of the permeability at the surface of silicon (*Q*0 *=* 2 J/cm2). 1 – no external emission, 2 – only photo-emission is taken into account, 3 – only thermionic emission is taken

Figure 9 shows dynamics of the real part of the permeability at the silicon surface for different values of the light flux density. If only two-photon photo-emission is taken into account, the surface acquires metal-like properties at the first femtoseconds of the pulse, and the permeability remains negative during the entire pulse. Increasing the radiation-flux density reduces the time to reach the steady-state level without changing the character of the dependence (Fig. 9a). When both photo-emission and thermionic emission are taken into account the picture changes (Fig. 9b). If *Q*<sup>0</sup> ≤ 1 J/cm2, Re*ε* smoothly returns to its initial value (curve 1). When the energy density increases the permeability oscillates during transition to its initial value (curves 2 and 3). These oscillations result from dependence of heat capacity of the electron gas on electron concentration. The higher light flux induces higher temperature of the electron gas and contribution of thermionic emission increases. As a result the permeability returns to its initial value faster. At the same time increasing flux density raises the electron concentration and, accordingly, the heat capacity of the electron gas increases thus reducing its temperature and the contribution of the thermal emission.

into account, 4 – combined action of photo-emission and thermal emission

This results in oscillating of dielectric permeability.

taken into account escape of the electrons speeds up (Fig. 8, curve 4).

**Figure 9.** Dynamics of the real part of the permeability on the surface of silicon for various values of the light-flux energy density. (a) Photoemission only, (b) Combined action of photoemission and thermal emission. 1 – *Q*0= 1 J/cm2, 2 – *Q*0 = 2 J/cm2, 3 – *Q*0 = 3 J/cm2

Figure 10 shows the calculated distribution of permeability Re*ε*(*z*), at the end of the pulse (*Q*<sup>0</sup> *=* 2 J/cm2). If there is no emission, a very thin metal-like layer appears on the surface (curve 1). If external photo-emission is taken into account, a metal-like layer (several tens of nanometers thick, curve 2) is formed. The loss of the electrons due to thermionic emission qualitatively changes the permeability distribution depth wards (curve 3). A layer with smaller permeability than its initial value is formed at a distance of tens of nanometers from the surface. Combined action of photo-emission and thermionic emission shifts this layer deeper (curve 4).

Let us consider again the conditions necessary for the SEW excitation in order to compare the results of numerical simulation with the above experimental data on femtosecond silicon microstructuring associated with the SEW excitation.

It is known that the formation of PSS oriented perpendicular to the polarization of laser radiation usually results from excitation of surface plasmon-polaritons under the laser pulse action on the metal. Surface plasmon-polaritons are partially longitudinal electromagnetic waves of the TM-type propagating along the interface between two media with the wave electro-magnetic field being localized near the interface. Excitation of plasmon-polaritons is

possible only if one of the media has positive dielectric permeability (*ε*1>0), while the real part of dielectric permeability of the other is negative (Re*ε*2<0). Also, the condition |Re*ε*2|>*ε*<sup>1</sup> must be satisfied. A negative dielectric permeability in a metal is determined by a high concentration of free electrons. In the case of the relatively long laser pulse action (tens of picoseconds or more) on the semiconductor the appearance of metal-like optical response is usually associated with the properties of the melt formed on the surface due to laser heating. Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 215

The refractive index of the waveguide layer in such optically layered structure (*n*2) exceeds

**Figure 11.** Schematic distribution of the electrons concentration in a semiconductor under the femtosecond laser pulse action (a) and the corresponding distribution of the dielectric permeability (b) for different emissions; *n*1, *n*2, *n*3, *n*4 - effective refractive indices of the layers as they are alternating

In addition, it is necessary to provide a certain minimum thickness of the waveguide layer at a given frequency. If action of laser radiation results in formation of a layer with a refractive index *n*3 and closer to surface a layer of thickness *h*<sup>1</sup> with a refractive index *n*2 provided

min 2 2 2 2

In particular, for silicon at a wavelength of 1.25 μm the minimum thickness of the

<sup>−</sup> <sup>≈</sup>

*n n <sup>h</sup>*

λ

2

waveguide layer is *h*min ~ 70 nm according to expression (22).

π

2 1 2 3

*n n n n*

arccos .

2 2 2 3

<sup>−</sup> <sup>−</sup> (22)

along the coordinate *z* 

*n*2>*n*3>*n*1, then [7]

the refractive index values of adjacent layers (*n*1*, n*3), (*n*2*>n*1*, n*2*>n*3).

**Figure 10.** Spatial distribution (into the depthward of the semiconductor) of the real part of the permeability at the end of the laser pulse. 1 – no external emission of electrons, 2 – photo-emission only, 3 – thermionic emission only, 4 – combined action of photo-emission and thermal emission

Semiconductors typically have positive value of their dielectric permeability in the visible and IR range. Under the femtosecond laser pulse action the concentration of nonequilibrium carriers in semiconductors can become so high that dielectric permeability ε would change its sign and creates conditions for excitation of surface plasmon-polaritons. In this case, the formation of PSS perpendicular to the polarization vector is experimentally observed.

The formation of structures parallel to the polarization vector is associated with the excitation of surface waveguide modes (TE-polaritons). It is necessary to create an optically layered structure for excitation of a waveguide mode at the semiconductor surface (Fig. 11). The refractive index of the waveguide layer in such optically layered structure (*n*2) exceeds the refractive index values of adjacent layers (*n*1*, n*3), (*n*2*>n*1*, n*2*>n*3).

214 Laser Pulses – Theory, Technology, and Applications

possible only if one of the media has positive dielectric permeability (*ε*1>0), while the real part of dielectric permeability of the other is negative (Re*ε*2<0). Also, the condition |Re*ε*2|>*ε*<sup>1</sup> must be satisfied. A negative dielectric permeability in a metal is determined by a high concentration of free electrons. In the case of the relatively long laser pulse action (tens of picoseconds or more) on the semiconductor the appearance of metal-like optical response is usually associated with the properties of the melt formed on the surface due to laser heating.

**Figure 10.** Spatial distribution (into the depthward of the semiconductor) of the real part of the permeability at the end of the laser pulse. 1 – no external emission of electrons, 2 – photo-emission only,

formation of PSS perpendicular to the polarization vector is experimentally observed.

Semiconductors typically have positive value of their dielectric permeability in the visible and IR range. Under the femtosecond laser pulse action the concentration of nonequilibrium carriers in semiconductors can become so high that dielectric permeability ε would change its sign and creates conditions for excitation of surface plasmon-polaritons. In this case, the

The formation of structures parallel to the polarization vector is associated with the excitation of surface waveguide modes (TE-polaritons). It is necessary to create an optically layered structure for excitation of a waveguide mode at the semiconductor surface (Fig. 11).

3 – thermionic emission only, 4 – combined action of photo-emission and thermal emission

**Figure 11.** Schematic distribution of the electrons concentration in a semiconductor under the femtosecond laser pulse action (a) and the corresponding distribution of the dielectric permeability (b) for different emissions; *n*1, *n*2, *n*3, *n*4 - effective refractive indices of the layers as they are alternating along the coordinate *z* 

In addition, it is necessary to provide a certain minimum thickness of the waveguide layer at a given frequency. If action of laser radiation results in formation of a layer with a refractive index *n*3 and closer to surface a layer of thickness *h*<sup>1</sup> with a refractive index *n*2 provided *n*2>*n*3>*n*1, then [7]

$$h\_{\min} = \frac{\mathcal{A}}{2\pi\sqrt{n\_2^2 - n\_3^2}} \arccos\sqrt{\frac{n\_2^2 - n\_3^2}{n\_2^2 - n\_1^2}}.\tag{22}$$

In particular, for silicon at a wavelength of 1.25 μm the minimum thickness of the waveguide layer is *h*min ~ 70 nm according to expression (22).

Let us consider the calculated spatial distribution of the silicon dielectric permeability under the femtosecond laser pulse action (see Fig. 10) in terms of the possible conversion of the incident light into surface plasmon-polaritons and waveguide modes. The surface polariton excitation requires transition of the semiconductor surface into a metal-like state, while formation of a dynamic optically-layered structure with a certain minimum thickness of the waveguide layer is necessary for the excitation and propagation of waveguide modes.

Interaction of Femtosecond Laser Pulses with Solids: Electron/Phonon/Plasmon Dynamics 217

σ

3=10-61cm3s2 forthree-photon absorption.

1=10-17cm2 for one-photon

The proposed method allowed to estimate the cross sections of multiphoton absorption in

Roman V. Dyukin, George A. Martsinovskiy, Olga N. Sergaeva, Galina D. Shandybina,

This work was supported by grants of RFBR 09-02-00932-а, 09-02-01065-а and State Contract

[1] Korte F., Nolte S., Chichkov B.N., Bauer T., Kamlage G., Wagner T., Fallnich C., Welling H. (1999) Far-Field and Near-Field Material Processing with Femtosecond Laser Pulses.

[2] Vorobyev D.Y., Guo C.L., (2010) Metallic Light Absorbers Produced by Femtosecond

[3] Chimmalgi A., Grigoropoulos C.P., Komvopoulos K. (2005) Surface Nanostructuring by Nano-Femtosecond Laser Assistant Force Microscopy. App. Phys. 97: 1043191 -

[4] Myers R.A., Farrell R., Karger A.M., Carey J.E., Mazur E. (2006) Enhancing Near-Infrared Avalanche Photodiode Performance by Femtosecond Laser Microstructuring.

[5] Etsion I. (2005) State of the Art in Laser Surface Texturing. J. of Tribology, 127(1): 248-

[6] Oktem B., Kalaycioglu H., Erdoğan M., Yavaş S., Mukhopadhyay P., Tazebay U. H., Aykaç Y., Eken K., Ilday F. Ö. (2010) Surface Texturing of Dental Implant Surfaces with an Ultrafast Fiber Laser. in Conference on Lasers and Electro-Optics, OSA Technical

[7] Libenson M.N. (2007) [Laser-Induced Optical and Thermal Processes in Solids and

[8] Libenson M.N. (2001) Non-Equilibrium Heating and Cooling of Metals Under Action of

[9] Carey J.E., Crouch C.H., Mazur E. (2003) Femtosecond-Laser-Assisted Microstructuring

[10] Shimosuma Y., Kazansky P.G., Qin J.R., Hirao K. (2003) Self-Organized Nanogratings in Glass Irradiated by Ultrashort Light Pulses. Phys. Rev. Lett. 91: 247405-247409.

Their Mutual Influence]. Saint-Petersburg: Nauka. 423 p. (in Russian)

*National Research University of Information Technologies, Mechanics and Optics,* 

Laser Pulses. Advances in Mechanical Engineering, 2010: 1-5.

Digest (CD) (Optical Society of America, 2010), JTuD15.

of Silicon Surfaces. Optics&Photonics News. 14(2): 32-36.

Super-Short Laser Pulse. Proc. SPIE. 4423: 1–7.

σ

metals. For example, for metals absorption cross section:

2=10 -28cm s fortwo-photon absorption,

*Department of Laser-Assisted Technologies and Applied Ecology,* 

absorption,

No P1134.

**5. References** 

104319112.

253.

σ

*Saint-Petersburg, Russia* 

**Acknowledgement** 

Vera V. Svirina and Eugeny B. Yakovlev

Appl. Phys. A, 69 [Suppl.]: S7–S11.

Appl Opt. 45(35): 8825-31.

**Author details** 

As follows from the numerical model (Fig. 10, curves 1-2) in case of relatively low emission a metal-like layer is formed at the surface. Within the layer thickness is about 50-60 nm dielectric permeability becomes negative. This provides the conditions for excitation of surface plasmon-polaritons (TM-type SEW), which is confirmed experimentally by formation of the microstructures perpendicularly to the polarization vector.

If emission rate is high as in case of thermo-emission (Fig. 10, curve 3) an optically layered structure is formed. Although the dielectric permeability does not change its sign, excitation of waveguide modes is possible.

Combination of both photo- and thermal-emission (Fig. 10, curve 4) results in formation of a dielectric layer of thickness ~ 60 nm and a metal-like layer of thickness about 40 nm. The refractive index of this dielectric layer is higher than both the refractive index of the air on one interface, and the refractive index of the metal-like layer on the other interface. Presence of such an optical structure enables excitation of waveguide mode, which results in the formation of periodic relief, which is parallel to the polarization vector on the incident radiation (see Fig. 6, right).

The above consideration showed that multiphoton emission and thermionic emission noticeably vary the optical properties of a semiconductor during a femtosecond pulse action. In particular, layer with different optical properties is formed at the surface and enables excitation of either surface polaritons or waveguide modes in the semiconductors. The considered model allows one to qualitatively and quantitatively interpret available experimental data. This approach allows one to use experimentally observed surface microstructures as relatively simple means of investigation of dynamics of the semiconductor surface properties under femtosecond action.

## **4. Conclusion**

The results of numerical simulation have shown that influence of emission processes on the electron gas temperature and lattice temperature of metals is negligible. Therefore, the emission can be neglected when assessing the parameters of metals processing by femtosecond laser pulse, which simplifies the numerical calculations.

The Coulomb explosion occurrence in metals requires high-power incident radiation, which is impossible for the real exposure conditions.

However, in semiconductors both types of extrinsic emission noticeably change distribution of dielectric permeability near surface providing conditions for excitation of surface polaritons or waveguide modes depending on laser power magnitude.

The proposed method allowed to estimate the cross sections of multiphoton absorption in metals. For example, for metals absorption cross section: σ1=10-17cm2 for one-photon absorption, σ2=10 -28cm s fortwo-photon absorption, σ3=10-61cm3s2 forthree-photon absorption.

## **Author details**

216 Laser Pulses – Theory, Technology, and Applications

of waveguide modes is possible.

radiation (see Fig. 6, right).

**4. Conclusion** 

Let us consider the calculated spatial distribution of the silicon dielectric permeability under the femtosecond laser pulse action (see Fig. 10) in terms of the possible conversion of the incident light into surface plasmon-polaritons and waveguide modes. The surface polariton excitation requires transition of the semiconductor surface into a metal-like state, while formation of a dynamic optically-layered structure with a certain minimum thickness of the waveguide layer is necessary for the excitation and propagation of waveguide modes.

As follows from the numerical model (Fig. 10, curves 1-2) in case of relatively low emission a metal-like layer is formed at the surface. Within the layer thickness is about 50-60 nm dielectric permeability becomes negative. This provides the conditions for excitation of surface plasmon-polaritons (TM-type SEW), which is confirmed experimentally by

If emission rate is high as in case of thermo-emission (Fig. 10, curve 3) an optically layered structure is formed. Although the dielectric permeability does not change its sign, excitation

Combination of both photo- and thermal-emission (Fig. 10, curve 4) results in formation of a dielectric layer of thickness ~ 60 nm and a metal-like layer of thickness about 40 nm. The refractive index of this dielectric layer is higher than both the refractive index of the air on one interface, and the refractive index of the metal-like layer on the other interface. Presence of such an optical structure enables excitation of waveguide mode, which results in the formation of periodic relief, which is parallel to the polarization vector on the incident

The above consideration showed that multiphoton emission and thermionic emission noticeably vary the optical properties of a semiconductor during a femtosecond pulse action. In particular, layer with different optical properties is formed at the surface and enables excitation of either surface polaritons or waveguide modes in the semiconductors. The considered model allows one to qualitatively and quantitatively interpret available experimental data. This approach allows one to use experimentally observed surface microstructures as relatively simple means of investigation of dynamics of the

The results of numerical simulation have shown that influence of emission processes on the electron gas temperature and lattice temperature of metals is negligible. Therefore, the emission can be neglected when assessing the parameters of metals processing by

The Coulomb explosion occurrence in metals requires high-power incident radiation, which

However, in semiconductors both types of extrinsic emission noticeably change distribution of dielectric permeability near surface providing conditions for excitation of surface

formation of the microstructures perpendicularly to the polarization vector.

semiconductor surface properties under femtosecond action.

is impossible for the real exposure conditions.

femtosecond laser pulse, which simplifies the numerical calculations.

polaritons or waveguide modes depending on laser power magnitude.

Roman V. Dyukin, George A. Martsinovskiy, Olga N. Sergaeva, Galina D. Shandybina, Vera V. Svirina and Eugeny B. Yakovlev *Department of Laser-Assisted Technologies and Applied Ecology, National Research University of Information Technologies, Mechanics and Optics, Saint-Petersburg, Russia* 

## **Acknowledgement**

This work was supported by grants of RFBR 09-02-00932-а, 09-02-01065-а and State Contract No P1134.

## **5. References**

	- [11] Kudryashov S.I., Emel'yanov V.I. (2001) Electron gas compression and Coulomb explosion in the surface layer of a conductor heated by femtosecond laser pulse. JETP Lett. 73(12): 666-670.

**Chapter 8** 

© 2012 Mazhukin, licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2012 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,

**Kinetics and Dynamics of Phase** 

Additional information is available at the end of the chapter

V.I. Mazhukin

**1. Introduction** 

[8], etc.

equilibrium.

http://dx.doi.org/10.5772/50731

**Transformations in Metals Under Action** 

Action of super-power and ultra-short laser pulses on highly absorbing condensed media is investigated in the past two decades [1] - [5]. The urgency of this problem is primarily determined by a variety of practical applications of pulsed laser irradiation. In a short period of time scientists have mastered such operations as ablation of elemental materials by femtosecond lasers (100 fs) [6], femtosecond nanostructuring [7], generation of metallic nanoparticles and nanostructures by laser ablation of massive targets by 40 fs and 2ps pulses

The use of super-power G~1012÷1015 W/cm2 and ultra-short τL≈10-12÷10-15 s laser pulses for dimensional processing of materials, such as cutting of various materials by picofemtosecond laser pulses [9], micro-drilling by a femtosecond laser [10 ], surface etching of metals (Al, Cu, Mo, Ni) and semiconductors (Si) [11], is accompanied by realization of unique physical conditions. In particular, the duration of action becomes comparable with the characteristic times of thermalization and phase transitions in matter. This leads to the need to address complex fundamental problems, including the heating of the material and the kinetics of phase transitions in a strong deviation from local thermodynamic

It should be noted that the majority of laser technologies is associated with the beginning of phase transformations in the material. In particular, the action of pico - femtosecond laser pulses of high intensity on solid targets is one of the ways to create individual particles with unique characteristics, or to form their streams, consisting of a cluster, liquid or solid fragments of the target. The formation of the particle flux in the pulsed laser ablation is observed for a wide range of materials: metals, semiconductors and insulators. The

and reproduction in any medium, provided the original work is properly cited.

**of Ultra-Short High-Power Laser Pulses** 


**Chapter 8** 
