Brillouin Zone integration over

The matrices A, B, Q are of infinite dimension since there are an infinite number (2N+1, N — oo) of k-values and thus an infinite number of k-states in each band. Moreover, there is an equation for each triplet formed by a k-value and two band indices. This triplet represents a particle-hole excitation that is vertical in order to preserve the momentum. As is the case in many polymeric techniques, the infinite sum over k is transformed into an integration in the first Brillouin zone ... 

In the approximation of a circular Brillouin zone, the desired sum over k in Eq. (2.2.15), if substituted by the integral over the angle between k and, acquires the form66... 

For the value FA to be calculated, the integration over the whole Brillouin zone should be performed which yields, in the limiting case of =90°, the value FA -5.44/Qijp (in the framework of the isotropic approximation). 

Fig. 25. Structure factor integrated over 2.3% of the surface Brillouin zone (radius of 5 mesh lengths in Fig. 24) vs. TtoclesL plotted with the rescaled energy (crosses) for the J i x overlayer on the triangular lattice. Rescaling involves multiphc-ation by a negative number and shifting by a constant. Temperature is measured in units of

The key features of this integral are that it is defined in reciprocal space and that it integrates only over the possible values of k in the Brillouin zone. Rather than examining in detail where integrals such as these come from, let us consider instead how we can evaluate them numerically. 

In Chapter 3, we discussed choosing k points for integrating over the Brillouin zone. The Monkhorst-Pack method for choosing k points may be used when... 

The total number of transitions per unit volume and time from Eq. (17) is obtained by dividing by t and summing over j (occupied bands) and / (unoccupied bands). Integrating over the Brillouin zone gives ... 

The first step is to calculate over the energy bands a function referred to as the energy distribution of the joint density of states. For polymers it amounts to a one-dimensional integration over the Brillouin Zone (BZ) of a rather complicated integrand such as found in the theoretical evaluation of optical spectra. For sufficiently high photon energies the final states are rather unstructured and the function to be compared with experiment may be simplified into occ... 

As an example. Fig. 3 plots the phonon dispersion curves for three highly S5mimetric directions in the Brillouin zone of the perfect ZnO crystal. Comparison of the theoretical and experimental frequencies shows good agreement for the acoustic branches. The densities of phonon states of the perfect ZnO crystal calculated by integrating over the Brillouin zone are displayed in Fig. 4. Comparison of the results of our calculation and a calcu-... 

Z) = 0.142 nm is the interatomic distance in graphene. Integration and summation in Eq. (3) is taken over the first Brillouin zone of the CNT which is the set of discrete lines due to the quantization of the transversal quasi-momentum [4]. 

The use of twisted boundary conditions is commonplace for the solution of the band structure problem for a periodic solid, particularly for metals. In order to calculate properties of an infinite periodic solid, properties must be averaged by integrating over the first Brillouin zone. 

The dynamic problem of vibrational spectroscopy must be solved to find the normal coordinates as linear combinations of the basis Bloch functions, together with the amplitudes and frequencies of these normal vibrations. These depend on k, and therefore the problem must be solved for a number of k-points to ensure an adequate sampling of the Brillouin zone. Vibrational frequencies spread in k-space, just as the Bloch treatment of electronic energy gave a dispersion of electronic energies in k-space. The number of vibrational levels whose energy lies between E and fc +d E is called the vibrational density of states. Vibrational contributions to the heat capacity and to the crystal entropy can be calculated by appropriate integrations over the vibrational density of states, just like molecular heat capacities and entropies are obtained by summation over molecular vibration frequencies. 

For this PP calculation the exact Ex has been combined with the LDA for Ec. In spite of the use of the exact exchange functional, FeO is predicted to be a metal, in contrast to experiment. At present, it is not clear whether this failure to reproduce the insulating ground state of FeO originates from the use of the LDA for Ec or from the technical limitations of the PP calculation (KLI approximation for vx, only three special it-points for the integration over the Brillouin zone, 3s elections in the core). It must be emphasized, however, that the band structure shown in Figure 4.6 is rather different from its LDA counterpart (Dufek et al. 1994), which emphasizes the importance of the exact vx for this system. 

If N -> oo the matrices , , X, Y and ft are of infinite dimension. Taking into account the periodic symmetry of the polymer the infinite sums over (crystal momentum) can be transformed into an integration over in the first Brillouin zone giving... 

Here the sum over k points has become an integral over the first Brillouin zone (with volume Vbz), because it has already been shown that k can be considered as a continuous variable. By Hmiting the integration to states with energy below Ep, a Heaviside step function 0 permits us to exclude the eigenvectors relative to empty states from the sum. The reason why this cannot be achieved by simply truncating the sum over the eigenvectors, like in the molecular case, will be clear when the main features of band structure are illustrated later on in this chapter. 

Solving the HF or KS equations in the present (CRYSTAL) scheme requires numerical integration over the first Brillouin zone because, in general, we do not possess an analytic expression for the eigenvalues and eigenvectors, as is the case of Hiickel s approximation. The question then becomes How many points need to be sampled, that is, in how many points must Eq. [25] be solved to get sufficiently accurate values of the observables of interest ... 

Linearized Augmented Plane Wave Method 1127 The integrals over the Brillouin zone are given by Vo... 

---