Pseudopotential techniques pseudopotentials

The pseudopotential density-functional technique is used to calculate total energies, forces on atoms and stress tensors as described in Ref. 13 and implemented in the computer code CASTEP. CASTEP uses a plane-wave basis set to expand wave-functions and a preconditioned conjugate gradient scheme to solve the density-functional theory (DFT) equations iteratively. Brillouin zone integration is carried out via the special points scheme by Monkhorst and Pack. The nonlocal pseudopotentials in Kleynman-Bylander form were optimized in order to achieve the best convergence with respect to the basis set size. 5... 

Table 2 Lattice parameter ao (in A) and elastic constants B, Cn, C12, and C44 for CoSi2 (in GPa), calculated using all-electron (FLAPW) and pseudopotential (VASP) techniques in the LDA and using GGC corrections.

The systems discussed in this chapter give some examples using different theoretical models for the interpretation of, primarily, UPS valence band data, both for pristine and doped systems as well as for the initial stages of interface formation between metals and conjugated systems. Among the various methods used in the examples are the following semiempirical Hartree-Fock methods such as the Modified Neglect of Diatomic Overlap (MNDO) [31, 32) and Austin Model 1 (AMI) [33] the non-empirical Valence Effective Hamiltonian (VEH) pseudopotential method [3, 34J and ab initio Hartree-Fock techniques. 

In molecular DFT calculations, it is natural to include all electrons in the calculations and hence no further subtleties than the ones described arise in the calculation of the isomer shift. However, there are situations where other approaches are advantageous. The most prominent situation is met in the case of solids. Here, it is difficult to capture the effects of an infinite system with a finite size cluster model and one should resort to dedicated solid state techniques. It appears that very efficient solid state DFT implementations are possible on the basis of plane wave basis sets. However, it is difficult to describe the core region with plane wave basis sets. Hence, the core electrons need to be replaced by pseudopotentials, which precludes a direct calculation of the electron density at the Mossbauer absorber atom. However, there are workarounds and the subtleties involved in this subject are discussed in a complementary chapter by Blaha (see CD-ROM, Part HI). 

Recent calculations of hyperfine parameters using pseudopotential-density-functional theory, when combined with the ability to generate accurate total-energy surfaces, establish this technique as a powerful tool for the study of defects in semiconductors. One area in which theory is not yet able to make accurate predictions is for positions of defect levels in the band structure. Methods that go beyond the one-particle description are available but presently too computationally demanding. Increasing computer power and/or the development of simplified schemes will hopefully... 

Some of the major areas of activity in this field have been the application of the method to more complex materials, molecular dynamics, [28] and the treatment of excited states. [29] We will deal with some of the new materials in the next section. Two major goals of the molecular dynamics calculations are to determine crystal structures from first principles and to include finite temperature effects. By combining molecular dynamics techniques and ah initio pseudopotentials within the local density approximation, it becomes possible to consider complex, large, and disordered solids. 

There have been a number of improvements in techniques, and more convenient models have been formulated however, the basic approach of the pseudopotential total energy method has not changed. This general approach or standard modd is applicable to a broad spectrum of solid state problems and materials when the dec-trons are not too localized. Highly correlated electronic materials require more attention, and this is an area of active current research. However, considering the extent of the accomplishments and die range of applications (see Table 14.3) to solids, dusters, and molecules, this approach has had a major impact on condensed matter physics and stands as one of the pillars of the fidd. 

The question of methanol protonation was revisited by Shah et al. (237, 238), who used first-principles calculations to study the adsorption of methanol in chabazite and sodalite. The computational demands of this technique are such that only the most symmetrical zeolite lattices are accessible at present, but this limitation is sure to change in the future. Pseudopotentials were used to model the core electrons, verified by reproduction of the lattice parameter of a-quartz and the gas-phase geometry of methanol. In chabazite, methanol was found to be adsorbed in the 8-ring channel of the structure. The optimized structure corresponds to the ion-paired complex, previously designated as a saddle point on the basis of cluster calculations. No stable minimum was found corresponding to the neutral complex. Shah et al. (237) concluded that any barrier to protonation is more than compensated for by the electrostatic potential within the 8-ring. 

The first explanation and use of such a pseudopotential is due to Heilman5 (1935) who used it in atomic calculations. More recently the pseudopotential concept was reformulated by Phillips and Kleinman7 who were interested in its application to the solid state.8-10 Research in both solid- and liquid-state physics with pseudopotentials was reviewed by Ziman,11 and work in the fields of atomic spectroscopy and scattering has been discussed by Bardsley.12 For an earlier review on applications to the molecular environment the reader is referred to Weeks et a/.13 In this article we shall concentrate on molecular calculations, specifically those of an ab initio nature. Our objective in Section 2 has been to outline the theoretical origins of the pseudopotential approximation, and in Section 3 we have described some of the techniques which have been used in actual calculations. Section 4 attempts to present results from a representative sample of pseudopotential calculations, and our emphasis has been to concentrate on particular molecules which have been the subjects of investigation by the various approaches, rather than to catalogue every available calculation. Finally, in Section 5, we have drawn some conclusions on the relative merits of the different methods and implementations of pseudopotentials. Some of the possible future developments are outlined in the context of the likely progress in quantum chemistry. 

However, this is definitely the technique for future calculations involving a large number of metal atoms. Furthermore, the idea behind the pseudopotential method is also applied in other types of Hamiltonians described below, e.g., valence effective Hamiltonian and semi-empirical methods. 

Further uses of pseudopotentials are numerous. The most obvious (and rather widely known ones) are to continue with the PP or MP Hamiltonians for a widely understood combination of the core and valence shells and to apply standard ab initio techniques to electrons in the valence subspace only. We do no elaborate further on this as the hybrid nature of the pseudopotential methods is rather obvious from the above and its more specific applications in a narrower QM/MM hybrid context will be described later. 

Electronic structure methods for studies of nanostructures can be divided broadly into supercell methods and real-space methods. Supercell methods use standard k-space electronic structure techniques separating periodically repeated nanostructures by distances large enough to neglect their interactions. Direct space methods do not need to use periodic boundary conditions. Various electronic structure methods are developed and applied using both approaches. In this section we will shortly discuss few popular but powerful electronic structure methods the pseudopotential method, linear muffin-tin orbital and related methods, and tight-binding methods. 

An increase in speed can be achieved by using pseudopotential calculations122151 >. In these type of calculations the inner shell electrons are approximated with a potential and the problem is reduced to a valence electron problem. This technique is very powerful for heavy atoms but the time saving is not more than roughly 50 % for molecules containing first-row atoms only 152). Ion-ligand interactions have been studied with pseudopotential calculations in several cases 153 157). 

