A simple analytic model of point-ion electrostatics has been previously proposed [R. Magri, S. -H. Wei, and A. Zunger, Phys. Rev. B 42, 11 388 (1990)] in which the magnitude of the net charge ${\mathit{q}}_{\mathit{i}}$ on each atom in an ordered or random alloy depends linearly on the number ${\mathit{N}}_{\mathit{i}}^{\mathrm{}\left(1\right)\mathrm{}}$ of unlike neighbors in its first coordination shell. Point charges extracted from recent large supercell (256–432 atom) local density approximation (LDA) calculations of ${\mathrm{Cu}}_{1\mathrm{-}\mathit{x}}$${\mathrm{Zn}}_{\mathit{x}}$ random alloys now enable an assessment of the physical validity and accuracy of the simple model. We find that this model accurately describes (i) the trends in ${\mathit{q}}_{\mathit{i}}$ vs ${\mathit{N}}_{\mathit{i}}^{\mathrm{}\left(1\right)\mathrm{}}$, particularly for fcc alloys, (ii) the magnitudes of total electrostatic energies in random alloys, (iii) the relationships between constant-occupation-averaged charges 〈${\mathit{q}}_{\mathit{i}}$〉 and Coulomb shifts 〈${\mathit{V}}_{\mathit{i}}$〉 (i.e., the average over all sites occupied by either A or B atoms) in the random alloy, and (iv) the linear relation between the site charge ${\mathit{q}}_{\mathit{i}}$ and the constant-charge-averaged Coulomb shift ${\mathit{V}}_{\mathit{i}}$ (i.e., the average over all sites with the same charge) for fcc alloys. However, for bcc alloys the fluctuations predicted by the model in the ${\mathit{q}}_{\mathit{i}}$ vs ${\mathit{V}}_{\mathit{i}}$ relation exceed those found in the LDA supercell calculations. We find that (a) the fluctuations present in the model have a vanishing contribution to the electrostatic energy. (b) Generalizing the model to include a dependence of the charge on the atoms in the first three (two) shells in bcc (fcc) — rather than the first shell only — removes the fluctuations, in complete agreement with the LDA data. We also demonstrate an efficient way to extract charge transfer parameters of the generalized model from LDA calculations on small unit cells. ©1996 The American Physical Society.

• Received 17 January 1996

DOI:https://doi.org/10.1103/PhysRevB.54.7843