Advances in Condensed Matter Physics

Volume 2013 (2013), Article ID 528208, 11 pages

http://dx.doi.org/10.1155/2013/528208

## Size Effects on Surface Elastic Waves in a Semi-Infinite Medium with Atomic Defect Generation

Institute on Laser and Information Technologies, Russian Academy of Sciences, Moscow 140700, Russia

Received 3 May 2013; Accepted 10 October 2013

Academic Editor: Michael C. Tringides

Copyright © 2013 F. Mirzade. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The paper investigates small-scale effects on the Rayleigh-type surface wave propagation in an isotopic elastic half-space upon laser irradiation. Based on Eringen’s theory of nonlocal continuum mechanics, the basic equations of wave motion and laser-induced atomic defect dynamics are derived. Dispersion equation that governs the Rayleigh surface waves in the considered medium is derived and analyzed. Explicit expressions for phase velocity and attenuation (amplification) coefficients which characterize surface waves are obtained. It is shown that if the generation rate is above the critical value, due to concentration-elastic instability, nanometer sized ordered concentration-strain structures on the surface or volume of solids arise. The spatial scale of these structures is proportional to the characteristic length of defect-atom interaction and increases with the increase of the temperature of the medium. The critical value of the pump parameter is directly proportional to recombination rate and inversely proportional to deformational potentials of defects.

#### 1. Introduction

During the last decades elastic wave propagation in solids (metals and semiconductors) with nonequilibrium atomic defects has received a lot of attention [1–4]. Various types of lattice defects (vacancies and interstitials), produced from the lattice site atoms due to intense external energy fluxes (laser and corpuscular radiations), introduce a significant strain of the medium as a result of the difference between the radii of lattice atoms and defects [5] and play an important role in surface modification of solids exposed to laser radiation [1, 5, 6]. The formation of atomic defects may occur also in a number of other technologies processes: in the laser fast recrystallization, laser annealing, multipulse laser etching, pulsed laser-assisted thin-film deposition, and so forth. Strains in an elastic wave cause a strain-induced drift of defects, whereas the strains and a variation in the temperature in the wave modulate the rates of generation and recombination of defects of the thermal-fluctuation origin (via variations in the energies of the defect formation and migration) [5, 6].

Several mathematical models have been considered to study the self-organization of two-dimensional (2D) ordered microscale concentration-strain (CS) structures (at high concentrations of atomic defects) on the surface of the solid half-space [7, 8] and in an isotropic solid elastic layer [9, 10] under the action of laser irradiation. In these studies, attention was also focused on the study of an influence of the strain-induced diffusion, generation, and annihilation of defects on the propagation characteristics of strain waves. In both [1–4] and [7–10] the analysis was based on coupled evolution equations for the atomic defect density fields and the classical local elasticity equations for the elastic displacement of the medium.

Classical continuum theory is characterized by the local character of stress (stress at a point depends only on the strain at that point) and does not contain internal length scale. The absence of the length scale creates several discrepancies in the predictions of mechanical responses, for example, nondispersive wave behavior (constant wave velocity, independent of frequency). For example, according to the classical elasticity, Rayleigh waves propagating on the surface of a semi-infinite isotropic elastic space are not dispersive at any frequency, whereas experiments and the atomic theory of lattice predict otherwise.

Limitations of the classical elasticity theory are also demonstrated in the study of the formation of coupled strain-defect nanometer sized ordered structures (short-range CS structures) on the surface of laser-irradiated solids and mechanical behavior of microstructured materials, because their behavior is characterized by nonlocal stresses and the existence of an internal length scale [11, 12]. The accurate analysis of dynamic behaviors of these structures cannot be correctly described by classical local elasticity theory. Conducting experiments with nanoscale size materials is found to be difficult and expensive. Therefore, development of appropriate nonlocal elasticity mathematical models for nanostructure formation and nanomaterials is of great importance.

