Interplay of order-disorder phenomena and diffusion in rigid binary alloys in the presence of vacancies: Monte Carlo simulations

A. DE VIRGILIIS; K. BINDER

PHYSICAL REVIEW B - CONDENSED MATTER AND MATERIALS PHYSICS

Año: 2006 vol. 73 p. 1 - 15

0163-1829

Transport phenomena are studied for a binary (AB) alloy on a rigid square lattice with nearest-neighbor attraction between unlike particles, assuming a small concentration cv of vacancies V being present, to which A (B) particles can jump with rates A (B) in the case where the nearest-neighbor attractive energy AB is negligible in comparison with the thermal energy kBT in the system. This model exhibits a continuous order-disorder transition for concentrations cA,cB=1–cA–cV in the range ccAc, with c=(1–m*–cV)/2, c=(1+m*–cV)/2, m*0.25, the maximum critical temperature occurring for c*=cA=cB=(1–cV)/2—i.e., m*=0. This phase transition belongs to the d=2 Ising universality class, demonstrated by a finite-size scaling analysis. From a study of mean-square displacements of tagged particles, self-diffusion coefficients are deduced, while applying chemical potential gradients allows the estimation of Onsager coefficients. Analyzing finally the decay with time of sinusoidal concentration variations that were prepared as initial condition, also the interdiffusion coefficient is obtained as function of concentration and temperature. As in the random alloy case (i.e., a noninteracting ABV model) no simple relation between self-diffusion and interdiffusion is found. Unlike this model mean-field theory cannot describe interdiffusion, however, even if the necessary Onsager coefficients are estimated via simulation.