Активная гравитационная масса этого этого облака оказывается
Неверно. Переход к формуле (21) неправильно обоснован - к облаку тех самых частиц эта формула отношения не имеет.
По сути вы правы. Здесь были допущены даже 2 ошибки. 1я в формуле (13). Она должна читаться как

(13).
Другая ошибка состоит в том, что расстояние (6) является расстоянием в системе отсчета, связанной с частицой облака, а при подсчете силы, действующей на тестовую частицу, надо перейти к расстоянию в системе отсчета, неподвижной относительно нее, K'=(t', x', y', z'). Не знаю, эти ли ошибки вы имели в виду. В результате активная гравитацитонная масса облака, состоящего из n релятивистских частиц, будет
где

это отношение скорости частиц в облаке к скорости света и

полный эллиптический интеграл 2-го рода,
https://osf.io/preprints/osf/9bkcz_v2.
Метрика (7) при малых

аппроксимируется канонической формой метрики постньютоновского пространства-времени

,
M. H. Soffel, Wen-Biao Han, Applied General Relativity: Theory and Applications in Astronomy, Celestial Mechanics and Metrology, Springer Nature, Cham, Switzerland, 2019,
при векторном потенциале

, описывающем гравитомагнитные эффекты, возникающие при движении материи.
Однако если рассматривать пару тел, описываемую метрикой (11), то координатные ускорения


в случае, если пространственный радиус-вектор частицы перпендикулярен линии движения тел (x'=0), ввиду

соответствуют активной гравитационной массе

.
Это будет выполняться и если частица расположена на линии движения тел (y'=z'=0) ввиду перехода к системе отсчета K'. Коэффициент

связывает пассивную гравитационную массу частицы и ее массу покоя при движении в слабом поле Шварцшильда в случае, если потенциал поля мал по сравнению с ее кинетической энергией. Этот результат получен при применении механики Лагранжа для анализа динамики частиц,
Application of Lagrangian mechanics to the analysis of particle dynamics in a gravitational field, PIRT- 2025, (4.26) , или В. Б. Беляев, Динамика в общей теории относительности: Вариационные методы, УРСС, 2017. Заметим, что в данном случае этот результат был получен без применения преобразований Лоренца к метрике Шварцшильда. К нему будет приближаться и масса облака

при малых

.