Formation Mechanism of Micro- and Nanocrystalline Surface Layers in Titanium and Aluminum Alloys in Electron Beam Irradiation

: The reported study discusses the formation of micro- and nanocrystalline surface layers in alloys on the example of Ti-Y and Al-Si-Y systems irradiated by electron beams. The study has established a crystallization mechanism of molten layers in the micro-and nanodimensional range, which involves a variety of hydrodynamic instabilities developing on the plasma–melt interface. Assuggested, micro-andnanostructuresformdue tothecombination ofthermocapillary, concentration and capillary, evaporation and capillary and thermoelectric instabilities. This mechanism has provided the foundation for a mathematical model to describe the development of structures in focus in the electron beam irradiation. The study has pointed out that thermoelectric ﬁeld strength E ≥ 10 6 V / m is attributed to the occurring combination of instabilities in micro- and nanodimensional ranges. A full dispersion equation of perturbations on the melt surface was analyzed.


Introduction
The development of pulsed high intensity technologies with the purpose of improving strength and tribology characteristics of high-duty machine elements represents one of the challenging issues confronting the present-day material science [1]. The material fracture is known to arise in the surface and near-surface layers because long-ranging stresses in them exceed the material strength limit [2]. It is necessary therefore to generate such a structure in the material surface, which would ensure high strength properties. Micro-, submicro-and nanocrystalline structures with a high level of hardness and wear resistance meet these requirements [3]. To form these layers concentrated energy flows are applied i.e., the heterogeneous plasma flow processing in conjunction with the succeeding low energy high current electron beam irradiation [4][5][6]. The latest research in this field [7,8] has highlighted microand nanocrystalline layers with a columnar structure ( Figure 1) form when irradiating titanium and silumin processed prior by the electric explosion plasma of yttrium powder. A thickness of a columnar crystallization layer is within a range from 1 to 1.5 µm (Figure 1a) in the Ti-Y system. Columns are inclined to the irradiation direction. This inclination may be connected with the longitudinal melt velocity, which, in turn, causes the Kelvin-Helmholtz instability. Transverse dimensions of columns vary from 100 to 300 nm ( Figure 1b) and depend on the irradiation conditions. Above this layer there A layer with cell crystallization in the Al-Si-Y system was determined to be 80 μm thick ( Figure  1c) [8]. Cell dimensions vary from 0.8 μm to 1.3 μm. (Figure 1d) [8]. An amorphous structure was detected in several zones of the surface layer, which is demonstrated by the halo effect on electron diffraction micro-patterns [8]. For silumin, cooling curves and isothermal curves of phase transition intersect; as a result, the volume of part of the crystal phase grows faster according to the Kolmogorov-Avrami formulation than the amorphous one. This approach including the idea that new phase nuclei arise sporadically owing to thermal fluctuations, makes it possible to follow the crystallization kinetics and estimate critical nuclei dimensions. However, it leaves open the questions of why a structure with unimodal and bimodal grain distribution (subgrains, second phase particles, etc.) forms and doping elements spread over grain boundaries and junctions. In our opinion, these issues can be resolved by the hydrodynamic approach [10][11][12]. It claims a crystal structure forms affected by a variety of hydrodynamic instabilities on the plasma-melt and meltcrystal interfaces. Small perturbations tend to occur on these boundaries because of thermal fluctuations [13,14]. To begin with, we assume these perturbations are harmonic. These perturbations increase under the instability initiating conditions, i.e., a positive real part of the A layer with cell crystallization in the Al-Si-Y system was determined to be 80 µm thick (Figure 1c) [8]. Cell dimensions vary from 0.8 µm to 1.3 µm. (Figure 1d) [8]. An amorphous structure was detected in several zones of the surface layer, which is demonstrated by the halo effect on electron diffraction micro-patterns [8]. For silumin, cooling curves and isothermal curves of phase transition intersect; as a result, the volume of part of the crystal phase grows faster according to the Kolmogorov-Avrami formulation than the amorphous one. This approach including the idea that new phase nuclei arise sporadically owing to thermal fluctuations, makes it possible to follow the crystallization kinetics and estimate critical nuclei dimensions. However, it leaves open the questions of why a structure with unimodal and bimodal grain distribution (subgrains, second phase particles, etc.) forms and doping elements spread over grain boundaries and junctions. In our opinion, these issues can be resolved by the hydrodynamic approach [10][11][12]. It claims a crystal structure forms affected by a variety of hydrodynamic instabilities on the plasma-melt and melt-crystal interfaces. Small perturbations tend to occur on these boundaries because of thermal fluctuations [13,14]. To begin with, we assume these perturbations are harmonic. These perturbations increase under the instability initiating conditions, i.e., a positive real part of the complex frequency, referred to as a growth rate α. A wavelength, when α = 0, is called a critical λ cr and instability starts developing from this point. It influences an initial size of structure elements. If a growth rate is maximal for a λ, this wavelength governs the most probable size of structure elements. It is a well-known fact that the gradient between temperature and concentration of doping elements arises in a molten material when processed by concentrated energy flows [15,16], resulting in thermal concentration and capillary instability. This instability is supposed to be one of the reasons for short wave micro-and nanometer periodic structures. If a material is irradiated in high vacuum conditions, the vapor return pressure is of significant importance, causing evaporation and capillary instability. The study [17] says a contribution of the evaporation pressure is decisive to the instability of capillary waves if a constraint k p ≥ − p v σ T is true, where p v and σ T -thermal coefficients of evaporation pressure and surface tension, respectively.
Moreover, thermal and concentration and capillary instabilities thermoelectric effects caused by a forming electric field are important to high temperature gradients~10 6 -10 12 K/m [18][19][20]. Once a temperature fluctuation arises, fluctuations of the electric field and its charge are registered. The effect of a constant electric field on this charge is critical for the power and instability provided so that the constant temperature gradient is high. The movement amplitude depends on the physical parameters of a liquid in focus and an active temperature gradient. The data on a molten layer thickness, when the thermal electric convection is a key mechanism responsible for the formation of cell structures, have demonstrated this effect is relevant for h~10-100 µm. As reported in studies [7,8,15,16], the thickness of molten layers ranges from~10 to 80 µm in electron beam irradiation. This fact suggests thermoelectric effects represent one of the basic formation mechanisms of micro-and nanodimensional structure elements in electron beam processing.
On the other hand, a periodical arrangement of atoms of the base material and alloying elements in high-speed crystallization cells and columns (Figure 1a,d) might stem from the arising and developing Mullins-Sekerka instability [21][22][23], which is a result of the inhomogeneous temperature and concentration of an alloying element on the crystallization front. However, this mechanism fails to explain the importance of substance evaporation and thermoelectrical effects in electron-beam treatment. Further research should be done to study this issue.
Thus, micro-and nanodimensional structures may form in the process of electron beam irradiation owing to the emergent combination of thermal, concentration, evaporation and capillary, and thermoelectric instabilities. Within this context this work seeks to investigate laws and conditions supporting the formation of the above structures when processing by low energy high current millisecond electron beams.