In contrast to local approach of zero-range internal interactions, the nonlocal elasticity theory, originated and developed in the last four decades, postulates that the stress at a point in a body depends not only on the strain at point but also on those at all other points of the body. Various nonlocal theories of linear elasticity have been proposed to describe the scale effects on the characteristics of the vibration and elastic wave propagation in the above-mentioned submicro- or nanosized structures. The basic constitutive equations and governing equations of linear nonlocal elasticity were derived by Eringen [11] and Kunin [12]. Later the nonlocal theories have been applied for the analysis of micro- and nanoscale plate-like structures, in which the small-scale effects become significant. In [13] Eringen considered vibration behavior of a nanoplate by using linear theory of nonlocal continuum mechanics. By using [11], Reddy [14] reformulated the classical and shear deformation beam and plate theories. The nonlocal scale influences on the wave dispersion properties of the nanoplates are discussed in detail in [15]. Vibration characteristics of nanoplates, based on three-dimensional theory of elasticity employing nonlocal continuum mechanics, have been discussed in [15]. A review of some other applications of nonlocal elasticity theories for nanostructures can be found in [16]. Propagation of longitudinal elastic and thermoelastic waves in an isotropic, homogeneous infinite medium with long-rang interactions has been studied by Nowinski [17, 18]. Other advances have been made by the application of nonlocal elasticity to such fields as the solid defects [19, 20] and fracture mechanics [21].

The present paper is concerned with the nonlocal elasticity theory of laser-excited solid half-space with mobile atomic point defects. We summarise the theory formulated in [8] to analyze the effects of nonlocal atom-atom and atom-defect interactions on the surface wave propagation in solids with defect generation. The changes in defect concentration are controlled by the (i) generation of defects by laser irradiation, (ii) their diffusion, and (iii) strain-induced diffusion. Dispersion equation that governs the propagation of elastic-concentration waves has been derived by solving a system of coupled partial differential equations. Some limiting cases of the dispersion equations are considered. For a small value of coupled parameter obvious expressions for the phase velocity and the attenuation (amplification) factor of surface elastic-concentration waves are determined. It is found that both phase velocity and attenuation coefficient are to be influenced by the presence of defect generation. We also obtain that at certain conditions concentration-elastic instabilities with the formation of regular nanosized structures in a system of atomic defects on the surface of the solids can be developed. We demonstrate that due to the nonlocal character of the defect-atom and atom-atom interactions the dispersive curve of the instability has two maxima. As a result, size distribution functions of the surface nanoscale nonhomogeneities having two maxima can be formed. The results of some earlier works are also deduced from the present formulation. To our knowledge, the problem in its present form was not investigated before.

#### 2. Governing Equations

In this section we give basic 2D equations governing the defect density and elastic displacement fields based on the nonlocal constitutive relations of the elasticity theory. Consider an isotopic elastic semi-infinite medium occupying a region , in a rectangular Cartesian coordinate system , where the origin is situated at any point on the plane boundary and points vertically downwards, that is, towards the bulk of the medium. The surface is supposed to be free from stresses. Let a plane elastic wave propagate along the -axis. Denoting and as the nonzero components of the displacement vector , we set , , and and consequently write for the strain tensor ():

Let us assume that an external energy flux (e.g., laser radiation) creates mobile atomic defects in a surface layer. The corresponding defect density profile results in a force that may induce strain field in medium. Let be the density of these defects of the th-type ( for vacancies (-defects) and for interstitials (-defects)). We limit our consideration to the case of only one type of atomic defects (for definiteness, -type defects).

The concentration field of atomic defects is dependent on temperature of the medium. One thus needs to know how the laser irradiation affects the local temperature of the surface at the laser spot. We will consider here situations where the laser irradiation only heats the solid (the light energy absorbed by the medium is transformed into heat) and that an equilibrium between laser radiation and the temperature field () is reached on time scales much shorter than the characteristic time scale of defect density evolution. Typically, the time scale for equilibration between photon absorption and defect generation is on the order of picoseconds, while that for defect diffusion is of the order of microseconds. We also assume that the contribution of thermal strains to deformation fields is negligible compared to lattice dilatation due to atomic defects and the phase changes and chemical reactions in the medium are absent.

