AM  Vol.11 No.11 , November 2020
Plasma Model of Generation and Slip of Linear Defects in Crystalline Materials
Abstract: In dielectrics and semiconductors, a plasma model of the generation and slip of dislocations is considered, where under shock loads in a generalized space of rectangular pulses an alternating field forms a distribution of pairs of photoelectrons and cations; these electrons with velocities Ve create δ-collisions with cold plasma from free electrons and holes with masses me and mh (mh ≫ me), they emit and absorb longitudinal electron plasma waves whose phase velocities wpw / kpw are close to or are equal to the velocities Ve, while the frequencies wpw and wave numbers kpw of the wave packet of plasma waves are complex, the short-wave components of this wave packet at kpw ⋅ ae ≫ 1 (ae -Debye screening radius) decay in the core linear defect, and its long-wavelength components propagate in the region of the medium surrounding the core of the defect at kpw ⋅ ae <≅ 1. When a defect is generated, the distribution of cations under the influence of the internal Coulomb field shifts to the region of the first peak (protrusion) of the electron plasma wave, thereby forming a vacancy valley. When sliding under the influence of an external electric field, a cationic plasma wave consisting of a vacancy valley and two cationic protrusions moves against the background of an additional potential relief created by an electron plasma wave near the core of the defect. It has been shown that δ-collisions create flows of dynamic large-scale correlations of plasma fluctuations in the form of asymptotics of different-time correlators of density and potential fluctuations as t → +∞.

1. Introduction

At present, the fundamental theory of the nuclei of linear defects in semiconductors is the Reed discrete static model [1] [2], where it was assumed that an additional half-plane of the edge dislocation introduces a chain of free or unsaturated bonds into the crystal. Due to the attraction of free electrons to these bonds in germanium, as well as in most cases, a negative linear charge of a dislocation core arises in silicon, which forms a cylindrical space charge in the environment surrounding the core with a radius R from 1 to 10 μm. Using the Poisson equation and the electrostatic energy of trapped electrons, Reed found the work of space charge formation for one additional electron in the chain

A p r = f E 0 ( 3 2 ln f f c 0.866 ) (1)

where f is the degree of filling of the chain is determined by the Fermi distribution function

f = [ 1 + exp ( E d E f ) k T ] 1 (2)

where E d is the dislocation level within the band gap; E f is the Fermi level; T is the absolute temperature; E 0 is the interaction energy of neighboring electrons in the chain, E 0 = e 2 / ϵ a c h ; a c h is the distance between neighboring trapped electrons; 1 / a c h is the density of additional electrons per unit length of the dislocation; ϵ —dielectric constant; f c is distribution function f = f c at R = a c h ; c is the distance between adjacent unsaturated bonds, the ratio a c h / c is of the order of 0.1. Reed also found the potential of the field ϕ 0 created by all other electrons of the linear charge and positive space charge at the location of an arbitrary injected electron

ϕ 0 = e ϵ c f ( 3 ln f f c 1.232 ) (3)

The numerical value of was confirmed experimentally by scanning with an electron beam except for its sign ( [1], p. 112). Of particular interest is the dependence of the electrostatic potential of the barrier ϕ b on the distance r from the dislocation line ( [1], Figure 8.25, p. 178), where ϕ b ( r ) obtained at an average concentration of conduction electrons n e s = 1 × 10 15 cm 3 ; a c h = 2.85 nm ; c = 0.4 nm ; c / a c h = 0.14 ; E d = 0.4 eV below the boundary of the conduction band, and consists of two branches: ascending at r = 0 from 0.18 V to the maximum of the curve 0.27 V at r m = 6 ÷ 10 nm and descending from r m to r f 100 nm . It should be noted that the distance r m is associated with the wavelength λ i e of the injected electron

r m = e 2 2 π ϵ k T ln R λ i e (4)

Here the question arises: How does the generation and slip of linear defect nuclei occur in the Reed model?

On the other hand, in the limiting case, near the electric breakdown and shock compression of crystalline dielectrics and semiconductors, there occurs a solid-plasma phase transition [3], where the combined influence of electric fields is E o u t = 10 9 - 10 10 V / m and hydrostatic pressures p = 10 9 - 10 11 Pa lead to two effects: 1) reduction of the band gap E g * ( p ) according to the linear law

E g * ( p ) = E g 0 ( p ) + a p p , (5)

where a p is the baric coefficient expressed in terms of the compressibility coefficient K and has a negative value

