Journal of MetallurgyVolume 2012 (2012), Article ID 517230, 6 pagesdoi:10.1155/2012/517230

Research Article

Study of Tracer Diffusion Mechanism in Amorphous Metal

P. H. Kien,¹ H. V. Hue,² and P. K. Hung³

¹Department of Physics, Thai Nguyen University of Education, Luong Ngoc Quyen, Thai Nguyen, Vietnam²Faculty of General Science, Ho Chi Minh City University of Food Industry, 140 Le Trong Tan, Ho Chi Minh, Vietnam³Department of Computational Physics, Ha Noi University of Science and Technology, 1 Dai Co Viet, Ha Noi, Vietnam

Received 11 March 2012; Revised 27 May 2012; Accepted 12 June 2012

Academic Editor: E. Mittemeijer

Copyright © 2012 P. H. Kien et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The 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.

1. Introduction

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 Ti₆₀N₄₀ and Fe₄₀Ni₄₀B₂₀ which show the cooperative activated movement more like diffusivity in solid state than in liquid [2, 17–20]. In [21], 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 × 10⁵ atoms in a cubic box with periodic boundary conditions. We use the Pak-Doyama potential [16], 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 10⁶ 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.

Figure 1: The schematic illustration of simplexes: (a) 4-simplex; (b) 6-simplex; (c) 5-simplex and neighboring atoms; (d) the new 5-simplex and neighboring atoms after RDA. The dash circle presents the CSP of four atoms, and the arrow is a path along which the MA jumps into CSP.

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 [22]. 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.

Table 1: Structure characteristics of obtained samples. Here, is the position of th peak and is the height of first peak of the RDF.

Systems

/
/
/
/

Sample A
2.51
1.67
1.99
2.52
3.39
3.64

Sample B
2.52
1.67
1.99
2.51
3.38
3.75

Sample C
2.50
1.66
2.01
2.54
3.44
3.78

aCo
2.49
1.69
1.93
2.49
3.35
3.20

bCo
—
1.65
1.90
2.57
3.46
—

xperimental data in [22–24].

Figure 2: The radial distribution function of the obtained samples.

Figure 3: The radial distribution function of sample A and the experimental data in [18].

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.

Table 2: The energy per atom () and the fraction of -simplexes in the obtained samples. Here, is number of atom in one simplex.

Systems
(eV)
4 (×10⁰)
5 (×10⁻¹)
6 (×10⁻²)
7 (×10⁻³)

Sample A
−0.9336
0.961
0. 374
0.173
0.02

Sample B
−0.9462
0.963
0. 358
0.167
0.00

Sample C
−0.9534
0.964
0.339
0.164
0.00

Figure 4: The number of simplexes per atom as function of radius .

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.

Table 3: The number of BB (); the pre-exponential factor () at 560 K in obtained samples.

Systems
(10⁻³)
(m²/s)

Samples A
8.39
0⁻⁶

Samples B
1.34
0⁻⁷

Samples C
0.31
0⁻⁷

cCo₈₉Zr₁₁
—
0⁻⁷

dCo₇₉Nb₁₄B₇
—
0⁻⁶

xperimental data in [25, 26].

(a)

(b)

Figure 5: Potential energy profiles for atom moving into simplex.

Figure 6: The snapshot of typical bubbles detected in the obtained sample.

Figure 7: The distribution of barrier height of diffusing atom.

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 [3]. 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].

Figure 8: The distribution of averaged square displacement of moving atoms.

Figure 9: The temperature dependence of diffusion coefficient; a line is drawn to show the Arrhenius behavior.

4. Conclusions

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.

References

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.

J. Pavlovsky, W. Ulfert, and W. Frank, “Self-diffusion of ⁵⁸Co in amorphous Co₇₉Nb₁₄B₇ during isothermal crystallization,” Materials Chemistry and Physics, vol. 36, no. 3-4, pp. 383–388, 1994.

W. Frank, A. Horner, P. Scharwaechter, and H. Kronmiiler, “Diffusion mechanisms in amorphous alloys,” Materials Science and Engineering, vol. 97, pp. 415–418, 1988.

A. K. Tyagi, M. P. Macht, and V. Naundorf, “Diffusion coefficients of ⁶³Ni in Fe₄₀Ni₄₀B₂₀ metallic glass,” Acta Metallurgica et Materialia, vol. 39, no. 4, pp. 609–617, 1991.

T. Schuler, J. Pavlovskýa, P. Scharwaechter, W. Ulfert, and W. Frank, “Change of the self-diffusion of ⁵⁸Co in Co₇₉Nb₁₄B₇ during transition from the amorphous to the nanocrystalline phase,” Nanostructured Materials, vol. 6, no. 5–8, pp. 863–867, 1995.

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.

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.

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.

S. K. Sharma, S. Banerjee, Kuldeep, and A. K. Jain, “Diffusion measurements in the Fe₈₂B₁₈ amorphous alloy by rutherford backscattering spectrometry,” Acta Metallurgica, vol. 36, no. 7, pp. 1683–1690, 1988.

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.

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.

G. Ruitenberg, P. de Hey, F. Sommerb, and J. Sietsma, “Pressure dependence of the free volume in amorphous Pd₄₀Ni₄₀P₂₀ and its implications for the diffusion process,” Materials Science and Engineering A, vol. 226-228, pp. 397–400, 1997.

Y. Limoge, “Activation volume for diffusion in a metallic glass,” Acta Metallurgica et Materialia, vol. 38, no. 9, pp. 1733–1742, 1990.

R. S. Averback, “Defects and diffusion in amorphous alloys,” MRS Bulletin, vol. 16, pp. 47–52, 1991.

Y. Limoge, “Microscopic and macroscopic properties of diffusion in metallic glasses,” Materials Science and Engineering A, vol. 226-228, pp. 228–236, 1997.

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.

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.

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.

V. Naundorf, M. P. Macht, A. S. Bakai, and N. Lazarev, “The pre-factor, D₀, of the diffusion coefficient in amorphous alloys and in grain boundaries,” Journal of Non-Crystalline Solids, vol. 224, no. 2, pp. 122–134, 1998.

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.

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.

Y. Waseda and S. Tamaki, “Structures of 3D-transition metals in the liquid state,” Philosophical Magazine, vol. 32, no. 2, pp. 273–281, 1975.

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.

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.

V. Naundorf, M.-P. Machta, A. S. Bakaib, and N. Lazarevb, “The pre-factor, D₀, of the diffusion coefficient in amorphous alloys and in grain boundaries,” Journal of Non-Crystalline Solids, vol. 224, no. 2, pp. 122–134, 1998.

J. Pavlovsky and J. Cermak, “Diffusion of ⁵⁹Fe in amorphous Co₇₉Nb₁₄B₇ alloy,” Scripta Metallurgica et Materialia, vol. 30, no. 10, pp. 1317–1320, 1994.