In this paper, we will consider the problem of the wave propagation in an elastic solid irradiated over a large area by CW or pulsed lasers. Furthermore, we will assume that the temperature field has reached its equilibrium value. Its evolution is sufficiently slow compared to atomic defect generation and can be considered as quasistationary. The solution of the heat conduction equation for this case is given by Duley [22].

Taking into account the defect generation, the constitutive equations of an isotropic nonlocal elastic solid are where is the nonlocal stress tensor, is the classical stress tensor, and are Lamé’s constants [23], is the strain tensor, is the dilatation, and . The functions and are known as atom-atom (short-range) and atom-defect (long-range) interaction kernels or moduli of nonlocality, which decay smoothly with distance. They reflect the influence of an independent constitutive variable at a point on a dependent constitutive variable at . Nonlocal moduli are physical properties of materials like other physical constants and need to be determined experimentally.

The basic difference between classical and nonlocal elasticity is in the presence of the volume integrals in (2) which indicates that the stress at depends on the strain and defect density at all other points of the body, at time . This signifies that the distant neighbors of a point have a role to play in the propagation of waves.

In absence of body forces the equations of motion and the equations of defect density dynamics have the following form: where is the density of the medium, is the thermal-fluctuation generation rate of atomic defects at sites ( is the constant and is the Boltzmann constant), and is the relaxation time of defects in the absence; and are the formation and migration energies for the defects in crystals.

The first term in the right-hand side of (4) takes into account laser-induced generation of defects, the second term represents diffusion with a coefficient , the third term corresponds to the drift of defects under the influence of the force resulting from the nonlocal interaction of defects with an inhomogeneous strain field, and the fourth term describes the rate of their disappearance due to recombination processes. It is assumed that the generation rate () is spatially uniform.

The expression for interaction energy of a defect with the strain field in a nonlocal elastic medium is given by the formula In the limit, when the kernel functions and become Dirac-delta functions, (2) reduces to the classical constitutive equations of local elasticity theory.

Thus the final dynamic field equations obtained are integrodifferential equations for the functions and . It seems to be obvious that a rigorous solution of such equations encounters serious but not insurmountable mathematical difficulties. However, these equations can be reduced to the partial differential forms under certain conditions with physical admissible kernels. We consider that the long-range internal influences of particles of the body are rather rapidly with increasing distance from the particle. Thus the 2D-kernel functions that characterize the nonlocal interaction in the -direction may be approximated in terms of delta-like functions as

Besides the nonlocal kernels and are symmetric with respect to the interchange of the arguments . Then expanding in (3) the functions and in Taylor series around and neglecting the terms higher than second-order gradients, we can obtain the following expression for the nonlocal constitutive relation with strain and defect concentration gradients: where and are the intrinsic characteristic length scale parameters, characterizing the defect-atom and atom-atom interactions, respectively; ; is deformation potential of the defect ( is the change of the volume of the medium under formation of a single defect and is the bulk modulus). The gradient coefficients and are dependent on the lattice or atomic chain model used and the interatomic potentials assumed. The values of these coefficients are taken to be order of interatomic distance. For the interaction energy we have from (5) .

To obtain the nonzero stress components in terms of the displacement components we substitute (1) into (7):

Then from (4), two nonzero equations of motion can be obtained:

Suppose that the stress components and defect density field satisfy the boundary conditions

We can express the defect density field as ( is a spatially homogeneous solution; is small nonhomogeneous perturbations). Inserting in (4) the nonlinear terms, we get the linearised equation as

The system of (9) and (11) is closely coupled. in (9) depends on the defect density field () and in (11) depends on elastic displacement field (). The system of equations thus becomes highly nonlinear.

#### 3. Solution of the Problem

To solve the above problem, we apply the Fourier integral transform in the following form: Substitution of this into (9) gives where where is the frequency of wave propagation, is the wave number, and ; the phase velocity is given by and attenuation constant by , where and mean, respectively, the real and imaginary parts of . and are the nonlocal elastic moduli; is the nonlocal constant characterizing lattice deformation due to atomic defects; , and are unknown functions (amplitude functions).

