В статье исследуется вклад поляризации пространственного заряда в комплексный нелинейный спектр диэлектрической проницаемости нематического жидкого кристалла (ЖК). Используя модель, основанную на системе уравнений Пуассона - Нернста - Планка (ПНП), был численно исследован процесс электромиграции примесных ионов под действием электрического поля. Результаты моделирования показывают, что при низких частотах (<10 Гц) и высоких приложенных электрических напряжениях (превышающих тепловое электрическое напряжение) происходит перераспределение ионов в области около электродов, что приводит к появлению нечетных гармоник в комплексном диэлектрическом спектре ЖК ячейки. Полученные данные могут быть использованы для идентификации низкочастотных релаксационных процессов в ЖК.