VASP计算振动能