a p = γ K (6)

where 9 × 10 12 K 200 × 10 12 ; γ ( 2 ÷ 3 ) × 10 19 J the generation of electron-hole pairs, the probability w of occurrence of which has the form [4] [5]

w = n e h / N = ( e E o u t a 0 * ) 2 2 π E g * exp ( E g * e a 0 * E o u t ln 1 / α ) (7)

where n e h and N are the average densities of electron-hole pairs formed per unit time and valence electrons, respectively; a 0 * is the effective lattice parameter; α is the ratio of the width of the valence band to the width of the conduction band. It was shown in [3] that when injected electrons move in the volume of a shock wave in the velocity range V i e = 4 × 10 3 ÷ 1 × 10 6 m / s and pressures p = 1.2 × 10 10 ÷ 1 × 10 11 Pa for alkali halide crystals of the NaCl type, the value of E g * ( p ) decreases from 8.6 to 1.5 eV.

It is known [6] that the energy of the ground state of a solid-state plasma is mainly determined by electrostatic energy, and as a second-order correction it contains the correlation energy W, which is divided into two components: W 1 , due to large-scale Coulomb correlations in the relative motion of electrons, and W s , due to the Pauli principle in the form of small-scale correlations with an effective radius of 0.1 - 0.2 nm. An analysis of the works [7] [8], which consider the processes of damage formation in ionic crystals under the influence of intense laser radiation, shows that here the decisive role is played by the influence of plasma arising both at the boundaries of blocks, grains, cleavage, and in the volumes of pores and cavities.

The well-known Reed model does not reflect the combined dynamics of conduction electrons captured by unsaturated bonds and the atomic collective near the dislocation line during generation and slip of the dislocation core. The question naturally arises: How is the generation and sliding of linear defects in a solid-state plasma under impact loads?

The aim of this work is to build a plasma model of the processes of generation and slip of linear defects in dielectrics and semiconductors.

2. The Theoretical Model

We find the Debye radius of screening of a solid-state plasma ( [9], p. 145)

a e = { k T 4 π n e e 2 } 1 / 2 (8)

where e is the electron charge; n e n e m , n e s , n e d are the average electron densities in metals, semiconductors, and dielectrics, respectively. At room temperature and n e m = 1 × 10 22 - 6 × 10 23 cm 3 ( [10], p. 380) the a e m radius in metals is 25 ÷ 30 nm, at low temperatures 4 - 6˚ K a e m = 4 - 5 nm . In semiconductors doped with impurity atoms, n e s 10 16 - 10 18 cm 3 , the radius a e s , according to (8), takes values from several mm to several cm. In pure Ge, n e s 10 13 - 10 14 cm 3 is comparable to or slightly larger than n e d alkali halide crystals.

At the first stage of the theory, we assume that under the influence of impurity ions, external electric fields and shock waves, the material already has an equilibrium neutral distribution of free electrons and holes with masses m e and m h for which: 1) the effective mass of electrons is m e * m h * -effective mass of holes; 2) the relaxation time of free electrons is τ r e τ r e c -the recombination time of the injected electrons and holes. We also emphasize that upon impact of the surface of a Ge, Si type semiconductor, an intermittent field injects valence electrons from atoms at lattice sites, turning them into cations; in the case of an insulator of the type Na+Cl injection of bond electrons and the formation of light cations Na+ take place.

Let us consider in the Cartesian coordinate system a certain volume V of crystalline material, where an alternating field E x i f acts in the generalized space of rectangular pulses along the x axis in the standing wave regime E 0 x i f δ ( t ) and E 0 x i f δ ( t + τ ) ( τ is the pulse duration) [11] in the slip plane x0z. The field E x i f creates in the right ( x > 0 ) and left ( x < 0 ) half-spaces of volume V a pair of photoelectrons and cations forming two flat beams with velocities V e and V e , while the cation velocities v c a t in both half-spaces will be consider small: v c a t V e . Each beam in its half-space creates δ-ollisions, through which it interacts with the plasma, while the chains of photoelectrons emit and absorb longitudinal plasma waves—harmonics exp [ i ( k r ω t ) ] of only certain frequencies ω p w and wave vector k p w , selectable using equations

ω p w k p w V e = 0 (9)

ε l ( k p w , ω p w ) = 0 (10)

