The computer simulation of vacancy annihilation to dislocations has been made by the Monte Carlo method to interpret the annealing kinetics of vacancies in deformed gold, which is described by the Harper form as C=C0exp(−K·tm), where C is the vacancy concentration, t is time, K and m are constants, as reported in the preceding paper. Three factors most important for the kinetics have been taken into account in the calculation; (a) vacancy flow due to the concentration gradient, (b) vacancy flow due to the gradient of the vacancy-dislocation interaction (elastic) potential, and (c) sink efficiency of dislocation. In the simulation, the sink efficiency has been treated in such a way that a vacancy at a nearest-neighbour site of dislocation must overcome an energy barrier ΔE on annihilation. The results are summarized as follows:(1) In the case where dislocation acts as a perfect sink (ΔE=0), the value of exponent m in the early stage of annealing indeed ranges from 0.55 to 0.69, depending on the dislocation density, in accord with the experiment, but in the later stage approaches 1.0 in marked contrast to the experiment. The similar tendency is found even when the interaction potential is assumed to be ten times as large as the theoretically reasonable value. The present results, therefore, do not support the Harper’s analysis that the value of m is uniquely defined by the form of the vacancy-dislocation interaction potential.(2) When a fixed sink efficiency (ΔE=const.>0) of dislocation is assumed, bearing in mind the extended dislocation, the obtained value of m becomes much higher than the experimental value, increasing with the increase in the energy barrier ΔE.(3) A variable sink efficiency model has been found to be most appropriate, in which the energy barrier ΔE is assumed to increase with the degree of bowing-out of dislocation as a function of the number of vacancies absorbed. The exponent m in this case varies with the dislocation density in a similar way as the experimental results; for the low dislocation density the whole process of annealing can be described by a single value of m, while for the high dislocation density the annealing process consists of two stages, i.e., the early stage with the exponent m′=0.7 and the later stage with the exponent m=0.5.
Read full abstract