We present pseudo-potential coefficients for the first two rows of the periodic table. The pseudo potential is of a novel analytic form, that gives optimal efficiency in numerical calculations using plane waves as basis set. At most 7 coefficients are necessary to specify its analytic form. It is separable and has optimal decay properties in both real and Fourier space. Because of this property, the application of the nonlocal part of the pseudopotential to a wave-function can be done in an efficient way on a grid in real space. Real space integration is much faster for large systems than ordinary multiplication in Fourier space since it shows only quadratic scaling with respect to the size of the system. We systematically verify the high accuracy of these pseudo-potentials by extensive atomic and molecular test calculations.

We generalize the concept of separable dual-space Gaussian pseudopotentials to the relativistic case. This allows us to construct this type of pseudopotential for the whole periodic table and we present a complete table of pseudopotential parameters for all the elements from H to Rn. The relativistic version of this pseudopotential retains all the advantages of its nonrelativistic version. It is separable by construction, it is optimal for integration on a real space grid, it is highly accurate and due to its analytic form it can be specified by a very small number of parameters. The accuracy of the pseudopotential is illustrated by an extensive series of molecular calculations. 71.15.Hx

No abstract

Methods exhibiting linear scaling with respect to the size of the system, so called O(N) methods, are an essential tool for the calculation of the electronic structure of large systems containing many atoms. They are based on algorithms which take advantage of the decay properties of the density matrix. In this article the physical decay properties of the density matrix will first be studied for both metals and insulators. Several strategies to construct O(N) algorithms will then be presented and critically examined. Some issues which are relevant only for self-consistent O(N) methods, such as the calculation of the Hartree potential and mixing issues, will also be discussed. Finally some typical applications of O(N) methods are briefly described.

A method is presented that can find the global minimum of very complex condensed matter systems. It is based on the simple principle of exploring the configurational space as fast as possible and of avoiding revisiting known parts of this space. Even though it is not a genetic algorithm, it is not based on thermodynamics. The efficiency of the method depends strongly on the type of moves that are used to hop into new local minima. Moves that find low-barrier escape-paths out of the present minimum generally lead into low energy minima.