А какие именно ограничения применимости здесь возникают?
Ну давайте посчитаем. В рамках классической теории при квазистационарном сжатии облака (если оно в каждый момент сжатия вириализуется) полная энергия равна половине потенциальной, а изменение потенциальной можно оценить как итоговую потенциальную после сжатия (начальная по модулю существенно меньше, ее отличием от нуля можно просто пренебречь). Итого с точностью до коэффициента порядка единицы (та самая половина и коэффициент, описывающий распределение плотности вдоль радиуса в облаке)

, где

- финальный радиус после сжатия. Приравниваем

, тогда

и итоговый радиус

. Иначе говоря, это радиус Шварцшильда (справа он и написан, только без коэффициента

).
Вывод - энергия порядка энергии покоя выделится как раз примерно тогда, когда объект сожмется в черную дыру. Очевидно, что на этой стадии классики недостаточно.