The vibrational spectrum of ice XI at thermal wavelengths using the CASTEP code, a first-principles simulation method, is investigated. A dual-track approach is constructed to verify the validity for the computational phonon spectrum: collate the simulated spectrum with inelastic neutron scattering experiments and assign the photon scattering peaks according to the calculated normal vibration frequencies. The 33 optical normal vibrations at the Brillouin center are illustrated definitely from the ab initio outcomes. The depolarizing field effect of the hydrogen bond vibrations at frequencies of 229 cm(-1) and 310 cm(-1) is found to agree well with the LST relationship. It is a convincing evidence to manifest the LO-TO splitting of hydrogen bonds in ice crystal. We attribute the two hydrogen bond peaks to the depolarization effect and apply this viewpoint to ordinary ice phase, ice Ih, which is difficult to analyse their vibration modes due to proton disorder.