Infrasonic Nanocrystal Formation in Amorphous NiTi Film: Physical Mechanism, Reasons and Conditions

: The physical mechanism, reasons and conditions of nanocrystal formation in an amorphous NiTi metal ﬁlm, stimulated by infrasonic action, are formulated. Nanostructural elements of an amorphous medium (relaxation centers) containing disordered nanoregions with two-level systems are considered to be responsible for this process. When exposed to infrasound, a large number of two-level systems are excited, signiﬁcantly contributing to inelastic deformation and the formation of nanocrystals. The physical mechanism of the nanocrystallization of metallic glass under mechanical action includes both local thermal ﬂuctuations and the additional quantum tunneling of atoms stimulated by shear deformation. A crystalline nanocluster appears as a result of local atomic rearrangement growing increasingly exposed to infrasound. It is possibly unstable in the absence of infrasound. When the radius of the nanocluster reaches a critical value, a potential well appears, in which a conducting electron is localized to form a phason (stable nanocrystal). Estimated values of the phason’s radius and the depth of the nanometer potential well is about 0.5 nm and 1 eV, respectively. It forms a condition of stable phason formation. The occurrence of the instability of the amorphous state and following transformation to the nanostructured state is based on the accumulation of the potential energy of inelastic deformation to a critical value equal to the latent heat of the transformation of the amorphous state into the nanostructured state.


Introduction
Amorphous metal alloys with a nanostructure are possessed by unique mechanical properties occurring due to intensive mechanical processing [1]. For nanostructuring an amorphous NiTi film, low-frequency mechanical vibrations are used at a temperature much lower than the glass transition temperature used for the remaining sample deformation within the elastic region [2][3][4]. Data obtained by the methods of X-ray diffraction and highresolution transmission electron microscopy indicate that the initial alloy is homogeneously amorphous, while the structure of the sample subjected to a mechanical vibration amplitude of 1 mkm is close to the structure of the initially untreated sample (the short-range order changes and free volume decreases). In the structure of the sample subjected to mechanical vibrations with an amplitude of 4 mkm, nonspherical clusters with correct atomic positions and sizes of 3-5 nm appear in an amorphous matrix.
It has been established [2] that, under mechanical action, the latent heat of crystallization decreases slightly. It means that the free energy of the state with nanocrystals becomes less in comparison with the free energy of the amorphous state wherein the chemical and phase composition of alloys with and without nanocrystals is the same. Therefore, nanocrystallization caused by infrasonic stimulation is not related to diffusion, but provoked by collective atomic rearrangements during structural relaxation.
Further experimental studies have shown that [3] an increase in the duration of oscillations up to 2 h enlarges the size of clusters and their growth in different directions.
Running on mechanical vibrations for 4 h, the film acquires a crystalline structure. It consists of irregular shape grains of different crystallographic orientations with amorphous islands embedded.
In [4], the influence of temperature, frequency and amplitude of vibrations (the duration of vibrations is 10 min) on the structure of amorphous alloys based on NiTi was investigated. It was found that, at room temperature, a vibration frequency of 20 Hz, a vibration amplitude of 4 µm (10 −6 m), and crystalline clusters 4-6 nm (10 −9 m) in size appear in the amorphous matrix. With an increase in the initial temperature of the sample (from 25 to 200 • C), the distribution and size of crystalline clusters practically do not change. The crystallization temperature practically does not change, and the latent heat of crystallization decreases with an increase in the vibration amplitude and the initial temperature of the sample. When the frequency (below 5 Hz) or the vibration amplitude (down to 1 µm) decreases, crystalline clusters do not appear in the amorphous matrix.
Despite the fact that the phenomenon of nanocrystal formation in amorphous NiTi film under the action of mechanical vibrations is described in the literature, the physical mechanisms, reasons and conditions of its occurrence remain unclear.
In accordance with classical concepts [5], the nucleation of crystals in amorphous alloys can occur due to either a homogeneous or heterogeneous mechanism. Approaching the glass transition temperature, the homogeneous nucleation occurs through the fluctuation nucleus formation with a radius larger than the critical one [5]. Heterogeneous nucleation occurs mainly at temperatures below the glass transition temperature placed on defects, boundaries, etc., [5,6] and due to "frozen in" crystallization centers [7]. The density of heterogeneous nuclei in metal melts is about 10 11 − 10 12 m −3 . Amorphous alloys, without nanosized crystallites, have a density noticeably exceeding these values (for example, in the Fe 80 B 20 alloy, the density of crystallites is 10 18 m −3 [8]). In nanostructured systems (for example, in Cu-Ti alloys [9]), the density of crystallites can reach 10 25 m -3 ; therefore, most of the nuclei arise through homogeneous nucleation during annealing.
As analysis shows, the mechanisms of the nanocrystal formation in amorphous NiTi film under stationary annealing and mechanical action differ significantly. During thermal annealing, the sample temperature is much higher than room temperature. The nanocrystallization process is slow and homogeneous. Nanocrystals appear throughout the entire volume by the thermal fluctuation mechanism. Under mechanical action, the medium temperature is room temperature; therefore, the probability of nanocrystal formation by the usual thermal fluctuation mechanism is much less than at the annealing temperature. However, during mechanical vibrations with an amplitude of 4 mkm, significant shear stresses arise in the film. At a critical value of shear stress, the medium loses shear stability, and inelastic deformation occurs. Therefore, the mechanism of nanocrystallization of metallic glass under mechanical action should include quantum tunneling of atoms or atomic groups [10][11][12][13][14] stimulated by shear deformation in addition to local thermal fluctuations.
The ability of amorphous metal alloys to cause irreversible deformation and nanocrystal formation is associated with the collectivized character of the metallic interatomic bond [15]. Irreversible processes of local collective atomic rearrangements can occur much more easily in the presence of conducting electrons. Therefore, when constructing the mechanism of the irreversible deformation and nanocrystallization of amorphous alloys, it is necessary to take into account the influence of the electronic subsystem [16][17][18]. Specifically, for metallic glass NiTi, it is necessary to consider the localization of a conducting electron in a nanometer potential well and the formation of a phason [18].
The study of nanocrystal formation under a different external surrounding is of interest to the modern science of solid state matter. At the same time, as we found in the literature, there is no consideration of the effect of infrasound nanocrystal formation, and its physical mechanisms and possible conditions are not established.
The present work aims to consider the physical mechanisms of nanocrystal formation in amorphous NiTi under the infrasound action. The mechanism of influence of infrasound is not known and has not been previously considered. Herein, we propose a model that takes into account the inhomogeneity of the amorphous medium, the so-called relaxation centers, capable of rearranging according to the mechanism of athermal shear deformations stimulated by infrasound. As a result of the atomic rearrangement of the relaxation center, a more ordered structure is formed. As a consequence, a potential well arises, in which a conductivity electron is localized.
As we supposed, the analytical formulation of the problem is of two objectives: (1) in the framework of the model to formulate the equations of the density operator, taking into account the relaxation term and the field of static shear deformation; (2) to determine the probabilities of transitions of an atom in a double-well potential by thermal and athermal, i.e., quantum, mechanisms.
At the same time, numerical estimations might be provided in order to make the process of crystallization predictable. (3) It is likely to know the value of the shear strain, epsilon, at which quantum tunneling (disordered nanocluster transforms to a nanocrystal) occurs according to the athermal type; (4) It is necessary to determine conditions for the appearance of a localized state of conductivity electron at the potential well of the nanocrystal. There is a need in quantitative estimations of the critical values of the potential well depth and the radius of a stable nanocrystal.