As an example we consider the system Si/Si(00 1). Over the years there have been several first principles DFT calculations aimed at establishing the relative stability of these structures [75, 77-80] the results are collected in Table 1. All of these DFT calculations make use of a supercell technique, pseudopotentials to represent the Si ion cores, and a basis set that consists of plane waves to expand the valence wave functions. Car-Parrinello techniques are used to simultaneously optimize the electronic... 

Except for occasional discussions of the basis set dependence of the results, the numerical implementation issues such as grid integration techniques, electron-density fitting, frozen-cores, pseudopotentials, and linear-scaling techniques, are omitted. 

Hutter J, Tuckerman ME, Parrinello M, Integrating the Car-Parrinello equations III Techniques for ultrasoft pseudopotentials, J Chem Phys, 102, 859 (1995)... 

Other complications are associated with the partitioning of the core and valence space, which is a fundamental assumption of effective potential approximations. For instance, for the transition elements, in addition to the outermost s and d subshells, the next inner s and p subshells must also be included in the valence space in order to accurately compute certain properties (54). A related problem occurs in the alkali and alkaline earth elements, involving the outer s and next inner s and p subshells. In this case, however, the difficulties are related to core-valence correlation. Muller et al. (55) have developed semiempirical core polarization treatments for dealing with intershell correlation. Similar techniques have been used in pseudopotential calculations (56). These approaches assume that intershell correlation can be represented by a simple polarization of one shell (core) relative to the electrons in another (valence) and, therefore, the correlation energy adjustment will be... 

---