where k ≡ ≡ ( k x , k y , k z ) ; V ≡ ≡ L 3 is the size of the system; and the integral is (as the summation) over all possible modes, which is assumed to be a finite region (bounded by the cut-off frequency).