where ε l ( k p w , ω p w ) is the longitudinal dielectric constant of the plasma, and obtained from the expression of the Landau collision integral ( [9], p.236). In our case, we will consider the intrinsic plasma of the material to be cold: v e d , v e s V e , where v e d , v e s are the average thermal velocities of plasma electrons in dielectrics (d) and semiconductors (s). Plasma electrons will be elastically scattered with respect to the direction of the beam only forward at large angles χ ( [12], p. 64). The roots of the dispersion law (10) are complex ω p w = ω p w + i ω p w , γ = ω p w , γ is the plasma wave attenuation decrement according to the law exp ( γ t )

γ = π 8 Ω e ( k p w a e ) 3 exp [ 1 2 ( k p w a e ) 2 3 2 ] (11)

where Ω e is the plasma frequency of electrons, Ω e = { 4 π n e e 2 / m e } 1 / 2 ; a e a e d , a e s . It should be noted here that the real part ω p w при k p w a e 1

ω p w = Ω e ( 1 + 3 2 ( k p w a e ) 2 ) (12)

and the imaginary part

ω p w = 1 2 Ω e ε l ( k p w , ω p w ) (13)

exponentially small with ε l . Hence, a plasma wave in the form of a long-wave harmonic is propagating if γ is exponentially small for k p w a e 1 . In contrast, short-wavelength harmonics are non-propagating; they make the main contribution to the converging Balescu-Lenard collision integral S t f B L and essentially form the core of a linear defect. Indeed, as shown in ( [9], p. 238), it is the region k p w a e 1 that makes the main contribution to the integral S t f B L . It should be noted here that the short-wave harmonics of the electron plasma wave create an internal Coulomb attraction of the cation distribution.

On the other hand, plasma waves can be considered as fluctuation waves. For times Δ t that are small compared with the period of the rectangular pulses ω r e p 1 , the plasma can be considered collisionless. If, within the standard procedure, the distribution function of plasma electrons f e and the scalar potential φ e are divided into the regular f ¯ e , φ ¯ e and fluctuation parts δ f e , δ φ e : f e = f ¯ e + δ f e and φ e = φ ¯ e + δ φ e , then the system of equations describing such a plasma has the form ( [9], p. 256):

f ¯ e t + v e f ¯ e r e φ ¯ e r f ¯ e x = e δ φ e r δ f e p (14)

Δ φ ¯ e = 4 π e f ¯ e d 3 p (15)

where the right side of (14) is the collision integral. It is clear that the terms quadratic in δ φ e and δ f e in these equations describe the effect of collisions on fluctuations. Neglecting these terms and taking f ¯ e = f ¯ e ( p ) ; φ ¯ e = 0 , according to [9], we obtain the equations for fluctuations δ f e and δ φ e for Δ t ω r e p 1

δ f e t + v e j δ f e x j e δ φ e x j f ¯ e p j = 0 (16)

Δ δ φ e = 4 π e δ f e d 3 p (17)

It should be noted here that f ¯ e ( p ) reflects the nonequilibrium stationary state of a collisionless plasma, in the absence of an external field it depends only on the momenta of plasma electrons and is a stationary solution (14), (15). Equations (16) and (17) allow us to construct the asymptotics of different-time correlators as t + δ f e ( t 1 , r 1 , p 1 ) δ f e ( t 2 , r 2 , p 2 ) and δ f e ( t 1 , r 1 , p 1 ) δ φ e ( t 2 , r 2 , p 2 ) . In ( [9], p. 260), Fourier images or spectral correlators of these asymptotics were found

( δ φ 2 ) ω k = 32 π 3 e 2 k 4 | ε l ( ω , k ) | 2 f ¯ e ( p ) δ ( ω k V e ) d 3 p (18)

( δ φ e δ f e ) ω k = e k ω k V e + i 0 f ¯ e p ( δ φ 2 ) ω k + 8 π 2 e k 2 ε l ( ω , k ) f ¯ e ( p ) δ ( ω k V e ) (19)

All three components of the bulky expression ( δ f e δ f e ) ω k (19) also contain the product f ¯ e ( p ) δ ( ω k V e ) , and the first (…)1-initially, and two the rest (…)2, (…)3 – after finding f ¯ e p , while (…)1 represents in the momentum space a δ-peak with an amplitude of ~ne, and (…)2, (…)3 in k -space the line spectrum with envelopes k 2 and k 1 , respectively.