The general solutions to (13) are where and are the roots of the equation where , are arbitrary constants, is the coupling constant of nonlocal defect-strain interaction, and

Then and are defined from (17) as follows:

Stress components are obtained by substituting (15) into (8):

#### 4. Dispersion Equations of the Wave in an Infinitive Medium

In this section as a special case we consider the generation of plane harmonic structures in unbounded nonlocal medium. Setting in (16) and (17) we obtain

Suppose now that one ignores the influence of the non-local properties of the medium and sets , , and . In this case (21) and (22) reduce to frequency equations ( and are the classical velocities of the longitudinal and shear waves, resp.), coinciding with the equations derived in the conventional local theory (see, e.g., [8]). In the nonlocal case, the frequency equation has, of course, a more complex structure.

To explore and delineate the strain and defect generation effects, we will seek solutions of (22) for small values of . For , (22) admits the following solutions (acoustical mode) and (diffusion mode). Now, for small , we may write Substituting (24) into (22), and, equating the coefficients of like powers of , we finally arrive to From (25) we obtain

If and the viscosity is taken into account [by adding in (9) the terms and , where ; ; and are the first and second viscosity coefficients] dispersion (28) describes attenuation of the amplitude of acoustic waves For the frequency spectrum of acoustic wave we have

Note that the frequency of acoustic waves is hardly changed. But the additional contribution to attenuation coefficient of waves arises. As attenuation decrement of acoustic waves is proportional to , according to (29), short waves fade fast.

If In this case there is a softening of frequencies of acoustic waves (instability of frequencies ) and this is related to taking into account the generation of atomic defects. It is necessary to notice that reduction of frequency occurs not up to zero but up to value . Equation (32) describes attenuation of the amplitudes of acoustic waves.

At from (22) we obtain dispersion relation (growth rate) of the volume concentration-elastic instability: or where .

Let us introduce the dimensionless growth rate , the wave number , the life-time of defects , and parameters . Dependence (34) in dimensionless variables takes the form

We consider two situations that differ from each other by the relation () of the squares of the characteristic lengths and . Figure 1 demonstrates the graphs of the growth rate calculated using (35) for two values of the ratio parameter and 10, respectively. The following typical values of parameters were used in the calculations: nm, , and cm^{2} s^{−1}. Comparison of the plots presented in Figure 1 shows that the dependences of the growth rates of the wave number for the two cases qualitatively differ from each other. It follows from Figure 1 that the growth rate exhibits a single maximum at for small values of and an additional maximum in the short-wavelength region at cm^{−1}. Therefore, depending on the relation between the characteristic lengths and , elastic-concentration structures have one or two maximum growth rates.

Figure 2 shows the dependences of the increment on the wave number for two values of the control parameter and 15 when the ratio parameter is . It follows from this figure that the growth rate has a single maximum not far above the threshold. However, at sufficiently high concentrations of the defects exceeding the threshold appears an additional maximum in the short-wavelength range. Therefore, as in the previous case, two gratings have the maximum growth rates.

The analytical expression for the maximum of the growth rates can be obtained provided that the nonlocal parameters are of the same order in (34) (i.e., ). Then, using (34), we have
where
is the critical value of the control parameter (defect generation rate). So, putting values into (37) for molybdenum: eV, erg cm^{−3}, K, and s, one can get the estimate s^{−1} cm^{−3}.

The dispersion curve has a maximum at , which determines the spatial scale of the dominant CS structures:

The maximum growth rate is

It is seen from (39) that the development of the elastic-concentration instability with the self-organization of CS structures in the volume of the solids occurs in a threshold manner when the defect generation rate exceeds the critical value: . At the elastic-concentration instability threshold (), the period tends to infinity . At high values of generation rates , the period asymptotically tends to its minimal value .

Near the threshold the time of formation of the CS structures that is critical slowing down takes place, which is characteristic of the phase transitions of the second kind.

The critical value of the generation rate derived here is directly proportional to recombination rate () and inversely proportional to deformational potentials () of defects. Putting into (37) for the generation rate, we get One notes that this dependence is very rapid, which is largely due to variation of the defect relaxation time with temperature.

