where N is the number of vibrational modes (equal to three times the number of atoms), ω the vibrational frequency, T the temperature, and h and k B Planck's and Boltzmann constants, respectively. To evaluate the heat capacity, the lattice vibrational frequencies can be extracted from the Hessian matrix, that is, the second-order partial derivatives of the lattice energy with respect to atom displacements.