Study of Tracer Diffusion Mechanism in Amorphous Metal
The statistical relaxation (SR) simulation has been conducted to study the behavior of simplexes and bubbles (BB) in amorphous Co metal containing atoms. The simulation reveals that the fraction of 4-simplex increases and -simplex decreases depending upon relaxation degree. The simulation found that a large number of BB vary upon relaxation degree, which could play a role of diffusion vehicle for Co atoms in amorphous matrix. The idea of the diffusion mechanism in amorphous metal is described as follows: the elemental atomic movement includes a jump of neighboring atom into the BB and then a collective displacement of a large number of atoms around BB.
The diffusion behavior in amorphous metals (AMs) has been intensively studied by both experiment and simulation for long times [1–8]. It is found that there are many specific properties of diffusion in AM compared to crystal matters. For example, the tracer diffusivity in a well-relaxed sample is much slower than the one in an as-quenched sample [9, 10]. This relaxation effect is interpreted by the reduction of vacancies in supersaturation until the relaxation is over. In the well-relaxed sample, conversely, the tracer atoms diffuse via collective movement of a group of neighboring atoms. The experimental studies in [11–13] on isotope effect, pressure dependence, and irradiation-enhanced diffusivity are sometimes in contradiction with the diffusion described above. Furthermore, there is not a clear definition of vacancy. Computer simulation, on the other hand, reveals unstableness of vacancies in amorphous matrix. Several works found a continuous spectrum of spherical voids in amorphous matter, but their size is less than atomic radius [14–16]. The free volume and two-level state theories are also employed to interpret the diffusion behavior of amorphous matter, but they cannot properly describe the diffusivity in some amorphous matters such as Ti60N40 and Fe40Ni40B20 which show the cooperative activated movement more like diffusivity in solid state than in liquid [2, 17–20]. In , Sietsma and Thijsse analyzed different types of holes in amorphous matters and found that the number of holes surrounded by ten or more atoms decreases strongly in a well-relaxed sample. Furthermore, they argue the importance of big holes for atomic diffusivity.
Our previous study shows that amorphous matters suffer from a number of spherical voids, the size of which is closed to atomic size. Because the concentration of those voids weakly depends on temperature, it is very sensitive to the relaxation degree. We call them “native vacancy,” which fully disappears until the amorphous matters transform into crystalline solid. However, the analysis of this study is based only on the geometrical consideration. Therefore, the account of potential barrier could provide more accurate elucidation of the voids for diffusion in amorphous matters. The aim of the present work is the study of simplexes and their role of them in diffusion in Co amorphous metals by using SR simulation method as a way to clarify diffusion mechanism in amorphous matters.
2. Calculation Method
The simulation has been conducted for the sample consisting of 2 × 105 atoms in a cubic box with periodic boundary conditions. We use the Pak-Doyama potential , and the density is taken from real amorphous metal. The inter-atomic potential has the following form: where is the distance between and atoms. Molecular dynamic (MD) simulation could reproduce the structure of AM; however, this method is time consuming. Thus, we employ the SR method, which provides the structure of AM similar to that of the low-temperature MD model. Accordingly, each atom moves by a length in the direction of the force acting on it from all remaining atoms. The length is equal to 0.01–0.4 . This movement is repeated many times until the system reaches an equilibrium state. More details about the SR method can be found in [16, 21].
Initial configuration is generated by randomly placing all atoms in a simulation box. Then, the sample is relaxed over several thousand steps until the system attains the equilibrium, for example, the energy of system fluctuates around a constant value and the pressure is equal to zero.The model A is constructed by relaxing with over 200 SR steps and then is treated with by 106 SR steps. To investigate the relaxation effect, two additional samples (samples B and C) are prepared with the same density as sample A, but their potential energy is lower (more stable state). Note that the lower the energy of the system the bigger its relaxation degree. More stable samples B and C can be constructed many times with relaxing the model A that has , which likes the shaking of atomic structure. Then, they are again relaxed with until the system reaches a new equilibrium.
Four neighboring atoms form a tetrahedron and have a circumsphere (CSP) whose surface passes through the vertices of this tetrahedron. Consider only such tetrahedron with CSP that does not contain any atom inside. For example, those atoms are the nearest neighboring. Let and be the radius of CSP and the number of atoms on the surface of CSP, respectively. The atom on surface of CSP is determined as one that is located from the center of CSP at a distance in range of . hereafter we call it -simplex (see Figures 1(a) and 1(b)). If and is large enough, then atoms form an atomic cage like a bubble (BB), for example, a large group of atoms gathered around a large void. The BB is unstable and it may break up leading to diffusion. Therefore, it is interesting to clarify which one among -simplexes is BB and how BB breaks up. Hence, the next calculation is performed as follows: first for every simplex, we test all atoms to perform jump inside the CSP. Second, the atom selected to move is inserted inside CSP and the system is relaxed until it reaches equilibrium.
3. Results and Discussion
3.1. The Sample
Figure 2 shows the radial distribution function (RDF) of the amorphous samples. As shown in Figure 2, although the energy per atom of the considered samples changes from −0.9336 to −0.9534 eV, the RDF for the considered samples is identical. This indicates the nonsensitivity of the function to the change in local microstructure of AMs. This result further indicates that the relaxation degree almost does not affect RDF, but it is reflected by the concentration of simplex. To test the validity of the constructed samples, we have compared our obtained RDF with the experimental data. As shown from Figure 3 the structural characteristics of our samples are in accordance with the experiment data . In addition, the function has a splitting second peak, which is often thought to be related to the icosahedrons in amorphous systems [17, 18, 20]. The main characteristics of the samples are listed in Table 1. The comparison with experiment shows agreement in the position and height of the first peak of RDF, which indicates the reality of the constructed sample.
3.2. Simplex and Role of Simplex
Table 2 presents the energy per atom and fraction of simplex found in the obtained samples. One can see that the energy per atom decreases from −0.9336 to −0.9534 eV, upon relaxation degree. The fraction of 4-simplex increases as the sample potential energy decreases, that is, the more stable state is, bigger the number of 4-simplexes is. For other kinds of simplex, we observe the opposite trend. Thus, it means that due to fast quenching from liquid the parking of atoms in AM is not efficient. Therefore, there is always an amount of structural defects like large void (free volume). The relaxation is accompanied by annihilation of those defects. As mentioned above, the bigger number is, the bigger the size of simplex and void inside it is. Therefore, Table 2 specifies that the monotony decrease in the number of -simplex with from less relaxed (sample A) to well relaxed state (sample C) indicates the annihilation of structural defects in amorphous matrix. This result also can be seen in the radius distribution of simplex as in Figure 4. From Figure 4, one can see clearly the significant decrease in the number of large simplex () in the sample A as compared to the sample C. We proposed that the large simplex () is broken up into 4-simplex under the relaxation degree. So as the concentration of 4-simplex increases, the large simplex decreases.
We now turn to discuss the role of large simplex for diffusivity in AM. For each simplex, we examine the potential energy profile (PEP) for neighboring atoms moving into CSP of the simplex. The PEP is determined as follows: each of atoms of the simplex is taken and then we move this atom step by step towards the center of CSP. The step length is set to be . At each step, the potential energy of the moving atom (MA) is recorded. Finally, we obtain the PEP for MA. Figure 5 displays the typical PEPs for several simplexes detected in our samples. Curve (b) (Figure 5(b)) shows the monotonous increase in the MA energy indicating that the MA could not jump into the CSP due to very high potential barrier. Curve (a) (Figure 5(a)), in converse, have a pronounced maximum. Here, the barrier height is determined by the maximum and initial points in PEP. They are like the PEP for tracer atom jumping in vacancy incrystalline lattice. Therefore, the MA can jump into CSP, and the simplex plays a role of the diffusion vehicle assisting tracer atom to move in amorphous matrix. The simplex having PEP like curves a, with corresponding MA, is called the bubbles (BB) and the diffusing atom (DA), respectively. Figure 6 shows some typical -bubbles in the obtained samples. The number of BB () is found in Table 3. As shown in Table 3, the number of BB decreases depending upon relaxation degree. The dependence of the number of DA on their corresponding barrier height is shown in Figure 7. Most frequent barrier height lies in the interval of 0.5–1.9 eV. Because the diffusion coefficient in general is proportional to , it is given as follows: Here, is the correlation factor, is the attempt frequency, is the jump length, are the effective migration entropy and energy for diffusion in disordered media, respectively, is the Boltzmann constant. Note that the parameter represents the mean square displacement of different jump lengths. Because the SR model is in fact the molecular dynamic (MD) model in the limit of zero temperature, the number of simplexes found in the samples is independent of temperature and it varies only with the relaxation degree. As a result, the activation energy is the migration energy (), and the preexponential factor is given as Furthermore, the decrease in diffusion coefficient under relaxation is ascribed by the annihilation of the bubbles. Therefore, the decrease in Co diffusivity is found to be 27 times for samples A and C.
Several experimental findings shows that (i) sudden temperature change during diffusion annealing results in instantaneous change of diffusion coefficient () and (ii) The self-diffusion enthalpy of transition metal for certain amorphous alloy seems to be the migration enthalpy, but not a sum of migration and formation enthalpy of diffusion vehicles as it is the case of self-diffusion in crystal . These experimental observations are interpreted as a result of direct diffusion mechanism occurring in relaxed amorphous sample. Our simulation result shows possible diffusivity via the bubble. It is consistent with the experiment because the concentration of bubbles is independent of temperature, and they vary only with relaxation degree. So, the activation energy for diffusion via the bubble is equal to the migration energy ().
To estimate the parameter in (2), we replace DA into CSP and then relax the system until a new equilibrium is attained. This procedure for convenience is called the remove diffusion atom (RDA). The distribution of mean square displacements for 100 RDAs is displayed in Figure 8. For sample C, the parameter of most RDAs is less than 10 (92 atoms). Because the jump length of DA lies in the interval of , the contribution of DA to will be essential and the RDA locates only in the small region nearby the bubble. In the case of sample A, we observe a very large value of (more than 10 ), which represents a collective movement of a large number of atoms. Obviously, the atomic movement is spread over a large volume inside the amorphous matrix as shown in Figure 1(d). Therefore, one may see two distinct diffusion mechanisms occurring in the as-quenched (sample A) and well-relaxed metastable states (sample C). The first one is like the hoping mechanism via bubble. The second one relates to the collective mechanism involving a large number of atoms in each elemental diffusion movement. Assume that , , , for sample C and 100 for sample A. The diffusion coefficient of Co atom can be estimated on the base of (2). Although the barrier heights for different jumps vary from 0.5 to 1.9 eV, the diffusion coefficient is found to obey the Arrehenius behavior as shown in Figure 9. According to (3), the pre-exponential factor is calculated as shown in Table 3. The calculation data is consistent with the experimental data for the amorphous samples [25, 26].
A comprehensive study on both the number of simplexes and a bubbles diffusion mechanism in AM has been carried out using SR simulation with the pair Pak-Doyama inter-atomic potential upon the change of relaxation degree. Conclusions are as followed: (i)The microstructure of the cobalt samples containing atoms is in agreement with the experimental data. The fraction of 4-simplex increases and -simplex () decreases depending upon relaxation degree. The simulation reveals that a large number of bubbles vary with relaxation degree in AM. (ii)A bubbles diffusion mechanism in AM is proposed for the elemental atomic movement, which includes a jump of neighboring atom into the bubbles and then collective displacement of a number of atoms. It may be seen two distinct diffusion mechanisms occurring in the as-quenched and well-relaxed samples. The decrease in pre-exponential factor upon relaxation is ascribed to the partial annihilation of bubbles. The obtained pre-exponential factor of Co atom is found to be in accordance with the experiment data.
J. Horváth, J. Ott, K. Pfahler, and W. Ulfert, “Tracer diffusion in amorphous alloys,” Materials Science and Engineering, vol. 97, pp. 409–413, 1988.View at: Google Scholar
J. Pavlovsky, W. Ulfert, and W. Frank, “Self-diffusion of 58Co in amorphous Co79Nb14B7 during isothermal crystallization,” Materials Chemistry and Physics, vol. 36, no. 3-4, pp. 383–388, 1994.View at: Publisher Site | Google Scholar
W. Frank, A. Horner, P. Scharwaechter, and H. Kronmiiler, “Diffusion mechanisms in amorphous alloys,” Materials Science and Engineering, vol. 97, pp. 415–418, 1988.View at: Publisher Site | Google Scholar
A. K. Tyagi, M. P. Macht, and V. Naundorf, “Diffusion coefficients of 63Ni in Fe40Ni40B20 metallic glass,” Acta Metallurgica et Materialia, vol. 39, no. 4, pp. 609–617, 1991.View at: Google Scholar
T. Schuler, J. Pavlovskýa, P. Scharwaechter, W. Ulfert, and W. Frank, “Change of the self-diffusion of 58Co in Co79Nb14B7 during transition from the amorphous to the nanocrystalline phase,” Nanostructured Materials, vol. 6, no. 5–8, pp. 863–867, 1995.View at: Publisher Site | Google Scholar
S. Flege, U. Fecher, and H. Hahn, “Diffusion in amorphous NiZrAl alloys,” Journal of Non-Crystalline Solids, vol. 270, no. 1–3, pp. 123–128, 2000.View at: Google Scholar
A. Grieschea, T. Zumkleya, M.-P. Machta, S. Suzukib, and G. Frohbergb, “Diffusion in PdCuNiP alloys from the amorphous to the liquid state,” Materials Science and Engineering A, vol. 375–377, pp. 285–287, 2004.View at: Publisher Site | Google Scholar
S. Chakravarty, M. Gupta, A. Gupta et al., “Fe and N self-diffusion in amorphous FeN: a SIMS and neutron reflectivity study,” Acta Materialia, vol. 57, no. 4, pp. 1263–1271, 2009.View at: Publisher Site | Google Scholar
S. K. Sharma, S. Banerjee, Kuldeep, and A. K. Jain, “Diffusion measurements in the Fe82B18 amorphous alloy by rutherford backscattering spectrometry,” Acta Metallurgica, vol. 36, no. 7, pp. 1683–1690, 1988.View at: Google Scholar
W. D. Roos, J. D. Plessis, and G. N. van Wyk, “Diffusion of silicon in Fe-based amorphous and crystalline alloys,” Applied Surface Science, vol. 40, no. 4, pp. 303–307, 1990.View at: Google Scholar
Y. Limoge, “Role of energetic disorder on diffusion in amorphous alloys,” Journal of Non-Crystalline Solids, vol. 117-118, no. 2, pp. 605–608, 1990.View at: Google Scholar
G. Ruitenberg, P. de Hey, F. Sommerb, and J. Sietsma, “Pressure dependence of the free volume in amorphous Pd40Ni40P20 and its implications for the diffusion process,” Materials Science and Engineering A, vol. 226-228, pp. 397–400, 1997.View at: Publisher Site | Google Scholar
Y. Limoge, “Activation volume for diffusion in a metallic glass,” Acta Metallurgica et Materialia, vol. 38, no. 9, pp. 1733–1742, 1990.View at: Google Scholar
R. S. Averback, “Defects and diffusion in amorphous alloys,” MRS Bulletin, vol. 16, pp. 47–52, 1991.View at: Google Scholar
Y. Limoge, “Microscopic and macroscopic properties of diffusion in metallic glasses,” Materials Science and Engineering A, vol. 226-228, pp. 228–236, 1997.View at: Google Scholar
P. K. Hung, H. V. Hue, and L. T. Vinh, “Simulation study of pores and pore clusters in amorphous alloys Co100-xBx and Fe100-yPy,” Journal of Non-Crystalline Solids, vol. 352, no. 30-31, pp. 3332–3338, 2006.View at: Publisher Site | Google Scholar
A. Zhu, G. J. Shiflet, and S. J. Poon, “Diffusion in metallic glasses: analysis from the atomic bond defect perspective,” Acta Materialia, vol. 56, no. 14, pp. 3550–3557, 2008.View at: Publisher Site | Google Scholar
P. K. Leung and J. G. Wright, “Structural investigations of amorphous transition element films: II. Chromium, iron, manganese and nickel,” Philosophical Magazine, vol. 30, no. 5, pp. 995–1008, 1974.View at: Publisher Site | Google Scholar
V. Naundorf, M. P. Macht, A. S. Bakai, and N. Lazarev, “The pre-factor, D0, of the diffusion coefficient in amorphous alloys and in grain boundaries,” Journal of Non-Crystalline Solids, vol. 224, no. 2, pp. 122–134, 1998.View at: Google Scholar
P. K. Hung, P. H. Kien, and L. T. Vinh, “Evidence of ‘microscopic bubbles’ and a new diffusion mechanism for amorphous alloys,” Journal of Physics, vol. 22, Article ID 035401, pp. 1–5, 2010.View at: Publisher Site | Google Scholar
J. Sietsma and B. J. Thijsse, “Characterization of free volume in atomic models of metallic glasses,” Physical Review B, vol. 52, no. 5, pp. 3248–3255, 1995.View at: Publisher Site | Google Scholar
Y. Waseda and S. Tamaki, “Structures of 3D-transition metals in the liquid state,” Philosophical Magazine, vol. 32, no. 2, pp. 273–281, 1975.View at: Google Scholar
P. K. Leung and J. G. Wright, “Structural investigations of amorphous transition element films I. Scanning electron diffraction study of cobalt,” Philosophical Magazine, vol. 30, no. 1, pp. 185–194, 1974.View at: Publisher Site | Google Scholar
P. K. Leung and J. G. Wright, “Structural investigations of amorphous transition element films—2. Chromium, iron, manganese and nickel,” Philosophical Magazine, vol. 30, no. 5, pp. 995–1108, 1974.View at: Google Scholar
V. Naundorf, M.-P. Machta, A. S. Bakaib, and N. Lazarevb, “The pre-factor, D0, of the diffusion coefficient in amorphous alloys and in grain boundaries,” Journal of Non-Crystalline Solids, vol. 224, no. 2, pp. 122–134, 1998.View at: Publisher Site | Google Scholar
J. Pavlovsky and J. Cermak, “Diffusion of 59Fe in amorphous Co79Nb14B7 alloy,” Scripta Metallurgica et Materialia, vol. 30, no. 10, pp. 1317–1320, 1994.View at: Publisher Site | Google Scholar