From (38), it follows that the period of the volume ordered CS structures is proportional to the characteristic length of the defect interaction with the crystal-lattice atoms, which lies in the obtained nanometer range. As temperature is increased, the period grows for a particular value of . For a particular value of , the period decreases with an increasing generation rate ().

A periodic structure of vacancy clusters in the bulk of a solid irradiated with nitrogen ions was first observed in 1971 by Evans in pure molybdenum [24]. The medium temperature was 1040°K. Subsequently, the same effect had been observed in other metals such as aluminum, tungsten, nickel, and niobium [25, 26]. Vacancy cluster sizes in the lattice were 2–4 nm, and values of the lattice parameter change in the range of (20–60) nm. At nm (for ) from (38) we obtain an estimate for the period of the volume structure nm. This is in qualitative agreement with the above-mentioned experimentally observed interval.

Note that the linear theory considered in this section describes early stage of the development of an elastic-concentration instability only. However, the nature of generated ordered volume structures (due to an instability) and the amplitudes of these structures as functions on material and irradiation conditions can only be determined by considering the influence of elastic nonlinear effects in the model.

#### 5. Dispersion Equations of the Surface Concentration-Strain Structures

Substituting expressions (20) into boundary conditions (10) and taking into account expressions (18), we find a set of three linear algebraic equations:

The condition of existence of a nontrivial solution of this system yields where , and are given by (16) and (19).

Equation (42) shows dispersive character of the Rayleigh waves propagating in an elastic half-space. This dispersive character of these waves propagation arises due to atomic defect generation in the medium with long-range interactions.

Note that the form of this dispersion equation is identical to the corresponding equation obtained in the context of concentration-elastic coupled theory, using classical constitutive equations of motion, though each of , and represents different expressions for classical elasticity and nonlocal elasticity.

Equation (21) is generalized form of the classical local Raleigh equations. The assumption of vanishing defect generation in this equation gives us exactly the same results as those obtained by Nowiński [17] for solids with long-range interactions; that is,

#### 6. Solution of Dispersion Equation

A general discussion of (42) presents considerable difficulties. The structure of (42) is identical to the structure of the dispersion equation for DEI considered in [8]. This equation can be readily solved numerically or graphically, but in some cases of interest it can also be solved analytically. We shall now consider the results for some limiting cases. We can rewrite (42) as where and .

Equation (45) has solutions describing qualitatively different types of instability: instability of frequencies of surface acoustic waves and generation of ordered surface (static) concentration-strain structures.

If the coupling term is small () the expressions that appear in the dispersion equation (25), to a first approximation, then read Then we may express the dispersion equation in the following form: Here For we obtain the nonlocal elastic solution given by Nowiński [17]. Introducing the dimensionless variable we may write the dispersion equation (47) in the following form: Now we may consider that is the increment of the value of due to . Then (49) can be rewritten in the form of where is a root of the equation .

Assuming that , we may expand both sides of (50) into a Taylor series in the vicinity of the point . Retaining only the first two terms, we obtain where and .

The real part of characterizes the change of the phase velocity and its imaginary part () defines the attenuation constant. We have This equation gives the following expressions:where , , and .

Then for the phase velocity we have Since , in this case there is a softening of frequencies of acoustic waves (instability of frequencies and ) and this is related to taking into account the generation of atomic defects. Equation (53b) describes attenuation of the amplitudes of acoustic waves.

Expanding and in powers of the small parameter (), we obtain from (21) the expression for the growth rate of the surface periodical static structures in the form where .

Introducing the dimensionless growth rate , the wave number , the life-time of defects , and parameters , , dependence (34) can be represented as

The analysis shows that the growth rate exhibits a single maximum at for small values of ratio parameter () and an additional maximum () in the short-wavelength region for large values of . Thus depending on the relation between the characteristic lengths and , elastic-concentration structures can have one or two (with and ) maximum growth rates.

Besides the growth rate has a single maximum not far above the threshold. However, at sufficiently high concentrations of the defects exceeding the threshold appears an additional maximum in the short-wavelength range. Therefore, as in the previous case, two gratings have the maximum growth rates.

