Приводятся результаты математического моделирования процесса фотоизомеризации многоатомных молекул. Динамика заселённости состояний описана на основе численного решения уравнения Лиувилля для оператора плотности молекулы при феноменологическом учёте её возбужденных состояний и на основе решений управляющего уравнения общей теории релаксации в форме Линдблада. Приводятся примеры эволюции заселённостей состояний, которая характеризуется много меньшей скоростью рассеяния светового импульса и послесвечения по сравнению со скоростью рассеяния и послесвечения соответствующей двухуровневой молекулы, то есть примеры динамики заселённости так называемых "тёмных" состояний.