Physical Mechanism of Inelastic Deformation of Metal Glasses
Inelastic deformation occurring as a response to the applied shear stress [19][20][21][22][23][24] is explained by local static displacements of a group of atoms (or an atom) from the initial equilibrium positions to the new ones at a distance less than interatomic spacing. The first microscopic model of elementary cooperative inelastic rearrangements in metallic glasses is based on the concept of a uniformly distributed free volume. It was proposed by Argon in 1979 [19]. In the model [20], it is assumed that relaxation centers appear due to the inhomogeneity of the glass structure, particularly in the regions with an excess free volume relative to the "ideal structure". At present, the model of the activation energies spectrum [21], as well as the model of directed structural relaxation oriented by an external force [22], is often used.
The model [23] proposes a method of description of the local structure of amorphous alloys on the basis of the concept of n, p, and τ defects. In [24], a polycluster model of the structure of an amorphous state was proposed. The elementary rearrangements of atomic configurations occur at intercluster boundaries under the action of external forces. At high temperatures, the rearrangements of configurations are, possibly, a mechanism of thermal fluctuations, whereas at low temperatures quantum tunneling is more physically suitable.
The structural elements of an amorphous medium containing free volume (relaxation centers [20]) rearranging not only as a result of thermal fluctuations but due to quantum tunneling are of particular interest. Atom moving in a double-well potential ( Figure 1) or a group of atoms with two configurations possessing different energies considered as twolevel systems [10][11][12][13][14]24]. The atom transition from one well to another (or the changing of one configuration of atoms to another), which means a transition from one energy level to another, can occur due to quantum tunneling. These effects are more probable when a two-level system is close to degenerating. Specifically, the levels of zero-point vibrations of an atom, in each of the neighboring wells, should differ slightly without taking into account tunneling, (i.e., there is a resonance detuning ∆ 0 ). A two-level system is characterized by an initial resonance detuning ∆ 0 << 1 eV and an initial frequency of tunnel junctions ω 0 K = ω 0 exp (− a a dB ). Here, a dB is the de Broglie wavelength and a is the width of the potential barrier in the double-well potential, ω 0 ∼ 10 13 Hz. The frequency, ω 0 K , on the one hand, must be much lower ∆ 0 in order for the initial splitting of the levels E 0 = ∆ 2 0 + 2 (ω 0 K ) 2 to be small~∆ 0 . On the other hand, ω 0 K should not be too small for the tunneling to manifest itself during experiment.
ting of the levels 0 = √∆ 0 2 + ℏ 2 ( 0 ) 2 to be small ~Δ0. On the other hand, 0 should not be too small for the tunneling to manifest itself during experiment. In the framework of the one-particle model of the interaction of phonons with two-level systems [10][11][12][13][14]24], the initial resonance detuning Δ0 and the initial frequency of tunneling transitions 0 are renormalized. During the passage of a sound wave, the atoms that form the potential curve of the two-level system are displaced. As a consequence, individual wells of the curve are distorted, which leads to a change in the resonance detuning ∆. The shape of the potential barrier also changes and shifts to the frequency of tunnel junctions . The distance between the two levels of the two-level system becomes to be equal to = √∆ 2 + ℏ 2 2 . For the linear approximation of deformation, we have Here, , are tensors of strain potentials, which depend on the local structure of the glass in the place where the two-level system is located.
-strain tensor, which is expressed through the components of the displacement vector The energy per atom per one deformation cycle is ≅ 10 −6 eV; therefore, ∆≅ ≅ 10 −6 eV. Since at the initial moment of time Δ 0 < 0, it is possible to estimate |Δ 0 | ≅ ∆≅ 10 −6 eV and the maximum value of ≅ 10 9 s −1 . To estimate the characteristic parameters of the atomic tunneling, we will perform a modeling calculation. Let us suppose an atom in the relaxation center is in a two-well potential or there is a group of atoms in which there are two configurations with slightly different energies [12][13][14]24]. The potential energy of such a subsystem can be represented as the sum of two single-well potentials and . At the initial moment of time (in the absence of mechanical load), the left and right potential holes have different depths-the left hole is deeper than the right one. The action of an external mechanical force on the system leads to a change in the distances between the atoms of the medium (deformation), which in turn changes the shape of the double-well potential.
Let the wave functions , and energy eigenvalues 0 , 0 be known for the  In the framework of the one-particle model of the interaction of phonons with twolevel systems [10][11][12][13][14]24], the initial resonance detuning ∆ 0 and the initial frequency of tunneling transitions ω 0 K are renormalized. During the passage of a sound wave, the atoms that form the potential curve of the two-level system are displaced. As a consequence, individual wells of the curve are distorted, which leads to a change in the resonance detuning ∆. The shape of the potential barrier also changes and shifts to the frequency of tunnel junctions ω K . The distance between the two levels of the two-level system becomes to be equal to E = ∆ 2 + 2 ω 2 K . For the linear approximation of deformation, we have Here, B ij , D ij are tensors of strain potentials, which depend on the local structure of the glass in the place where the two-level system is located. ε ij -strain tensor, which is expressed through the components of the displacement vector u i The energy per atom per one deformation cycle is e a ∼ = 10 −6 eV; therefore, δ∆ ∼ = e a ∼ = 10 −6 eV. Since at the initial moment of time ∆ 0 < 0, it is possible to estimate |∆ 0 | ∼ = δ∆ ∼ = 10 −6 eV and the maximum value of ω K ∼ = 10 9 s −1 .
To estimate the characteristic parameters of the atomic tunneling, we will perform a modeling calculation. Let us suppose an atom in the relaxation center is in a two-well potential or there is a group of atoms in which there are two configurations with slightly different energies [12][13][14]24]. The potential energy of such a subsystem can be represented as the sum of two single-well potentials U L and U R . At the initial moment of time (in the absence of mechanical load), the left and right potential holes have different depths-the left hole is deeper than the right one. The action of an external mechanical force on the system leads to a change in the distances between the atoms of the medium (deformation), which in turn changes the shape of the double-well potential.
Let the wave functions Ψ L , Ψ R and energy eigenvalues E 0 L , E 0 R be known for the Schrödinger equation with single-well potentials U L , U R . For simplicity, when superimposing the field of the static shear component of deformation ε, only the energy eigenvalues The complete Hamiltonian of a one-dimensional subsystem with a two-well potential U L (ε) + U R (ε), U L (ε) = U L + Bε, U R (ε) = U R − Bε in the field of a static shear strain component ε has the form. (2) Due to the assumed weak overlap of wave functions Ψ L , Ψ R , (S = dxΨ L Ψ R 1), solving the Schrödinger equation with Hamiltonian (2), we find the stationary levels E ± ( Figure 2) and stationary states Ψ ± , in which an atom is delocalized between two wells: Due to the assumed weak overlap of wave functions , , ( = ∫ ≪ 1) solving the Schrödinger equation with Hamiltonian (2), we find the stationary levels ± (Figure 2) and stationary states ± , in which an atom is delocalized between two wells: Let us investigate the tunneling transition of an atom from one potential well to an other in the field of the static shear component of deformation at a finite temperature using the formalism of the density operator ρ [25]. In this case, the atom making transi tions between the wells interacts with the thermostat. Using eigenvectors , as ba sis vectors, leaving only the first-order terms in the parameter , we obtain differentia equations for the matrix elements of the density operator ρ: The diagonal matrix element has the meaning of the probability of finding an atom in the left well. For simplicity, the system under study is characterized by two relaxation times: the longitudinal relaxation time of the system = with an atom in the lef (right) potential well and = the transverse relaxation time of off-diagonal matrix elements.
When the time of the tunneling transition of an atom significantly exceeds the time of longitudinal or transverse relaxation, quantum coherence is violated. The probability of detecting a particle in a certain state can be found from the solution of the classica kinetic equations obtained by reducing the system Equations (5)-(7). The relaxation time has an exponential dependence on temperature 0 −1 exp ( ), where is the height of the interwell barrier. When the temperature is high (the relaxation time is short) and the shear deformation is less than the critical one < * , * = − 0 − 0 2 , the two-level system Let us investigate the tunneling transition of an atom from one potential well to another in the field of the static shear component of deformation ε at a finite temperature using the formalism of the density operator ρ [25]. In this case, the atom making transitions between the wells interacts with the thermostat. Using eigenvectors Ψ L , Ψ R as basis vectors, leaving only the first-order terms in the parameter S, we obtain differential equations for the matrix elements of the density operator ρ: The diagonal matrix element ρ LL has the meaning of the probability of finding an atom in the left well. For simplicity, the system under study is characterized by two relaxation times: the longitudinal relaxation time of the system τ LL = τ RR with an atom in the left (right) potential well and τ LR = τ RL the transverse relaxation time of off-diagonal matrix elements.
When the time of the tunneling transition of an atom significantly exceeds the time of longitudinal or transverse relaxation, quantum coherence is violated. The probability of detecting a particle in a certain state can be found from the solution of the classical kinetic equations obtained by reducing the system Equations (5)-(7). The relaxation time has an exponential dependence on temperature ω −1 0 exp ( V kT ), where V is the height of the interwell barrier. When the temperature is high (the relaxation time is short) and the shear deformation is less than the critical one ε < ε * , ε * = − E 0 L −E 0 R 2B , the two-level system relaxes into a thermodynamically equilibrium state with ρ 0 LL > 1 2 , ρ 0 RR < 1 2 , in which the atom is in the left well with a higher probability than in the right well. When the deformation is equal to the critical ε = ε * , the two-level system relaxes into the equilibrium state with ρ 0 LL = 1 2 , ρ 0 RR = 1 2 , in which the atom is with equal probability in the left and right wells. When the deformation is greater than the critical one ε > ε * , the two-level system relaxes to the equilibrium state with ρ 0 LL < 1 2 , ρ 0 RR > 1 2 , in which the atom is in the right well with a greater probability than in the left well.
If the period of quantum oscillations is much shorter than the time of longitudinal and transverse relaxation, then the process of tunneling transition is coherent. The tunnel transition probability is an oscillating function of time, i.e., described by quantum kinetics. A linear combination of stationary wave functions Ψ(x, t), satisfying the time-dependent Schrödinger equation with Hamiltonian (2) and the initial condition Ψ(t = 0) = Ψ L , makes it possible to find the probability of finding an atom in the left well: The function W(t) oscillates from unity to the minimum value, which is determined by the difference in energies in two states ∆ and the magnitude of the product U 2 . If ∆ = 0, then W(t) = 1 − sin 2 ( U t). The probability of finding an atom in the left well varies from one to zero with a period τ = π U , and the atom spends the same time in both wells. If ∆ 2 4U 2 , The atom is in the left potential well almost all of the time. When the shear strain ε is zero, the case ∆ = ∆ 0 = E 0 L − E 0 R < 0, ∆ 2 4U 2 is realized, and the atom is only in the left well. At a critical value of deformation = ε * , the case ∆ = 0 is reached and the atom tunnels into the right well. Further increase in deformation leads to the case ∆ > 0 ∆ 2 4U 2 , and the atom is only in the right well. The static displacement of an atom is an elementary act of structural relaxation of a medium under load and is accompanied by stress relaxation. This leads to a change in the positions of atoms from the nearest environment and "freezing" of the changed shape of the double-well potential. This process is reversible; after a slow removal of the load, the atom returns to the left well again. This is how the elementary inelastic reversible deformation proceeds in metallic glass at T = 0.
Metallic glass contains a large number of relaxation centers containing two-level systems, which are described by a double-well potential, resonance detuning ∆, and tunnel transition frequency ω K . Two groups of such double-well potentials can be distinguished-soft (with an interwell barrier V ∼ 10 −4 − 10 −2 eV) and hard (with V ∼ 10 −1 − 1 eV) [12][13][14]24]. In the temperature range T ∼ = 1 − 300 K, thermal jumps between the states of a two-level system with a barrier V ∼ 10 −4 − 2.4·10 −2 eV are activated. Therefore, the quantum tunneling mechanism is more probable for a two-level system with a hard potential.
Let us estimate the temperature T * at which the frequency of the quantum sub-barrier tunneling of an atom ω K becomes equal to the frequency of the above-barrier fluctuation hopping ω T . Using the Gamow formula and Arrhenius's law, equating the exponents, we get From (9), it can be seen that T * increases √ V proportionally and decreases 1 m , 1 a proportionally. For a two-level system with a hard potential, the characteristic values are V ∼ = 1 eV, m ∼ = 10 −25 kg, and a dB ∼ = 0.56·10 −12 m. With a barrier width a = 3·10 −11 , 1.5·10 −11 , 0.9·10 −11 m (0.1, 0.05, 0.03 interatomic distance), we obtain, respectively, T * ∼ = 65, 216, 721 K and the quantum tunneling frequency ω K ∼ = 10 −10 20, 100, 000 Hz. For local structural rearrangements with small displacements of atoms (much less than the interatomic distance), the tunneling effect is significant for room temperatures, and becomes the main at lower temperatures. Therefore, a possible physical mechanism of the inelastic deformation of metallic glasses (at room temperatures and below) is the quantum tunneling of some atoms in a double-well potential, or within an atomic group.
Under mechanical action on an amorphous material, a large number of two-level systems are excited, making a significant contribution to inelastic deformation and the formation of nanocrystals. The ability of amorphous metal alloys to deform irreversibly and to form nanocrystals is associated with the collectivized metallic character of the interatomic bond. In the presence of conduction electrons, irreversible processes of local collective atomic rearrangements can be carried out much more easily.

