Vacancy Migration

The atomistic process of vacancy migration consists of one atom next to the vacant site jumping into this site and leaving behind another vacant site. The jump is thermally activated, and transition state theory predicts a diffusion coefficient for vacancy migration in cubic crystals of the form

Dv = Wvd02 exp(sm)exp(—EV/kgT)

= DVexp(—EV1 / kg T) [13]

Here, vLV is an average frequency for lattice vibrations, d0 is the nearest neighbor distance between atoms, S^ is the vacancy migration entropy, and E^ is the energy for vacancy migration. It is in fact the energy of an activation barrier that the jumping atom must over­come, and when it temporarily occupies a position at the height of this barrier, the atomic configuration is referred to as the saddle point of the vacancy. It will be considered in greater detail momentarily.

Values obtained for Щ from experimental mea­surements are shown in Figure 8 as a function of the melting point. While we notice again a trend similar to that for the vacancy formation energy, we find that Em for fcc and bcc metals apparently follow different correlations. However, the correlation for bcc metals is rather poor, and it indicates that E^ may be related to fundamental properties of the metals other than the melting point.

The saddle point configuration of the vacancy involves not just the displacement of the jumping atom but also the coordinated motion of other atoms that are nearest neighbors of the vacancy and of the jumping atom. These nearest neighbor atoms

Table 3

Vacancy relaxation volumes for metals

Metal

go(Jm 2)

m (GPa)

n

HV (eV)

V Vі/V (model)

V Vі1 /V (experiment)

Ag

1.19

33.38

0.354

1.11 ± 0.05

—0.247 ± 0.005

Al

1.1

26.18

0.347

0.67 ± 0.03

—0.311 ± 0.003

—0.05, —0.38

Au

1.45

31.18

0.412

0.93 ± 0.04

—0.262 ± 0.003

—0.15 to —0.5

Cu

1.71

54.7

0.324

1.28 ± 0.05

—0.259 ± 0.005

—0.25

Ni

2.28

94.6

0.276

1.79 ± 0.05

—0.236 ± 0.004

—0.2

Pb

0.57

10.38

0.387

0.58 ± 0.04

—0.282 ± 0.005

Pd

1.91

53.02

0.374

1.7, 1.85

—0.239, —0.225

Pt

2.40

65.1

0.393

1.35 ± 0.05

—0.260 ± 0.003

—0.24, —0.42

Cr

2.23

117.0

0.209

2.0 ± 0.3

—0.218 ± 0.02

a-Fe

2.31

90.4

0.278

1.4, 1.89

—0.278, —0.245

—0.05

Mo

2.77

125.8

0.293

3.2 ± 0.09

—0.191 ± 0.004

—0.1

Nb

2.54

39.6

0.397

2.6, 3.07

—0.284, —0.258

Ta

2.76

89.9

0.324

2.2, 3.1

—0.264, —0.228

V

2.51

47.9

0.361

2.2 ± 0.4

—0.298 ± 0.028

W

3.09

160.2

0.280

3.1, 4.1

—0.201, —0.161

lie at the corners of a rectangular plane as shown in
Figure 9. As the jumping atom crosses this plane,
they are displaced such as to open the channel. This
coordinated motion can be viewed as a particular
strain fluctuation and described in terms of phonon
excitations. In this manner, Flynn11 has derived the
following formula for the energy of vacancy migra-
tion in cubic crystals.

Em = ______ 15CH C44(C11 ~ C12 VX 14i

V 2[C11(C11 — C12) + C44(5C11 — 3 C12)]

Here, a is the lattice parameter, Ch1, C12, and C44 are
elastic moduli, and X is an empirical parameter that
characterizes the shape of the activation barrier and
can be determined by comparing experimental
vacancy migration energies with values predicted by
eqn [14]. Ehrhart eta/.7,12 recommend that X = 0.022
for fcc metals and X = 0.020 for bcc metals.

2.5

 

In the derivation of Flynn,11 only the four nearest neighbor atoms are supposed to move, while all other atoms are assumed to remain in their normal lattice positions. On the other hand, Kornblit eta/.13 treat the expansion of the diffusion channel as a quasistatic elastic deformation of the entire surrounding mate­rial. The extent of the expansion is such that the opened channel is equal to the cross-section of the jumping atom, and a linear anisotropic elasticity cal­culation is carried out by a variational method to determine the energy involved in the channel expan­sion. A vacancy migration energy is obtained for fcc metals of

 

image024

[15]

 

image025
image026
image027

Figure 9 Second nearest neighbor atom (blue) jumping
through the ring of four next-nearest atoms (green) into
adjacent vacancy in a fcc structure.

 

image028image029image030

image031

("£max _ e™1). As a result, Kornblit14 assumes that the vacancy migration energy for bcc metals is given by

 

To determine the preexponential factor for self­diffusion

 

if E^ _ £yin < 2kT if EVlax _ E™11 > 2kT

 

max EV ’ max min 2EV EV ;

 

DSd — VLV*2exp( (4 + sm) /k) [21]

requires the values for the entropy 4 + and for the attempt frequency nLV. Based on theoretical esti­mates, Seeger and Mehrer5 recommend a value of 2.5 k for the former. The atomic vibration of nearest neighbor atoms to the vacancy is treated within a sinusoidal potential energy profile that has a maxi­mum height of Ef. For small-amplitude vibrations, the attempt frequency is then given by

 

m

EV —

 

[18]

 

Using the formulae of Flynn and Kornblit, we compute the vacancy migration energies and compare them with experimental values in Figure 10.

With a few exceptions, both the Flynn and the Kornblit values are in good agreement with the exper­imental results.

The self-diffusion coefficient determines the transport of atoms through the crystal under condi­tions near the thermodynamic equilibrium, and it is defined as

 

image032

Vlv

 

Dsd — D? CVq — VLV«2exp( (sV + Sm)/k) exp(_Qso/kT) — D0Dexp(_QsD/kT) [19]

where the activation energy for self-diffusion is

Qsd — 4 + Em [20]

The most accurate measurements of diffusion coef­ficients are done with a radioactive tracer isotope of the metal under investigation, and in this case one obtains values for the tracer self-diffusion coeffi­cient dTd — fDsD that involves the correlation factor f For pure elemental metals of cubic structure, f is a constant and can be determined exactly by com — putation.15 For fcc crystals, f— 0.78145, and for bcc crystals, f— 0.72149.

 

[22]

 

Vlv

 

crystals where M is the atomic mass.

In contrast, Flynn11 assumes that the atomic vibra­tions can be derived from the Debye model for which the average vibration frequency is

 

image033

2

 

• Evm, Flynn, fcc О Evm, Kornblit, fcc ■ Evm, Flynn, bcc □ Evm, Kornblit, bcc

 

■SI

>,

O)

® 1.5

ф

c

q

оз

D)

1

>s

о

c

я

о

03

>

оз

0.5

ф

 

О

г

 

Q •

 

image034