No abstract

No abstract

We have constructed a large database of 171 first-principles total energy calculations of aluminum structures with coordination number ranging from 0 to 12 and nearest neighbor distance from 2.0 to 5.7 A. These calculations throw quantitative light on the nature of metallic bonding. We have tested a large range of "'glue" models with this database: none gives a rms error less than 0,10 eV per atom. We attribute this minimum error to bonding effects beyond the glue formalism, and as such it sets the ultimate reliability of any glue scheme.

Ab initio. calculations are presented of the cohesive energies of aluminium in a number of diverse hypothetical structures which span a wide range of the coordination number, C, from C = 0 to C = 12. The calculations have been performed to investigate the nature of multi-atom bonding, its dependence on C and to form a database for testing and developing empirical and semi-empirical models. The results support the saturation of cohesive energy for large C predicted by several simple theoretical models. Calculations on the same structures using semi-empirical schemes suggest that these methods might have a greater degree of accuracy than had previously been believed.

A fast method is presented for the calculation of total energies within the density functional formalism for systems that require a large number of k-points. Approximate solutions at a large number of k-points are obtained from exact self-consistent solutions at a smaller number of k-points using an extension of the k.p method. The method is demonstrated by a calculation of the total energy of FCC aluminium at a series of lattice parameters. An analysis is presented that partitions the errors in the calculations into those inherent in any finite sampling procedure and those that are specific to the k.p method. The extra error due to the k.p method is shown to be around 10% of energy differences for these calculations but a simple modification to the calculations requiring no more computational time reduces this error to 2%. For these simple calculations the time saved by using the k.p method is fairly small but for larger calculations the method is several orders of magnitude quicker. Errors from either source are shown to be far less important in the electronic potentials than in the eigenvalue sums. It is concluded that in any total energy calculation fewer k-points are required to describe the potential than to calculate the eigenvalue sum. This could be exploited to gain an order of magnitude saving in computational time in any calculation.