Electron Localized States in Potential Nanometer Well and Phason Formation
Let us assume that the electron moves in a three-dimensional spherically symmetric potential well having experienced the action of central force. In the simplest case, the potential energy U(r) depends only on the distance between the electron and the force center r. Let the three-dimensional potential well have the form U(r) = −U 0 at r ≤ a and be equal to zero at > a. In a spherical coordinate system, the Schrödinger equation for stationary states has the form [26] 1 As solutions to Equation (10), which are finite at r = 0 and vanish at r → ∞ , we can take: Since we consider the particle inside a potential well, it is necessary to count < 0. Thus, the problem of Equations (10) and (11) has been reduced to the problem of the motion of an electron in a one-dimensional potential well. Therefore, the energy levels are determined in the same way. The only difference now is in the necessity to discard states with even wave functions and leave only states with odd wave functions. The fundamental difference between a one-dimensional potential well and a threedimensional one is that for a one-dimensional well there is at least one eigenvalue of energy with even wave function. In the case of a spherically symmetric rectangular well, this may not be the case. It can be seen from the relation k 2 + α 2 = 2mU 0 2 if 2mU 0 2 a 2 < ( π 2 ) 2 , i.e., U 0 < π 2 2 8ma 2 (12) that the curve given by the equation tg( ka 2 ) = − k α will never intersect the circle given by the equation k 2 + α 2 = 2mU 0 2 . This means that when condition (12) is satisfied, not a single level of the discrete energy spectrum will appear in the potential well (due to the fact/the depth of the potential is not enough/the power of the well is too low).
Substituting the values of the fundamental constants ∼ = 1·10 −34 J·s, the electron mass m ∼ = 9·10 −31 kg, for the radii of the first, second, and third coordination sphere a 1 = 0.35·10 −9 m, a 2 = 0.7·10 −9 m, and a 3 = 1.05·10 −9 m, we obtain U 0 (a 1 ) = 0.7 eV, U 0 (a 2 ) = 0.18 eV, and U 0 (a 3 ) = 0.08 eV, respectively. From physical considerations, it follows that the perfect local order can propagate starting from the first coordination sphere. The minimum radius of the potential well is equal to the interatomic distance a 1 = 0.35·10 −9 m. For a localized state of an electron to appear in a potential well with a minimum radius, the depth of the potential well of the resulting nanocrystal should be approximately 1 eV.
At a high interfacial surface energy inherent in solid-phase transformation, the fluctuation formation of a region of a new crystalline phase leads to a significant change in the thermodynamic potential and may turn out to be unstable (nonequilibrium). However, in the presence of the attraction of a conduction electron in the formed region of a new crystalline phase, the electron can be localized in there. In this case, during the formation of a new region, if the decrease in the electron energy exceeds the growth of the thermodynamic potential energy of a new crystalline phase, then the fluctuation is stabilized. The arising thermodynamically stable formation of the new phase region with a localized electron in it (nonequilibrium in the absence of an electron) is called a phason [18].
In Section 1 of this work, it is shown that the formation of the nanoregion of the crystalline phase in an amorphous metal system can also occur by quantum tunneling stimulated by shear deformation. Therefore, the process of the collective atomic rearrangement of the relaxation center into a more ordered configuration (a transition between states in a two-level system, stimulated by shear deformation) can be accompanied by the localization of a conduction electron in the resulting potential well.
Let us consider a center of relaxation of radius R, which has experienced a shear deformation, to which the potential well with depth −U corresponds to (with respect to the average potential energy of the amorphous state). During the transition between states in a two-level system, the change in the thermodynamic potential ∆Φ is expressed by the formula where E e (R) is the energy of an electron in a potential well with a depth −U and m * is the effective mass of an electron. Here, ϕ < 0 is the difference between the densities of thermodynamic potentials of the new ordered state of atoms in the relaxation center and the initial disordered amorphous state, and σ s is the interphase surface energy. Using the condition 2R = λ = 2π k , where λ is the wavelength of the electron, and substituting k into the expression for the kinetic energy of the electron (13), we find 2 k 2 2m * = 2 π 2 2m * R 2 . ∆Φ becomes negative (the deformed relaxation center is stable) at U, greater than a certain value U 0 . The minimum value ∆Φ corresponds to the critical value R 0 , which determines the radius of the stably existing deformed relaxation center. Let us determine the value of R 0 and U 0 . Equating

