Calculation of point-defect entropy in metals

We investigate previously used methods and propose a new method for atomist ic calculations of point-defect entropies in metals within the harmonic app roximation to lattice vibrations. The key problem is to predict accurately the defect formation entropy in a macroscopic crystal from atomistic calcul ations performed on a small system containing relatively few atoms. The res ults of atomistic calculations may depend significantly on the system size, geometry and boundary conditions. Two previously used methods, which we ca ll the supercell and embedded-cluster methods, are analysed in two ways: fi rstly, within a linear elasticity model of a point defect and, secondly, by atomistic calculations for a vacancy and an interstitial in copper using a n embedded-atom potential. The results of atomistic calculations confirm th e linear elasticity analysis and show that the supercell method is much mor e accurate than the embedded-cluster method. However, the latter is useful for computing the defect core entropy, which turns out to be a well-defined physical quantity. We propose a new method of defect entropy calculations that combines the embedded-cluster method with a quasicontinuum approximati on outside the cluster. This method, which we call an elastically corrected embedded-cluster method, has an accuracy comparable with that of the super cell method and extends defect entropy calculations towards larger system s izes.