Site / Method  
AG Paulus Archive
AG Paulus Website

Method of increments  a wavefunctionbased ab initio correlation method for solids
The field of manybody theory has been developing at high speed in recent years. A significant amount of effort is directed toward affordable descriptions of electronic correlation in infinite periodic systems. The most widely used approach here is densityfunctional theory (DFT). However, the present DFTbased approaches have in practice many shortcomings. When applying quantumchemical correlation methods to extended systems one can make use of the property that the correlation hole is fairly local. For a solid it is useful, therefore, to switch from the description with delocalized Bloch states to localized Wannier orbitals. This is the startingpoint for the socalled local correlation methods. One method of this type, the method of increments [1], combines HF calculations for periodic systems with correlation calculations for finite embedded clusters, and the total correlation energy per unit cell of a solid is written as a cumulant expansion in terms of contributions from localized orbital groups of increasing size. A schematic overview of the procedure is given in Fig. 1. Since 1992, when H. Stoll has applied the idea of local excitations in the electronic theory of solids for considering correlation energy of diamond [1], silicon [2] and graphite [3] firstly, the method of increments has been successfully applied for various insulators and semiconductors ranging from strongly bound covalent compounds to weakly bound van der Waals raregas crystals (for a review see [4]). The problem to extend the method to metallic systems is twofold. The first problem consists in the modeling of a metallic solid by a finite cluster and connected with this the question about the localisation of orbitals in a metal arises. As a solution we suggest an embedding scheme [5] which has itself no metallic character but can mimic the metal in the internal region, where the atoms are correlated. Within the incremental scheme we allow for a delocalisation of the orbitals and therefore account for the metallic character of the systems. With this approach the convergence rate, especially with the order of increments, is not as good as for insulators, but still three and higherbody terms account for only about 15 % of the correlation contribution to the binding. Up to now we have carefully described closed shell metals [613]. However, it is possible to extend this approach also to openedshell cases. So, the correlation energy of a metallic Liring has been considered applying a projection technique for nonorthogonal occupied and virtual orbitals within the incremental scheme [14,15]. Even an embedding scheme was suggested for 3dimensional lithium relying on pairs of atoms [16]. In the case of lithium the s band is crossing the Fermi level and is half occupied. In this regard the second problem occurring in metals has to be mentioned, namely the quasidegeneracy of the orbitals around the Fermi surface. Then a singlereference correlation method can fail and the incremental scheme has to be extended to the use of multireference correlation treatments. This can be the case not only in openedshell systems, but also, for example, in barium [17], where the dbands are crossing the Fermi level. Thus, with the method of increments we have a tool, which allows the evaluation the groundstate properties of metals. Starting from the HF solution for the extended systems, which is a routine calculation nowadays, we can localize the orbitals and perform an incremental expansion for the correlation energy. An asset of the incremental scheme is that applying it one gets a possibility not only to obtain systematically improvable result, but also to investigate influence of individual contributions on the correlation energy. This opens up a large field of applications. With such a method at hand we can try to understand unusual structures of some metals. Preliminary work has been done in our laboratory for gallium [18], which forms Gadimers within the metallic phase. One more class of the object affecting the great interest is intermetallic compounds, especially Laves phases [19]. Here our approach could be useful for understanding their stability and phase formation as a function of components. An advantage of the method of increments in comparison with other local correlation methods is that the translational symmetry of the system can be used to reduce the effort of the calculation. Translational symmetry is, however, not a prerequisite of the method. Therefore it is possible to extend the method to systems with reduced symmetry like surfaces and molecules on the surface [20]. Impurities and functional groups in large biological molecules can be described as well. If the description for metals and the one for impurities are combined, the whole field of Kondo physics opens up. But one has to be aware of the fact that for the small energy scales involved in the Kondo effect require highly accurate quantumchemical methods, and very extended basis sets are necessary to achieve reliable results. [1] H. Stoll, Phys. Rev. B 46, 6700 (1992). 
↑ Top
View
Edit
Attributes
History
Attach
Print
This page was last modified on April 27, 2009, at 01:14 PM 