Quantum mechanical calculation of the OH vibrational frequency in crystalline solids

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 omega(h). (ii) at the anharmonic level, by decoupling the OH stretching mode from the bulk phonons and by numerically solving the one-dimensional Schrodinger equation associated with the OH potential energy to yield the fundamental omega(01) and the first overtone omega(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, omega(e)x(e) (2 omega(01) - omega(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.