Here, R Iα is the α = x, y, z Cartesian component of the position of the I-th nucleus and m I is its atomic mass. By diagonalizing the mass-weighted Hessian, one obtains the eigenvalues ω 2 p , which correspond to the squared angular frequencies of the harmonic vibrations, as well as the corresponding normal modes Q p .