The process of generating nuclei of linear defects is divided into several stages: the formation of short-wavelength harmonics of the plasma wave δ f e and their stabilization by moving the distribution of photoelectrons to the region of the first trough δ f e , and the distribution of cations to the region of the first peak (protrusion) δ f e and the formation of a vacancy valley. Indeed, using the diffusion equations, or rather the mass transfer of cations and photoelectrons ( [9], p. 118; 2, p. 157)

n c a t t = d i v [ μ c a t k T x n c a t + n c a t μ c a t z e E i n s ( δ f e ( + ) 1 ) ] (20)

n p h e t = d i v [ μ p h e k T x n p h e + n p h e μ p h e e E i n s ( δ f e ( ) 1 ) ] (21)

where n c a t and n p h e , μ c a t and μ p h e are the average distribution densities of cations and photoelectrons and their mobility; E i n s ( δ f e ( + ) 1 ) and E i n s ( δ f e ( ) 1 ) are the internal Coulomb fields caused by the distribution of plasma electrons in the regions of the first protrusion and the first trough δ f e , we will determine

the conditions for stabilization of the defect nucleus. At n c a t t = n p h e t = 0 , the

drift currents caused by the fields E i n s ( δ f e ( + ) 1 ) and E i n s ( δ f e ( ) 1 ) are balanced by mass transfer currents from the gradients densities n c a t and n p h e , while the photoelectrons are redistributed to the potential valley of the first trough δ f e , and the cations to the potential valleys of the first protrusion δ f e both on the left and on the right, provided that such stabilization time is shorter than the recombination time at the holes, and this is always true [2] [9]. In other words, the distributions of photoelectrons and cations are in quasistable states in an additional potential relief created by the short-wavelength harmonic δ f e . In this case, the variation of the force acting on the photoelectron (cation) from the side of external and internal electric fields

F j = F j o u t ( E j o u t ) n e p x j [ e E j o w n + e φ e p ( x j x j ) f e p ( x j x j ) d 3 x j ] (22)

