The OH vibrational frequency of four crystalline compounds ranging from ionic (brucite, Mg(OH) 2 , and portlandite, Ca(OH) 2 ) to semi-covalent (edingtonite, as representative of free surface OH groups in silica, and acid chabazite, as representative of acid zeolites) has been investigated at quantum mechanical level with the CRYSTAL program using the B3LYP hybrid functional. The OH vibration is calculated in two ways: (i) in the harmonic approximation, by diagonalizing the fully coupled dynamical matrix to yield the harmonic frequency ! h . (ii) at the anharmonic level, by decoupling the OH stretching mode from the bulk phonons and by numerically solving the one-dimensional Schro¨dinger equation associated with the OH potential energy to yield the fundamental ! 01 and the first overtone ! 02 frequencies. The harmonic and anharmonic frequencies differ by more than 150 cm À1 . In the cases where direct comparison is possible (brucite, portlandite and edingtonite), the experimental and calculated frequencies differ by less than 10 cm À1 ; the calculated anharmonicity constant, ! e x e ¼ (2! 01 À ! 02 )/2, is systematically smaller than the experimental value by about 10 cm À1 . The effect of the computational parameters on the computed frequencies is explored, with particular attention to the grid used for the construction of the DFT exchange and correlation contribution to the Hamiltonian and the accuracy in the geometry optimisation.