The significance of multipole interactions for the stability of regular structures composed from charged particles.

Identifying the forces responsible for stabilising binary particle lattices is key to the controlled fabrication of many new materials. Experiments have shown that the presence of charge can be integral to the formation of ordered arrays; however, a complete analysis of the forces responsible has not included many of the significant lattice types that may form during fabrication. A theory of many-body electrostatic interactions has been applied to six lattice stoichiometries, AB, AB2, AB3, AB4, AB5 and AB6, to show that induced multipole interactions can make a very significant (>80 %) contribution to the total lattice energy of arrays of charged particles. Particle radii ratios which favour global minima in electrostatic energy are found to be the same or a close match to those observed by experiment. Although certain lattice types exhibit local energy minima, the calculations show that many-body rather than two-body interactions are ultimately responsible for the structures observed by experiment. For a lattice isostructural with CFe4, a particle size ratio not previously observed is found to be particularly stable due to many-body effects

An integral equation approach to calculate electrostatic interactions in many-body dielectric systems

In this article, a numerical method to compute the electrostatic interaction energy and forces between many dielectric particles is presented. The computational method is based on a Galerkin approximation of an integral equation formulation, which is sufficiently general, as it is able to treat systems embedded in a homogeneous dielectric medium containing an arbitrary number of spherical particles of arbitrary size, charge, dielectric constant and position in the three-dimensional space. The algorithmic complexity is linear scaling with respect to the number of particles for the computation of the energy which has been achieved through the use of a modified fast multipole method. The method scales with the third power of the degree of spherical harmonics used in the underlying expansions, for general three-dimensional particle configurations. Several simple numerical examples illustrate the capabilities of the model, and the influence of mutual polarization between particles in an electrostatic interaction is discussed