Intermolecular vibrational modes and H-bond interactions in crystalline urea investigated by terahertz spectroscopy and theoretical calculation.

The characteristic absorption spectra of crystalline urea in 0.6-1.8 THz region have been measured by terahertz time-domain spectroscopy at room temperature experimentally. Five broad absorption peaks were observed at 0.69, 1.08, 1.27, 1.47 and 1.64 THz respectively. Moreover, density functional theory (DFT) calculation has been performed for the isolated urea molecule, and there is no infrared intensity in the region below 1.8 THz. This means that single molecule calculations are failure to predict the experimental spectra of urea crystals. To simulate these spectra, calculations on a cluster of seven urea molecules using M06-2X and B3LYP-D3 are performed, and we found that M06-2X perform better. The observed THz vibrational modes are assigned to bending and torsional modes related to the intermolecular H-bond interactions with the help of potential energy distribution (PED) method. Using the reduced-density-gradient (RDG) analysis, the positions and types of intermolecular H-bond interactions in urea crystals are visualized. Therefore, we can confirm that terahertz spectroscopy can be used as an effective means to detect intermolecular H-bond interactions in molecular crystals.