where p = d , s ; n e p = f e p is the average density of plasma electrons; The φ e p is scalar potential of photoelectrons (cations) in the Coulomb field of plasma waves and the distributions of cations (photoelectrons) reflect either oscillations at the bottom of the corresponding potential valleys or tunneling through barriers when the threshold beam velocities V e t h r are reached. Note that photoelectrons and cations under the influence of external fields move through a periodic sequence of potential barriers + U ( x ) and wells U ( x ) , while the photoelectron distribution is Boltzmann’s: f ¯ p h e ~ n p h e exp ( ε / k T e ) , where ε = p x 2 / 2 m e ± U ( x ) , U ( x ) = e φ e . If we put ω r e p δ f e on the right side of equation (16), and then multiply the left and right sides of (16) by the variation δ F x and average in direct space over coordinates and time, then the correlator obtained at different times can be reduced to correlators (18), (19), (19'). In what follows, we restrict ourselves to the case of low temperatures, when such nonlinear effects as scattering by lattice vibrations and ionized impurities do not play a significant role [2]. Here, in δ-collisions, the equilibrium distribution f 0 is displaced in the momentum space along the p x axis abruptly by | P x t h r | both in the positive and negative directions of this axis without changing the shape f 0 ( p ) . Suppose that, under the influence of an external uniform electric field E o u t , the distribution of photoelectrons in the defect nucleus shifts in the positive direction of the + p x axis. According to (21), the second peak (protrusion) δ f e ( ) 2 serves as a barrier to the passage of the drift current of these electrons. In turn, the cations in δ f e ( + ) 1 in one direction + p x are affected by the Coulomb attraction of the photoelectron distribution, which plays a decisive role, and the less significant Coulomb attraction of the hole distribution due to the small density, and also the field E o u t in the opposite direction. Here, the barrier for cations is the trough f e ( ) 1 . In other words, when sliding, when the photoelectron velocity reaches the threshold values V s l t h r > V g e n t h r , the ion plasma wave, consisting of a vacancy valley and two cationic protrusions, moves against the background of an additional potential relief created by the short-wavelength harmonic δ f e .

Let us consider the interaction of two dislocations d1 and d2 with arbitrary edge and screw components of the Burgers vector. In the framework of the plasma model, the long-wave propagating harmonics of plasma waves d1 (d2) interact with the oscillation currents J o s c e and rotations J t u r n e in the additional potential relief of the dislocation core d2 (d1) created by short-wave harmonics that do not propagate in this core i.e. J o s c e , J t u r n e ( d 1 ) J o s c e , J t u r n e ( d 2 ) . The smoothness of electric fields from J o s c e and magnetic fields from axial components of spiral currents from d1 (d2) to d2 (d1) core allows us to use quasiclassical dynamics in the form of Hamilton equations

r ˙ e = H s p e (23)

p ˙ e = H s r e (24)

where instead of the Hamiltonian of the system H s of two subsystems of pairs of photoelectrons and cations d1 and d2 we substitute the interaction Hamiltonian H 1 in the interaction representation [11]

H 1 = μ ν J μ ν ( r ) A μ ν ( r r ) d r , ν = e , c a t ; μ = o s c , t u r n (25)

where in our case the vector potential A μ ν depends on the difference between the coordinates of the charged particles of the source d1 (d2) and the test charge d2 (d1), which transformes (25) into the convolution integral, maximally simplifying the description of the force interaction d1 and d2 both in direct and in the reverse k -space by means of the Borel convolution theorem, transforming (25) upon its Fourier transform into the usual product of the Fourier components of currents. Here H 1 is directly proportional to the scalar product k 1 k 2 , where k 1 and k 2 are the wave vectors of charged particles by d1 and d2, and the nature of the interaction is determined by the law cos α . For 0 < α < π / 2 d1 and d2 are attracted, and for π / 2 < α < π they repel. Note that in the equilibrium state d1 and d2 in k -space, the right-hand side of (24) contains δ ( k 1 + k 2 ) . When considering accumulations of dislocations, a statistical description is applicable using binary correlation functions of currents and stresses σ k l i f , σ m n e l in the corresponding spaces.

Of particular interest is the experimental confirmation of the plasma model. In Figure 1, the AFM image of the source of Frank-Read dislocations located near the surface of a single-crystal KCl sample is shown. On the installation SPM Femtoscan, the change in the topography of the crystal surface under the influence of a pulsed magnetic field in the contact regime was studied by the method of scanning atomic force microscopy (AFM). Cantilevers Micromasch CSC-12 with a nominal tip radius of no more than 10 nm were used for scanning. Samples of a single-crystal KCl with dimensions 5 × 5 × 2 mm were obtained by cleavage along the cleavage planes, and were subjected to treatment with a weak pulsed magnetic field with induction B = 0.12T at room temperature in the time interval ~30 s. The maximum result of this action was found after the after effect time 48 hours and corresponded to the appearance of linear defects in the form of Frank-Read sources near the cleavage surface. To obtain an AFM image by atomic force microscopy, we used the contact quasistatic regime, where the tip of the probe was at a distance of 0.1 - 0.3 nm from the surface. Here, the source dislocation segment occupies a common field of view of 5 × 103 nm × 5 × 103 nm. AFM image includes three cross-sections of a common field: one central in the middle of the segment, two others near the attachment points. We give a qualitative analysis of the AFM image:

1) The surface relief on the AFM image consists of valleys and protrusions of various depths and widths. Cross-sections of the general field reflect the redistribution of the substance at the time of passage of the probe;

2) An insulated dislocation core, in the environment of which there are no linear (dislocations) and flat (interface, ledges, steps of terraces) defects, has a central trough relief with a depth of 4 - 6 nm and a width of 100 - 140 nm, which

Figure 1. AFM-image of the Frank-Reed source. The general field of view and three cross-sections: the central and two extreme near the fixation points.

corresponds to a vacancy valley, and two protrusions, the dislocation cores corresponding to the cationic protrusions are of approximately the same parameters, while the distribution of the protrusions and troughs of relief is symmetrical with respect to the axial plane passing through the minimum of the central relief trough (dislocation line) and has an approximately oscillatory character; a line approximately parallel to the abscissa axis, relative to which the sum of the areas above and below it is equal to zero, will be called the average level of the profile L ¯ p r ;

3) The approach of a dislocation to a deep step leads to a decrease in L ¯ p r from the dislocation line to this step and to an increase in the wavelength and

profile amplitude in this interval;

4) The convergence of two dislocations of the same sign near the source fixation points leads to a significant deepening: up to 9 - 10 nm, and expansion: up to 200 - 300 nm, the troughs of the vacancy valleys of the nuclei of the dislocations, separated by a wide protrusion, while L ¯ p r monotonically rises upward from both the left and on the right.

3. Discussion of the Results