Problem Formulation
Generally, conditions initiating the instability of small perturbations on the melt-plasma interface are specified via searching for nontrivial solutions to homogenous differential equations with homogenous kinematic and dynamic boundary conditions. The existence of a parameter range providing such solutions means a spontaneous transition of the substance into a new state is possible [24]. Similarly to studies [11,12,17], we focus on a viscous heat conducting and evaporating liquid, which occupies a layer −h < z < 0 and −∞ < x< +∞ on a free surface z = η(x, t). If affected by an electron beam, a temperature profile in the liquid layer is T 0 (z), where T 0 -nonperturbed temperature ( Figure 2). These profiles change over time due to thermal conductivity and diffusion, but if a characteristic time of evolving perturbations is shorter than their transformation period, a function T 0 (z) can be considered constant in the instability analysis [10,17]. A melt temperature is a sum of a nonperturbed T 0 (z) and perturbed T(x,z,t) temperature. Assume a wavenumber vector of perturbed temperature and melt velocity is directed similarly to the X axis, these perturbations depend therefore on coordinates x and time t varying according to a function exp(ωt − ikx), where k-wavenumber, ω-complex cyclic frequency calculated as ω = α + iΩ, α-growth rate of perturbations, Ω-cyclic frequency. For a temperature gradient, as mentioned in Introduction, a thermoelectric filed arises in a liquid  [18][19][20] and contributes to the convective melt flow, it is determined from E = γ∇T, where γ -thermoelectric coefficient, ∇T-temperature gradient.
Metals 2020, 10, x FOR PEER REVIEW 4 of 15 thermoelectric filed arises in a liquid metal [18][19][20] and contributes to the convective melt flow, it is determined from A system of the Navier-Stokes equation and convective thermal conductivity equation is linearized and written as follows: where u, w -disturbance vector components of velocity, ρ-density, ν-kinematic viscosity, χ-temperature conductivity, G0-temperature gradient, p, T, η-pressure and temperature perturbations. No elements in the Equation (1) are related to the electric field, which is acceptable for conductive liquids [25]. Here, an electric field effect is included into boundary conditions on the surface z = 0: where η -perturbation of the surface along the axis z, -surface tension at a melting temperature of the material Tm, T σ -temperature coefficient of the surface tension, pe-electric field-melt pressure, el xz τ -tangential component of Maxwell stress tensor, v p′ -temperature coefficient of the evaporation pressure.
The kinematic boundary conditions and those for the thermal conductivity equation are formulated as follows: A system of the Navier-Stokes equation and convective thermal conductivity equation is linearized and written as follows: where u, w -disturbance vector components of velocity, ρ-density, ν-kinematic viscosity, χ-temperature conductivity, G 0 -temperature gradient, p, T, ηpressure and temperature perturbations. No elements in the Equation (1) are related to the electric field, which is acceptable for conductive liquids [25]. Here, an electric field effect is included into boundary conditions on the surface z = 0: where η -perturbation of the surface along the axis z, σ = σ 0 + σ T (T − T m )-surface tension, σ 0 -surface tension at a melting temperature of the material T m , σ T -temperature coefficient of the surface tension, p e -electric field-melt pressure, τ el xz -tangential component of Maxwell stress tensor, p v -temperature coefficient of the evaporation pressure.
The kinematic boundary conditions and those for the thermal conductivity equation are formulated as follows: Electric elements of the boundary conditions (2) are determined on the following grounds [25,26]: electric field strength is E 0 on the stable horizontal surface of a liquid, the electric potential therefore ϕ 0 = −E 0 z. If a perturbation arises, the electric potential takes the form: ϕ = ϕ 0 + ϕ 1 , where ϕ 1 -small potential perturbation satisfying the Laplace equation ∆ϕ 1 = 0 on the condition of z → +∞, ϕ 1 → 0. This perturbation has a form of the normal mode with an amplitude damping on the axis z . Given that the electric potential along the wavelength is 0, we obtain −E 0 η + φ 1 (z = η) = 0, i.e., φ 1 (z = 0) = E 0 η, and the electric field pressure at E z = E 0 + kφ 1 takes the form: p e = εε 0 2 E 2 0 + 2kE 2 0 η + k 2 E 2 0 η 2 , where ε-dielectric permittivity, ε 0 -electric constant. Omitting elements of the second order of smallness, p e = εε 0 kE 2 0 η is obtained. A tangential component in the Maxwell stress tensor for z = η is written τ el In this case boundary conditions (2) are formulated: (1) is sought in the form: . Substituting (5) into boundary conditions (3) and (4), a system of linear algebraic equations in constants A 1 and A 2 is obtained: where The temperature gradient is calculated as follows: where q = E s t 0 -power density of an electron beam, E s -surface power density of an electron beam, q out -evaporation-conditioned surface power density, κ-thermal conductivity of a liquid metal. The role of evaporation for the power density [27] is calculated as follows: where A out -work function of an atom escaping from the melt, k B -Boltzmann constant, m-atom mass, T V -evaporation temperature.
Metals 2020, 10, 1399 6 of 14 The existence of a nontrivial solution to the system (6) needs its determinant to be zero, which in turn is a dispersion equation: The Prandtl number for liquid metals takes a value Pr << 1. In this case the dispersion Equation (9) is formulated: Substituting z = k 1 k and ω = ω ν (z 2 − 1), an Equation (10) is written in the form: where C = ω c ω ν , C 2 = ω T Pr 2ω ν , C 3 = ω p Pr 2ω ν , C 4 = ω E 2ω ν . Instable solutions comply with constraints Re(ω) > 0 and Re(z) > 0. With the help of these solutions a relationship between the growth rate of perturbations α = Re(ω) and wavenumber (wavelength) may be established. Using this correlation, we determine a wavelength, when the instability (α = 0) arises, and a value λ for α = α m . The instability of the plasma-melt interface is explored on the example of titanium and aluminum alloys. The data for calculations are provided in Table 1. The thermoelectric coefficient for a temperature above the Debye one is estimated by the formula γ = f k B T D eT [31], where T D -Debye temperature, e-charge carrier, f -constant with values ensuring γ 100 µV/K.
The correlation between the growth rate of perturbations on the interface and wavelength is estimated for an electron beam power density of 35 J/cm 2 and a pulse time of 150 µs.

