A semiempirical molecular dynamics model is developed to study the vibrational frequencies of uracil at very low kinetic temperature by using the Fourier transform of velocity autocorrelation function of trajectories of molecular dynamics simulations. The finite difference harmonic method is used to assign the vibrational frequency of each mode. The calculated frequencies are found to be in good agreement with experimental measurements. Moreover, we make up for the lost vibrational modes in experiments self-consistently. A total of 30 vibrational modes and their corresponding frequencies are reported.