∂∆Φ(R) ∂R
to zero, for a large interfacial surface energy |ϕ|R 0 2σ s , we obtain For the value of σ s = 20 erg/cm 2 (1 erg = 10 −7 J; 1 cm = 10 −2 m), which is characteristic of metals [18], and the condition |ϕ|R 0 2σ s , the critical radius of the deformed relaxation center is R 0 ∼ = 0.5 nm, and the value U 0 is approximately 1 eV. It is clear from physical considerations that the minimum radius of a stable deformed relaxation center is approximately equal to the radius of the second coordination sphere. Since the theoretically estimated radius of a phason coincides in order of magnitude with the experimentally founded value of the radius of the nanocrystal [2][3][4]15], it is reasonable to assume that the possible mechanism for the nanocrystal formation in an amorphous NiTi metal film is the phason formation.

Discussion: Physical Representation, Reasons and Conditions for Nanocrystal Formation, and Model Parameters Estimation
In the framework of the model, the equations of the density operator are formulated, taking into account the relaxation term and the field of static shear deformation. The probabilities of the transitions of an atom in a double-well potential by thermal and athermal, i.e., quantum, mechanisms are determined. The obtained results allow us to formulate a physical idea of the nanocrystal formation and make conclusions about proper conditions for the crystallization.
At the initial moment of time, an amorphous metallic medium contains a large number of nanoregions (relaxation centers) as well as free volume and has a disordered atomic structure. The single action of infrasound (shear deformation) on the relaxation center, an atom (located in a double-well potential) or a group of atoms (which have two configurations with slightly different energies) provokes both thermal fluctuation (thermal mechanism) and quantum tunneling to new equilibrium positions (athermal mechanism). At a critical value of deformation ε * (corresponding to the material proportionality limit), the resonance detuning becomes zero and the atom tunnels. At room temperature, local structural rearrangement with small static displacements of atoms (much less than the interatomic distance) occurs with a tunneling frequency higher than the thermal hopping frequency. Thus, at room temperature and lower temperatures, the athermal (quantum) mechanism of inelastic deformation prevails. This conclusion is consistent with the experimental result that with an increase in the initial temperature of the sample (from 25 to 200 • C), the distribution and size of crystalline clusters practically does not change [4]. As a result of inelastic deformation, a nanocluster with a crystalline atomic order appears at the relaxation center, but it may turn out to be unstable in the absence of infrasound.
Since the relaxation time and time of tunnel transition depend on temperature differently, both mechanisms (thermal and athermal) are realized upon the deformation of metallic glasses. With the decreasing of temperature, the thermal motion of atoms is getting to be frozen; as a consequence, the relaxation rate decreases exponentially accordingly to the Arrhenius law. Hence, at low temperatures, the rate of the tunneling transition exceeds the relaxation rate and coherent tunneling will occur. As the temperature rises, the relaxation rate increases faster than the tunneling rate. These values are initially compared, and then the relaxation rate can significantly exceed the tunneling rate. The transition process becomes thermofluctuational, and an ordered nanocluster arises due to thermal fluctuations. Thus, by changing the temperature, one can change the kinetics of the process of inelastic deformation of metallic glasses from quantum to classical.
The conditions for the appearance of a localized state of a conduction electron at the potential well of a nanocrystal are determined. The critical values of the potential well depth and the radius of a stable nanocrystal are calculated. Analysis of the obtained results allows us to formulate the following conclusions.
Under repeated exposure to infrasound, a nanocluster with a crystalline atomic order grows and reaches a critical radius. The nanocluster achieved the value of the critical radius and creates a local potential well in which the conduction electron appears in a bounded state. If the electron is attracted to the region of the unstable crystalline phase and is localized in there, then a decrease in the electron energy can be compensated for the surface energy of the interface. Thus, stable formation arises, also known as a phason (stable nanocrystal). The physical reason for the appearance of the stable nanocrystal is in the conduction electron localization by the potential well created by the nanocrystal.
The radius of the phason is determined by the value of the interfacial surface energy σ s , and an increase σ s leads to a decrease in the radius of the phason. Physically, the radius of the phason cannot be less than the radius of the first coordination sphere. For the value σ s = 20 erg/cm 2 , which is characteristic of metals [18], the critical radius of the phason is R 0 ∼ = 0.5 nm. The critical value of the depth of the potential well U 0 , i.e., the condition under which a phason with R 0 ∼ = 0.5 nm is formed, is approximately 1 eV.
Since the theoretical estimation of the critical radius of a phason coincides in order of magnitude with the experimentally founded value of the radius of a nanocrystal [2][3][4]15], it is reasonable to assume that the mechanism of the formation of a nanocrystal in a metallic glass is the mechanism of the phason formation. When conduction electrons are localized on nanocrystals formed as a result of inelastic deformation, they lower the energy of the system [2,4]. As a result, the state with nanocrystals becomes stable and a locally stable nanostructure appears (with respect to the globally stable crystal structure [3]).
It has been experimentally established [2][3][4]15] that the morphology of nanocrystals depends on the method of nanocrystallization. The difference in the morphology of nanocrystals, isotropic upon annealing [15] and anisotropic upon exposure to infrasound [2][3][4], is associated with different mechanisms of nucleation and growth of nanocrystals. During nanocrystallization by high-temperature annealing, the nucleation and growth of nanocrystals are controlled by thermal fluctuation processes. However, at room temperature, the contribution of thermal fluctuation processes is small, inelastic deformation and nanocrystallization occur due to quantum transitions in nanoscale two-level systems, stimulated by infrasonic vibrations. The anisotropy of the nanocrystal shape is determined by the strain potential tensor, which depends on the local glass structure in the place where the two-level system is located.
Let us estimate the density of the stored potential energy of static displacements of atoms (inelastic deformation) under infrasonic action on the NiTi amorphous film. The dimensions of the film are known from the experiment [2][3][4]: length l = 1·10 −2 m, width r = 1.6·10 −3 m, and thickness d = 4·10 −5 m. The time of the cyclic mechanical action is 10 min, and the frequency of the action is 20 Hz. The amplitude of the tensile displacement of the infrasound is A = 4·10 −6 m, and the relative longitudinal elongation of the film is = A l = 4·10 −4 . It is known that the potential energy density of homogeneous longitudinal elastic deformation is equal to = E 2 ·ε 2 , where E is the modulus of longitudinal elasticity (Young's modulus). The NiTi amorphous film is an alloy in which the main elements are nickel (45 percent) and titanium (41 percent). It is known that E Ni ≈ 220 GPa, E Ti ≈ 110 GPa, and their average value is equal to E ≈ 165 GPa. Taking into account the average value E ≈ 165 GPa and the relative elongation ε = 4·10 −4 , the potential energy density is equal to e ≈ 1.3·10 4 J/m 3 .
Further, let us estimate the numerical value of the potential energy of uniform longitudinal elastic deformation per atom of the medium. Since the concentration of atoms in the medium is approximately equal to 0.9·10 29 atom/m 3 , the energy per atom per one deformation cycle is equal to e a ≈ 1·10 −6 eV/atom. During the time of exposure to infrasound, n = 12, 000 cycles of mechanical deformation occur; therefore, the energy supplied to the system during the entire exposure time is n·e a ≈ 1.2·10 −2 eV/atom. It is more probable that the small part (estimations are about 10 percent) of this energy can be converted into potential energy of inelastic deformation. Therefore, the value of the stored potential energy of inelastic deformation is 0.1·n·e a ≈ 1.2·10 −3 eV/atom, which is much less than the kinetic energy of an atom at room temperature kT 0 ≈ 0.024 eV/atom.
The latent heat of amorphous-crystalline transformation, determined through the temperature difference between amorphous and crystalline states ∆T ≈ 150 K in the TiCu alloy [27], is L a = k∆T ≈ 0.012 eV/atom, which we can consider as an example. For the NiTi alloy, the latent heat of transformation is L a ≈ 0.01 eV/atom [4]. It has been experimentally established [2][3][4] that the latent heat of amorphous-crystalline transformation slightly decreases after 12,000 cycles of infrasonic exposure. Therefore, the free energy of the amorphous state with nanocrystals slightly decreases in comparison with the free energy of the amorphous state. It is reasonable to assume that the latent heat of transformation of an amorphous state into a nanostructured state (nanocrystals in an amorphous matrix) is much less than L a ≈ 0.01 eV/atom, for example, 0.1·L a ≈ 0.001 eV/atom. Thus, the latent heat of transformation of an amorphous state into a nanostructural state 0.1·L a ≈ 0.001 eV/atom coincides in order of magnitude with the value of the stored potential energy of inelastic deformation 0.1·n·e a ≈ 1.2·10 −3 eV/atom.
As was experimentally established in [4], when the frequency is below 5 Hz or the vibration amplitude is down to 1 µm, crystalline clusters do not appear in the amorphous matrix. Therefore, the following conclusion can be formulated: the occurrence of the instability of the amorphous state and following transformation to the nanostructured state is based on the accumulation of the potential energy of inelastic deformation to a critical value equal to the latent heat of the transformation of the amorphous state into the nanostructured state.

Conclusions
The relaxation centers including disordered nanoregions with two-level systems are the causes giving rise to the nanocrystal formation. When exposed to infrasound, the two-level systems, of which there are many, are excited, significantly contributing to inelastic deformation and the formation of nanocrystals. The physical mechanism of the nanocrystallization of metallic glass under mechanical action includes both local thermal fluctuations and the additional quantum tunneling of atoms stimulated by shear deformation. A crystalline nanocluster appears as a result of local atomic rearrangement growing increasingly exposed to infrasound. It is possibly unstable in the absence of infrasound. When the radius of the nanocluster reaches a critical value, the potential well appears. Bounding the conducting electron, it forms the phason. The estimated values (in case of the NiTi film) of the phason's radius and the depth of the nanometer potential well are 0.5 nm and 1 eV, respectively.
We consider the fact that the proposed microscopic mechanism makes it possible to formulate physical reasons and conditions for the appearance of nanocrystals, and we describe results and regularities obtained experimentally both quantitatively and qualitatively.