The obvious expression for the maximum of the growth rates can be obtained if the scale coefficients () are of the same order in (34):

The growth rate reaches maximum at where is the critical value of the control parameter, at exceeding of which the self-organization of surface-ordered CS structures becomes possible due to development of the elastic-concentration instability. As a result, the surface deformation lattice and the corresponding relief modulation with the period grow in time on the surface, as well as the defect concentration lattice (defect accumulations at the extrema of the surface relief).

The spatial scale of surface SC structures with wave number is given by

The maximum growth rate is Near the threshold the time of formation of the CS structures that is critical slowing down takes place, which is characteristic of the phase transitions of the second kind.

The period of the ordered CS structures formed on the surface is proportional to the characteristic length of the nonlocal defect-atom interaction, which lies in the obtained nanometre range. As temperature is increased, the period grows for a particular value of . For a particular value of , the period decreases with an increasing generation rate ().

Thus at certain conditions concentration-elastic instabilities in a system of atomic defects in unbounded condensed media as well as on the free surface of a solid half-space can be developed. We have observed that if the pump parameter is above the critical value, due to concentration-elastic instability, coupled strain-defect periodic nanometer sized structures on the surface of solids arise. A mechanism on the development of the instability is due to the coupling between defect dynamics and the elastic field of the solids. Laser radiation (or, in general, a flux of particles) generates high concentrations of mobile atomic defects in the surface layer of the irradiated material. When a fluctuation harmonic of the elastic deformation field appears in a medium because of the generation of atomic defects, the strain-induced drift of atomic defects occurs. This is a consequence of nonlocal defect-strain interaction. The strain-induced flux of defects gives rise to periodic spatial-temporal fields of the defect concentration. The redistribution of defects creates forces proportional to their gradients. These forces lead in turn to additional growth of strain fluctuations. When a certain critical value of the defect density or rate of defect generation is exceeded, diffusion-elastic instabilities develop as a result of positive feedback, which result in the formation of ordered CS structures.

As an example let us consider the formation of periodic surface CS structures in laser-irradiated semiconductors (in particular, *CdTe*). To evaluate the concentration of generated lattice defects (), we consider here conditions when the duration of a laser pulse () exceeds the defect-relaxation time (). In this case the density of defects on the surface of the solid reaches a steady-state value
where
is the defect generation rate, is the steady-state value of the temperature field on the surface, and and are constants. If (uniform irradiation) and the optical absorption length () is sufficiently less than heat diffusion length , the maximum temperature rise at the substrate surface owing to the laser pulse action may be evaluated [22] as
where is the reflectivity coefficient, is the thermal conductivity coefficient, and is the thermal diffusion coefficient. So putting values into (64) for *CdTe* ( W/cm°K, , W/cm^{2}, pulse duration s, and ) one can get °K. Then, taking cm^{3}, erg/cm^{3}, , and eV, a value of 1/cm^{3} may be estimated for the critical density of defect concentration (), which is several orders of magnitude less than the concentration of the host atoms and shows that this mechanism of the formation of ordered structures may be realized on practice. For the period () of the resultant surface CS structure we have an estimate of nm, which follows from expression (60) (for typical values of parameters g/cm^{3}, = 30 eV, cm/s, , nm, cm^{2}/s, and s). The maximum growth increment of the instability is ; the instability increment exceeds the reciprocal of the duration of a laser pulse acting on the surface of a solid ().

The appearance of periodic nanostructures of the surface relief was detected experimentally [27] when the surface of the p-*CdTe* semiconductor crystal was irradiated by single pulses of the multimode ruby laser (pulse duration s) with a homogeneous distribution of intensity. The period of the relief structure (with lateral size of the nanoclusters (19–23) nm) was 74 nm, which is in agreement with the predictions of our theory.

#### 7. Conclusions

