Vibrational corrections to the Verdet constants of nine molecules (H-2, N-2, CO, H2O, CH4, benzene, toluene, p-xylene, and o-xylene) were calculated with pure density functional theory (DFT), hybrid DFT, and an approximate coupled-cluster theory. Comparisons are made for the accuracy of the vibrational averages among different methods and with respect to experimental data where available. It is found that vibrational corrections to magneto-optical rotation can be as large as 10% of the equilibrium value. Hybrid DFT with the B3LYP hybrid functional offers reasonable accuracy at a relatively inexpensive computational cost for accurate calculations of vibrationally averaged Verdet constants.