Results and Discussion
First, we consider a situation, when thermoelectric effects and evaporation pressure may be omitted. Here, a gradient of nonperturbed temperature G 0 = 3.66 ·10 7 K/m-start of the instability-occurs in titanium alloys when the wavelength is 38 µm, and a maximal growth rate is registered for a wavelength of 113 µm (Figure 3 curve 1). A value of the thermoelectric field for γ~100 µV/K ranges up to E 0~1 0 3 V/m, in case of such a field value the thermoelectric effect is almost irrelevant to a growth rate of perturbations on the interface. The reason is that the thermoelectric coefficient is calculated for a solid in the pre-melting state. In a liquid substance charge, carrying is highly dependent on the convective flow, which intensifies the thermoelectric effect. Therefore, the value γ is~10 4 -10 5 µV/K; and electric field strength E 0~1 0 5 -10 6 V/m, respectively. For such field values the growth rate of perturbations increases significantly (Figure 3 curve 2), whereas the critical wavelength is two times shorter (17 µm). A maximal growth rate is registered for λ = 50 µm. For an evaporation pressure of 2 ·10 5 Pa the critical wavelength increases to 22 µm with a maximal growth rate for a wavelength of 92 µm (Figure 3 curve 3). In this case the thermoelectric effect makes no contribution to instability on the interface (Figure 3 curve 4). In aluminum alloys the same trend is observed (Figure 3). Ignoring evaporation pressure and thermoelectric effects for a nonperturbed temperature gradient G 0 = 1.81 × 10 7 K/m, the critical wavelength is 39 µm, and a maximal growth rate of perturbations is attributable to a wavelength of 155 µm (Figure 4 curve 1). With thermoelectric effects, we obtain λ cr = 13 µm and λ m = 53 µm (Figure 4 curve 2). For an evaporation pressure of 2 × 10 5 Pa-λ cr = 17 µm and λ m = 87 µm. The thermoelectric effect shifts the maximal growth rate to the wavelength λ m = 76 µm, and the critical wavelength is λ m = 18 µm in these conditions (Figure 4 curve 3 and 4). From the experimental data it can be seen [7,8] that cell structures develop due to the electron beam irradiation varying in a range from 100 to 600 nm, depending on processing conditions. The low frequency approximation (10) has shown this range of wavelengths is possible for electric field strength E 0~1 0 6 -10 7 V/m (Figure 5 curve 1 and 2). In titanium alloys λ cr = 0.17 µm and λ m = 0.3 µm, and in aluminum alloys-λ cr = 0.13 µm and λ m = 0.26 µm.
In the calculations presented the effect of a doping element on the melt flow is ignored. In fact, the occurrence of doping elements is essential to the melt surface tension and to the thermocapillary instability as well [32]. Depending on the nature of a base material, a doping element may be a surface-active substance, which reduces (surface-inactive substance) or increases the surface tension of the melt or have no effect on the surface tension of a material. This study investigates the influence of an yttrium concentration gradient on the thermocapillary flow of a titanium or aluminum melt. From the experimental data on the surface tension it is apparent [33] that yttrium is a surface-active substance as related to titanium, whereas it is a surface-inactive substance to aluminum. Including the influence of a doping element concentration, the surface tension vs. temperature and concentration relationship is formulated as follows: where σ C -concentration coefficient of the surface tension, C-surface concentration. A convective diffusion equation needs to be added to the system (1): where D-diffusion coefficient, G 1 -concentration gradient. The boundary conditions for tangential stresses are written: Metals 2020, 10, x FOR PEER REVIEW 8 of 15     In the calculations presented the effect of a doping element on the melt flow is ignored. In fact, the occurrence of doping elements is essential to the melt surface tension and to the thermocapillary instability as well [32]. Depending on the nature of a base material, a doping element may be a surface-active substance, which reduces (surface-inactive substance) or increases the surface tension of the melt or have no effect on the surface tension of a material. This study investigates the influence of an yttrium concentration gradient on the thermocapillary flow of a titanium or aluminum melt. From the experimental data on the surface tension it is apparent [33] that yttrium is a surface-active substance as related to titanium, whereas it is a surface-inactive substance to aluminum. Including the influence of a doping element concentration, the surface tension vs. temperature and concentration relationship is formulated as follows: where σC-concentration coefficient of the surface tension, C-surface concentration. A convective diffusion equation needs to be added to the system (1): where D-diffusion coefficient, G1-concentration gradient. The boundary conditions for tangential stresses are written: To consider the effect of a doping element concentration in the low frequency approximation, a product Pr . The use of the low frequency approximation is reasonable provided that Pr << 1 and Sc << 1. For liquid metals the Schmidt number is ~10 [34]. It was assessed under the assumption of a stable liquid. In practice, the convective flow intensifies mass transfer processes; thus, the diffusion coefficient is accepted to be one order of magnitude higher than the kinematic viscosity coefficient. Figure 5 gives curves of a growth rate on the melt interface of yttrium-doped titanium for a surface concentration gradient of G1 = 10 6 m −2 /m. To consider the effect of a doping element concentration in the low frequency approximation, a product ω T Pr is to be replaced by ω T Pr + ω C Sc in the dispersion Equation (10), where Sc = ν/D -the Schmidt number, ω C = σ C G 1 ρν . The use of the low frequency approximation is reasonable provided that Pr << 1 and Sc << 1. For liquid metals the Schmidt number is~10 [34]. It was assessed under the assumption of a stable liquid. In practice, the convective flow intensifies mass transfer processes; thus, the diffusion coefficient is accepted to be one order of magnitude higher than the kinematic viscosity coefficient. Figure 5 gives curves of a growth rate on the melt interface of yttrium-doped titanium for a surface concentration gradient of G 1 = 10 6 m −2 /m.
As can be seen, if thermoelectric and evaporation effects are ignored, the critical wavelength is 26 µm, and a maximal growth rate aligns with a wavelength of 74 µm ( Figure 6, curve 1). The consideration of thermoelectric effects results in a critical wavelength of 17 µm with a maximal growth rate conforming to a wavelength of 48 µm ( Figure 6, curve 2). For an evaporation pressure of 2 × 10 5 Pa λ cr = 11 µm and λ m = 36 µm thermoelectric effects cause a slightly increasing maximal growth rate of perturbations ( Figure 6, curves 3 and 4). In aluminum alloys another phenomenon is observed for a concentration gradient of G 1 = 10 6 m −2 /m. If ω E =0 and ω p =0, the critical wavelength is 7 µm, and the maximal growth rate corresponds to 25 µm (Figure 7, curve 1). For ω E 0 and ω p =0 the critical wavelength is 7 µm, and the growth rate maximum is detected for wavelengths of 10 and 44 µm (Figure 7, curve 2). For ω E = 0 and ω p 0 the critical wavelength is 7 µm, and a maximal growth rate is attained for a wavelength of 27 µm (Figure 7, curve 3). For ω E 0 and ω p 0 the growth rate function also displays two maximums (Figure 7, curve 4). The first maximum corresponds with a wavelength of 10 µm and the second one-110 µm. The occurrence of the first maximum might result from the interaction of concentration and capillary and thermoelectric instabilities, the second maximum might be caused by a temperature gradient. The results obtained explain adequately the formation of surface-periodic structures in a micrometer range. Within the low frequency approximation, similarly to the above considered situation, submicro-and nanodimensional cell structures may develop for a thermoelectric coefficient~10 −1 V/K or an evaporation pressure~10 11 Pa. Thus, the inclusion of the yttrium concentration makes a spectrum of capillary waves in the Al-Si-Y system more complex, whereas this effect is not detected in the Ti-Y system. maximum might be caused by a temperature gradient. The results obtained explain adequately the formation of surface-periodic structures in a micrometer range. Within the low frequency approximation, similarly to the above considered situation, submicro-and nanodimensional cell structures may develop for a thermoelectric coefficient ~10 −1 V/K or an evaporation pressure ~10 11 Pa. Thus, the inclusion of the yttrium concentration makes a spectrum of capillary waves in the Al-Si-Y system more complex, whereas this effect is not detected in the Ti-Y system.  An analysis of the dispersion Equation (10) has revealed that the low frequency approximation even with included thermoelectric and evaporation capillary effects cannot provide an adequate explanation why a micro-and nanodimensional cell crystallization structure forms; we seek therefore for functions of the perturbations growth rate on the base of solutions to the dispersion Equation (9). For this purpose, we use the substitution k z k k z k  An analysis of the dispersion Equation (10) has revealed that the low frequency approximation even with included thermoelectric and evaporation capillary effects cannot provide an adequate explanation why a micro-and nanodimensional cell crystallization structure forms; we seek therefore for functions of the perturbations growth rate on the base of solutions to the dispersion Equation (9). For this purpose, we use the substitution k 1 = z 1 k, k 2 = z 2 k and ω = ω ν (z 2 − 1), z 2 = (1 + Pr(z 2 1 − 1)) 1/2 . This substitution transforms the Equation (9) into an algebraic equation of the 16th degree in z 1 , which is too lengthy to be given here. Roots of this equation satisfying the constraints Re(ω) > 0, Re(z 1 ) > 0, Re(z 2 ) > 0 are instable. Figure 8 provides the data of how the growth rate of perturbations depends on a wavelength in the Ti-Y system. It is apparent that without thermoelectric and evaporation capillary effects there are two maximums on this curve. The first maximum is found for a wavelength of 18 μm, and the second one appears if the wavelength is 294 μm (Figure 8a). If a thermoelectric field ~10 5 V/m arises, the first maximum is detected for a wavelength of 8 μm, and the second one 43 μm (Figure 8b). An increase 6 Figure 8. A growth rate characteristic of perturbations on the plasma-melt interface for the Ti-Y system found when solving the Equation (9): (a) ignoring thermoelectric and evaporation and capillary effects, (b) for thermoelectric field strength of 10 5 V/m, (c) for thermoelectric field strength of 10 6 V/m. It is apparent that without thermoelectric and evaporation capillary effects there are two maximums on this curve. The first maximum is found for a wavelength of 18 µm, and the second one appears if the wavelength is 294 µm (Figure 8a). If a thermoelectric field~10 5 V/m arises, the first maximum is detected for a wavelength of 8 µm, and the second one 43 µm (Figure 8b). An increase in thermoelectric field strength up to 10 6 V/m eliminates the second maximum, and makes the first one conform to a wavelength of 0.31 µm (Figure 8c). In the Al-Si-Y system we find similar transformations in the growth rate vs. wavelength relationship (Figure 9). Without thermoelectric effects the growth rate maximum is for a wavelength of 182 μm (Figure  9a), and two maximums are detected for wavelengths of 4 μm and 48 μm, given that field strength is 10 5 V/m (Figure 9b). For E0 ~ 10 6 V/m-λm = 0.26 μm (Figure 9c). In the Ti-Y system the maximum is shifted to a wavelength of 9.6 μm for an evaporation pressure of 2 × 10 5 Pa, ignoring thermoelectric Without thermoelectric effects the growth rate maximum is for a wavelength of 182 µm (Figure 9a), and two maximums are detected for wavelengths of 4 µm and 48 µm, given that field strength is 10 5 V/m (Figure 9b). For E 0~1 0 6 V/m-λ m = 0.26 µm (Figure 9c). In the Ti-Y system the maximum is shifted to a wavelength of 9.6 µm for an evaporation pressure of 2 × 10 5 Pa, ignoring thermoelectric effects. Given that ω E 0 and ω p 0, a maximal growth rate is registered for wavelengths of 4.2 µm and 19.6 µm. In the Al-Si-Y system similar transformations are observed.

Conclusions
The detected correlations between the thermoelectric effect and early instability stage of the thermocapillary melt flow in the electron beam irradiation suggest that surface micro-and nanodimensional structures form due to the combination of thermal, concentration, evaporation and capillary and thermoelectric instabilities. The findings provide some support for the significant importance of thermoelectric and concentration and capillary phenomena to the growth rate of perturbations on the plasma-melt interface. The low frequency approximation in the Al-Si-Y system (9) demonstrates two maximums of this correlation, whereas only one maximum is detected in the Ti-Y system. A growth rate maximum is attributed to submicro-and nanometer wavelengths for a thermoelectric coefficient~10 −1 V/K both in Ti-Y and Al-Si-Y system. The results obtained may be useful when searching for appropriate conditions of the electron beam irradiation to further the formation of micro-and nanostructures.