文摘
We present an analysis of the vibrational modes of a model of hexagonal ice, ice Ih, comprised of 192 monomers with a core region of 105 monomers, using the ab initio WHBB potential energy surface [Wang, Y.; Shepler, B.; Braams, B.; Bowman, J. M. J. Chem. Phys. 2011, 134, 094509]. A standard normal-mode analysis and a local-monomer normal-mode analysis of 105 core monomers are performed to obtain harmonic frequencies and state densities of the 鈥減seudo-translation鈥?(0鈥?00 cm鈥?), 鈥渓ibration鈥?(500鈥?100 cm鈥?), monomer bend fundamental (1600 cm), and O鈥揌 stretch (3000鈥?700 cm鈥?) bands. In addition, the coupled local-monomer model is used to obtain the vibrational density of states in the bend fundamental and O鈥揌 stretch regions. The harmonic and local-monomer vibrational density of states obtained from core monomers are in good agreement with those of inelastic neutron scattering spectra, especially the latter, which accounts for anharmonic coupling of monomer modes. Full deuteration is also considered, and the vibrational density of states is again compared to experiment, where good agreement is found.