## 1. Introduction

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.

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.

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 ultrashort-pulse 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*≈10^{22} 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.

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 *k*_{B}
*-* 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}*/m*_{e})^{1/2}=1,5 10^{8} cm/s is also high compared to the speed of sound in metals, that is typically: *υ*_{F} ~ 10^{5} 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 (*δ*_{s} ~ 10^{-6} cm), but in the deeper layer of *l*_{d}=(*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 *l*_{D}, is further transferred depthward into the material through heat conduction. The characteristic time of the absorbed energy transfer in metal is *τ*_{ei} = 1/*ν*_{ei} = 10^{-12}-10^{-11} s.

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

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

(*ν*_{ii} frequency of ion-ion collisions) the electron gas and the lattice can be described separately with electron *T*_{e} and lattice *T*_{i} 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:

where *α* is the absorption coefficient in metals: *α* ~ 10^{5} 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:

where *υ*_{F} the electron velocity on the Fermi surface ~ 10^{8} cm/s; *σ*_{ee} electron - electron interaction cross-section, *σ*_{ee} ≈ 510^{16} cm^{2}; *k*_{B}*T*_{e} Fermi smearing region. Then, at *T*_{e} ~ 10^{3} K, the value of *ν*_{ee} 10^{14} s^{1}, and the time to establish an equilibrium distribution of the electron gas *τ*_{ee} ~ 1/*ν*_{ee} ~ 10^{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}(*T*_{e}-*T*_{i}), *β*_{ei} ~ 10^{10} W/(cm^{3}K).

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

where *ρ*_{i}*c*_{i} the volumetric heat capacity of the lattice [10^{7} Ws/(cm^{3}K)]. Substituting numerical values one obtains *ν*_{ei}~10^{11} s^{1} and *τ*_{ei}~10^{11} s.

Comparing the expressions (3) and (4), one can show that (1) is always satisfied with the flux densities *q*_{0} ≤ 10^{9} W/cm^{2} and the electron gas in a metal is described by temperature *T*_{e}. Similarly, the condition (2) is also satisfied, therefore temperature *T*_{i} 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.

**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:

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ν*=*ΔE*_{mn}, 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ν*≠*ΔE*_{mn}, the system transits into a virtual state, where the lifetime is determined by the energy-time uncertainty relation: *δt*≥*ћ*/*δE*, *δE*=│*hν*-*ΔE*_{mn}│.

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}=_{1}^{2}n_{0}. Similar consideration provides _{3}=_{1}^{3}n_{0}^{2}. Consequently for m+1-photon absorption cross section is _{m+1}=_{1}^{m+1}n_{0}^{m}. For example, for =10^{5}cm^{-1} and n=10^{22}cm^{-3} one can obtain _{1}=10^{-17}cm^{2}, _{2}=10^{-28}cm s, _{3}=10^{-61}cm^{3}s^{2}.

## 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 T_{e} and electron concentration n, and the equation describing the temporal evolution of the electron concentration.

where λ_{e}, λ_{i}, c_{e}, c_{i} are the electron and lattice heat conductivity and heat capacity, β_{ei} – electron-ion energy transfer coefficient, q_{v} – 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: _{e}=_{e}^{2}_{ee}c_{e}/3, electron-ion energy transfer coefficient: β_{ei}=c_{e}/t_{ei}, electron velocity: _{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:

where j_{e} – is the heat flow carried away by the emitted electrons, T_{n} – the initial temperature. Evolution of the electron density distribution is described by the following diffusion equation:

with initial and boundary conditions:

where n_{0}=10^{22} cm^{-3} is the initial concentration of electrons evenly distributed over the volume, F_{e} is the flow of the electrons resulted from thermionic emission and photo-emission F_{e}=F_{t}+F_{mpho}.

The thermionic emission is determined by the law of Richardson:

here B is the Richardson coefficient, φ_{e} is the work function, q_{e} 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:

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]:

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]:

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.

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=q_{m}exp(-(t-t_{m})^{2}/t_{m1}^{2}), t_{m}=100 fs, t_{m1}=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.

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.

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 > 10^{15} W/cm^{2}.

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.

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

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).

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/cm^{2}) or parallel to the direction of polarization (Q_{0} ~ 2 J/cm^{2}).

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.

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_{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^{-10}s.

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, σ_{m} - 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".

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

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:

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 σ_{1} ≈ α_{max}/n_{max} ~ 10^{–17} cm^{2} (n_{max} 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}≈σ_{1}^{2}n(z, t)τ_{0} ~ 10^{–28} cm s. Expression (14) describes external thermo-emission (Richardson low).

Temperature of electrons gas (T_{e}) is determined by heat conduction equations (15), where heat capacity of electron gas is a function of temperature and concentration.

where F_{2pho} is determined from (13) for two-photon photoeffect, F_{t} - 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.

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

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

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

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 cm^{2}/s, τ_{е} =10^{-14}s, μ = 2 10^{8} cm/s, Q_{0} = 2 J/cm^{2}, 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:

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

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, γ = 10^{14} 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).

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_{0} ≤ 1 J/cm^{2}, 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.

Figure 10 shows the calculated distribution of permeability Reε(z), at the end of the pulse (Q_{0} = 2 J/cm^{2}). 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}|>ε_{1} 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.

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}).

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_{1} with a refractive index n_{2} provided n_{2}>n_{3}>n_{1}, then [7]

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.

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^{-17}cm^{2} for one-photon absorption, _{2}=10^{ -28}cm s for two-photon absorption, _{3}=10^{-61}cm^{3}s^{2} for three-photon absorption.