Based on the nonlocal elastic theory, which accounts for the small-scale effect by assuming the stress at a reference point as a function of the strain at every point in the body, the governing equations for the elastic displacement vector and defect density fields in an isotropic laser-excited solid have been derived. The defect-density dynamics are controlled by the generation and annihilation processes and also the strain-stimulated transport of defects. We have provided exact formal solutions for the displacements and defect density fields in an elastic solid half-space. A linear analysis is used to obtain dispersion equations of surface waves through the solution of three coupled integropartial differential equations, by assuming that the surface of the solids is stress-free. A procedure for determining the phase velocity and the attenuation (amplification) factor has been discussed. The proposed analysis is applied for the special case of very short wavelengths. The phase velocity and attenuation constants of the waves get modified due to the defect subsystem of the solid. We have shown that the dispersion equations predict two types of the instability: instability of frequencies of acoustic waves and generation of ordered (surface or volume) structures.

We observed that the small-scale effect plays an important role in the dispersion behaviors for larger wave numbers. We have demonstrated that due to the nonlocal character of the defect-atom and atom-atom interactions the dispersion curves of the instability have two maxima sufficiently far above the instability threshold. When the regular surface relief is generated, two maxima of the growth rate must give rise to two scales of the surface relief. These two scales of the relief are the characteristic of the surface nanostructure self-organization upon both laser and ion-beam irradiation of semiconductors.

In the case of the laser molecular-beam epitaxy (MBE) and the laser-controlled deposition process of atoms on substrates, the extrema (minima or maxima) of the surface relief formed by the CS structure serve as the nucleation and growth centers of the nanoparticles. Similarly, in the processes of multipulse laser etching, the rapid material removal also takes place on the extremes of the CS field, where the removal rate takes maximum values. As a result, on the surface of substrate, size distribution functions of the nanoscale nonhomogeneities having two maxima can be formed. A transformation of the unimodal distribution to the bimodal (with two maxima) distribution upon a variation in the irradiation regime (control parameter) is typical of both the MBE and the laser-induced surface nanostructuring of metals and semiconductors. Formation of bimodal metal (Au) nanoparticle size distribution function upon the pulsed laser deposition was detected experimentally in [28].

The results obtained in this work are expected to be helpful in designing the nanostructures in small-scale physical devices.

#### References

