Thermally stimulated currents in ice. Correction for orientational defects

The Bernal and Fowler model of ice [1] has been supplemented with modern ideas about relaxing defects arising as a result of migration of a proton along hydrogen bonds [2, 3]. We develop a theory for thermally stimulated currents (TSC) in hexagonal ice. We propose a kinetic equation taking into account relaxation of the structural defects supporting proton conductivity in the crystals. We obtain formulas qualitatively describing thermally stimulated depolarization currents and also a formula for calculating the activation energy for structural defects, allowing us to compare with experiment. Due to the unwieldiness of the calculations, we present only the derivation of the kinetic equations and analysis of the solution.