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