Part II – February 1969 - Papers - Elastic Calculation of the Entropy and Energy of Formation of Monovacancies in Metals

The formation of a monovacancy in a metal is simulated in an elastic model by the displacement of the surface of a small spherical cavity in a large elastic continuum. The application of linear elasticity to this distortion results in a well- known formula for the energy and an expression for the concomitant entropy change due both to the shear strain in the continuum and also to the dilation of the solid resulting from the boundary conditions at the surface of the solid. Elastic data (the sliear modulus and its temperature coelficient) are used to calculate the entropy and energy of formation for many metals. Despite the simplicity of the assumptions involved, the agreement between the calculated entropies and energies and experimental values is remarkably good. In recent years there has been a large increase in measurements of the absolute concentration of mono-vacancies in metals as a function of temperature. Hence new data for both the energy and the noncon-figurational entropy of formation of monovacancies has become available. Recent measurements&apos; of the anomalous (non-Arrhenius) self-diffusion in many bcc metals has also focused interest on the prediction of the thermodynamic parameters of mono- and multi-vacancies in those metals for which no data are available. Damask and Dienes&apos; have discussed the various theoretical calculations of the energy of formation EL, of a monovacancy. These include simple models involving the breaking of atomic bonds on moving atoms from the interior of a crystal to the surface, models combining elastic calculations with surface-energy terms and detailed quantum mechanical calculations. The simler models give the correct order of magnitude of &, but tend to overestimate it by a factor of about two. The quantum mechanical calculations4"7 have been carried out for the noble and alkali metals with generally reasonably good agreement with the available Ef data. The calculation of entropy of formation Sfv14 lnvolves a fundamental calculation of the perturbation of the phonon spectrum caused by the creation of a vacancy. Huntington, Shirn. and wajda8 have given an approximate evaluation of sJV by considering an Einstein model for the localized vibrations in the immediate neighborhood of the defect and then using elastic theory to calculate the entropy associated with the shear stress field in the distorted crystal (as originally proposed by Zenerg). They also included a term due to the dilation of the crystal. They obtained a value of 1.47k for copper, in good agreement with the experimental value (1.50k). However, Nardelli and Tetta- manzi1° have recently shown that neglecting the coupling between atoms (Einstein Model) may lead to a serious error so the agreement may be somewhat fortuitous. In this work simple linear elastic theory is used to calculate the entropy and energy of formation of mono-vacancies. Despite the simplicity of some of the assumptions involved, the agreement with the available experimental data is remarkable. However. the reasonable degree of success in the application of linear elastic calculations to the excess entropy of a solute atom in a dilute solid solution1&apos; indicates that the application of elastic theory to vacancies. where the interaction of different atomic species is not involved, may not be inappropriate. THE ELASTIC MODEL The metal is assumed to be a spherical elastic continuum. A small spherical cavity of volume V = 4i;v:&apos;/3 is cut from the center. removed. and dissolved rever-sibly in the bulk of the material. TO a good approximation no net atomic bonds are broken and the material does not undergo a volume change although the externally measured volume of the body would increase by V. The radius of the sphere of metal is much larger than r Next a negative pressure is applied to the cavity causing its surface to be displaced inward by an amount simulating the relaxation of the lattice around a monovacancy. In this model the energy and entropy accompanying the distortion are taken as 4, and <. As a first approximation the equation of state for the solid is taken as: r = ro(i + *~D LiJ where K is the bulk modulus. P the hydrostatic pressure. Vo the volume of the material at 0°K and zero pressure. and d+/dT = 30. where 0 is the linear thermal expansion coefficient. The variation of entropy with hydrostatic pressure is given by the Maxwell equation: These equations give the entropy change resulting from increasing the hydrostatic pressure from 0 to P as: and since • we have: This is the entropy arising from the dilation resulting