To obtain the transition rules in Sec. 3.2, the gradient approximation Eq.3.12 and the wavefunction Eq. 3.22 are used. The momentum matrix elements Eq. 3.23 are obtained as
| (A.1) |
where ⟨Ap|H|Bq⟩ = ⟨Bp|H|Aq⟩ = t for p = q and p = q ± 1, otherwise the matrix elements are zero. Therefore, Eq. A.1 can be written as
| (A.2) |
Here the relation sin + sin = 2 sin cos is employed. Using Eq. 3.21, Eq. A.2 can be written as
| (A.3) |
The summation over the sine functions in Eq. A.3 determines the transition rules. Using some trigonometric identities one can write this summation as
| (A.4) |
If n ± m = 2k + 1, where k is a non-zero integer, both terms in the bracket of Eq. A.4 will be zero. In the case of n±m = 2k, both terms in the bracket will be equal to -1, therefore, the summation will be again zero. However, if n = m, the fist term in will be equal to N and the second term will be equal to -1. Therefore, only transitions between valence and conduction subbands with the same band-index are allowed
| (A.5) |