- F. Mirzade, “Nonlinear strain waves interacting with laser-induced carries of the local disorder,” in
*Laser Technologies of Materials Treatment*, V. Y. Panchenko, Ed., pp. 220–277, Fizmatlit, Moscow, Russia, 2009, Russian. View at Google Scholar - F. Mirzade, “Elastic wave propagation in a solid layer with laser-induced point defects,”
*Journal of Applied Physics*, vol. 110, no. 6, Article ID 064906, 2011. View at Publisher · View at Google Scholar · View at Scopus - F. Mirzade, “A model for the propagation of strain solitary waves in solids with relaxing atomic defects,”
*Journal of Applied Physics*, vol. 103, no. 4, Article ID 044904, 2008. View at Publisher · View at Google Scholar · View at Scopus - F. Mirzade, “Influence of atomic defect generation on the propagation of elastic waves in laser-excited solid layers,”
*Physica B*, vol. 406, no. 24, pp. 4644–4651, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. M. Kosevich,
*Physical Mechanics of Real Crystals*, Nauka, Moscow, Russia, 1986. - F. Mirzade, V. Y. Panchenko, and L. A. Shelepin, “Laser control of processes in solids,”
*Uspechi Phyzicheskih Nauk*, vol. 39, no. 1, pp. 3–30, 1996. View at Google Scholar - F. Mirzade, “Concentration-elastic instabilities in a solid half-space,”
*Physica Status Solidi B*, vol. 246, no. 7, pp. 1597–1603, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Mirzade, “On diffusion-elastic instabilities in a solid half-space,”
*Physica B*, vol. 406, no. 1, pp. 119–124, 2011. View at Publisher · View at Google Scholar · View at Scopus - D. Walgraef and N. Ghoniem, “Modeling laser-induced deformation patterns: nonlinear effects and numerical analysis,”
*Journal of Computer-Aided Materials Design*, vol. 6, no. 2, pp. 323–335, 1999. View at Google Scholar · View at Scopus - F. Mirzade, “Finite-amplitude strain waves in laser-excited plates,”
*Journal of Physics: Condensed Matter*, vol. 20, no. 27, Article ID 275202, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. C. Eringen,
*Nonlocal Continuum Field Theories*, Springer, New York, NY, USA, 2002. - A. I. Kunin,
*Elastic Media with Microstructure*, Springer, New York, NY, USA, 1982. - A. C. Eringen, “On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves,”
*Journal of Applied Physics*, vol. 54, no. 9, pp. 4703–4710, 1983. View at Publisher · View at Google Scholar · View at Scopus - J. N. Reddy, “Non-local nonlinear formulations for bending of classical and shear deformation theories of beams and plates,”
*International Journal of Engineering Science*, vol. 48, no. 11, pp. 1507–1518, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Alibeigloo, “Free vibration analysis of nano-plate using three-dimensional theory of elasticity,”
*Acta Mechanica*, vol. 222, no. 1-2, pp. 149–159, 2011. View at Publisher · View at Google Scholar · View at Scopus - B. Arash and Q. Wang, “A review on the application of nonlocal elastic models in modeling of carbon nanotubes and graphenes,”
*Computational Materials Science*, vol. 51, no. 1, pp. 303–313, 2012. View at Publisher · View at Google Scholar · View at Scopus - J. L. Nowiński, “On non-local theory of wave propagation in elastic plates,”
*ASME Journal of Applied Mechanics*, vol. 51, no. 3, pp. 608–613, 1984. View at Publisher · View at Google Scholar - J. L. Nowinski, “On the propagation of thermoelastic waves in media with long-range cohesive forces,”
*Journal of Thermal Stresses*, vol. 10, no. 1, pp. 17–27, 1987. View at Google Scholar · View at Scopus - K. L. Pan and J. Fang, “Interaction energy of dislocation and point defect in bcc iron,”
*Radiation Effects and Defects in Solids*, vol. 139, no. 2, pp. 147–154, 1996. View at Google Scholar · View at Scopus - K. L. Pan, “An image force theorem for a screw dislocation near a crack in non-local elasticity,”
*Journal of Physics D*, vol. 27, no. 2, pp. 344–346, 1994. View at Publisher · View at Google Scholar · View at Scopus - K.-L. Pan, “Interaction of a dislocation with a surface crack in nonlocal elasticity,”
*International Journal of Fracture*, vol. 69, no. 4, pp. 307–318, 1995. View at Publisher · View at Google Scholar · View at Scopus - W. W. Duley,
*Laser Processing and Analysis of Materials*, Plenum Press, New York, NY, USA, 1983. - L. D. Landau and E. M. Lifshitz,
*Theory of Elasticity*, Pergamon Press, Oxford, UK, 3rd edition, 1986. - J. H. Evans, A. W. MacLean, A. A. Ismail, and D. Love, “Concentrations of plasma testosterone in normal men during sleep,”
*Nature*, vol. 229, pp. 261–262, 1971. View at Publisher · View at Google Scholar - V. K. Sikka and J. Moteff, “Damage in neutron-irradiated molybdenum,”
*Journal of Nuclear Materials*, vol. 54, no. 2, pp. 325–345, 1974. View at Google Scholar · View at Scopus - V. F. Zelensky, I. M. Nekludov, and L. S. Ozhigov,
*Some Problems in the Physics of Radiation Damage in Materials*, Naukova Dumka, Kiev, Ukraine, 1979. - A. Baidullaeva, A. B. Vlasenko, L. F. Cuzan et al., “Formation of nanosize structures on a surface of
*p*-CdTe crystals due to the ruby laser single pulse irradiation,”*Semoconductors*, vol. 39, pp. 1064–1067, 2005. View at Google Scholar - J. Gonzalo, A. Perea, D. Babonneau et al., “Competing processes during the production of metal nanoparticles by pulsed laser deposition,”
*Physical Review B*, vol. 71, no. 12, Article ID 125420, 8 pages, 2005. View at Publisher · View at Google Scholar · View at Scopus