Спектральные наблюдения рентгеновских пульсаров, получаемые на обсерваториях NuSTAR и Integral, по-прежнему требуют теоретической интерпретации. Существующие несамосогласованные модели оказываются недостаточными как по причине отсутствия хорошего согласия с наблюдениями, так и вследствие достижения такого согласия с использованием неоправданно большого числа модельных параметров. Поэтому настоящий проект направлен на создание самосогласованной модели радиационно-доминированной аккреционной колонки, в которой формируются спектры рентгеновских пульсаров при темпах аккреции, превышающих $\\sim 10^{17}$~г$/$c.\r\nНовизна исследования заключается в решении кинетического уравнения для излучения одновременно с уравнением Эйлера, уравнением неразрывности и уравнением для электронной температуры. Такой подход позволяет проводить вычисления с использованием сечений рассеяния и поглощения, зависящих от частоты, и учитывать поляризационные эффекты, возникающие в сверхсильном магнитном поле.\r\nСопоставление ожидаемых результатов проекта с современными наблюдательными данными позволит сделать важные выводы о физической природе излучающих областей рентгеновских пульсаров.