Plane-wave basis

There are a variety of other approaches to understanding the electronic structure of crystals. Most of them rely on a density functional approach, with or without the pseudopotential, and use different bases. For example, instead of a plane wave basis, one might write a basis composed of atomic-like orbitals ... 

The higher-order bulk contribution to the nonlmear response arises, as just mentioned, from a spatially nonlocal response in which the induced nonlinear polarization does not depend solely on the value of the fiindamental electric field at the same point. To leading order, we may represent these non-local tenns as bemg proportional to a nonlinear response incorporating a first spatial derivative of the fiindamental electric field. Such tenns conespond in the microscopic theory to the inclusion of electric-quadnipole and magnetic-dipole contributions. The fonn of these bulk contributions may be derived on the basis of synnnetry considerations. As an example of a frequently encountered situation, we indicate here the non-local polarization for SFIG in a cubic material excited by a plane wave (co) ... 

Jordan oompared the use of plane wave and oonventional Gaussian basis orbitals within density funotional oaloulations in ... 

Figure B3.2.4. A schematic illustration of an energy-independent augmented plane wave basis fimction used in the LAPW method. The black sine fimction represents the plane wave, the localized oscillations represent the augmentation of the fimction inside the atomic spheres used for the solution of the Sclirodinger equation. The nuclei are represented by filled black circles. In the lower part of the picture, the crystal potential is sketched.

The projector augmented-wave (PAW) DFT method was invented by Blochl to generalize both the pseudopotential and the LAPW DFT teclmiques [M]- PAW, however, provides all-electron one-particle wavefiinctions not accessible with the pseudopotential approach. The central idea of the PAW is to express the all-electron quantities in tenns of a pseudo-wavefiinction (easily expanded in plane waves) tenn that describes mterstitial contributions well, and one-centre corrections expanded in tenns of atom-centred fiinctions, that allow for the recovery of the all-electron quantities. The LAPW method is a special case of the PAW method and the pseudopotential fonnalism is obtained by an approximation. Comparisons of the PAW method to other all-electron methods show an accuracy similar to the FLAPW results and an efficiency comparable to plane wave pseudopotential calculations [, ]. PAW is also fonnulated to carry out DFT dynamics, where the forces on nuclei and wavefiinctions are calculated from the PAW wavefiinctions. (Another all-electron DFT molecular dynamics teclmique using a mixed-basis approach is applied in [84].)... 

Plane waves are often considered the most obvious basis set to use for calculations on periodic sy stems, not least because this representation is equivalent to a Fourier series, which itself is the natural language of periodic fimctions. Each orbital wavefimction is expressed as a linear combination of plane waves which differ by reciprocal lattice vectors ... 

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... 

It is not easy to see why the authors believe that the success of orbital calculations should lead one to think that the most profound characterization of the properties of atoms implies such an importance to quantum numbers as they are claiming. As is well known in quantum chemistry, successful mathematical modeling may be achieved via any number of types of basis functions such as plane waves. Similarly, it would be a mistake to infer that the terms characterizing such plane wave expansions are of crucial importance in characterizing the behavior of atoms. 

The KS equation (Eq. 23) when expressed in terms of a plane-wave basis set takes a very simple form ... 

In order to determine the relations among orbital exponents in a basis which will follow these guidelines, we look at the matrix elements contributing to on(q). To that end we consider eq. 17 for the plane wave operator (eq. 21) which involves evaluation of terms of the sort (p e i /i). We wish to determine how these matrix elements behave as a function of orbital exponent and momentum transfer, and we then propose a scheme for choice of orbital exponents that will keep the BSR satisfied to as high momentum transfer as possible. 

Kresse G, Furthmuller J. 1996a. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci 6 15-50. 

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). 

We should also mention that basis sets which do not actually comply with the LCAO scheme are employed under certain circumstances in density functional calculations, i. e., plane waves. These are the solutions of the Schrodinger equation of a free particle and are simple exponential functions of the general form... 

A substantial improvement in the theoretical modeling of this substrate, able to overcome the drawbacks of both former approaches (standard GGA and localized atomic basis set B3LYP) is the implementation of the HSE03 [36] hybrid functionals in a plane-wave formalism. 

For the conduction electrons, it is reasonable to consider that the inner-shell electrons are all localized on individual nuclei, in wave functions very much like those they occupy in the free atoms. The potential V should then include the potential due to the positively charged ions, each consisting of a nucleus plus filled inner shells of electrons, and the self-consistent potential (coulomb plus exchange) of the conduction electrons. However, the potential of an ion core must include the effect of exchange or antisymmetry with the inner-shell or core electrons, which means that the conduction-band wave functions must be orthogonal to the core-electron wave functions. This is the basis of the orthogonalized-plane-wave method, which has been successfully used to calculate band structures for many metals.41... 

Formally, each orthogonalized-plane-wave basis function may be written as (1 - P), where ijjk is a plane wave and P is the projection operator such that Pif/k gives the core-state component of Il>k ... 

---