Determining point charge arrays that produce accurate ionic crystal fieldsfor atomic cluster calculations

In performing atomic cluster calculations of local electronic structure def ects in ionic crystals, the crystal is often modeled as a central cluster o f 5-50 ions embedded in an array of point charges. For most crystals, howev er, a finite three-dimensional repeated array of unit cells generates elect rostatic potentials that are in significant disagreement with the Madelung (infinite crystal) potentials computed by the Ewald method. This is illustr ated for the cubic crystal CaF2. We present a novel algorithm for solving t his problem for any crystal whose unit cell information is known: (1) the u nit cell is used to generate a neutral array containing typically 10 000 po int charges at their normal crystallographic positions; (2) the array is di vided into zone 1 (a volume defined by the atomic cluster of interest), zon e 2 (several hundred additional point charges that together with zone 1 fil l a spherical volume), and zone 3 (all other point charges); (3) the Ewald formula is used to compute the site potentials at all point charges in zone s 1 and 2; (4) a system of simultaneous linear equations is solved to find the zone 3 charge values that make the zone 1 and zone 2 site potentials ex actly equal to their Ewald values and the total charge and dipole moments e qual to zero, and (5) the solution is checked at 1000 additional points ran domly chosen in zone 1. The method is applied to 33 different crystal types with 50-71 ions in zone 1. In all cases the accuracy determined in step 5 steadily improves as the sizes of zones 2 and 3 are increased, reaching a t ypical rms error of 1 mu V in zone 1 for 500 point charges in zone 2 and 10 000 in zone 3. (C) 2000 American Institute of Physics. [S0021-9606(00)3010 3-9].