1) Here a natural question arises: What is the reason for the appearance of the plasma model of linear defect nuclei? First, the impossibility of sufficiently reliably representing the essence and mechanisms of electroplastic and magnetoplastic effects [13]. Secondly, within the framework of the theory of atom-vacancy states, using the continuum model of concentration waves, it was shown that a dislocation is born as a solution of a nonlinear wave equation ( [14], p. 23) without revealing the structure of the defect nucleus. Third, the well-known atomic models [15] are based on the consideration of displacement fields caused by the introduction or removal of excess atoms into the crystal lattice without taking into account the density of charged particles, and the discrete model of the Reed dislocation core [1] [2] in semiconductors does not allow displaying generation and slip defect nuclei without corresponding displacements of excess atoms. Fourth, over the past 30 years, a number of new experimental methods have been created and tested, including atomic force microscopy, which allow us to take a fresh look at the problem of generation and slip of linear defect nuclei. It is clear that for a reliable experimental substantiation of the plasma model, a sufficiently large data set will be required.

2) There is a fundamental difference between the Reed model and the plasma model. In the Reed model, conduction electrons captured by unsaturated bonds of a chain of atoms along the edge of the half-plane inserted into the crystal form a linear charge. In the plasma model, as a result of the internal photoelectric effect, pairs of cations and valence electrons knocked out from semiconductor atoms are formed, as well as pairs of cations of alkali halide crystals, other dielectrics and bond electrons under the influence of an alternating field. In this case, the distribution of conduction electrons forms short-wave non-propagating harmonics of plasma waves, forming the structure of the defect nucleus and, at the same time, an additional potential relief of the crystal in which the distribution of the above pairs moves.

Cite this paper: Busov, V. and Vladimirovna, G. (2020) Plasma Model of Generation and Slip of Linear Defects in Crystalline Materials. Applied Mathematics, 11, 1167-1177. doi: 10.4236/am.2020.1111079.

[1]   Matare, H.F. (1971) Electronics of Defects in Semiconductors. Wiley Interscience. A Division of J. Wiley and Sons, Inc. New York, London, Sydney, Toronto.

[2]   Seeger, K.H. (1973) Semiconductors Physics. Springer-Verlag, Wien, New York.

[3]   Vershinin, Yu.N., Il’ichev, D.S. and Morozov, P.A. (2000) The Effect of Shock Compression of Solid Dielectrics on the Injection of Valence Electrons in Strong Electric Fields. Journal of Technical Physics, 70, 84-87.

[4]   Fever, P. (1952) Electronic States in Crystals under Large Over-All Perturbations. Journal Physic Review, 88, 92-101.

[5]   Keldysh, L.V. (1958) Behavior of Non-Metallic Crystals in Strong Electric Fields. Journal of Experimental and Technical Physics, 33, 994-1003.

[6]   Reims, S. (1972) Theory of Many-Electron Systems. North-Holland Publishing Company, Amsterdam, London.

[7]   Gorshkov, B.G. (1982) Investigation of the Mechanisms of Destruction of Ionic Crystals under the Action of Pulsed Laser Radiation of the Nanosecond Range. Transactions of Phys. Inst. Acad. Nauk, 136, 81-136.

[8]   Kazantsev, S.G. (2003) Laser Resistance of Promising Materials of IR-Optics. Mechanisms of Damage to Materials of Optical Elements of Power Optics. Optics of the Atmosphere and Ocean, 6, 390-401.

[9]   Lifshits, E.M. and Pitaevsky, L.P. (1979) Physical Kinetics. Nauka, Moscow.

[10]   Kittel, C. (1976) Introduction to Solid State Physics. 4th Edition, J. Wiley and Sons, New York, London, Sydney, Toronto.

[11]   Busov, V.L. (2019) Dynamic Equations of Evolution in the Nuclei of Linear Defects of Crystalline Materials during the Collision of Solids. PhysicalMesomechanic, 22, 91-96.

[12]   Landau, L.D. and Lifshits, E.M. (1973) Mechanics. Nauka, Moscow.

[13]   Golovin, Yu.I. (2004) Magnetoplasticity of Solids (Review). Journal of Solid State Physics, 46, 769-803.

[14]   Panin, V.E., Likhachev, V.A. and Grinyaev, Yu.V. (1985) Structural Levels of Deformation of Solids. Nauka, Siberian Branch, Novosibirsk.

[15]   Mirkin, L.I. (1968) The Physical Basis of Strength and Plastisity. Publishing House of Moscow State University, Moscow.